- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
控制系统数字仿真课程设计
1 ?课程设计应达到的目的
1、 通过Matlab仿真熟悉课程设计的基本流程;
2、 掌握控制系统的数学建模及传递函数的构造;
3、 掌握控制系统性能的根轨迹分析;
4、 学会分析系统的性能指标;
2 ?课程设计题目及要求
设计要求
1、 进行系统总体设计,画出原理框图。(按给出的形式,自行构造数学模型,构造成 1 个零点,三个极点的三阶系统,主导极点是一对共轭复根)
G (s) =10(s+2)/(s+1)(s 2+2s+6)
2、 构造系统传递函数,利用 MATLAB绘画系统的开环和闭环零极点图; 份别得
至到闭环和开环的零极点图)参考课本P149页例题4-30
clear;
num = [10,20];
den =[1 3 8 6];
pzmap( nu m,de n)
3、 利用MATLAB绘画根轨迹图,分析系统随着根轨迹增益变化的性能。并估算超 调量=16.3%时的K值(计算得到)。参考课本P149页例题4-31
clear
num=[10,20];
den=[1 3 8 6];
sys=tf( nu m,de n);
rlocus(sys)
hold on
jjx(sys);
s=jjx(sys);
[k,Wcg]=imwk(sys)
set(fi ndobj(marker,x),markersize,8,l in ewidth,1.5,Color,k); set(fi ndobj(marker,o),markersize,8,l in ewidth,1.5,Color,k);
fun cti on s=jjx(sys) sys=tf(sys);
num=sys .nu m{1};
den=sys.de n{1};
p=roots(de n);
z=roots( nu m);
n=len gth(p);
m=le ngth(z);
if nm
s=(sum(p)-sum(z))/( n-m)
sd=[];
if n argout1
for i=1: n-m
sd=[sd,s]
end
sysa=zpk([],sd,1);
hold on;
[r,k]=rlocus(sysa);
for i=1: n-m
plot(real(r(i,:)),imag(r(i,:)),k:);
end
end
else
disp;
s=[];
end
fun cti on [k,wcg]=imwk(sys) sys=tf(sys)
num=sys .nu m{1}
den=sys.de n{1};
asys=allmargi n(sys);
wcg=asys.GMFreque ncy;
k=asys. GainMargin;
Rc ct Lc cus-1Real Axis3一 xy
Rc ct Lc cus
-1
Real Axis
3一 xy
4、 编写M文件作出单位阶跃输入下的系统响应、分析系统单位阶跃响应的性能指 标;(以阶跃信号作为输入,利用函数得到系统的阶跃响应,并求出对应的性能指标, 上升时间,峰值时间,超调量,调节时间, 稳态误差(计算得到))参考自控课本
P105页例题3-22
clear
num=[10,20];
den=co nv([1 0],[1 2 6]);
sys=tf( nu m,de n);
sys = feedback(sys,1);
[y,t] = step(sys);
Plot(t,y)
ytr=fi nd(y=1)
rise_time=t(ytr(1))
[ymax,tp]=max(y);
peak_time=t(tp)
max_overshoot=ymax-1
s=le ngth(t);
while y(s)0.98 y(s)1.02
s=s-1;
end
setti ng_time=t(s+1)
plot(t,y,k,t,o nes(le ngth(t),1),k-.)
axis([0 10 0 2.2])
title(Plot of Unit-Step Response Curves,Position,[5 2.22],FontSize,8)
x-abe-(-Time(sec)-zposmorrl9.86」5L-F0nfsizQ.oo) y-abe-(-Response-zposmorr」o.25 二-Fonfsize-oo)
riuE—time u
F燼g
peakktime H
口 ■舟一
maxIOTerwhoot =
p言g
mEr+r+inmltimE H
魚
华侨大学厦门工学院
课程设计报告
题 目
专业、班级:
学生姓名:
学 号:
指导教师:
分 数:
2011年 月 日
课程设计报告要求:
一、 目录(小四,宋体)
二、 内容
文档评论(0)