2

在下面编写的代码图中，我应该看到 75 到 200 uV 值之间的幅度，但是我的图中有些奇怪，你能帮我并建议，我的代码有什么问题，我该如何纠正？

谢谢你们

这是代码：

clear all;
close all;
Fs = 200
t= 0:1/Fs:180
y_in=zeros (1, length(t));
  for  i = 1:18
F = randi ([4 7], 1);% frequency
A = randi ([75 200],1);% amplitude A=75~200 uV.
y_tmp =A*sin (2*pi*F*t);
 y_in=y_in+ y_tmp;
  end;
L=length (y_in);
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y_new = fft(y_in,NFFT)/L;
f_new = Fs/2*linspace(0,1,NFFT/2+1);

figure(1)
plot(f_new,2*abs(Y_new(1:NFFT/2+1))) ;
title('Single-Sided Amplitude Spectrum of y_in(t)')
xlabel('Frequency (Hz)')
ylabel('|y_in(f)|')
4

1 回答 1

1

由于您循环了 18 次并使用

F = randi([4 7],1);

建立最终信号的某些频率会多次出现，因此它们的幅度会相互叠加。

因此，您在绘图中的 fft 之后看到的是频率幅度的总和

F = 4, 5, 6, 7 Hz
于 2012-06-11T21:08:12.273 回答