数字信号处理实验二MATLAB运行代码

clear all; close all; x1=[1,1,1,1]; x2=[1,2,3,4,4,3,2,1]; x3=[4,3,2,1,1,2,3,4]; x6=x6(16,32); n1=8; n=0:n1-1; x4=cos(0.25*pi*n); n2=16; n=0:n2-1; x5=sin(0.125*pi*n); figure(1); [x1a,k1,d1]=dft1(x1,8); [x1b,k2,d2]=dft1(x1,64); subplot(3,1,1); stem([0:3],x1,'.') ylabel('x1(n)'); subplot(3,1,2); stem(k1*d1,abs(x1a),'.'); xlabel('w'); ylabel('|x1a|(N=8)'); subplot(3,1,3); stem(k2*d2,abs(x1b),'.'); xlabel('w'); ylabel('|x1b|(N=64)'); figure(2); [x2a,k1,d1]=dft1(x2,8); [x2b,k2,d2]=dft1(x2,64); subplot(3,1,1); stem([0:7],x2,'.') ylabel('x2(n)'); subplot(3,1,2); stem(k1*d1,abs(x2a),'.'); xlabel('w'); ylabel('|x2a|(N=8)'); subplot(3,1,3); stem(k2*d2,abs(x2b),'.'); xlabel('w'); ylabel('|x2b|(N=16)'); figure(3); [x3a,k1,d1]=dft1(x3,8); [x3b,k2,d2]=dft1(x3,64); subplot(3,1,1); stem([0:7],x3,'.') ylabel('x3(n)'); subplot(3,1,2); stem(k1*d1,abs(x3a),'.'); xlabel('w'); ylabel('|x3a|(N=8)'); subplot(3,1,3); stem(k2*d2,abs(x3b),'.'); xlabel('w'); ylabel('|x3b|(N=16)'); figure(4); [x4,k1,d1]=dft1(x4,8); subplot(3,1,1); stem([0:7],x4,'.') ylabel('x4(n)'); subplot(3,1,2); stem(k1*d1,abs(x4),'.'); xlabel('w'); ylabel('|x4|(N=8)'); figure(5); [x5,k1,d1]=dft1(x5,16); subplot(3,1,1); stem([0:15],x5,'.') ylabel('x5(n)'); subplot(3,1,2); stem(k1*d1,abs(x5),'.'); xlabel('w'); ylabel('|x5|(N=16)'); figure(6); [x6,k1,d1]=dft1(x6,16); subplot(3,1,1); stem([0:15],x6,'.') ylabel('x6(n)'); subplot(3,1,2); stem(k1*d1,abs(x6),'.'); xlabel('w'); ylabel('|x6|(N=16)'); n3=8; n=0:n3-1; x4=cos(0.25*pi*n); x5=sin(0.125*pi*n); x7=x4+x5; x8=x4+x5*i; figure(7); [x7a,k1,d1]=dft1(x7,8); [x7b,k2,d2]=dft1(x7,16); subplot(3,1,1); stem([0:7],x7,'.') ylabel('x7(n)'); subplot(3,1,2); stem(k1*d1,abs(x7a),'.'); xlabel('w'); ylabel('|x7a|(N=8)'); subplot(3,1,3); stem(k2*d2,abs(x7b),'.'); xlabel('w'); ylabel('|x7b|(N=16)'); figure(8); [x8a,k1,d1]=dft1(x8,8); [x8b,k2,d2]=dft1(x8,16); subplot(4,1,1); stem([0:7],real(x8),'.') ylabel('x8(n)实部'); subplot(4,1,2); stem([0:7],imag(x8),'.') ylabel('x8(n)虚部'); subplot(4,1,3); stem(k1*d1,abs(x8a),'.'); xlabel('w'); ylabel('|x8a|(N=8)'); subplot(4,1,4); stem(k2*d2,abs(x8b),'.'); xlabel('w'); ylabel('|x8b|(N=16)');