example4_8.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 19 行

%例4-8;example4_8
figure(1)
subplot(2,2,1)
N=64;n=0:N-1;t=0.01*n;q=n*2*pi/N;
x=2*sin(4*pi*t)+5*cos(8*pi*t);
plot(x)
title('原始信号')
subplot(2,2,3)
Y=fft(x,N);
plot(q,abs(Y))
title('FFT N=64')
x1=2*sin(4*pi*t)+5*cos(8*pi*t)+2.5*randn(1,N);
subplot(2,2,2)
plot(x1)
title('加噪信号')
subplot(2,2,4)
Y1=fft(x1,N);
plot(q,abs(Y1))
title('加噪FFT  N=64')