设信号为x(t)=sin(2πf1t)+sin(2πf2 t)+随机噪声,f1=50Hz,f2=120Hz,以取样频率200hz和1000hz对x(t)进行取样,并绘制x(t)和采样后的信号幅x(t)和采样后的信号幅频特性
代码供参考:
Fs = 200; % 采样频率 1000Hz自己改Fs=1000
T = 1/Fs; % 采样时间间隔
L = 1500; % 信号长度
t = (0:L-1)*T; % 时间向量
f1 = 50;
f2 = 120;
X = sin(2*pi*f1*t) + sin(2*pi*f2*t) + randn(1,L)*0.3;%加上随机噪声
Y = fft(X);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L;
figure(1);clf%画图
subplot(2,1,1)
plot(t,X,'r-')
title('X(t)幅值图')
xlabel('时间t')
ylabel('幅值x(t)')
subplot(2,1,2)
plot(f,P1,'b-')
title('x(t)的频谱图')
xlabel('f (Hz)')
ylabel('|P1(f)|')
200Hz采样频率
1000Hz采样频率
可见200Hz采样频率的频域图不能正确反映120Hz的频率峰值,1000Hz可以都反映出来