[工学]MATLAB大作业09机械5班 韦靖 0901901012.doc

梧州学院 MATLAB与系统仿真 专 业：机械设计制造及其自动化 班 级：09机械5班 学 号：0901901012 姓 名：韦 靖 邮 箱：wei_jing1988@163.com 指导老师: 姚 金 2012年6月27日 一、春天是鲜花的季节，水仙花就是其中最迷人的代表，数学上有个水仙花数，它是这样定义的：“水仙花数”是指一个三位数，它的各位数字的立方和等于其本身，比如：153=1^3+5^3+3^3。要求编m文件输出所有存在的水仙花数。（本题10分） 解： 题意中“水仙花数”是指一个三位数说明100—999之间的数，在MATLAB编程M文件，运行可得。 程序如下： for i=100:999; a=floor(i/100); b=floor(i/10)-a*10; c=i-(a*100+b*10); if i==a^3+b^3+c^3; i end end 运行结果： 本程序源程序附件名为： 二、带饱和输出特性PID控制系统的被控对象为=，其PID控制器为===5+。以与为前向通道的单位负反馈系统。试对于PID控制器分别采用以下3种办法：PID控制3个分量叠加、传递函数方框图与PID控制子系统绘制其Simulink仿真模型图，并对其进行阶跃仿真。然后对PID控制子系统进行封装。（本题20分） 解： 1.分析：根据题意，已知PID控制器为=5+，即=5，=0.048，=0。PID的三个分量。仿真时可采用下列处理措施：①PID控制3个分量叠加法②传递函数方框图法③PID 控制子系统及其Simulink仿真模型图法。 2. PID的三个分量叠加法则。 本图为 PID控制三个分量叠加的Simulink模型syb1.mdl 3.将PID的三个分量叠加作用直接绘制早系统中Simulink模型syb1.mdl如图1-1所示。图中的阶跃信号设为“1”，即单位阶跃作用，三信号叠加模块“Add”，混路器模块Mux,饱和特性为“+”1.1或“-1.1”，时间延迟“60”s，仿真时间为“2000”。模块的设置参数如图所示 STEP设置参数如下 Integrator设置参数如下 Add设置参数如下 Saturation参数设置如下 Transfer Fcn设置参数如下 Transport Delay设置参数 Gain参数设置 Gain1参数设置 Gain2参数设置 Mux参数设置 Scope运行结果如下： 4.传递函数的方框图法： 将PID的三个分量叠加作用（求和）通分后求得一个传递函数，即 ===5+=，将传递函数绘制在系统的Simulink模型syb2.mdl如图1-7所示。对模型syb2.mdl进行仿真，结果同样得到1-6所示的系统的单位阶跃给定响应曲线。 图为 PID控制传递函数方框的Simulink模型sybl2.md 根据PID控制传递函数方框的Simulink模型syb2.mdl进行仿真得到如图1-8所示系统单位阶跃给定响应曲线。 运行结果如下： 5.PID系统控制法。 ①定义子系统。 先到Simulink库浏览器的“Commonly Used Blocks”模块组中把标准功能模块“Subsystem”复制到空白模型设计区域（“untitled”）内，还复制其他模块，将其连成系统如图1-9所示，模型命名为“syb32.mdl”后存盘。 PID控制子系统的Simulink模型syb32.mdl 用鼠标左键双击“Subsystem”，将会弹出“syb32/Subsystem”的设计区。 在“syb32/Subsystem”的设计区域中进行如图1-10所示的PID控制器的设计（复制各个标准功能模块，连线构成系统模型），或者将设计好的PID控制器模型复制到其中，并在输入与输出端假如“In”和“Out”两个功能模块（见图1-10）。模型设计完后，要注意存盘。至此，子系统创建完成。 PID控制子系统Subsystem的Simulinkmsyb32/Subsystem 6.用PID控制子系统封装的步骤。 在syb32.mdl中的Subsystem模块更换为一个没有涉及的孔Subsystem模块。 在“syb32/Subsystem”的设计区域中进行如图1-10所示的PID控制器的设计，输入数据并存盘。 再在模型syb32.mdl中，选中“Subsystem”，选择主菜单【Edit】，执行其下的【Mask Subsystem...】或【Edit Mask】菜单项命令，打开【Mask edit】即子系统的封装编辑器。封装编辑器中有四张标签页，所谓完成封装操作，就是按要求填写接下来的四张标签页， 填写第1张标签页“Icon”。此时仅在“Drawing comm