I want to do a Fourier transformation of an interpolation function that is the result of NDsolve along one variable. Currently I am using the following approach, which numerically integrates the kernel
fourier[f_,t_,t1_,t2_,Nt_,w0_,T0_,w_]:=Module[{h=(t2-t1)/Nt,f2=f*Exp[I*(w-w0)*t*T0]},
1/Sqrt[2*\[Pi]]* h/2*(N[f2/.t->t1]+N[f2/.t->t2]+2*Sum[N[f2/.t->t1+j*h],{j,1,Nt-1}])
]
In the end, this is trapezoidal rule (example code is attached). Is there any approach which can do faster numerical integration on interpolation functions?
Attachments: