![1how to do Fourier transform with current data - 2019 12 28.png](https://www.mathworks.com/matlabcentral/answers/uploaded_files/257529/1how%20to%20do%20Fourier%20transform%20with%20current%20data%20-%202019%2012%2028.png)
how to do Fourier transform with current data
5 views (last 30 days)
Show older comments
dang khoa tran
on 29 Dec 2019
Commented: Star Strider
on 29 Dec 2019
Hello,
I have a set of data as shown in the attactment, in which the "t_spectrum" shows the time and "y_spectrum" shows the electric field amplitude. And from this I want to do Fourier transform from these data. I have tried fft function but it seems that the spectrum is not reasonable.
So could anyone please recommend some approaches to that?
Thank you very much
0 Comments
Accepted Answer
Star Strider
on 29 Dec 2019
Try this:
E = dlmread('y_spectrum.txt',' ',1,0);
t = dlmread('t_spectrum.txt',' ',1,0);
Ts = mean(diff(t));
Fs = 1/Ts;
Tstd = std(diff(t));
Fs = 12E+15; % Sampling Frequency
[Er,tr] = resample(E,t,Fs); % Resample To Regularly-Spaced Intervals
figure
plot(t,E,'.b')
hold on
plot(tr,Er,'-r')
hold off
grid
legend('Original E', 'Resampled E')
Fn = Fs/2; % Nyquist Frequency
L = numel(tr); % Vector Lengths
Erm = Er - mean(Er);
FT_Erm = fft(Er)/L;
Fv = linspace(0, 1, fix(L/2)+1)*Fn;
Iv = 1:numel(Fv);
[pks,locs] = findpeaks(abs(FT_Erm(Iv))*2, 'MinPeakHeight',3E+10);
figure
plot(Fv, abs(FT_Erm(Iv))*2)
hold on
plot(Fv(locs), pks, '.r')
hold off
grid
text(1.5E+15, 7.5E+10, sprintf('Peak %10.3E Frequency %10.3E\n', [pks Fv(locs).'].'))
producing:
![1how to do Fourier transform with current data - 2019 12 28.png](https://www.mathworks.com/matlabcentral/answers/uploaded_files/257529/1how%20to%20do%20Fourier%20transform%20with%20current%20data%20-%202019%2012%2028.png)
I am not certain what does not appear to be ‘reasonable’ about it. It looks good to me.
0 Comments
More Answers (1)
See Also
Categories
Find more on Transforms in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!