函数零极点发布判断系统的稳定性响应的方法.doc

正文课 程 设 计 内 容 编制Matlab程序，完成以下功能，根据系统函数求出系统的零极点分布图并求解系统的单位脉冲响应；根据零极点分布图判断系统的稳定性；比较不同零极点发布对系统频率响应特性的影响；绘制相关信号的波形。具体要求如下： 下面四种二阶网络的系统函数具有相同的极点发布： 分别画出各系统的零极点分布图，并判断系统的稳定性； 分别画出系统的幅频特性和相频特性曲线； 分别求出系统的单位脉冲响应，并画出其波形。分别画出各系统的零极点分布图，并判断系统的稳定性； MATLAB源程序： B=[1,0,0]; A=[1,-1.6,0.9425]; %设置系统函数系数向量B和A zplane(B,A),grid on %绘制零极点图 legend(零点,极点) title(零极点分布图) 图（1）系统零极点分布图 程序运行结果如图所示。由图可知，该因果系统的极点全部在单位圆内，故系统是稳定的。 由程序判断稳定性，程序如下： function stab(A) A=[1,-1.6,0.9425] %设置系统函数系数向量A disp(‘系统极点为：’) p=roots(A) %求极点，并显示 disp(‘系统极点模的最大值’) M=max(aps(p)) %求模的最大值并显示 if M1 disp(‘系统稳定’),else,disp(‘系统不稳定’),end %判断系统稳定性 图（2）程序运行结果 MATLAB源程序： B=[1,-0.3,0]; A=[1,-1.6,0.9425]; %设置系统函数系数向量B和A zplane(B,A),grid on %绘制零极点图 legend(零点,极点) title(零极点分布图) 图（3）零极点分布图 程序运行结果如图所示。由图可知，该因果系统的极点全部在单位圆内，故系统是稳定的。 由程序判断稳定性，程序如下： function stab(A) A=[1,-1.6,0.9425] %设置系统函数系数向量A disp(‘系统极点为：’) p=roots(A) %求极点，并显示 disp(‘系统极点模的最大值’) M=max(aps(p)) %求模的最大值并显示 if M1 disp(‘系统稳定’),else,disp(‘系统不稳定’),end %判断系统稳定性 图（4）程序运行结果 MATLAB源程序： B=[1,-0.8,0]; A=[1,-1.6,0.9425]; %设置系统函数系数向量B和A zplane(B,A),grid on %绘制零极点图 legend(零点,极点) title(零极点分布图) 图（5）零极点分布图 程序运行结果如图所示。由图可知，该因果系统的极点全部在单位圆内，故系统是稳定的。 由程序判断稳定性，程序如下： function stab(A) A=[1,-1.6,0.9425] %设置系统函数系数向量A disp(‘系统极点为：’) p=roots(A) %求极点，并显示 disp(‘系统极点模的最大值’) M=max(aps(p)) %求出模的最大值并显示 if M1 disp(‘系统稳定’),else,disp(‘系统不稳定’),end %判断系统稳定性 图（6）程序运行结果 MATLAB源程序： B=[1,-1.6,0.8]; A=[1,-1.6,0.9425]; %设置系统函数系数向量B和A zplane(B,A),grid on %绘制零极点图 legend(零点,极点) title(零极点分布图) 图（7）零极点分布图 程序运行结果如图所示。可见，该因果系统的极点全部在单位圆内，