I want to find the peak values of the out put signal.

Question

rohail khan 2018-4-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/393546-i-want-to-find-the-peak-values-of-the-out-put-signal

回答： Birdman 2018-4-9

I have generated a waveform ,added some noise to it and then used fft to recover the signal. The code is below

% Generate example time series
t = 0:1/1000:0.2-1/1000; 
% Sampled at 1 ms for .2 sec. Make the sample nuber even. 
% This works better for FFT
% 3 freqeuncies
w1 = 2*pi*60;
w2 = 2*pi*71;
w3 = 2*pi*89;
% 3 amplitudes
A1 = 1;
A2 = 1/2;
A3 = 1/3;
% generate the three sinusoides
f1 = A1*sin(w1*t);
f2 = A2*sin(w2*t);
f3 = A3*sin(w3*t);
f=f1+f2+f3;
Then, add noise to the amplitude of each sinusoid and shift its phase randomly. I will leave the variance of the noise to your decision. Plot the resulting sum of the three waves.
% 3 seperate frequency curve
plot(t,f1,t,f2,t,f3)
title ('3 seperate curves')
xlabel('sec')
hold off
% 3 combined  frequency plot
plot(t,f)
title('Combined frequency plot')
xlabel('sec') % f is the combined wave 
hold off
Use fft to find the frequencies, 60 Hz, 71Hz, and 89 Hz.
F = fft(f);
W = (0:length(F)-1)*1000/length(F); 
% Complet the rest
len=length(t);
fshift=(-len/2:len/2-1)*(1000/len)
yshift=fftshift(F)
plot(fshift,abs(yshift)./100)
title('recovered')
xlabel('Hz')
ylabel('amplitude')
hold off

I get the waveform

I can see the peak values from the plot but I want to use findpeaks command to get them.Each time I try to use findpeaks, I get an error message. How can I find the peak values using findpeaks in this question