MATLAB_6 SIMULINK动态系统仿真
- 格式:ppt
- 大小:863.50 KB
- 文档页数:17
Simulink动态系统仿真入门Simulink是基于MA TLAB的图形化仿真设计环境,是MATLAB 提供的进行动态系统建模、仿真和综合分析的集成软件包。
它使用图形化的系统模块对动态系统进行描述,并在此基础上采用MATLAB 的计算引擎对动态系统在时域内进行求解。
它可以处理的系统包括:线性、非线性、离散、连续及混合、单任务、多任务离散事件等。
在MATLAB7.X版本中,可以直接在Simulink环境中运作的工具箱和模型库很多,已经覆盖了航天、航空、通信、控制、信号处理等等诸多领域,涉及内容专业性很强。
1、Simulink系统的启动由于Simulink和MATLAB是高度集成在一起的,因此启动Simulink必须先启动MA TLAB。
在MA TLAB启动Simulink可以通过在命令窗口输入Simulink,或者点击MATLAB工具栏的Simulink 快速启动图标。
启动Simulink后,出现Simulink的主窗口,选择主菜单File中的New\model,即可以打开系统模型编辑器。
下图依次是MATLAB 主窗口、Simulink主窗口和系统模型编辑窗口,图中的箭头表示了操作顺序。
在打开一个新的系统模型文件以后,用户可以从Simulink模块库中选择适合的系统模块或自定义模块来建立系统模型。
我们通过一个简单的例子来分步说明Simulink建模和仿真的能力。
1)在MATLAB 窗口运行Simulink。
打开Simulink模块库浏览器。
2)点击Source子库前的“+”展开库,可以看到各种信源模块。
3)点击新建图标,打开一个空白型的模型窗口。
4)用鼠标选中需要的信源模块,把它拖入新建的空白模型编辑窗口,生成一个正弦波的复制品。
5)同样将信宿库Sinks中的示波器Scope拷贝到模型窗口。
6)利用鼠标完成两个模块的连线操作,完成一个简单的模型。
7)为进行仿真,双击示波器模块,打开示波器显示屏。
目录MATLAB/Simulink 与控制系统仿真课程设计 (1)一、课题设计目的 (3)二、课题设计要求 (3)1.实现单回路控制系统的设计及仿真 (3)2.实现串联控制系统的设计与仿真 (3)3.实现反馈前馈控制系统设计及仿真 (3)三、课题设计内容与步骤 (3)1.实现单回路控制系统的设计及仿真 (3)(1)原始单回路的单位阶跃曲线: (4)(2) P 调节 (4)(3) PI 调节 (5)(4) PID 调节 (5)(5) 总结: (6)2. 实现串联控制系统的设计与仿真 (6)(1).建立开环传递函数 (6)(2).设计串联校正器的滞后环节 (8)(3).设计串联校正器的超前环节 (8)(4).对照校正先后的系统频率响应 (9)(5).系统校正先后的阶跃响应曲线 (10)(6)结论: (12)3.实现反馈前馈控制系统设计及仿真 (12)(1).开环传递函数模型 (12)(2).分析系统的频率响应特性 (13)(3).设计反馈校正器环节 (13)(4).设计前馈校正器环节 (14)(5).设计前向通道 (15)(6).设计前向通道与前馈通道并联连接 (16)(7).设计反馈前馈校正器环节 (16)(8).对照校正先后的系统频率响应 (17)(9).系统校正先后的阶跃响应曲线 (18)(10).总结: (20)四、心得体味 (20)进行PID 控制器设计与应用1.实现单回路控制系统的设计及仿真2.实现串联控制系统的设计与仿真3.实现反馈前馈控制系统设计及仿真P 调节器公式Wc (s) =20.经P 控制后的单回路的单位阶跃曲线:P 控制只改变系统的增益而不影响相位,它对系统的影响主要反映在系统的稳态误差和稳定性上,增大比例系数可提高系统的开环增益、减小系统的稳态误差,从而提高系统的控制精度,但这会降低系统的相对稳定性。
PI 调节器公式Wc (s) =20+3/s.经PI 控制后的单回路的单位阶跃曲线:PI 控制器可以使系统在进入稳态后无稳态误差。
二并联杆数控螺旋面钻头尖刃磨机的机构仿真一、仿真原理图1二并联杆数控螺旋面钻头尖刃磨机床示意图图2 二并联杆数控螺旋面钻头尖刃磨机床刃磨原理图重要假设条件:1、二并联杆数控螺旋面钻头尖刃磨机床是通过两组并联杆(2,a和3,b)保证动平台4只在空间中做水平运动,而没有翻转运动。
每一组并联杆是由空间相互平行的4根杆件组成,由于组内各杆件受力相同,所以将其简化成平面机构如图2。
构件a,b是保证动平台4只做水平运动的辅助平行杆,所以可以假设将机构中杆件a,b省略,而动平台4只做水平移动,没有翻转运动,也就是4相对于地面的夹角θ4恒等于0。
2、直线电机的次子有两个(1和5)但是在加工过程中并不是两者同时运动,所以假设5与导轨固联。
3、假设机床在工作过程中动平台4只受到树直向上的恒力作用,且作用在其中心位置。
基于以上假设机床平面结构示意图如图3。
图3二并联杆数控螺旋面钻头尖刃磨机床简化机构平面结构示意图二、建立仿真方程C2=cos(θ2) S2=sin(θ2) C3=cos(θ3) S3=sin(θ3)一)力方程(分别对各个杆件进行受力分析)对动平台4:受力分析如图4F24x+F43x=m4*Ac4x (1)F24y+F43y=m4*Ac4y (2)F24y*rc4-F43y*rc4=0 (3)图4动平台4的受力分析对并联杆2:受力分析如图5F12x+F24x=-m2*Ac2x (4)F12y+F24y=-m2*Ac2y (5)F12x*rc2*S2+F12y*rc2*C2-F24x*rc2*S2-F24y*rc2*C2=I2*α2 (6)图5并联杆2的受力分析对直线电机滑块1:受力分析如图6Fm+F12x=m1*r1_dot_dot (7)Fy=F12y (8)图6直线电机滑块1的受力分析对并联杆3:受力分析如图7图7并联杆3的受力分析二)闭环矢量运动方程(矢量图如图8)图8 闭环矢量图矢量方程为:R1+R2=R3+R4将上述矢量方程分解为x和y方向,并分别对方程两边对时间t求两次导数得:r1_dot_dot+r2*α2*S2+r2*w2^2*C2=r3*α3*S3+r3*w3^2*C3 (12)r2*α2*C2-r2*w2^2*S2=r3*α3*C3-r3*w3^2*S3 (13)三)质心加速度的矢量方程F13x+F43x=-m3*Ac3x (9)F13y+F43y=-m3*Ac3y (10)F43x*r3*C3+F43y*r3*S3= I3*α3 (11)图9质心加速度的矢量示意图矢量关系:Ac3=Rc3_dot_dotAc4=R3_dot_dot+ Rc4_dot_dotAc2=R3_dot_dot+ R4_dot_dot+ Rc2_dot_dot (_dot_dot 表示对时间求两次导数)将上述三个矢量方程分别分解为x 和y 方向,则它们等效为以下六个方程;Ac3x=-rc3*w3^2*C3-rc3*α3*S3 (14) Ac3y=-rc3*w3^2*S3+rc3*α3*C3 (15) Ac4x=-r3*w3^2*C3-r3*α3*S3 (16) Ac4y=-r3*w3^2*S3+r3*α3*C3 (17) Ac2x=-r3*w3^2*C3-r3*α3*S3-rc2*w2^2*C2-rc2*α2*S2 (18) Ac2y=-r3*w3^2*S3+r3*α3*C3-rc2*w2^2*S2+rc2*α2*C2 (19) 力未知量为:F12x,F12y,F24x,F24y,F43x,F43y,F13x,F13y,Fy,Fm 引入的加速度有:α2,α3,r1_dot_dot ,Ac3x ,Ac3y ,Ac4x ,Ac4y ,Ac2x ,Ac2y三、系统方程的组装将所有19个方程组装成矩阵形式1010000000000200000010100000000002002222222200000020000000000101000000000000400001010000000000004000101000000000000010000000010010000000100000010000000000000010100000000m m rc S rc C rc S rc C I m m m ⋅⋅-⋅-⋅-----300000000101000000003000000333300003000000000000000002233100000000000000002233000000000000000002233010000000000000002233001000000000000000330001000000000000m m r C r S I r S r S r C r C rc S r S rc C r C rc S ⋅⋅-⋅-⋅⋅-⋅⋅⋅⋅⋅-⋅003300001000000000000033000001000000000000330001rc C r S r C ⎛⎫⎪ ⎪⎪⎪⎪ ⎪⎪⎪ ⎪⎪⎪ ⎪ ⎪ ⎪⎪⎪ ⎪⎪⎪ ⎪⎪⎪ ⎪⎪⎪⋅-⎪⋅ ⎪⎪ ⎪⋅-⎝⎭120120240240434301301300020333^2322^2222^2233^2322^222233441F x F y F x F y F x p F y F x F y Fy Fm r w C r w C r w S r w S rc w C Ac x Ac y Ac x Ac y Ac x Ac y r αα••⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⨯= ⎪ ⎪ ⎪⋅⋅-⋅⋅ ⎪ ⎪⋅⋅-⋅⋅ ⎪-⋅⋅ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭33^2322^2233^2333^2333^2333^2333^23r w S rc w S r w S rc w C rc w S r w C r w S ⎛⎫ ⎪⎪⎪ ⎪⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎪ ⎪⎪⎪ ⎪ ⎪ ⎪-⋅⋅ ⎪⋅⋅+⋅⋅ ⎪ ⎪-⋅⋅ ⎪ ⎪⋅⋅ ⎪-⋅⋅ ⎪⎪ ⎪⋅⋅⎝⎭ 四、初始条件的设定假设图3位置就是初始位置。