sy12.m

数字信号处理工具箱

N=64;f1=1/16;
n=[0:1:N-1];
xn=sin(2*pi*0.15*n)+cos(2*pi*(0.15+f1)*n)+cos(2*pi*(0.15+2*f1)*n)
subplot(2,2,1); stem(n,xn);
title('x(n)');xlabel('n');
axis([0,N,-2.5,5.5]); set(gca,'fontsize',10)
X=fft(xn,N)'
k=[0:1:N/2];w=2*pi/N*k;
magX=abs(X(1:1:N/2+1));
subplot(2,2,2);plot(w/pi,magX);title('FFT Magnitude');
xlabel('频率(单位:pi)')
N=128;f1=1/16;
n=[0:1:N-1];
xn=sin(2*pi*0.15*n)+cos(2*pi*(0.15+f1)*n)+cos(2*pi*(0.15+2*f1)*n);
subplot(2,2,3); stem(n,xn);
title('x(n), 0 <= n <= 15');xlabel('n');
axis([0,N,-2.5,5.5]); set(gca,'fontsize',10)
X1=fft(xn,N)'
k1=[0:1:N/2];w1=2*pi/N*k1;
magX1=abs(X1(1:1:N/2+1));
subplot(2,2,4);plot(w1/pi,magX1);title('FFT Magnitude');
xlabel('频率(单位:pi)')