1-1 clc; 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-2 b=[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 的波形'); 2 clear; 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)) 3 b=[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('系统的冲激响应'); 4 b=[1 1]; a=[1 -0.9 0.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; 5 x11=[1 1 1 1]; 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-y12 x21=[1 -1 1 -1]; x22=[1 0 -1 0]; y21=circonvt(x21,x22,5) y22=conv(x21,x22) y23=[y21(1:1:5),zeros(1,2)] e2=y23-y22 n=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-y33 n=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-y42 6 clear; clc; x1=[2 1 1 2]; x2=[1 -1 -1 1]; n=[0:8-1]; y11=circonvt(x1,x2,4) y12=...
