零极点对系统的性能影响分析
零极点对系统性能的影响分析
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 得到图
图4
1()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 得到图
图5
2()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 得到图
图6
3()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 得到图
图7
4()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 得到图
图8
5()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);
得到图：
图15
1()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); 得到图：
图16
2()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);
图17
3()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);
图18
4()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);
图19
5()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的系统性能参数。如以下表2所示：
表 2
根据表1可画出lgp 与各个指标的关系曲线，如以下图20，图21，图22，图23和图24。因为原函数中的lga 的值为负无穷，所以无法在图中直接反映，
所以图20，图21，图22，图23和图24五个图反映的是，极点距离原点的远近对系统性能的影响。
图20 曲线峰值Mr与lg（p）的关系
图21 上升时间tr与lg（p）的关系
系统稳定性意义以及稳定性的几种定义.
系统稳定性意义以及稳定性的几种定义 一、引言： 研究系统的稳定性之前，我们首先要对系统的概念有初步的认识。 在数字信号处理的理论中，人们把能加工、变换数字信号的实体称作系统。由于处理数字信号的系统是在指定的时刻或时序对信号进行加工运算，所以这种系统被看作是离散时间的，也可以用基于时间的语言、表格、公式、波形等四种方法来描述。从抽象的意义来说，系统和信号都可以看作是序列。但是，系统是加工信号的机构，这点与信号是不同的。人们研究系统还要设计系统，利用系统加工信号、服务人类，系统还需要其它方法进一步描述。描述系统的方法还有符号、单位脉冲响应、差分方程和图形。 电路系统的稳定性是电路系统的一个重要问题，稳定是控制系统提出的基本要求，也保证电路工作的基本条件；不稳定系统不具备调节能力，也不能正常工作，稳定性是系统自身性之一，系统是否稳定与激励信号的情况无关。对于线性系统来说可以用几点分布来判断，也可以用劳斯稳定性判据分析。对于非线性系统的分析则比较复杂，劳斯稳定性判据和奈奎斯特稳定性判据受到一定的局限性。 二、稳定性定义： 1、是指系统受到扰动作用偏离平衡状态后，当扰动消失，系统经过自身调节能否以一定的准确度恢复到原平衡状态的性能。若当扰动消失后，系统能逐渐恢复到原来的平衡状态，则称系统是稳定的，否则称系统为不稳定。 稳定性又分为绝对稳定性和相对稳定性。 绝对稳定性。如果控制系统没有受到任何扰动，同时也没有输入信号的作用，系统的输出量保持在某一状态上，则控制系统处于平衡状态。 （1）如果线性系统在初始条件的作用下，其输出量最终返回它的平衡状态，那么这种系统是稳定的。 （2）如果线性系统的输出量呈现持续不断的等幅振荡过程，则称其为临界稳定。（临界稳定状态按李雅普洛夫的定义属于稳定的状态，但由于系统参数变化等原因，实际上等幅振荡不能维持，系统总会由于某些因素导致不稳定。因此从工程应用的角度来看，临界稳定属于不稳定系统，或称工程意义上的不稳定。） （3）如果系统在初始条件作用下，其输出量无限制地偏离其平衡状态，这称系统是不稳定的。 实际上，物理系统的输出量只能增大到一定范围，此后或者受到机械制动装置的限制，或者系统遭到破坏，也可以当输出量超过一定数值后，系统变成非线性的，从而使线性微分方程不再适用。因此，绝对稳定性是系统能够正常工作的前提。
零极点对系统的性能影响分析
零极点对系统性能的影响分析 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%
极点及系统稳定性
极点对系统性能影响 一．控制系统与极点 自动控制系统根据控制作用可分为：连续控制系统和采样控制系统，采样系统又叫离散控制系统。通常把系统中的离散信号是脉冲序列形成的离散系统,称为采样控制系统。连续控制系统即指控制量为连续的模拟量如时变系统。 系统的数学模型一般由系统传递函数表达。传递函数为零初始条件下线性系统响应（即输出）量的拉普拉斯变换（或z 变换）与激励（即输入）量的拉普拉斯变换之比。记作Φ（s ）=Xo （s ）/Xi （s ），其中Xo （s ）、Xi （s ）分别为输出量和输入量的拉普拉斯变换。 特征方程的根称为极点。如试Φ﹙S ﹚= C [∏（S-Pi ）/∏(S-Qi) ]中Q1 Q2 Q3 …… Qi ……即为系统的极点。 二．极点对系统的影响 极点--确定了系统的运动模态；决定了系统的稳定性。下面对连续系统与离散系统分别进行分析： ⑴连续系统 理论分析：连续系统的零极点分布有如下几种形式 设系统函数为： 将H(S)进行部分分式展开： 1n a s -+++
系统冲激响应H(S)的时域特性h(t)随时间衰减的信号分量完全由系统函数H(S)的极点位置决定。每一个极点将决定h(t)的一项时间函数。 稳定性：由上述得知Y(S)= C [∏（S-Pi ）/(S-Qi) ]可分解为Y(S)=C1/(S-τ1)+ C2/(S-τ2)+ C3/(S-τ3)+……+ Ci/(S-τi)+…… 则时间响应为 …… 由于特征方程的根不止一个，这时，应把系统的运动看成是多个运动分量的合成。只要有一个运动分量是发散的，则系统是不稳定的。因此，特征方程所有根的实部都必须是负数，亦即所有的根都在复平面的左半平面。 通过复变函数幅角定理将S 由G 平面映射到GH 平面。 如果封闭曲线 F 内有Z 个F(s)的零点，有P 个F(s)的极点，则s 沿 F 顺时针转一圈时，在F(s)平面上，F(s)曲线绕原点顺时针转的圈数R 为z 和p 之差，即R ＝z －p 。 若R 为负,表示F(s)曲线绕原点逆时针转过的圈数。 F(s)的分母是G0(s)的分母，其极点是G0(s)的极点；其分子是?(s)的分母，即?(s)的特征多项式，其零点是?(s)的极点。 取D 形曲线（D 围线）如图所示，是整个右半复平面。 且设D 曲线不经过F(s)的任一极点或零点。 s 沿D 曲线顺时针变化一周，F(s)顺时针包围原点的周数为： n=z-p=F(s)在右半复平面的零点数(闭环传函在右半复平面极点数) -F(s)在右半复平面的极点数(开环传函在右半复平面极点数) 所以闭环系统稳定的充分必要条件是： n=- p =-开环传函在右半复平面的极点数 1212()n s t s t s t n y t C e C e C e =+++0()0()0()0()t s y t y t Ce y t y t t ααααα=<→?? ===??>→∞? →∞（1）只有一个实根：时，时，恒量时，()()121()0cos()00j t j t t s j y t C e C e C e t t αωαωααωαω?αα+-=±=+? →∞（2）有一对复根：时，收敛时，等幅振荡时，发散
零极点分布对系统频率响应的影响
备注：(1)、按照要求独立完成实验内容。 (2)、实验结束后，把电子版实验报告按 要求格式改名(例：09 号_张三 _实验七.doc)后，实验室统一刻 盘留档。 实验三零极点分布对系统频 率响应的影响 一、实验目的 1. 掌握系统差分方程得到系统函数的方法； 2. 掌握系统单位脉冲响应获取系统函数的方法； 3. 掌握用系统函数零级点分布的几何方法分析研究系统的频率响应 二、实验原理 在MA TLAB 中，可以用函数[z，p，K]=tf2zp ( num ，den)求得有理分式形式的系统转移函数的零、极点，用函数zplane( z，p)绘出 零、极点分布图；也可以用函数 zplane( num，den)直接绘出有理分式形式的系统转移函数的零、极点分布图。 另外，在MA TLAB 中，可以用函数[r，p，k]=residuez(num，den)完成部分分式展开计算；可以用函数sos=zp2sos( z，p，K )完成三、实验内容(包括代码与产生的图形) 1. 假设系统用下面差分方程描述： y(n)=x(n)+ay(n-1) 假设a=0.7, 0.8, 0.9 ，分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。 B=1; A=[1,-0.7]; subplot(3,3,1);zplane(B,A); xlabel(' 实部Re'); ylabel(' 虚部Im'); title('y(n)=x(n)+0.7y(n-1) 传输函数零、极点分布'); grid on [H,w]=freqz(B,A,'whole'); subplot(3,3,4); 将高阶系统分解为 2 阶系统的串联。plot(w/pi,abs(H),'linewidth',2);
绘制离散系统零极点图.
绘制离散系统零极点图：zplane() 滤波器 绘制离散系统零极点图：zplane() zplane(Z,P) 以单位圆为基准绘制零极点图，在图中以'o'表示零点，以'x'表示极点，如果存在重零极点，则在它们的右上方显示其数目。如果零极点是用矩阵来表示，在不同行内的零极点用不同的颜 色来表示。 zplane(B, A) 输入的是传递函数模型，则函数将首先调用root 函数以求出它们的零极点。 [H1, H2, H3]=zplane(Z,P) 函数返回图形对象的句柄。其中，H1返回的是零点线的句柄；H2返回的是极点线的句柄；H3返回的是轴和单位圆线条句柄。如果有重零极点，它还包括显示在其右上方 的文本句柄。 例：设计一个数字椭圆带阻滤波器，具体要求是：通带截止频率是 wp1=1500Hz，wp2=2500Hz，阻带截止频率是ws1=1000Hz，ws2=3000Hz，在通带内的最大衰减为0.5dB，在阻带内的最小衰减 为60dB 程序设计如下： wp1=1500; wp2=2500; ws1=1000; ws2=3000; Fs=100 00Hz; rp=0.5; rs=60; wp=[wp1,wp2]; ws=[ws1,ws2]; [n,wn]=ellipord(wp/(Fs/2), ws/(Fs/2), rp, rs); [num,den]=ellip(n, rp, rs, wn, 'stop'); [H, W]=freqz(num, den); figure; plot(W*Fs/(2*pi), abs(H)); grid; xlabel('频率/Hz'); ylabel('幅值'); figure; impz(num, den); figure; grpdelay(num, den); figure; zplane(num, den); FREQZ 是计算数字滤波器的频率响应的函数
零极点分布对系统频率响应的影响
备注：（1）、按照要求独立完成实验内容。 （2）、实验结束后，把电子版实验报告 按要求格式改名（例：09号_张 三_实验七.doc）后，实验室统一 刻盘留档。 实验三零极点分布对系统频 率响应的影响 一、实验目的 1.掌握系统差分方程得到系统函数的方法； 2.掌握系统单位脉冲响应获取系统函数的方法； 3.掌握用系统函数零级点分布的几何方法分析研究系统的频率响应 二、实验原理 在MA TLAB中，可以用函数[z，p，K]=tf2zp （num，den）求得有理分式形式的系统转移函数的零、极点，用函数zplane（z，p）绘出零、极点分布图；也可以用函数zplane（num，den）直接绘出有理分式形式的系统转移函数的零、极点分布图。 另外，在MA TLAB中，可以用函数 [r，p，k]=residuez（num，den）完成部分分式展开计算；可以用函数sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系统的串联。三、实验内容（包括代码与产生的图形） 1. 假设系统用下面差分方程描述： y(n)=x(n)+ay(n-1) 假设a=0.7, 0.8, 0.9 ，分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。 B=1; A=[1,-0.7]; subplot(3,3,1);zplane(B,A); xlabel('实部Re'); ylabel('虚部Im'); title('y(n)=x(n)+0.7y(n-1)传输函数零、极点分布'); grid on [H,w]=freqz(B,A,'whole'); subplot(3,3,4); plot(w/pi,abs(H),'linewidth',2);
判断系统稳定性
摘要 现今数字信号处理理论与应用已成为一门很重要的高新科学技术学科，通过功能强大的MATLAB软件与数字信号处理理论知识相互融合在一起，既使我们对数字信号处理的理论知识能够有更加深厚的解也提高了动手能力，实践并初步掌握了MATLAB 的使用。 根据本次课题要求，通过使用MATLAB，方便了对系统函数的繁琐的计算，并且直观形象的用计算机进行模拟仿真，通过观察图，由图像的特征从而进一步的对系统进行形象的分析。 本课题中给出了系统函数，对其稳定性进行分析我们可以通过MATLAB画零极图观察极点的分布，另外还可以通过MATLAB分析系统的单位阶跃响应、单位脉冲响应、幅频相频特性的图形更加具体的对系统进行分析。 关键字：离散系统函数、MATLAB、零极点分布、系统稳定性。
一、设计原理 1.设计要求 （1）：根据系统函数求出系统的零极点分布图并且判断系统的稳定性。 （2）：求解系统的单位阶跃响应，并判断系统的稳定性。 （3）：求系统的单位脉冲响应，并判断系统的稳定性 （4）：求出各系统频率响应，画出幅频特性和相频特性图（zp2tf,zplane，impz等） 2、系统稳定性、特性分析 进行系统分析时我主要利用MATLAB软件绘制出系统零极点的分布图、单位脉冲响应图、单位阶跃响应图等。采用MATLAB 软件进行设计时我调用了软件本身的一些函数来对课题进行绘图和分析。诸如zplane、impz、stepz、freqz等。 对系统函数的零极图而言：极点在单位圆内,则该系统稳定,极点在单位圆外，则该系统为非稳定系统。 当极点处于单位圆内，系统的冲激响应曲线随着频率的增大而收敛；当极点处于单位圆上，系统的冲激响应曲线为等幅振荡；当极点处于单位圆外，系统的冲激响应曲线随着频率的增大而发散。 系统的单位阶跃响应若为有界的则系统为稳定系统。由以上的判据配合图形对系统的稳定性进行分析，达到我们的课程要求。 系统函数H（z）的零极点分布完全决定了系统的特性，若某系统函数的零极点已知，则系统函数便可确定下来。 因此，系统函数的零极点分布对离散系统特性的分析具有非常重要意义。通过对系统函数零极点的分析，可以分析离散系统以下几个方面的特性： (1)系统单位样值响应h(n)的时域特性； (2)离散系统的稳定性； (3)离散系统的频率特性；
matlab实验四 系统的零极点分析
实验四连续时间系统复频域分析和离散时间系统z域分析 一.实验目的： 1.掌握连续信号拉氏变换和拉氏反变换的基本实现方法。 2.熟悉laplace函数求拉普拉斯变换，ilaplace函数求拉氏反变换 的使用。 3.掌握用ztrans函数，iztrans函数求离散时间信号z变换和逆z 变换的基本实现方法。 4.掌握用freqs函数，freqz函数由连续时间系统和离散时间系统 系统函数求频率响应。 5.掌握zplane零极点绘图函数的使用并了解使用零极点图判断系 统稳定性的原理。 二、实验原理： 1.拉氏变换和逆变换 原函数()() ?象函数 f t F s 记作：[()]() =→拉氏变换 L f t F s 1[()]() -=→拉氏反变换 L F s f t 涉及函数：laplace,ilapace. 例如：
syms t;laplace(cos(2*t)) 结果为：ans =s/(s^2+4) syms s;ilaplace(1./(s+1)) 结果为：ans = exp(-t) 2. 系统传递函数H(s)或H(z)。 12121212...()()()...m m m n n n b s b s b B s H s A s a s a s a ----+++==+++ 112112...()()()...m m m n n n b z b z b B z H z A z a z a z a --+--++++==+++ 其中，B 为分子多项式系数，A 为分母多项式系数。 涉及函数：freqz,freqs. 3. 系统零极点分布与稳定性的判定。 对于连续时间系统，系统极点位于s 域左半平面，系统稳定。 对于离散时间系统，系统极点位于z 域单位圆内部，系统稳定。 涉及函数：zplane. 三、 实验内容 1． 验证性实验 a) 系统零极点的求解和作图
二阶系统性能改善与稳定性
例1 系统结构图如图所示。求开环增益K分别为10，0.5，0.09时系统的动态性能指标。 计算过程及结果列表 K 计算 10 0.5 0.09 开环 传递 函数 )1 ( 10 ) ( 1+ = s s s G )1 ( 5.0 ) ( 2+ = s s s G )1 ( 09 .0 ) ( 3+ = s s s G 闭环 传递 函数10 10 ) ( 2 1+ + = Φ s s s 5.0 5.0 ) ( 2 2+ + = Φ s s s 09 .0 09 .0 ) ( 2 3+ + = Φ s s s 特征 参数 ? ? ? ?? ? ? ? = = = ? = = = 81 arccos 158 .0 16 .3 2 1 16 .3 10 ξ β ξ ω n ? ? ? ?? ? ? ? = = = ? = = = 45 arccos 707 .0 707 .0 2 1 707 .0 5.0 ξ β ξ ω n ?? ? ? ? = ? = = = 67 .1 3.0 2 1 3.0 09 .0 ξ ω n 特征 根 12 .3 5.0 2,1 j ± - = λ5.0 5.0 2,1 j ± - = λ ? ? ? - = - = 9.0 1.0 2 1 λ λ ? ? ? = = 11 .1 10 2 1 T T 动态 性能 指标 2 2 1 00 00 1.01 1 60.4 3.5 3.5 7 0.5 p n s n t e t ξπξ π ξω σ ξω -- ? == ? - ? ? == ? ? ?=== ? ? ? ? ? ? ? ? ? ? ? = = = = = - = - - 7 5.3 5 238 .6 1 1 2 2 n s n p t e t ξω σ ω ξ π ξ ξπ() 1221 11 9 31 ,0 s s p T T t t T T t λλ σ ?== ? =?= ? ?=∞= ?
信号与系统_——零极点及稳定性响应
实验七、系统极零点及其稳定性 三、已知下列传递函数H(s)或H(z)，求其极零点，并画出极零图。 1. b=[3 -9 6]; a=[1 3 2]; zplane(b,a) 2. b=[1]; a=[1 0]; zplane(b,a)
3. b=[1 0 1]; a=[1 2 5]; zplane(b,a)
4. b=[1.8 1.2 1.2 3]; a=[1 3 2 1]; zplane(b,a) 五、求出系统的极零点，判断系统的稳定性。 5、先求出分子分母多项式系数 >> syms s >> zs=100*s*(s+2)^2*(s^2+3*s+2)^2; >> expand(zs) ans = 100*s^7+1000*s^6+4100*s^5+8800*s^4+10400*s^3+6400*s^2+1600*s >> syms s >> ps=(s+1)*(s-1)*(s^3+3*s^2+5*s+2)*((s^2+1)^2+3)^2; >> expand(ps) ans = -32-80*s-48*s^2+8*s^4-16*s^3+28*s^6+20*s^5+44*s^7+30*s^8+s^13+8*s^11+23*s^9+3*s^12 +11*s^10 再求出极零点 b=[100 1000 4100 8800 10400 6400 1600 0]; a=[1 3 8 11 23 30 44 28 20 8 -16 -48 -80 -32];
[z,p]=tf2zp(b,a) 求解结果： z = -2.0005 + 0.0005i -2.0005 - 0.0005i -1.9995 + 0.0005i -1.9995 - 0.0005i -1.0000 + 0.0000i -1.0000 - 0.0000i p = 1.0000 0.7071 + 1.2247i 0.7071 - 1.2247i 0.7071 + 1.2247i 0.7071 - 1.2247i -1.2267 + 1.4677i -1.2267 - 1.4677i -0.7071 + 1.2247i -0.7071 - 1.2247i -0.7071 + 1.2247i -0.7071 - 1.2247i -1.0000 -0.5466 极点不是都在左半平面，因此系统不稳定。 6、clear all; clc; num=conv([1 -1.414 1],[1 1]); den=conv([1 0.9 0.81],[1 -0.3]); [z,p]=tf2zp(num,den) zplane(z,p); z = -1.0000 0.7070 + 0.7072i 0.7070 - 0.7072i
实验二：系统稳定性和稳态性能分析
实验二：系统稳定性和稳态性能分析 主要内容： 自动控制系统稳定性和稳态性能分析上机实验 目的与要求： 熟悉 MATLAB 软件对系统稳定性分析的基本命令语句 熟悉 MATLAB 软件对系统误差分析的 Simuink 仿真 通过编程或 Simuink 仿真完成系统稳定性和稳态性能分析 一 实验目的 1、研究高阶系统的稳定性，验证稳定判据的正确性； 2、了解系统增益变化对系统稳定性的影响； 3、观察系统结构和稳态误差之间的关系。 二 实验任务 1、稳定性分析 欲判断系统的稳定性，只要求出系统的闭环极点即可，而系统的闭环极点就是闭环传递函数的分母多项式的根，可以利用MATLAB 中的tf2zp 函数求出系统的零极点，或者利用root 函数求分母多项式的根来确定系统的闭环极点，从而判断系统的稳定性。 （1）已知单位负反馈控制系统的开环传递函数为0.2( 2.5)()(0.5)(0.7)(3)s G s s s s s +=+++，用 MA TLAB 编写程序来判断闭环系统的稳定性，并绘制闭环系统的零极点图。 （2）已知单位负反馈控制系统的开环传递函数为( 2.5)()(0.5)(0.7)(3)k s G s s s s s +=+++，当取k =1，10，100用MA TLAB 编写程序来判断闭环系统的稳定性。 只要将（1）代码中的k 值变为1，10，100，即可得到系统的闭环极点，从而判断系统的稳定性，并讨论系统增益k 变化对系统稳定性的影响。 2、稳态误差分析 （1）已知如图所示的控制系统。其中2(5)()(10) s G s s s +=+，试计算当输入为单位阶跃信号、单位斜坡信号和单位加速度信号时的稳态误差。 从 Simulink 图形库浏览器中拖曳Sum （求和模块）、Pole-Zero （零极点）模块、Scope （示波器）模块到仿真操作画面，连接成仿真框图如右上图所示： （2）若将系统变为I 型系统，5()(10) G s s s =+，在阶跃输入、斜坡输入和加速度信
实验Z变换离散系统零极点分布和频率分析
实验三 Z 变换、离散系统零极点分布和频率分析 一、 实验目的 ● 学会运用MATLAB 求离散时间信号的z 变换和z 反变换； ● 学会运用MATLAB 分析离散时间系统的系统函数的零极点； ● 学会运用MATLAB 分析系统函数的零极点分布与其时域特性的关系； ● 学会运用MATLAB 进行离散时间系统的频率特性分析。 二、 实验仪器：电脑一台，MATLAB6.5或更高级版本软件一套。 三、 实验原理及实例分析 （一）离散时间信号的Z 变换 1．利用MATLAB 实现z 域的部分分式展开式 MATLAB 的信号处理工具箱提供了一个对F(Z)进行部分分式展开的函数residuez()，其调用形式为： [r,p,k]=residuez(num,den) 式中，num 和den 分别为F(Z)的分子多项式和分母多项式的系数向量，r 为部分分式的系数向量，p 为极点向量，k 为多项式的系数向量。 【实例3-1】 利用MATLAB 计算3 21431818 ) (-----+z z z z F 的部分分式展开式。 解：利用MATLAB 计算部分分式展开式程序为 % 部分分式展开式的实现程序 num=[18]; den=[18 3 -4 -1]; [r,p,k]=residuez(num,den) 2．Z 变换和Z 反变换 MATLAB 的符号数学工具箱提供了计算Z 变换的函数ztrans()和Z 反变换的函数iztrans （），其调用形式为 )()(F iztrans f f ztrans F ==
上面两式中，右端的f 和F 分别为时域表示式和z 域表示式的符号表示，可应用函数sym 来实现，其调用格式为 ()A sym S = 式中，A 为待分析的表示式的字符串，S 为符号化的数字或变量。 【实例3-2】求（1）指数序列()n u a n 的Z 变换;(2)()() 2 a z az z F -= 的Z 反变换。 解 （1）Z 变换的MATLAB 程序 % Z 变换的程序实现 f=sym('a^n'); F=ztrans(f) 程序运行结果为： z/a/(z/a-1) 可以用simplify( )化简得到 : -z/(-z+a) （2）Z 反变换的MATLAB 程序 % Z 反变换实现程序 F=sym('a*z/(z-a)^2'); f=iztrans(F) 程序运行结果为 f = a^n*n （二）系统函数的零极点分析 1. 系统函数的零极点分布 离散时间系统的系统函数定义为系统零状态响应的z 变换与激励的z 变换之比，即 ) () ()(z X z Y z H = （3-1） 如果系统函数)(z H 的有理函数表示式为：
零极点对系统的影响
增加零极点以及零极点分布对系统的影响一般说来，系统的极点决定系统的固有特性，而零点对于系统的暂态响应 和频率响应会造成很大影响。以下对于零极点的分布研究均是对于开环传递函 数。 零点一般是使得稳定性增加，但是会使调节时间变长，极点会使调节时间变短，是系统反应更快，但是也会使系统的稳定性变差。在波特图上反应为，增加一个零点会在幅频特性曲线上增加一个+20db/10倍频的曲线，幅频曲线上移，增加一个极点，会在幅频特性曲线上增加一个-20db/10倍频的曲线，幅频曲线下移。 在s左半平面增加零点时，会增加系统响应的超调量，带宽增大，能够减小系统的调节时间，增快反应速度，当零点离虚轴越近，对系统影响越大，当零点实部远大于原二阶系统阻尼系数ξ时，附加零点对系统的影响减小，所以当零点远离虚轴时，可以忽略零点对系统的影响。从波特图上来看，增加一个零点相当于增加一个+20db/10倍频的斜率，可以使的系统的相角裕度变大，增强系统的稳定性。 在s右半平面增加零点，也就是非最小相位系统，非最小相位系统的相位变化范围较大，其过大的相位滞后使得输出响应变得缓慢。因此，若控制对象是非最小相位系统，其控制效果特别是快速性一般比较差，而且校正也困难。对于非最小相位系统而言，当频率从零变化到无穷大时，相位角的便变化范围总是大于最小相位系统的相角范围，当ω等于无穷大时，其相位角不等于-（n-m）×90o。非最小相位系统存在着过大的相位滞后，影响系统的稳定性和响应的快速性。 在s左半平面增加极点时，系统超调量%pσ减小，调整时间st(s)增大，从波特图上看，s左半平面增加一个极点时，会在幅频特性曲线上增加一个-20db/10倍频的曲线，也就意味着幅频特性曲线会整体下移，导致相角域度减小，从而使得稳定性下降。当极点离原点越近，就会增大系统的过渡时间，使得调节时间增加，稳定性下降，当系统影响越大当极点实部远大于原二阶系统阻尼系数ξ时，附加极点对系统的影响减小，所以当极点远离虚轴时可以忽略极点对系统的影响。 在s右半平面增加极点会导致系统不稳定。 最小相位系统 从传递函数角度看,如果说一个环节的传递函数的极点和零点的实部全都小于或等于零,则称这个环节是最小相位环节.如果传递函数中具有正实部的零点或极点,或有延迟环节,这个环节就是非最小相位环节. 对于闭环系统,如果它的开环传递函数极点或零点的实部小于或等于零,则称它是最小相位系统.如果开环传递函中有正实部的零点或极点,或有延迟环节,则称系统是非最小相位系统.因为若把延迟环节用零点和极点的形式近似表达时(泰勒级数展开),会发现它具有正实部零点. 最小相位系统具有如下性质: 1,最小相位系统传递函数可由其对应的开环对数频率特性唯一确定;反之亦然. 2,最小相位系统的相频特性可由其对应的开环频率特性唯返航一确定;反之亦然. 3,在具有相同幅频特性的系统中,最小相位系统的相角范围最小.
滤波器稳定性与极点
在数字信号处理中，系统的稳定性是一个很重要的问题，比如说在滤波器的设计中，都要求系统必须稳定，否则是无法使用的。那么，如何判断系统是否稳定呢？ 从定义上说，如果输入有界，则输出必定有界的系统是稳定的。从数学上可以推导出，因果系统冲击响应Z变换的收敛域包含单位圆的系统是稳定的。从零点极点的角度，则是系统函数的所有极点都在单位圆内的系统是稳定的。如何来理解呢？ 我们先以一个简单的单极点系统为例来理解系统的稳定性。比如有一个单极点系统： H(z)=1/(1-2z-1) 表示的是如下的如下的信号处理过程：系统当前输出是当前的输入加上2倍的系统上一时刻输出之和。这个系统是不稳定的，因为当前输出需要放大上一个时刻的输出，这也就是说，系统存在的自激的过程，直观上我们就可以很好地理解，自激系统是不稳定的。从分析极点的角度看，这个系统的极点为2，在单位圆外，与数学上的分析是一致的。极点在单位圆内的要求，对一阶极点而言，实际上也就是直观上要求系统不能自激。 对于高阶极点的情况，由代数学可知，高阶极点可进行分式的分解，也即是高阶极点可以分解成多个一阶极点并联而成的系统，在并联系统中，只要有一个系统不稳定，整个系统就是不稳定的。这与数学上要求的所有极点都在单位圆内是对应的。对于更一般的既包含零点又包含极点的系统，可以看成一个全零点系统和全极点系统串接而成，零点与系统的稳定性无关，分析和结论与高阶全极点系统完全一致。 在滤波器的设计中，可以很方便地通过调整极点改变滤波器的特性。而在许多设计精巧的滤波器中，极点往往在单位圆上或单位圆附近，在实际中还要考虑量化及数的精度等问题，确保系统的稳定性。
离散系统的频率响应分析和零极点分布
离散系统的频率响应分析和零极点分布 Company number：【0089WT-8898YT-W8CCB-BUUT-202108】
实验2 离散系统的频率响应分析和零、极点分布一、实验目的 通过MATLAB仿真简单的离散时间系统，研究其时域特性，加深对离散系统的冲激响应，频率响应分析和零、极点分布的概念的理解。 二、基本原理 离散系统的时域方程为 其变换域分析方法如下： 频域 ) ( ) ( ) ( ] [ ] [ ] [ ] [ ] [ω ω ωj j j m e H e X e Y m n h m x n h n x n y= ? - = * =∑∞ -∞ = 系统的频率响应为 ω ω ω ω ω ω ω jN N j jM M j j j j e d e d d e p e p p e D e p e H - - - - + + + + + + = = ... ... ) ( ) ( ) ( 1 1 Z域 ) ( ) ( ) ( ] [ ] [ ] [ ] [ ] [z H z X z Y m n h m x n h n x n y m = ? - = * =∑∞ -∞ = 系统的转移函数为 N N M M z d z d d z p z p p z D z p z H - - - - + + + + + + = = ... ... ) ( ) ( ) ( 1 1 1 1 分解因式 ∏- ∏- = ∑ ∑ = = - = - = - = - N i i M i i N i i k M i i k z z K z d z p z H 1 1 1 1 ) 1( ) 1( ) ( λ ξ ，其中i ξ 和i λ 称为零、极点。 在MATLAB中，可以用函数[z，p，K]=tf2zp（num，den）求得有理分式形式的系统转移函数的零、极点，用函数zplane（z，p）绘出零、极点分布图；也可以用函数zplane（num，den）直接绘出有理分式形式的系统转移函数的零、极点分布图。 另外，在MATLAB中，可以用函数 [r，p，k]=residuez（num，den）完成部分分式展开计算；可以用函数sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系统的串联。
极点与系统稳定性
极点对系统性能影响 一．控制系统与极点 自动控制系统根据控制作用可分为：连续控制系统和采样控制系统，采样系统又叫离散控制系统。通常把系统中的离散信号是脉冲序列形成的离散系统,称为采样控制系统。连续控制系统即指控制量为连续的模拟量如时变系统。 系统的数学模型一般由系统传递函数表达。传递函数为零初始条件下线性系统响应（即输出）量的拉普拉斯变换（或z变换）与激励（即输入）量的拉普拉斯变换之比。记作Φ（s）=Xo（s）/Xi（s），其中Xo（s）、Xi（s）分别为输出量和输入量的拉普拉斯变换。 特征方程的根称为极点。如试Φ﹙S﹚= C [∏（S-Pi）/∏(S-Qi) ]中Q1 Q2 Q3 ……Qi ……即为系统的极点。 二．极点对系统的影响 极点--确定了系统的运动模态；决定了系统的稳定性。下面对连续系统与离散系统分别进行分析： ⑴连续系统 理论分析：连续系统的零极点分布有如下几种形式 设系统函数为： 将H(S)进行部分分式展开：
系统冲激响应H(S)的时域特性h(t)随时间衰减的信号分量完全由系统函数H(S)的极点位置决定。每一个极点将决定h(t)的一项时间函数。 稳定性：由上述得知Y(S)= C [∏（S-Pi ）/(S-Qi) ]可分解为Y(S)=C1/(S-τ1)+ C2/(S-τ2)+ C3/(S-τ3)+……+ Ci/(S-τi)+…… 则时间响应为 …… 由于特征方程的根不止一个，这时，应把系统的运动看成是多个运动分量的合成。只要有一个运动分量是发散的，则系统是不稳定的。因此，特征方程所有根的实部都必须是负数，亦即所有的根都在复平面的左半平面。 通过复变函数幅角定理将S 由G 平面映射到GH 平面。 如果封闭曲线 F 内有Z 个F(s)的零点，有P 个F(s)的极点，则s 沿 F 顺时针转一圈时，在F(s)平面上，F(s)曲线绕原点顺时针转的圈数R 为z 和p 之差，即R ＝z －p 。 若R 为负,表示F(s)曲线绕原点逆时针转过的圈数。 F(s)的分母是G0(s)的分母，其极点是G0(s)的极点；其分子是?(s)的分母，即?(s)的特征多项式，其零点是?(s)的极点。 取D 形曲线（D 围线）如图所示，是整个右半复平面。 且设D 曲线不经过F(s)的任一极点或零点。 s 沿D 曲线顺时针变化一周，F(s)顺时针包围原点的周数为： n=z-p=F(s)在右半复平面的零点数(闭环传函在右半复平面极点数) -F(s)在右半复平面的极点数(开环传函在右半复平面极点数) 所以闭环系统稳定的充分必要条件是： n=- p =-开环传函在右半复平面的极点数 1212()n s t s t s t n y t C e C e C e =+++ 0()0()0()0()t s y t y t Ce y t y t t ααααα=<→?? ===??>→∞?→∞（1）只有一个实根：时，时，恒量时，()()121 ()0cos()00j t j t t s j y t C e C e C e t t αωαωααωαω?αα+-=±=+? →∞ （2）有一对复根：时，收敛时，等幅振荡时，发散
零极点对系统性能的影响分析
摘要 本次课程设计主要是分析零极点对系统性能的影响。首先从根轨迹、奈奎斯特 曲线、伯德图和阶跃响应四方面分析原开环传递函数时的系统性能，然后在原开环 传递函数基础上增加一个零点，并且让零点的位置不断变化，分析增加零点之后系 统的性能，同时与原系统进行分析比较，发现增加的零点与虚轴的距离决定了对系 统影响的大小；再在原开环传递函数基础上增加一个极点，并且令极点位置不断变 化，分析增加极点后系统的性能，同时与原系统进行分析比较，同样发现增加的极 点与虚轴的距离决定了对系统的影响大小。 关键词：零极点开环传递函数系统性能 MATLAB 谐振带宽 The curriculum design is mainly the analysis of effect of zero pole on the performance of the system. First from the root locus, Nyquist curve, Bode diagram and step response analysis of four aspects of the original open-loop transfer function of the system performance, and then in the original open-loop transfer function is added on the basis of a zero, and let the zero point position changes continuously, increase system performance analysis of zero, at the same time and the original system analysis that increase, the zeros and the imaginary axis distance determines the impact on the system size; adding a pole in the original open-loop transfer function based on pole position, and make the changes, analysis of increasing performance point system, at the same time and the analysis of the original system, also found that increasing pole and the imaginary axis distance determines the impact on the size of the system. Keywords: zero pole open loop transfer function of system performance of MATLAB resonant bandwidth
实验-Z变换、零极点分析
（一）离散时间信号的Z 变换 1．利用MATLAB 实现z 域的部分分式展开式 MATLAB 的信号处理工具箱提供了一个对F(Z)进行部分分式展开的函数residuez()，其调用形式为： [r,p,k]=residuez(num,den) 式中，num 和den 分别为F(Z)的分子多项式和分母多项式的系数向量，r 为部分分式的系数向量，p 为极点向量，k 为多项式的系数向量。 【实例1】 利用MATLAB 计算321431818)(-----+z z z z F 的部分分式展开式。 解：利用MATLAB 计算部分分式展开式程序为 % 部分分式展开式的实现程序 num=[18]; den=[18 3 -4 -1]; [r,p,k]=residuez(num,den) 2．Z 变换和Z 反变换 MATLAB 的符号数学工具箱提供了计算Z 变换的函数ztrans()和Z 反变换的函数iztrans （），其调用形式为
)()(F iztrans f f ztrans F == 上面两式中，右端的f 和F 分别为时域表示式和z 域表示式的符号表示，可应用函数sym 来实现，其调用格式为 ()A sym S = 式中，A 为待分析的表示式的字符串，S 为符号化的数字或变量。 【实例2】求（1）指数序列()n u a n 的Z 变换;(2)()()2a z az z F -= 的Z 反变换。 解 （1）Z 变换的MATLAB 程序 % Z 变换的程序实现 f=sym('a^n'); F=ztrans(f) 程序运行结果为： z/a/(z/a-1) 可以用simplify( )化简得到 : -z/(-z+a) （2）Z 反变换的MATLAB 程序 % Z 反变换实现程序 F=sym('a*z/(z-a)^2'); f=iztrans(F) 程序运行结果为 f = a^n*n （二）系统函数的零极点分析 1. 系统函数的零极点分布 离散时间系统的系统函数定义为系统零状态响应的z 变换与激励的z 变换之比，即 )()()(z X z Y z H = （3-1）
实验z变换、零极点分析
1. 学会运用MATLAB 求离散时间信号的z 变换和z 反变换； （一）离散时间信号的Z 变换 1．利用MATLAB 实现z 域的部分分式展开式 MATLAB 的信号处理工具箱提供了一个对F(Z)进行部分分式展开的函数residuez()，其调用形式为： [r,p,k]=residuez(num,den) 式中，num 和den 分别为F(Z)的分子多项式和分母多项式的系数向量，r 为部分分式的系数向量，p 为极点向量，k 为多项式的系数向量。 【实例1】 利用MATLAB 计算3 21431818 ) (-----+z z z z F 的部分分式展开式。 解：利用MATLAB 计算部分分式展开式程序为 % 部分分式展开式的实现程序 num=[18]; den=[18 3 -4 -1]; [r,p,k]=residuez(num,den)
2．Z 变换和Z 反变换 MATLAB 的符号数学工具箱提供了计算Z 变换的函数ztrans()和Z 反变换的函数iztrans （），其调用形式为 )() (F iztrans f f ztrans F == 上面两式中，右端的f 和F 分别为时域表示式和z 域表示式的符号表示，可应用函数sym 来实现，其调用格式为 ()A sym S = 式中，A 为待分析的表示式的字符串，S 为符号化的数字或变量。 【实例2】求（1）指数序列()n u a n 的Z 变换;(2)()() 2 a z az z F -= 的Z 反变换。 解 （1）Z 变换的MATLAB 程序 % Z 变换的程序实现 f=sym('a^n'); F=ztrans(f) 程序运行结果为： z/a/(z/a-1) 可以用simplify( )化简得到 : -z/(-z+a) （2）Z 反变换的MATLAB 程序 % Z 反变换实现程序 F=sym('a*z/(z-a)^2'); f=iztrans(F) 程序运行结果为 f = a^n*n （二）系统函数的零极点分析