数字信号处理MATLAB实验50题.doc

数字信号处理MATLAB实验50题————————————————————————————————作者：————————————————————————————————日期：个人收集整理勿做商业用途个人收集整理勿做商业用途个人收集整理勿做商业用途1—1clc；b=[1，1］;a=[1，—0.5］;subplot（3，1,1);zplane（b,a）;title（'因果系统零极图’）;n=0:50；x=3＊cos（pi*n/3)；y=filter(b,a,x）；subplot（3，1，2);stem(n,x,’.');title(’输入x的波形’）；subplot（3，1,3)；stem(n，y,’.'）;title（'输出y的波形’)；1-2b=[1，1，1]；a=［1,0。5,-0。25］;subplot(3,1,1）；zplane(b，a）;title（'因果系统零极图’）；n=0：50；x=3*cos(pi＊n/3);y=filter(b,a，x)；subplot(3,1，2);stem（n,x，’。’）；title('输入x的波形')；subplot（3，1,3)；stem（n,y,'。’);title('输出y的波形’)；2clear；clc；b=[0,1］;a=［1，-1,—1];x=impseq（0，-5，50）;n=-5：50；h=filter(b，a，x）;stem（n，h,'。');title('单位脉冲响应’）sum(abs(h）)3b=［2];a=[1—0.8-0。5］；subplot(4，1，1)；zplane(b,a)；title('系统的零极图'）；［H,w］=freqz(b，a,100,’whole’）;magH=abs（H);phaH=angle（H)；subplot（4,1，2）；plot(w/pi，magH);title（’系统的幅频响应’)；subplot(4，1，3)；plot(w/pi,phaH/pi）;title(’系统的相频响应');n=0：100；x=impseq（0，0,100）;h=filter(b,a,x）；subplot(4，1，4)；stem（n,h,'.'）;title（'系统的冲激响应’）;4b=[11]；a=［1—0。90。81］；[H，w]=freqz（b,a,400，'whole’）；magH=abs(H)；phaH=angle（H);subplot(4,1,1）；plot(w/pi,magH）；title（'系统的幅频响应');subplot(4，1，2);plot（w/pi，phaH/pi)；title(’系统的相频响应’）；n=0:200;x=sin(pi＊n/3)+5＊cos(pi*n);y=filter(b,a,x）;subplot(4，1,3）；plot(n,x）；title(’输入信号X’);subplot（4，1,4);plot(n，y）；title（’输出信号Y’）;grid;5x11=［1111];n=0:5;x12=cos(pi＊n/4);y11=circonvt（x11,x12,8)y12=conv(x11，x12)y13=［y11（1：1：8),zeros（1,1）]e1=y13—y12x21=[1-11—1]；x22=[10—10］;y21=circonvt(x21,x22，5)y22=conv（x21，x22）y23=［y21（1:1:5)，zeros（1，2）］e2=y23-y22n=0:15；x31=cos(2＊pi*n/32)；x32=sin（2＊pi＊n/32);y31=circonvt(x31,x32，32）y32=conv(x31，x32)y33=［y32（1：1:31），zeros（1,1）]e3=y31-y33n=0:9；x41=（0.8)。^n；x42=（—0。8).^n；y41=circonvt（x41,x42，15）y42=conv（x41,x42)y43=[y41（1：1：15),zeros(1，4)］e4=y43-y426clear;clc;x1=［2112］；x2=[1-1—11]；n=[0:8-1];y11=circonvt（x1,x2，4）y12=circonvt(x1，x2，7）y13=circonvt（x1，x2,8）y2=conv(x1,x2)%N最小值77——1x=[2,2，2,2，2，2，2，2］；w=［0:1：500］*2*pi/500;［H]=freqz（x，1,w）;magH=abs(H);phaH=angle(H);subplot（2，2，1)；plot(w/pi,magH）;gridxlabel（'')；ylabel（'|x｜'）；title(’