兰理工 流控实验实验12仿真Word格式文档下载.docx
- 文档编号:22105312
- 上传时间:2023-02-02
- 格式:DOCX
- 页数:28
- 大小:1.30MB
兰理工 流控实验实验12仿真Word格式文档下载.docx
《兰理工 流控实验实验12仿真Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《兰理工 流控实验实验12仿真Word格式文档下载.docx(28页珍藏版)》请在冰豆网上搜索。
deng=[121]
numh=[1]
denh=[11]
[num,den]=feedback(numg,deng,numh,denh,-1);
结果为:
numg=
1
deng=
121
numh=
denh=
11
s+1
---------------------
s^3+3s^2+3s+2
5
程序
numg=[172424]
deng=[110355024]
numc=[105]
denc=[10]
denh=[0.011]
[num1,den1]=series(numc,denc,numg,deng)
[num,den]=feedback(num1,den1,numh,denh,-1)
结果
172424
110355024
numc=
105
denc=
10
0.01001.0000
num1=
01075275360120
den1=
1103550240
00.100010.750077.7500278.6000361.2000120.0000
0.01001.100020.3500110.5000325.2400384.0000120.0000
0.1s^5+10.75s^4+77.75s^3+278.6s^2+361.2s+120
---------------------------------------------------------------
0.01s^6+1.1s^5+20.35s^4+110.5s^3+325.24s^2+384s
+120
6
num=[6.8,61.2,95.2]
den=[1,7.5,22,19.5,0]
[z,p,k]=tf2zp(num,den
6.800061.200095.2000
1.00007.500022.000019.50000
z=
-7
-2
p=
0
-3.0000+2.0000i
-3.0000-2.0000i
-1.5000
k=
6.8000
8;
numg=[1]
deng=[10]
numc=[900]
denc=[19]
[num1,den1]=series(numg,deng,numc,denc)
[num,den]=cloop(num1,den1,-1)
900
19
00900
190
19900
---------------
s^2+9s+900
(8)
deng=[10]
numc=[900]
denc=[19]
[num1,den1]=series(numg,deng,numc,denc)
[num,den]=cloop(num1,den1,-1)
step(num,den)
matlab进行时与分析
p=[13211];
r=roots(p)
r=
-1
-2.2056
-1.0000
0.1028+0.6655i
0.1028-0.6655i
2
P=poly(y)
p=poly(r)
1.00003.00002.00001.00001.0000
4;
num=[2536]
den=[16116]
[r,p,k]=residue(num,den)
2536
16116
-6.0000
-4.0000
3.0000
-3.0000
-2.0000
-1.0000
2
第五5
num=[123]
den=[1331]
123
1331
1.0000
0.0000
2.0000
第6
num=[6.861.295.2]
den=[17.52219.50]
[z,p,k]=tf2zp(num,den)
6.8000
7
num=[3,2,5,4,6]
den=[134272]
[z,p]=tf2zp(num,den)
pzmap(num,den)
32546
134272
0.4019+1.1965i
0.4019-1.1965i
-0.7352+0.8455i
-0.7352-0.8455i
-1.7680+1.2673i
-1.7680-1.2673i
0.4176+1.1130i
0.4176-1.1130i
-0.2991
P31系统动态性能分析
num=[25]
den=[1425]
t=0:
0.01:
step(num,den,t)
gridon
xlabel('
t'
)
ylabel('
c(t)'
25
1425
3
x=0:
2*pi
y=sin(x)
t=0:
plot(x,y)
z=[]
p=[-1-3i-1+3i]
k=[3]
sys=zpk(z,p,k)
10
step(sys,t)
4线性系统的频域分析
num=[1]
den=[10.81]
nyquist(num,den)
title('
NyquistPlotofG(s)=1/(s^2+0.8s+1)'
1.00000.80001.0000
伯德图
num=[0,0,50];
den=[25,2,1]
bode(num,den)
grid
title('
BodeDiagramofG(s)=10(s+1)/s(s+7)'
1010
170
用matlab球稳定裕量
num=[7]
conv([2],[1232])
nyquist(num,den)
[Gm,Pm,wcg,wcp]=margin(num,den)
printsys(num,den)
grid
BodeDiagramofG(s)=7/2(s^3+2s^2+3s+2)'
(2)
num=100*conv([15],[15])
den=conv([11],[119])
G(s)=100*(s+5)^2/(s+1)(s^2+s+9)'
10010002500
12109
100s^2+1000s+2500
-----------------------
s^3+2s^2+10s+9
Gm=
Inf
Pm=
85.4365
wcg=
NaN
wcp=
100.3285
numg=[100]
deng=conv([10],[0.041])
[num,den]=cloop(numg,deng,-1)
printsys(num,den,'
-'
6时间延迟环节
den=[11]
G=tf(numden)
set(G'
Td'
1)
figure
(1)
nyquist(G)
figure
(2)
step(G,t)
最后
num=[100]
den=[0.0410]
bode(num,den)
L(w)'
100
-------------
0.04s^2+s
28.0243
46.9701
deng=[0.0410]
bode(numg,deng,'
step(numg,deng)
[Gm1,Pm1,wcg1,wcp1]=margin(numg,deng)
numc=[0.0251]
denc=[0.011]
[num,den]=series[numgdengnumcdenc]
bode(num,den,'
--'
figure(3)
figure(4)
[Gm2,Pm2,wcg2,wcp2]=margin(num,den)
0.04001.00000
Gm1=
Pm1=
wcg1=
wcp1=
0.02501.0000
002.5000100.0000
0.00040.05001.00000
Gm2=
Pm2=
48.0885
wcg2=
wcp2=
59.6177
[num,den]=series(numg,deng,numc,denc)
grid
grid
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 兰理工 流控实验实验12仿真 理工 实验 12 仿真
![提示](https://static.bdocx.com/images/bang_tan.gif)