零极点对系统性能的影响分析1任务步骤1.分析原开环传递函数G0（s）的性能，绘制系统的阶跃响应曲线得到系统的暂态性能（包括上升时间，超调时间，超调量，调节时间）；2.在G0（s）上增加零点，使开环传递函数为G1（s），绘制系统的根轨迹,分析系统的稳定性；3.取不同的开环传递函数G1（s）零点的值，绘制系统的阶跃响应曲线得到系统的暂态性能（包括上升时间，超调时间，超调量，调节时间）；4.综合数据，分析零点对系统性能的影响5.在G0（s）上增加极点，使开环传递函数为G2（s），绘制系统的根轨迹，分析系统的稳定性；6.取不同的开环传递函数G2（s）极点的值，绘制系统的阶跃响应曲线得到系统的暂态性能（包括上升时间，超调时间，超调量，调节时间）；7.综合数据，分析极点对系统性能的影响。

8.增加一对离原点近的偶极子和一对距离原点远的偶极子来验证偶极子对消的规律。

2原开环传递函数G0（s）的性能分析2.1 G0（s）的根轨迹取原开环传递函数为：Matlab指令：num=[1];den=[1,0.8,0.15];rlocus(num,den);得到图形：图1 原函数G0(s)的根轨迹根据原函数的根轨迹可得：系统的两个极点分别是-0.5和-0.3，分离点为-0.4，零点在无限远处，系统是稳定的。

2.2 G0（s）的阶跃响应Matlab指令：G=zpk([],[-0.3,-0.5],[1])sys=feedback(G,1)step(sys)得到图形：图2 原函数的阶跃响应曲线由阶跃响应曲线分析系统暂态性能： 曲线最大峰值为1.12，稳态值为0.87， 上升时间tr=1.97s 超调时间tp=3.15s 调节时间ts=9.95s ，2=∆超调量%p σ=28.3%3 增加零点后的开环传递函数G1（s）的性能分析为了分析开环传递函数的零点对系统性能的影响，现在在原开环传递函数的表达式上单独增加一个零点S=-a,并改变a值大小，即离虚轴的距离，分析比较系统性能的变化。

所以增加零点后的开环传递函数为：开环传递函数表达式：3.1 G1（s）的根轨迹因为后面利用阶跃响应来分析时将取的零点均在实轴的负半轴，那么只要了解其中一个开环传递函数稳定，那么其它的稳定也可以推知。

所以取a=1画出根轨迹来观察系统的稳定性。

当a=1时，开环传递函数的表达式为：Matlab指令：num=[1,1];den=[1,0.8,0.15];rlocus(num, den)得到图图3 G1（s ）的根轨迹曲线根据G1（s ）的根轨迹可得：根轨迹均在左半平面，只是多了一个零点，系统仍然是稳定的，并且可以推知，只要零点在实轴的负半轴上，系统都是稳定的。

3.2 增加不同零点时G1（s ）的阶跃响应3.2.1 当a=0.01的阶跃响应当a=0.01时，对应的闭环传递函数为Matlab 指令： num=[100,1]; den=[1,100.8,1.15]; step(num,den) grid on 得到图图41()s 的阶跃响应曲线由阶跃响应曲线分析系统暂态性能： 曲线最大峰值为0.992，稳态值为0.87， 上升时间tr=0.0434s超调时间tp=0.139s 调节时间ts=197s ，2=∆超调量%p σ=11.4%3.2.2 当a=0.1的阶跃响应当a=0.1时，对应的闭环传递函数为Matlab 指令： num=[10,1]; den=[1,10.8,1.15]; step(num,den) grid on 得到图图52()s φ的阶跃响应曲线由阶跃响应曲线分析系统暂态性能：由图可知，曲线最大峰值为0.931，稳态值为0.87， 上升时间tr=0.256s 超调时间tp=0.685s 调节时间ts=12.4s ，2=∆超调量%p σ=7.02%3.2.3当a=1的阶跃响应当a=1时，对应的闭环传递函数为Matlab 指令 num=[1,1]; den=[1,1.8,1.15]; step(num,den) grid on 得到图图63()s φ的阶跃响应曲线由阶跃响应曲线分析系统暂态性能：由图可知，曲线最大峰值为0.905，稳态值为0.87， 上升时间tr=2.04s 超调时间tp=2.97s 调节时间ts=4.43s ，2=∆超调量%p σ=4.03%3.2.4当a=10的阶跃响应当a=10时，对应的闭环传递函数为Matlab 指令 num=[0.1,1]; den=[1,0.9,1.15]; step(num,den) grid on 得到图图74()s φ的阶跃响应曲线由阶跃响应曲线分析系统暂态性能：由图可知，曲线最大峰值为1.07，稳态值为0.87， 上升时间tr=1.98s 超调时间tp=3.15s 调节时间ts=7.73s ，2=∆超调量%p σ=23.5%3.2.5当a=100的阶跃响应当a=10时，对应的闭环传递函数为Matlab 指令 num=[0.01,1]; den=[1,0.81,1.15]; step(num,den) grid on 得到图图85()s φ的阶跃响应曲线由阶跃响应曲线分析系统暂态性能：由图可知，曲线最大峰值为1.11，稳态值为0.87， 上升时间tr=1.96s 超调时间tp=3.11s 调节时间ts=9.84s ，2=∆超调量%p σ=27.7%3.3增加零点后对系统性能的影响分析根据图2，图4，图5，图6，图7，图8，可以得到原函数以及在原开环传递函数上增加一个零点s=a ，a 分别取0.01,0.1,1,10,100的系统性能参数。

