sy11.m

来自「数字信号处理工具箱」· M 代码 · 共 20 行

N=40;n=[0:N-1];
t=0.01*n;
x=2*sin(4*pi*t)+5*cos(8*pi*t);
k=[0:1:N/2];w=2*pi/N*k;
X=fft(x,N)'
magY=abs(X(1:1:N/2+1));
subplot(2,2,1);stem(n,x,'.');title('signal x(n)');grid;
subplot(2,2,2);plot(w/pi,magY);title('FFT N=40');
xlabel('频率(单位：pi)');ylabel('|X|');grid

N=128;n=[0:N-1];
t=0.01*n;
x=2*sin(4*pi*t)+5*cos(8*pi*t);
k=[0:1:N/2];w=2*pi/N*k;
X=fft(x,N)'
magY=abs(X(1:1:N/2+1));
subplot(2,2,3);stem(n,x,'.');title('signal x(n)');grid;
subplot(2,2,4);plot(w/pi,magY);title('FFT N=128');
xlabel('频率(单位：pi)');ylabel('|X|');grid