b=1;c1=095;c2=05;a1=1-c1;a2=1-c2;y前系数的值H1w1=freqzba11000whole;H2w2=freqzba21000whole;subplot211plotw1absH1;xlabel频率弧度;ylabel幅度;title 系统的幅频响应;grid;subplot212plotw2absH2;xlabel频率弧度;ylabel幅度;grid;在以上代码基础

n=0:99; x=n.*(n>=0); c1=0.95; c2=0.5; a1=[1,-c1]; a2=[1,-c2]; b=1; [H1,w1]=freqz(b,a1,1000,'whole'); [H2,w2]=freqz(b,a2,1000,'whole'); X=fft(x,1000); X_mag=abs(X); w=linspace(-pi,pi,1000); figure; subplot(2,1,1) plot(w,X_mag); xlabel('频率（弧度）');ylabel('幅度'); title ('信号频域特性'); grid; subplot(2,1,2) plot(w,abs(H1).*X_mag); hold on; plot(w,abs(H2).*X_mag); xlabel('频率（弧度）');ylabel('幅度'); legend('c=0.95','c=0.5'); title ('系统处理后的频域特性'); grid;