如以下表1所示：表1根据表1可画出lga 与各个指标的关系曲线，如以下图9，图10，图11，图12和图13。

因为原函数中的lga 的值为负无穷，所以无法在图中直接反映，所以图9，图10，图11，图12和图13五个图反映的是，零点距离原点的远近对系统性能的影响。

图9 曲线峰值Mr 与lg （a ）的关系图10 上升时间tr 与lg （a ）的关系图11 超调时间与lg（a）的关系图12调节时间与lg（a）的关系图13 超调量与lg（a）的关系结论：1.增加不同的零点对系统参数有不同的影响；2.曲线峰值与超调量受到影响后的值与原值没有重合，上升时间，超调时间与调节时间与原值有重合；3.随着a的增加（或者说随着零点渐渐远离零点），曲线峰值受到的影响（取绝对值来看）和超调量受到的影响均是先增后减；上升时间受到的影响，超调时间受到的影响，调节时间受到的影响均是先减后增再减；4.当a=100时，也就是零点距离原点最远时，增加的零点对系统的影响最小，可以预见，当零点与原点的距离趋近于无穷远时，系统性能受到的影响趋近于0。

4 增加极点后的开环传递函数G2（s）的性能分析为了分析开环传递函数的极点对系统性能的影响，现在在原开环传递函数的表达式上单独增加一个极点S=-p,并改变p值大小，即离原点的距离，分析比较系统性能的变化。

所以增加零点后的开环传递函数为：4.1 G2(s)的根轨迹因为后面利用阶跃响应来分析时将取的极点均在实轴的负半轴，那么只要了解其中一个开环传递函数稳定，那么其它的稳定也可以推知。

所以取p=1画出根轨迹来观察系统的稳定性。

当p=1时，开环传递函数G2(s)的表达式为Matlab指令：num=[1];den=[1,1.8,0.95,0.15];rlocus(num,den);h = findobj(gcf, 'Type','line');set(h, 'LineWidth', 3);得到图：图14 原函数G0(s)的根轨迹根据G(s)的根轨迹可得：根轨迹均在左半平面，只是多了一个极点，系统仍然是稳定的，并且可以推知，只要极点在实轴的负半轴上，系统都是稳定的。

4.2增加不同极点时G2（s）的阶跃响应4.2.1 当p=0.01的阶跃响应当p=0.01时，对应的闭环传递函数为：Matlab指令：num=[1];den=[100,81,15.8,1.15];step(num,den);h = findobj(gcf, 'Type','line');set(h, 'LineWidth', 3);得到图：图151()s φ的阶跃响应曲线由阶跃响应曲线分析系统暂态性能： 曲线最大峰值为0.875，稳态值为0.87， 上升时间tr=37.1s 超调时间tp=44.5s 调节时间ts=31.7s ，2=∆超调量%p σ=0.569%4.2.2 当p=0.1的阶跃响应当p=0.1时，对应的闭环传递函数为num=[1];den=[10,9,2.3,1.15]; step(num,den);h = findobj(gcf, 'Type','line'); set(h, 'LineWidth', 3); 得到图：图162()s φ的阶跃响应曲线由阶跃响应曲线分析系统暂态性能： 曲线最大峰值为1.37，稳态值为0.87， 上升时间tr=5.84s 超调时间tp=9.58s 调节时间ts=69.7s ，2=∆超调量%p σ=57.2%4.2.3 当p=1的阶跃响应当p=1时，对应的闭环传递函数为Matlab 指令： num=[1];den=[1,1.8,0.95,1.15]; step(num,den);h = findobj(gcf, 'Type','line'); set(h, 'LineWidth', 3);图173()s φ的阶跃响应曲线由阶跃响应曲线分析系统暂态性能： 曲线最大峰值为1.45，稳态值为0.87， 上升时间tr=2.59s 超调时间tp=4.38s 调节时间ts=50s ，2=∆超调量%p σ=66.4%4.2.4 当p=10的阶跃响应当p=10时，对应的闭环传递函数为Matlab 指令： num=[1];den=[0.1,1.08,0.815,1.15]; step(num,den);h = findobj(gcf, 'Type','line'); set(h, 'LineWidth', 3);图184()s φ的阶跃响应曲线由阶跃响应曲线分析系统暂态性能： 曲线最大峰值为1.16，稳态值为0.87， 上升时间tr=1.97s 超调时间tp=3.18s 调节时间ts=10.5s ，2=∆超调量%p σ=33.7%4.2.5 当p=100的阶跃响应当p=100时，对应的闭环传递函数为Matlab 指令： num=[1];den=[0.01,1.008,0.8015,1.15]; step(num,den);h = findobj(gcf, 'Type','line'); set(h, 'LineWidth', 3);图195()s φ的阶跃响应曲线由阶跃响应曲线分析系统暂态性能： 曲线最大峰值为1.12，稳态值为0.87， 上升时间tr=1.95s 超调时间tp=3.19s 调节时间ts=10s ，2=∆超调量%p σ=28.8%4.3增加极点后对系统性能的影响分析根据图2，图15，图16，图17，图18，图19，可以得到原函数以及在原开环传递函数上增加一个零点s=-p ，p 分别取0.01,0.1,1,10,100的系统性能参数。