i have a signal with many frequencies. how to find this frequencies and reconstruct the signal from these frequencies

Question

rohith bharadwaj 2017-8-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/352240-i-have-a-signal-with-many-frequencies-how-to-find-this-frequencies-and-reconstruct-the-signal-from

评论： rohith bharadwaj 2017-8-29

fivelamda.m

I've attached a MATLAB code, in that timeseries is the signal for which frequencies have to be found. I've tried using FFT but not getting the answer. my code for frequency and regeneration: t =0:1:2048; y=timeseries;

ffty = fft(y);

ffty = abs(ffty(1:ceil(length(y)/2))); % ffty(ffty<2)=0; [~,locs] = findpeaks(ffty); freqs = (locs-1)/t(end); signal_1=0; for n=1:1:length(freqs) signal_1 = (signal_1+sin(2*pi*freqs(n)*t)); end hold on plot(smooth(signal_1));

hold on plot(y,'b');

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Manan Mishra 2017-8-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/352240-i-have-a-signal-with-many-frequencies-how-to-find-this-frequencies-and-reconstruct-the-signal-from#answer_278258

在 MATLAB Online 中打开

I think this should work for you:

ffty = fft(y);
P2 = abs(ffty/length(ffty));
P1 = P2(1:ceil(length(y)/2)); 
P1(2:end-1) = 2*P1(2:end-1);
[amp,locs] = findpeaks(P1); 
freqs = (locs-1)/t(end); 
signal_1=0; 
for n=1:length(freqs) 
    signal_1 = (signal_1+amp(n)*sin(2*pi*freqs(n)*t)); 
end
hold on;
plot(smooth(signal_1));
plot(y,'b');

Here, P2 is the two-sided spectrum and P1 is the single-sided spectrum.

A two-sided power spectrum displays half the energy at the positive frequency and half the energy at the negative frequency. Therefore, to convert a two-sided spectrum to a single-sided spectrum, you discard the second half of the array and multiply every point except for DC by two.

Also, while regenerating the signal in the "for" loop, you need to take the amplitude of corresponding frequencies into consideration. Hence, you need the first output "amp" of "findpeaks".

You can also eliminate the "for" loop by using vectorization:

signal_1=sum(amp'.*(sin(2*pi*freqs'*t)));