自动控制原理图.docx
- 文档编号:24776314
- 上传时间:2023-06-01
- 格式:DOCX
- 页数:26
- 大小:382.65KB
自动控制原理图.docx
《自动控制原理图.docx》由会员分享,可在线阅读,更多相关《自动控制原理图.docx(26页珍藏版)》请在冰豆网上搜索。
自动控制原理图
Wn=0.4;
t=0:
0.1:
80;
num=[Wn^2];
zeta1=0;dem1=[12*zeta1*WnWn^2];
zeta2=0.5;dem2=[12*zeta2*WnWn^2];
zeta3=1.0;dem3=[12*zeta3*WnWn^2];
zeta4=1.5;dem4=[12*zeta4*WnWn^2];
zeta5=2.0;dem5=[12*zeta5*WnWn^2];
[y1,x,t]=impulse(num,[dem1,0],t);
[y2,x,t]=impulse(num,[dem2,0],t);
[y3,x,t]=impulse(num,[dem3,0],t);
[y4,x,t]=impulse(num,[dem4,0],t);
[y5,x,t]=impulse(num,[dem5,0],t);
plot(t,y1,t,y2,t,y3,t,y4,t,y5)
gridon;
t=0:
0.1:
80;
num=[Wn^2];
zeta1=0;dem1=[12*zeta1*WnWn^2];
zeta2=0.5;dem2=[12*zeta2*WnWn^2];
zeta3=1.0;dem3=[12*zeta3*WnWn^2];
zeta4=1.5;dem4=[12*zeta4*WnWn^2];
zeta5=2.0;dem5=[12*zeta5*WnWn^2];
[y1,x,t]=step(num,dem1,t);
[y2,x,t]=step(num,dem2,t);
[y3,x,t]=step(num,dem3,t);
[y4,x,t]=step(num,dem4,t);
[y5,x,t]=step(num,dem5,t);
plot(t,y1,t,y2,t,y3,t,y4,t,y5)
gridon;
Wn=0.4;
t=0:
0.1:
80;
num=[Wn^2];
zeta1=0;dem1=[12*zeta1*WnWn^2];
zeta2=0.5;dem2=[12*zeta2*WnWn^2];
zeta3=1.0;dem3=[12*zeta3*WnWn^2];
zeta4=1.5;dem4=[12*zeta4*WnWn^2];
zeta5=2.0;dem5=[12*zeta5*WnWn^2];
[y1,x,t]=step(num,[dem1,0],t);
[y2,x,t]=step(num,[dem2,0],t);
[y3,x,t]=step(num,[dem3,0],t);
[y4,x,t]=step(num,[dem4,0],t);
[y5,x,t]=step(num,[dem5,0],t);
plot(t,y1,t,y2,t,y3,t,y4,t,y5)
gridon;
num=[1];
den=[110];
pzmap(num,den)
sgrid
num=[1];
den=[1812];
pzmap(num,den)
sgrid
num=[12];
den=[17159];
rlocus(num,den)
sgrid
num=[13];
den=[112701000];
rlocus(num,den)
sgrid
num=[15.5117.5];
den=[13.556.250];
rlocus(num,den)
sgrid
num=[15.5117.5];
den=[13.556.250];
[k,p]=rlocfind(num,den)
rlocus(num,den)
sgrid
num=1;
den=1;
bode(num,den);
grid
num=1;
den=[1,0];
bode(num,den);
grid
num=[1,0];
den=1;
bode(num,den);
grid
num=1;
den=[1,1];
bode(num,den);
grid
num=[1,1];
den=1;
bode(num,den);
grid
num=1;
wn=10;
zeta=[0.1:
0.1:
1.2];
holdon
forz=zeta
den=[(1/wn)^22*z/wn1];
sys=tf(num,den);
bode(sys)
end
title
holdoff
num=[1];
zeta1=[0.1];den1=[0.010.2*zeta11];
zeta2=[0.3];den2=[0.010.2*zeta21];
zeta3=[0.5];den3=[0.010.2*zeta31];
zeta4=[0.7];den4=[0.010.2*zeta41];
zeta5=[0.8];den5=[0.010.2*zeta51];
zeta6=[0.9];den6=[0.010.2*zeta61];
zeta7=[1.0];den7=[0.010.2*zeta71];
zeta8=[1.1];den8=[0.010.2*zeta81];
zeta9=[1.2];den9=[0.010.2*zeta91];
[re1,im1]=nyquist(num,den1);
[re2,im2]=nyquist(num,den2);
[re3,im3]=nyquist(num,den3);
[re4,im4]=nyquist(num,den4);
[re5,im5]=nyquist(num,den5);
[re6,im6]=nyquist(num,den6);
[re7,im7]=nyquist(num,den7);
[re8,im8]=nyquist(num,den8);
[re9,im9]=nyquist(num,den9);
plot(re1,im1,re2,im2,re3,im3,re4,im4,re5,im5,re6,im6,re7,im7,re8,im8,re9,im9);
grid
title('振荡环节奈氏图')
num=[1,1];
den=[1,0];
bode(num,den);
grid
num=[1,1,1];
den=[1,0];
bode(num,den);
grid
sys=tf([250200],[14.2100.8200]);
subplot(2,1,1)
nyquist(sys)
subplot(2,1,2)
bode(sys)
h1=tf([250200],[10.20]);
h2=tf([1],[14100]);
h=h1*h2;
[num,den]=tfdata(h);
[mag,phase,w]=bode(num,den);
subplot(211);
semilogx(w,20*log10(mag));grid;
subplot(212);
semilogx(w,phase);grid;
[gm,pm,wcg,wcp]=margin(mag,phase,w)
G1=tf(35,conv([0.0020.2110],[0.005,1]));
[magphasew]=bode(G1);
magdb=20*log10(mag);
subplot(211);
bode(G1);
margin(mag,phase,w);
G2=tf(35*[0.2310],conv([0.0020.2110],[0.005,1,0]));
[magphasew]=bode(G2);
magdb=20*log10(mag);
subplot(212);
bode(G2,'--');
margin(mag,phase,w);
G1=tf(35,conv([0.0020.2110],[0.005,1]));
G2=tf([0.2,3,10],[10]);
figure;
G1_c=feedback(G1,1);
G2_c=feedback(G2,1);
step(G1_c);
hold;
step(G2_c,'--');
G1=tf(35,[0.000001,0.00305,0.11,1,0]);
G2=tf(35*[0.2,3,10],conv([0.000001,0.00305,0.11,1,0],[1,0]))
bode(G1)
hold
bode(G2,'--')
figure
G1_c=feedback(G1,1)
G2_c=feedback(G2,1)
step(G1_c)
hold
step(G2_c,'--')
G1=tf(100,[0.04,1,0]);
G2=tf(100*[0.025,1],conv([0.04,1,0],[0.01,1]))
bode(G1)
hold
bode(G2,'--')
figure
G1_c=feedback(G1,1)
G2_c=feedback(G2,1)
step(G1_c)
hold
step(G2_c,'--')
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 自动控制 原理图