数字信号处理实验(独家).docVIP

clc; clear all; x1=[ones(1,8),zeros(1,50)]; x2=ones(1,20); B=[0.05,0.05]; A=[1,-0.9]; y1=filter(B,A,x1); y2=filter(B,A,x2); y3=impz(B,A,40); subplot(221); stem(y1); xlabel(n); ylabel(response(y1)); title(R8(n)的滤波器系统响应y1); subplot(222); stem(y2); xlabel(n); ylabel(response(y2)); title(u(n)的滤波器系统响应y2); subplot(223); stem(y3); xlabel(n); ylabel(response(y3)); title(滤波器系统的单位脉冲响应y3); h1=ones(1,10); x4=[ones(1,8),zeros(1,8)]; y4=conv(h1,x4); subplot(224); stem(y4); xlabel(n); ylabel(response(y4)); title(R8(n)与系统h1的卷积y4); figure(2); h2=[1,2.5,2.5,1]; y5=conv(h2,x4); subplot(221); stem(y5); xlabel(n); ylabel(response(y5)); title(R8(n)与系统h2的卷积y5); A=[1,-1.8237,0.9801]; B=[1/100.49,-1/100.49]; x5=ones(1,150); y6=filter(B,A,x5); subplot(222); stem(y6); xlabel(n); ylabel(response(y6)); title(u(n)的谐振器系统响应y6); n=0:90; x6=sin(0.014*n)+sin(0.4*n); y7=filter(B,A,x6); subplot(223); stem(y7); xlabel(n); ylabel(response(y7)); title(x(n)的谐振器系统响应y7); clc; clear all; x1n=[ones(1,4)]; M=8;xa=1:(M/2); xb=(M/2):-1:1; x2n=[xa,xb];x3n=[xb,xa]; X1k8=fft(x1n,8);X1k16=fft(x1n,16); X2k8=fft(x2n,8); X2k16=fft(x2n,16); X3k8=fft(x3n,8);X3k16=fft(x3n,16); subplot(3,2,1);stem(abs(X1k8)); title(x1(n)的8点DFT);xlabel(ω/π);ylabel(幅度); axis([0,2,0,1.2*max(abs(X1k8))]) subplot(3,2,2);stem(abs(X1k16));box on title(x1(n)的16点DFT);xlabel(ω/π);ylabel(幅度); axis([0,2,0,1.2*max(abs(X1k16))]) subplot(3,2,3);stem(abs(X2k8)); box on title(x2(n)的8点DFT);xlabel(ω/π);ylabel(幅度); axis([0,2,0,1.2*max(abs(X2k8))]) subplot(3,2,4);stem(abs(X2k16));box on title(x2(n)的16点DFT);xlabel(ω/π);ylabel(幅度); axis([0,2,0,1.2*max(abs(X2k16))]) subplot(3,2,5);stem(abs(X3k8));box on title(x3(n)的8点DFT);xlabel(ω/π);ylabel(幅度); axis([0,2,0,1.2*max(abs(X3k8))]) subplot(3,2,6);stem(abs(X3k16)); box on title(x3(n)的16点DFT);xlabel(ω/π);ylabel(幅度); axis([0,2,0,1.2*max(abs(X3k16))]) N=8;n=0:N-1; x4n=cos(pi*n/4); x5n=cos(pi*n/4)+cos(pi*n/8); X4k8=fft(x4n); X5k8=fft(x5n); N=16;n=0:N-1; x4n=cos(pi*n/4); x5n=cos(pi*n/4)+cos(pi*n/8);