自控原理_.doc

附录
二阶系统传递函数程序
num=[10];
den=[0.1110];
t=0:0.1:20
step(num,den,t)
[y,x,t]=step(num,den,t)
maxy=max(y)
yss=y(length(t))
pos=100*(maxy-yss)/yss
fori=1:201
ify(i)==maxy
n=i;
end
end
tp=(n-1)*0.1
y1=1.05*yss
y2=0.95*yss
i=201
whilei>0
i=i-1
ify(i)>=y1y(i)<=y2
m=1
break
end
end
ts=(m-1)*0.1
title“单位阶跃响应”
gridon
程序1
二阶系统传递函数程序
num=[10];
den=[0.110];
G=tf(num,den);
[nc,dc]=cloop(num,den,-1)
sys=tf(nc,dc)
figure(1)
step(nc,dc);
gridon；
figure(2)
rlocus(G);
gridon；
程序2
二阶系统传递函数程序
num=[10];den=[0.1110];sys=tf(num,den);ltiview(sys);
程序3

三阶系统传递函数程序
num=[10];
den=[0.0510.61110];
t=0:0.1:20
step(num,den,t)
[y,x,t]=step(num,den,t)
maxy=max(y)
yss=y(length(t))
pos=100*(maxy-yss)/yss
fori=1:201
ify(i)==maxy
n=i;
end
end
tp=(n-1)*0.1
y1=1.05*yss
y2=0.95*yss
i=201
whilei>0
i=i-1
ify(i)>=y1y(i)<=y2
m=1
break
end
end
ts=(m-1)*0.1
title“单位阶跃响应”
gridon
程序4

三阶系统传递函数程序
num=[10];
den=[0.0510.6110];
G=tf(num,den);
[nc,dc]=cloop(num,den,-1)
sys=tf(nc,dc)
figure(1)
step(nc,dc);
gridon；
figure(2)
rlocus(G)
gridon；

程序5

三阶系统传递函数程序
num=[10];den=[0.0510.61110];sys=tf(num,den);ltiview(sys);
程序6

没增加零极点原程序
num=[10];
den=conv([10],[0.11]);
G=tf(num,den);
[nc,dc]=cloop(num,den,-1)
sys=tf(nc,dc)
figure(1)
step(nc,dc);
figure(2)
margin(G);
figure(3)
rlocus(G);
figure(4)
nyquist(feedback(G,1));
程序7

增加一个极点程序
num=[10];
den=conv([10],conv([0.11],[14]));
G=tf(num,den);
[nc,dc]=cloop(num,den,-1)
sys=tf(nc,dc)
figure(1)
step(nc,dc);
figure(2)
margin(G);
figure(3)
rlocus(G);
figure(4)
nyquist(feedback(G,1));
程序8
增加二个极点程序
num=[10];
den=conv([10],conv([0.11],conv([12-4i],[12+4i])));
G=tf(num,den);
[nc,dc]=cloop(num,den,-1)
sys=tf(nc,dc)
figure(1)
step(nc,dc);
figure(2)
margin(G);
figure(3)
rlocus(G);
figure(4)
nyquist(feedback(G,1));
程序9

增加一个零点程序
num=conv([010],[14]);
den=conv([10],[0.11]);
G=tf(num,den);
[nc,dc]=cloop(num,den,-1)
sys=tf(nc,dc)
figure(1)
step(nc,dc);
figure(2)
margin(G);
figure(3)
rlocus(G);
figure(4)
nyquist(feedback(G,1));
程序10

增加二个零点程序
num=conv([12+4i],[12-4i]);
den=conv([0.10],[0.11]);
G