在z变换中，对于系统响应函数H(z) 绘制频率响应曲线

注意点：

绘制零极点图的函数 zplane(B,A);%%B为分母的系数矩阵  A为分子的系数矩阵 第一个数为z^0 第二个数为z^-1的系数 以此类推

freqz(B,A) %%计算频率响应H(e^jw)

angle(H) %%  绘制相频响应

实现代码如下：

%%已知H(z) = 1-z^(-N)
%% N=8时 绘制频率响应曲线
B=[1 0 0 0 0 0 0 0 -1];A=-1;
subplot(2,2,1);
zplane(B,A);%%绘制零极点图
[H w] = freqz(B,A);%%计算频率响应H(ejw)
subplot(2,2,2);plot(w/pi,abs(H));%绘制频率响应曲线
xlabel('\omega^pi');ylabel('|H(e^j^w)|');axis([0 1 0 2.5])
subplot(2,2,4);
plot(w/pi,angle(H));%绘制相频响应
xlabel('\omega^pi');ylabel('phi(\omega)');

绘制图像如下： 

判断系统稳定的程序：

%%判断系统是否稳定
A=[1 3/4  1/8];
P=roots(A);
M=max(abs(P));
if(M<1) disp('系统稳定');else disp('系统不稳定');
end