实验一 控制系统建模与分析

c) 程序如下：

G=tf([100,500],conv([1,6,-16],[1,20]));

G_close=feedback(G,1);

step(G_close)

运行结果：

6、已知系统传递函数为：

G(S)?53S2?4S?6

（1）编写M文件求取系统的性能指标（超调量，上升时间，调节时间） 先建立一个stepchar.m的文件

程序为：

function [pos,tr,ts,tp]=stepchar(g0,delta)

[y,t]=step(g0); [mp,ind]=max(y); dimt=length(t); yss=y(dimt);

pos=100*(mp-yss)/yss;

tp=t(ind); for i=1:dimt

if y(i)>=0.8

tr=t(i); break; end

end;

for i=1:length(y)

if y(i)<=(1-delta)*yss || y(i)>=(1+delta)*yss ts=t(i); end

end

再建立一个M文件，程序：

sys=tf(5,[3,4,6]);

[pos,tr,ts,tp]= stepchar(sys,0.02)

两个m文件存在同一个文件夹中。 运行结果： pos =

18.4228 tr =

1.5895 ts =

5.8811 tp =

2.5432

（2）利用LTI Viewer对系统进行分析，在一个分析窗口中画出系统的单位阶跃响应、单位脉冲响应、伯德图、伯德幅频图、相频图、奈奎斯特曲线、尼克尔斯

曲线。并从单位阶跃响应曲线图上读取最大超调量，调节时间，上升时间，与（1）中的结果进行对比。 程序如下：

sys=tf([5],[3,4,6]);

ltiview

运行结果：

超调量为18.6%；调节时间为5.85s；上升时间为1.13s。结果与（1）的结果基本一样。