数字信号处理课程设计
a=[0.95,0.5];A=[1,-a(1)];
B=1;[H,w]=freqz(B,A,'whole');
subplot(1,2,1);plot(w/pi,abs(H));
xlabel('\omega/\pi');
ylabel('|H(e^j^\omega)|');
title('a=0.95 幅频响应特性');
axis([0,2,0,2.5]);grid on;
A=[1,-a(2)];
B=1;[H,w]=freqz(B,A,'whole');
subplot(1,2,2);plot(w/pi,abs(H));
xlabel('\omega/\pi');
ylabel('|H(e^j^\omega)|');
title('a=0.5 幅频响应特性');
axis([0,2,0,2.5]);grid on;