仿真试验2系统的频率响应和稳定性讲解.docx
《仿真试验2系统的频率响应和稳定性讲解.docx》由会员分享,可在线阅读,更多相关《仿真试验2系统的频率响应和稳定性讲解.docx(20页珍藏版)》请在冰豆网上搜索。
![仿真试验2系统的频率响应和稳定性讲解.docx](https://file1.bdocx.com/fileroot1/2022-10/20/ed64ceff-4863-4265-8bac-a74f9e4f524a/ed64ceff-4863-4265-8bac-a74f9e4f524a1.gif)
仿真试验2系统的频率响应和稳定性讲解
2系统的频率响应和稳定性研究
.实验目的
1.绘制并观察典型系统的开环幅频曲线。
2.绘制并观察典型系统的开环对数频率曲线。
3.运用恩奎斯特准则判断闭环系统的稳定性。
二.实验要求
1.根据所给开环传递函数的机构形式,绘制相应的开环幅频曲线和开环对数频率曲线。
2.如绘制的开环幅相曲线不封闭,或用文字说明所缺部分曲线的走向,或在图上添加所缺曲线;曲线与(-1,j0)点的几何关系应足够清晰,能够支持判断结论的导出。
3.对该开环传递函数构成的单位负反馈系统的稳定性做出判断,说明理由;假如闭环不稳定,则应指出不稳定极点的数目。
三.实验内容
1.根据所给开环传递函数的结构形式,首先绘制出相应的开环幅频曲线和开环对数频率曲线。
2.对于存在积分环节的开环传递来说,因为得到的开环幅相曲线不封闭,所以需在图上添加所缺曲线,以使曲线与(-1,j0)点的几何关系清晰,支持判断结论的准确导出。
3.最后,利用开环幅频稳定判据(恩奎斯特准则)或开环对数频率稳定判据对开环传递函数构成的单位负反馈系统的稳定性作出判断;假如闭环不稳定,则指出不稳定极点的数目
K
(1)开环传递函数的形式为G1(Ts1)(Ts1),其中K,T1,T2(T1s1)(T2s1)可取大于0的任意数。
1举例,如令T1=1,T2=2,K=1,则G1,此时的指令如下:
(s1)(2s1)零极点形式的传递函数指令:
G=zpk([],[-1,-1/2],1);得到开环幅频曲线(恩奎斯特曲线):
figure
(1);nyquist(G);得到开环对数频率曲线:
figure
(2);margin(G);可以利用零极点形式的时域指令进行验证结果,也就是看闭环实部根是否都<0,此时的指令如下:
由零极点形式转换为因子式形式:
[n1,d1]=zp2tf([],[-1,-1/2],1);G=n1+d1;时域闭环根:
roots(G);
BodeDiagram
Gm=InfdB(atInfrad/sec),Pm=93.3deg(at0.666rad/sec)20
)Bd(edutinga
-2-1012
1010101010
Frequency(rad/sec)
因子式形式的开环频域指令:
因子式形式的传递函数指令:
G=tf([0,0,1],[2,3,1])得到开环幅频曲线(恩奎斯特曲线):
figure(3);nyquist(G)得到开环对数频率曲线:
figure(4);margin(G)可以利用零极点形式的时域指令进行验证结果,也就是看闭环实部根是否都<0,此时的指令如下:
由零极点形式转换为因子式形式:
n1=[0,0,1],d1=[2,3,1];G=n1+d1;
时域闭环根:
roots(G);
0.8
2)
sixAyranigam
NyquistDiagram
0.6
0.4
0.2
0
-0.2
-0.4
-0.6
-1-0.8-0.6-0.4-0.200.20.40.60.8RealAxis
-0.8
BodeDiagram
)Bd(edutinga
-45
-90
-135
-180
-2
10
G2
-1
10
0
10
1
10
2
10
Frequency(rad/sec)
K
(T1s1)(T2s1)(T3s1),
其中K,T1,T2,T3
可取大于0
的任意数。
如令T1=1,T2=2,T3=3,K=1,则G2(s1)(2s1)(3s1)
NyquistDiagram
5
4
3
2
is1ixA
0
Im-1
-2
-3
-4-5
-2-10123456
RealAxis
BodeDiagram
Gm=4.44dB(at1rad/sec),Pm=17.7deg(at0.781rad/sec)
50
0
-50
-100
-2
10
-1
10
0
10
1
10
2
10
Frequency(rad/sec)
.2
0
sixAyranigam
)Bd(edutinga
)ged(esah
NyquistDiagram
-0.6
-0.4-0.2
RealAxis
00.20.4
51.000-
2
5050
-4-938
--11
BodeDiagram
Gm=InfdB(atInfrad/sec),Pm=Inf
-1
10
0
10
1
10
2
10
1s(Ts
Frequency(rad/sec)
其中K,T1可取大于0的任意
如令T1=1,K=1,则G3s(s1)
0
)Bd(edutinga
20
15
NyquistDiagram
10
5
-5
-10
-15
-1-0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.10RealAxis
-20
BodeDiagram
Gm=InfdB(atInfrad/sec),Pm=51.8deg(at0.786rad/sec)
50
-50
-100
-90
-180
-2
10
e-135
seahP-135
P
-1
10
0
10
1
10
2
10
Frequency(rad/sec)
sixAyranigam
-0.5
-1
NyquistDiagram
0.4
0.3
0.2
0.1
0
-0.1
-0.2
-0.3
-0.4
-0.8-0.6-0.4-0.200.20.40.60.81RealAxis
BodeDiagram
)Bd(edutinga
-90
-2
10
-45
-1
10
0
10
Frequency(rad/sec)
1
10
2
10
s(T1s1)(T2s1)
令T1=1,T2=2,K=1
NyquistDiagram
-0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.10RealAxis
sixAyranigam
BodeDiagram
Gm=InfdB(atInfrad/sec),Pm=51.8deg(at0.786rad/sec)
)Bd(edutinga
0
005
0-93
1-1
)ged(esah
-1
10
0
10
1
10
2
10
Frequency(rad/sec)
sixAyranigam
-0.5
-1
NyquistDiagram
0.4
0.3
0.2
0.1
-0.1
-0.2
-0.3
-0.4
-0.8-0.6-0.4-0.200.20.40.60.81RealAxis
BodeDiagram
)Bd(edutinga
-90
-2
10
5)G5
-45
-1
10
0
10
1
10
2
10
Frequency(rad/sec)
K(Tas1),s(T1s1)(T2s1),
其中
K可取大于0的任意数。
令Ta=1,T1=1,T2=2,K=1
10
NyquistDiagram
40
60
20isxAran0
Im
-20
-60
-4
-3.5
-3
-2.5
-2
RealAxis
-1.5
-1-0.50
)Bd(edutinga
00
BodeDiagram
Gm=InfdB(atInfrad/sec),Pm=28deg(at0.94rad/sec)60
)ged(esah
-2
10
-1
10
0
10
1
10
Frequency(rad/sec)
11
BodeDiagram
Gm=Inf,Pm=-180deg(at0rad/sec)0
-5
-10
-15
-20
-25
-300
-90
-2-101
10101010Frequency(rad/sec)
K
(6)G62,其中K,T1可取大于0的任意数。
s(T1s1)
令T1=1,K=1
12
NyquistDiagram
-0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.10RealAxis
sixAyranigam
BodeDiagram
Gm=InfdB(atInfrad/sec),Pm=51.8deg(at0.786rad/sec)
)Bd(edutinga
0
005
0-93
1-1
)ged(esah
-1
10
0
10
1
10
2
10
Frequency(rad/sec)
13
NyquistDiagram
15
10
5
0
-5
-10
-15
-1-0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.1RealAxis
sixAyranigam
-20
BodeDiagram
Gm=InfdB(atInfrad/sec),Pm=51.8deg(at0.786rad/sec)
Frequency(rad/sec)
K(Ts1)
7)G7s2(Tas1),TaT1,其中K可取大于0的任意数
令Ta=2,T1=1,K=1
14
sixAyranigam
)Bd(edutinga
)ged(esah
NyquistDiagram
-0.5
0
86
RealAxis
BodeDiagram
Gm=Inf,Pm=109deg(at0.707rad/sec)
-1
10
Frequency(rad/sec)
0
10
0.5
1
10
15
NyquistDiagram
0
0.2
0.4
0.6
15
10
5
-5
-10
-15
0.8
RealAxis
60
40
0
-75
-80
-85
-90
10
10
10