matlab课程设计报告书 2011
时间:2025-07-08
时间:2025-07-08
《计算机仿真及应用》 课程设计报告书
学号: 08057102, 08057127 班级: 自动化081 姓名 陈婷,万嘉
目 录
一、设计思想
二、设计步骤
三、调试过程
四、结果分析
五、心得体会
六、参考文献
选题一、
考虑如下图所示的电机拖动控制系统模型,该系统有双输入,给定输入R(t)和负载输入M(t)。
1、 编制MATLAB程序推导出该系统的传递函数矩阵。
2、 若常系数增益为:C1=Ka=Km=1,Kr=3,C2=0.8,Kb=1.5,时间常数T1=5,
T2=0.5,绘制该系统的根轨迹、求出闭环零极点,分析系统的稳定性。若R(t)和(要求C1,Ka,Km,Kr,M(t)分别为单位阶跃输入,绘制出该系统的阶跃响应图。C2,Kb, T1,T2所有参数都是可调的)
一.设计思想
题目分析:
系统为双输入单输出系统,采用分开计算,再叠加。
要求参数均为可调,而matlb中不能计算未赋值的函数,那么我们可以把参数设置为可输入变量,运行期间根据要求赋值。
设计思路:
使用append命令连接系统框图。
选择‘参数=input('inputanumber:')’实现参数可调。
采用的方案:
将结构框图每条支路稍作简化,建立各条支路连接关系构造函数,运行得出相应的传递函数。
在得出传递函数的基础上,使用相应的指令求出系统闭环零极点、画出其根轨迹。 通过判断极点是否在左半平面来编程判断其系统是否稳定。
二.设计步骤
(1)将各模块的通路排序编号
(2)使用append命令实现各模块未连接的系统矩阵 (3)指定连接关系
(4)使用connect命令构造整个系统的模型
三.调试过程
出现问题分析及解决办法:
在调试过程出现很多平时不注意且不易寻找的问题,例如输入的逗号和分号在 系统运行时不支持中文格式,这时需要将其全部换成英文格式,此类的程序错误需要细心。
在实现参数可调时初始是将其设为常量,再将其赋值进行系统运行,这样参数可调性差,后用 ‘参数=input('inputanumber:')’实现。 最后是在建立通路连接关系时需要细心。
四.结果分析
源代码:
Syms C1 C2 Ka Kr Km Kb T1 T2 C1=input('inputanumber:') C2=input('inputanumber:') Ka=input('inputanumber:') Kr=input('inputanumber:') Km=input('inputanumber:') Kb=input('inputanumber:') T1=input('inputanumber:') T2=input('inputanumber:') G1=tf(C1,[0 1]); G2=tf(Ka*Kr,[0 1]); G3=tf(Km,[T1 1]); G4=tf(1,[T2 1]); G5=tf(1,[1 0]); G6=tf(-C2,1); G7=tf(-Kb,1); G8=tf(-1,1);
Sys=append(G1,G2,G3,G4,G5,G6,G7,G8)
Q=[1 0 0;2 1 6;3 2 7;4 3 8;5 4 0;6 5 0;7 4 0;8 0 0;]; INPUTS1=1; OUTPUTS=5;
Ga=connect(Sys,Q,INPUTS1,OUTPUTS) INPUTS2=8; OUTPUTS=5;
Gb=connect(Sys,Q,INPUTS2,OUTPUTS) rlocus(Ga)
rlocus(Gb) step(Ga,Gb) Z1=tzero(Ga) Z2=tzero(Gb) p=pole(Ga) p=pole(Gb) if p(real)<0 'yes' else 'no' end
仿真框图:
函数Ga(R(t)为输入)的根轨迹
Ga传函及其表达式:
Ga=connect(Sys,Q,INPUTS1,OUTPUTS)
Transfer function: 1.2 ------------------------ s^3 + 2.2 s^2 + s + 0.96
函数Gb(M(t)为输入)的根轨迹
Gb 传函及其表达式: Gb=connect(Sys,Q,INPUTS2,OUTPUTS)
Transfer function: -2 s - 0.4 ------------------------ s^3 + 2.2 s^2 + s + 0.96
由图可知:闭环极点靠近虚轴,系统快速性不好
函数Ga(上),Gb(下)的阶跃响应曲线
3.仿真结果分析
阶跃响应是衡量系统控制性能的优劣和定义时域性能的指标,所以可由此图得出两个函数的性能指标。
① G(a) 函数
延迟时间:上升时间:
t
t
r
2.5s
1.5s
p
峰值时间:t 5s
超调量: 调节时间:稳态误差:
1.85 1.25
100% 48%
1.25
t
s
25s
1 h( ) 0
e
ss
② G(b)函数: 延迟时间:上升时间:
t
t
r
5.5s
2s
p
峰值时间:t 7.5s
超调量: 调节时间:稳态误差:
1.5 1
100% 50% 1
t
s
27.5s 0
e
ss
以上数据含有误差,其中,延迟时间,上升时间和峰值时间表征系统响应初始阶段的快慢,调节时间表示系统过渡过程持续的时间,是系统的快速性的指标,超调量反应系统响应过程的平稳性,稳态误差反应系统复现输入的最终精度。
选题二、
考虑下列非线性系统,其数学描述为非线性微分方程组:
dx/dt y(t) z(t)
dy/dt x(t) ay(t)
dz/dt b [x(t) c]z(t)
其中a=b=0.2,c=5.7,x(0) y(0) z(0) 0
利用Simulink仿真工具建模,并绘制出各个状态变量的时间响应曲线。
一. 设计思想
对分方程组在零初始条件下进行拉氏变换,并作出各元件的方框图。 按照系统中各变量的传递顺序,依次将各元件的子结构图连接起来。
二. 设计步骤
根据动态结构图用simulink中的元件画出系统的仿真框图。
三. 调试过程
详细调试过程:
建立好模型后,点‘三角形’运行键运行,
在commend window中输入plot(tout,yout),运行,得到仿真的输出波形。
出现问题:
三个输出波形都在一个figure里面,不知道哪个输出对应哪个状态变量。
分析及解决办法:
< …… 此处隐藏:2286字,全部文档内容请下载后查看。喜欢就下载吧 ……上一篇:酒店客房服务质量检查内容和标准
下一篇:整脏治乱实施方案