自动控制MATLAB仿真实验三四.docx
《自动控制MATLAB仿真实验三四.docx》由会员分享,可在线阅读,更多相关《自动控制MATLAB仿真实验三四.docx(14页珍藏版)》请在冰点文库上搜索。
![自动控制MATLAB仿真实验三四.docx](https://file1.bingdoc.com/fileroot1/2023-7/16/268fe9d3-c40f-443d-bae4-e2a5083b4d92/268fe9d3-c40f-443d-bae4-e2a5083b4d921.gif)
自动控制MATLAB仿真实验三四
num=[50]
den=[1,3,-10]
nyquist(num,den)
axis([-62-20])
title('Nyquistͼ')
w=[0,logspace(-2,2,200)]
wn=0.7
tou=[0.1,0.4,1,1.6,2]
forj=1:
5
sys=tf([wn*wn],[1,2*tou(j)*wn,wn*wn])
bode(sys,w)
holdon
end
gtext('tou=0.1')
gtext('tou=0.4')
gtext('tou=1.0')
gtext('tou=1.6')
gtext('tou=2.0')
num=5*[0.0167,1]
den=conv(conv([1,0],[0.03,1]),conv([0.0025,1],[0.001,1]))
sys=tf(num,den)
w=logspace(0,4,50)
bode(sys,w)
grid
[Gm,Pm,Wg,Wc]=margin(sys)
num=[100]
den=[1,15,50,0]
nyquist(num,den)
title('Nyquistͼ')
num=[800]
den=[1,15,50,0]
nyquist(num,den)
title('Nyquistͼ')
num=[2000]
den=[1,15,50,0]
nyquist(num,den)
title('Nyquistͼ')
num=100*[0.5,1]
den=conv(conv([1,0],[1,1]),conv([0.1,1],[0.05,1]))
sys=tf(num,den)
w=logspace(0,4,50)
bode(sys,w)
grid
[Gm,Pm,Wg,Wc]=margin(sys)
num=2
num=
2
>>den=conv([0.1,1,0],[0.3,1])
den=
0.03000.40001.00000
>>G=tf(num,den)
Transferfunction:
2
----------------------
0.03s^3+0.4s^2+s
>>sisotool
num=1.5
num=
1.5000
>>den=[1,14,40.02]
den=
1.000014.000040.0200
>>G=tf(num,den)
Transferfunction:
1.5
------------------
s^2+14s+40.02
>>sisotool