信号与系统实验7连续系统零极点分析.docx
《信号与系统实验7连续系统零极点分析.docx》由会员分享,可在线阅读,更多相关《信号与系统实验7连续系统零极点分析.docx(22页珍藏版)》请在冰豆网上搜索。
![信号与系统实验7连续系统零极点分析.docx](https://file1.bdocx.com/fileroot1/2022-12/14/f933a350-44dd-438c-a5ac-2046f21fa4bf/f933a350-44dd-438c-a5ac-2046f21fa4bf1.gif)
信号与系统实验7连续系统零极点分析
实验七连续时间系统S域零极点分析
一、目的
(1)掌握连续系统零极点分布与系统稳定性关系
(2)掌握零极点分布与系统冲激响应时域特性之间的关系
(3)掌握利用MATLAB进行S域分析的方法
二、零极点分布与系统稳定性
根据系统函数H(s)的零极点分布来分析连续系统的稳定性是零极点分析的重要应
用之一。
稳定性是系统固有的性质,与激励信号无关,由于系统函数H(s)包含了系统
的所有固有特性,显然它也能反映出系统是否稳定。
对任意有界信号f(t),若系统产生的零状态响应y(t)也是有界的,则称该系统为稳定系统,否则,则称为不稳定系统。
上述稳定性的定义可以等效为下列条件:
时域条件:
连续系统稳定充要条件为h(t)dt,即冲激响应绝对可积;
复频域条件:
连续系统稳定的充要条件为系统函数H(s)的所有极点位于S平面的左半平面。
系统稳定的时域条件和频域条件是等价的。
因此,只要考察系统函数H(s)的极点
分布,就可判断系统的稳定性。
对于三阶以下的低阶系统,可以利用求根公式方便地求出极点位置,从而判断系统稳定性,但对于告阶系统,手工求解极点位置则显得非常困难。
这时可利用MATLAB来实现这一过程。
例7-1:
已知某连续系统的系统函数为:
H(s)
试用MATLAB求出该系统的零极点,
s23s2
8s42s33s2s5
画出零极点图,并判断系统是否稳定
解:
调用实验六介绍的绘制连续系统零极点图函数sjdt即可解决此问题,对应的
MATLAB命令为:
a=[82315];b=[132];
[p,q]=sjdt(a,b)
运行结果为:
P=
-0.6155-0.6674i-0.6155+0.6674i0.4905-0.7196i0.4905+0.7196i
q=
-2-1
绘制的零极点图如图7-1所示。
由程序运行结果可以看出,该系统在S平面的右半平面有一对共轭极点,故该系统
是一个不稳定系统。
三、零极点分布与系统冲激响应时域特性
设连续系统的系统函数为H(s),冲激响应为h(t),贝U
H(s)0h(t)estdt
显然,H(s)必然包含了h(t)的本质特性。
对于集中参数的LTI连续系统,其系统函数可表示为关于s的两个多项式之比,即
B(s)(sq」)
(7-1)
H(s)少C誤
A(s)(sPi)
i1
连续系诜零极点图
实釉
u
1
J
图7-1例7-1的系统零极点图
其中qj(j1,2,,M)为H(s)的M个零点,Pi(i1,2,,N)为H(s)的N个极点。
若系统函数的N个极点是单极点,则可将H(s)进行部分分式展开为:
Nk・
H(s)丄(7-2)
i1sPi
从式(7-1)和(7-2)可以看出,系统冲激响应h(t)的时域特性完全由系统函数H(s)的极点位置决定。
H(s)的每一个极点将决定h(t)的一项时间函数。
显然,H(s)的极点位置不同,则h(t)的时域特性也完全不同。
下面利用例子说明H(s)的极点分布与h(t)时域特性之间的关系。
例7-2:
已知连续系统的零极点分布如图7-2所示,试用MATLAB分析系统冲激响应h(t)的时域特性。
解:
系统的零极点图已知,则系统的系统函数H(s)就可确定。
这样就可利用绘制连续
系统冲激响应曲线的MATLAB函数impulse。
,将系统冲激响应h(t)的时域波形绘制出来。
对于图7-2(a)所示的系统,系统函数为H(s)1,即系统的极点位于原点,绘制
s
冲激响应时域波形的MATLAB命令如下:
a=[10];
b=[1];
impulse(b,a)
绘制的冲激响应h(t)波形如图7-3(a)所示,此时h(t)为单位阶跃信号。
图7-2例7-2的系统零极点图
叶」址阳耳城n萌
"ineaecI
(c)
Q由一_百・—*
DD?
n,4OEOB『m(3tc.
(a)(b)
(e)
(d)
『」
图7-3例7-2的系统冲激响应时域波形图
对于图7-2(b)所示的系统,系统函数为H(s)丄,即系统的极点为位于S平
s
面左半平面的实极点,令
2,绘制冲激响应时域波形的MATLAB命令如下:
a=[12];
b=[1];
impulse(b,a)
绘制的冲激响应h(t)波形如图7-3(b)所示,此时h(t)为衰减指数信号
对于图7-2(c)所示的系统,系统函数为H(s),即系统的极点为位于S平
面右半平面的实极点,令a=[1-2];
b=[1];
s
2,绘制冲激响应时域波形的MATLAB命令如下:
impulse(b,a)
绘制的冲激响应h(t)波形如图7-3(c)所示,此时h(t)为随时间增长的指数信号
对于图7-2(d)所示的系统,系统函数为H(s)2,即系统的极点为位
(s)
于S平面左半平面的一对共轭极点,令0.5、4,绘制冲激响应时域波形的
MATLAB命令如下:
a=[1116.25];
b=[1];impulse(b,a,5)
绘制的冲激响应h(t)波形如图7-3(d)所示,此时h(t)为按指数衰减的正弦振荡信号。
对于图7-2(e)所示的系统,系统函数为H(s),即系统的极点为位于S
s
平面虚轴上的一对共轭极点,令4,绘制冲激响应时域波形的MATLAB命令如下:
a=[1016];
b=[1];
impulse(b,a,5)
绘制的冲激响应h(t)波形如图7-3(e)所示,此时h(t)为等幅正弦振荡信号。
对于图7-2(f)所示的系统,系统函数为H(s)2—2,即系统的极点为位
(s)
于S平面右半平面上的一对共轭极点,令0.5、4,绘制冲激响应时域波形的
MATLAB命令如下:
a=[1-116.25];
b=[1];
impulse(b,a,5)
绘制的冲激响应h(t)波形如图7-3(f)所示,此时h(t)为按指数增长的正弦振荡信号。
从上述程序运行结果和绘制的系统冲激响应曲线,可以总结出以下规律:
系统冲激响应h(t)的时域特性完全由系统函数H(s)的极点位置决定,H(s)位于S平面左半平面的极点决定了h(t)随时间衰减的信号分量,位于S平面虚轴上的极点决定了冲激响应的稳态信号分量,位于S平面右半平面的极点决定了冲激响应随时间增长的信号分量。
三、由连续系统零极点分布分析系统的频率特性
由前面分析可知,连续系统的零极点分布完全决定了系统的系统函数H(s),显然,
H(j)的方法-
系统的零极点分布也必然包含了系统的频率特性。
下面介绍如何通过系统的零极点分布来直接求出系统的频率响应
H(j)的一种直
几何矢量法,以及如何用MATLAB来实现这一过程。
几何矢量法是通过系统函数零极点分布来分析连续系统频率响应
观而又简便的方法。
该方法将系统函数的零极点是为S平面上的矢量,通过对这些矢量
的模和幅角的分析,即可快速确定出系统的幅频响应和相频响应。
其基本原理如下:
设某连续系统的系统函数为:
M
B(s)(sqj)
H(s)-(-)C4
A(s)(sP)
i1
其中qj(J
1,2,,M)为H(s)的M个零点,Pi(i1,2,,N)为H(s)的N个极点。
则频率
响应为:
M
(sqj)
H(j)H(s)SjC稈(7-3)
(spj
i1
现在从几何矢量空间的角度分析S平面,即将S平面的任一点看成是从原点到该点
的矢量,则j即是从S平面原点到虚轴上角频率为的点的矢量。
同理,qj(j1,2,,M)
和Pi(i1,2,,N)即是从S平面原点到系统函数各零点和极点的矢量。
现在考虑矢量jqj,由矢量运算可知,它实际上就是零点q到虚轴上角频率为
的点的矢量,如图7-3所示;而矢量j
Pi则是极点Pi到虚轴上角频率为的点的矢量
图7-3连续系统几何矢量法示意图
jqBjJ
jPiAej
其中,Bj为矢量j的幅角。
因此有:
qj的模,
j为该矢量的幅角;
A为矢量j
Pi的模,i为该矢量
Mjj
(7-4)
J)
H(j)C|H(j)ej()
(Aeji)
i1
则系统的幅频特性和相频特性为:
M
Bj
H(j
)
cj1
C—N
(7-5)
()
A
i1
MN
ji
j1i1
(7-6)
由上述分析可以得出如下结论:
连续系统的幅频响应H(j)等于系统函数所有零点到虚轴上角频率为的点
的距离之积与系统函数所有极点到虚轴上角频率为的点的距离之积的比值;
连续系统的相频响应()等于系统函数所有零点到虚轴上角频率为的点的
矢量相角之和与系统函数所有极点到虚轴上角频率为的点的矢量相角之和的
差值。
让矢量j沿着虚轴变化,即角频率由0~进行改变,便可直观地求出系统幅频
响应和相频响应随的变化,从而分析出系统的频率特性。
根据上述结论,若已知系统的零极点分布,即可直接由几何矢量法分析出系统的频率特性。
上述过程可用MATLAB快速实现。
用MATLAB实现已知系统零极点分布,求系统频率响应,并绘制其幅频特性和相频特性曲线的程序流程如下:
(1)定义包含系统所有零点和极点位置的行向量q和p;
(2)定义绘制系统频率响应曲线的频率范围向量fl和f2、频率取样间隔k,并产生频率等分点向量f;
(3)求出系统所有零点和极点到这些等分点的距离;
(4)求出系统所有零点和极点到这些等分点的矢量相角;
(5)根据式(7-5)和(7-6)求出fl到f2频率范围内各频率等分点的|H(j)和
();
(6)绘制f1~f2频率范围内系统的幅频特性曲线和相频特性曲线。
下面是完成上述分析过程的MATLAB实用函数splxy()。
functionsplxy(f1,f2,k,p,q)
%根据系统零极点分布绘制系统频率响应曲线程序
%f1、f2:
绘制频率响应曲线的频率范围(即频率起始和终止点,单位为赫兹)%p、q:
系统函数极点和零点位置行向量
%k:
绘制频率响应曲线的频率取样间隔
p=p';
q=q';
f=f1:
k:
f2;%定义绘制系统频率响应曲线的频率
范围
w=f*(2*pi);
y=i*w;
n=length(p);
m=length(q);
ifn==0%如果系统无极点
yq=ones(m,1)*y;
vq=yq-q*ones(1,length(w));
bj=abs(vq);
cosaij=angle(vq)./pi.*180;
ai=1;
%如果系统无零点
thetai=O;
elseifm==0
yp=ones(n,1)*y;
vp=yp-p*ones(1,length(w));ai=abs(vp);
thetai=angle(vp)./pi.*180;
bj=1;
cosaij=0;
else
yp=ones(n,1)*y;
yq=ones(m,1)*y;
vp=yp-p*ones(1,length(w));vq=yq-q*ones(1,length(w));ai=abs(vp);
thetai=angle(vp)./pi.*180;
bj=abs(vq);
cosaij=angle(vq)./pi.*180;endsubplot(121);Hw=prod(bj,1)./prod(ai,1);plot(f,Hw);
title('连续系统幅频响应曲线')xlabel(濒率w(单位:
赫兹)')ylabel('F(jw)')subplot(122);
Angw=sum(cosaij,1)-sum(thetai,1);plot(f,Angw);
title('连续系统相频响应曲线')xlabel(濒率w(单位:
赫兹)')ylabel('Angle(jw)')
下面举例说明如何调用该函数。
例7-3:
已知系统的系统函数为H(s),试用MATLAB绘制出该系统的
(s50)(s100)
频率特性曲线(幅频曲线和相频曲线)。
解:
通过调用上述实用函数求解,命令如下:
q=[0];
p=[-100-50];
f仁0;
f2=100;
k=0.01;
splxy(f1,f2,k,p,q);
运行结果如图7-4所示。
连擴乔统相频吋应脏线
100
00
B0
4J
-23
-40•
-60-
□1111
020^0ED00100
坝率创(莹位:
赫勒
2D4000EO100
频率卑程:
描船)
图7-4系统幅频特性与相频特性曲线
四、实验内容
已知连续系统的系统函数分别如下所示,试用MATLAB绘制系统的零极点图,并
根据零极点图判断系统的稳定性,同时绘制出系统的频率特性曲线(提示:
系统的零极
点可用MATLAB
(1)H(s)
(3)H(s)
函数roots求得)。
Ss2
32
3s5s4s6
2
3s(s9)
s420s264
(2)H(s)
(4)H(s)
2(^4s5)
2
s4s5
1
32
s2s2s1
function[p,q]=sjdt(A,B)
%绘制连续系统零极点图程序
%A:
系统函数分母多项式系数向量
%B:
系统函数分子多项式系数向量
%p:
函数返回的系统函数极点位置行向量%q:
函数返回的系统函数零点位置行向量p=roots(A);
q=roots(B);
p=p';
q=q';
x=max(abs([pq]));
x=x+0.1;
y=x;
elf
%求系统极点
%求系统零点
%将极点列向量转置为行向量
%将零点列向量转置为行向量%确定纵坐标范围
%确定横坐标范围
axis([-xx-yy]);axis('square')plot([-xx],[00])plot([00],[-yy])
holdon
%确定坐标轴显示范围
%画横坐标轴
%画纵坐标轴
%画极点
%画零点
%标注标题
plot(real(p),imag(p),'x')plot(real(q),imag(q),'o')title('连续系统零极点图')text(0.2,x-0.2,'虚轴')
text(y-0.2,0.2,'实轴')
functionsplxy(f1,f2,k,p,q)%根据系统零极点分布绘制系统频率响应曲线程序
%f1、f2:
绘制频率响应曲线的频率范围(即频率起始和终止点,单位为赫兹)%p、q:
系统函数极点和零点位置行向量
%k:
绘制频率响应曲线的频率取样间隔
p=p';
q=q';
f=f1:
k:
f2;%定义绘制系统频率响应曲线的频率范围
w=f*(2*pi);
y=i*w;
n=length(p);m=length(q);
ifn==0%如果系统无极点
yq=ones(m,1)*y;
vq=yq-q*ones(1,length(w));
bj=abs(vq);cosaij=angle(vq)./pi.*180;
ai=1;
thetai=0;
elseifm==0%如果系统无零点
yp=ones(n,1)*y;
vp=yp-p*ones(1,length(w));
ai=abs(vp);thetai=angle(vp)./pi.*180;
bj=1;cosaij=0;
else
yp=ones(n,1)*y;yq=ones(m,1)*y;
vp=yp-p*ones(1,length(w));
vq=yq-q*ones(1,length(w));ai=abs(vp);
thetai=angle(vp)./pi.*180;
bj=abs(vq);
cosaij=angle(vq)./pi.*180;
endsubplot(121);
Hw=prod(bj,1)./prod(ai,1);
plot(f,Hw);
title('连续系统幅频响应曲线')xlabel('频率w(单位:
赫兹)')ylabel('F(jw)')
subplot(122);
Angw=sum(cosaij,1)-sum(thetai,1);
plot(f,Angw);
title('连续系统相频响应曲线')xlabel('频率w(单位:
赫兹)')ylabel('Angle(jw)')
1.求零极点
A=[3546];
B=[112];
[p,q]=sjdt(A,B)
-1.6101
-0.0283-1.1142i-0.0283+1.1142i
q=
-0.5000-1.3229i-0.5000+1.3229i
连娃痕垛参磁点也
1j5・虚錮
0.5・
0—
£.弓-
-1.5-
Ii.1--・iI
■1.&-1C.5□0.611.5
系统为稳定系统
绘制出系统的频率特性曲线:
f1=0;
f2=5;
k=0.01;
splxy(f1,f2,k,p,q);
2.求零极点
A=[145];B=[2-810];[p,q]=sjdt(A,B)
-2.0000-I.OOOOi-2.0000+I.OOOOi
2.0000-1.0000i2.0000+1.0000i
系统为稳定系统
绘制出系统的频率特性曲线:
f1=0;f2=50;k=0.01;splxy(f1,f2,k,p,q);
3Q9rf1I1
IJB
an
sogI-
SB
TO
ioa丨-
叫
_.
14
>2
0--3Q
壬〒刀<50
類嚴”練It)
3.求零极点
A=[1020064];B=[30-270];
[p,q]=sjdt(A,B)
0.0000-4.0000i0.0000+4.0000i
0-2.0000i
0+2.0000i
-3
Eftt
系统为不稳定系统绘制出系统的频率特性曲线:
f1=0;
f2=5;
k=0.01;
T»
wfl
rm
4.求零极点:
A=[1221];
B=[1];
[p,q]=sjdt(A,B)
-1.0000
-0.5000-0.8660i-0.5000+0.8660i
Emptymatrix:
1-by-0
1
-
E
o.s
-
0.6
L
04
02
0
-02
-
44
JJ6
L
48
L
E
-1
-1
4JB-OS-D4-012
0204
0*6
系统为稳定系统
绘制出系统的频率特性曲线:
f1=0;
f2=5;
k=0.01;
splxy(f1,f2,k,p,q);
逋您墓电E峯射咱应由⑥曲戎
帅皇rl氯恆:
林U7(Affi:
炉IE,