ex_04_09_1.m

来自「清华大学胡广书教授< 数字信号处理理论算法实现>>第二版一书的」· M 代码 · 共 23 行

%  ex_04_09_1.m :  习题 4.9

clear all;

N=4096;
 f1=.1;
 fs=1;
 a1=5;
 w=2*pi/fs;
 x=a1*sin(w*f1*(0:N-1))+randn(1,N);
 
 % 应用FFT 求频谱；
 figure('color','w')
 subplot(3,1,1);
 plot(x(1:N/100));
% title('x(n)');
 f=-0.5:1/N:0.5-1/N;
 X=fft(x); 
 %
 subplot(3,1,2);
 plot(f,fftshift(abs(X)));
 %title('|X(k)|');