matlab的应用实例.doc

南京工程学院 期末考查报告 课 题 名 称： 控制系统仿真 学 生 姓 名： 学 号： 指 导 教 师： 所 在 院 系： 专 业 名 称： 2012年 12月 03 日 目录 一、题目及要求…………………………………………………………………….2 二、基本题 ……………………………………………………………………………3 三、弹簧阻尼振荡的simulink仿真……………………………………7 四、心得体会 ………………………………………………………………………12 题目及要求 一、基本题 要求：1、求解上述多项式等于10的解，并将结果输出保存到硬盘中； 2、将多项式以子程序的方式描述出来； 3、求当 时的取值； 4、采用绘图工具将的函数曲线描绘出来，并且在曲线图上标注； 二、找一个具体的对象进行建模，然后运用所学知识（子函数调用和自定义Simulink模块）进行仿真,将m语言编程和模块仿真结合起来。 三、整个论文采用小四字体，行间距为20磅，字体为宋体，图表必须统一标注，标注字体大为小采用五号，公式采用必须采用公式编辑器，且需要对公式标注。文章段落清晰，章节之间需插入分页符。必须标明文献引用。 二、基本题 要求：1、求解上述多项式等于10的解，并将结果输出保存到硬盘中； 2、将多项式以子程序的方式描述出来； 3、求当 时的取值； 4、采用绘图工具将的函数曲线描绘出来，并且在曲线图上标注； 解答过程： 要求1：求解上述多项式等于10的解，并将结果输出保存到硬盘中 分析： 先清屏从内存中清除变量和函数，题目要求多项式等于10，即求解的解，用roots（）实现，结果存入result中，save为 .mat文件。计算结果见图1-1 图1-1 要求2：将多项式以子程序的方式描述出来 主程序： clear %清除所有变量 x=input(请输入x的值) %提示输入 if x==0 disp(y=0) else y=jisuan(x) end 子程序： function result=jisuan(x) result=x^5+x^4+x^3+x^2+x 执行主程序计算x=4的值结果如图1-2所示 图1-2 要求3：求当 时的取值； 方法一 分析：使用循环法 程序如下 for x=0:pi/100:10 disp(x=) disp(x) y=x^5+x^4+x^3+x^2+x disp(...............................................) end 结果如图1-3所示 图1-3 方法二 利用polyval语句 x=[0:pi/100:10]; p=[1 1 1 1 1 ]; polyval(p,x) 图1-4和图1-5中列出了部分值 图1-4 图1-5 要求4：采用绘图工具将的函数曲线描绘出来，并且在曲线图上标注 程序 clear %清除变量 x=[0:pi/100:10] %定义x范围 p=[1 1 1 1 1] y=polyval(p,x) plot(x,y) %plot绘图命令的使用 xlable(x) ylabel(y) %设置x，y坐标轴的名称 title(y=x^5+x^4+x^3+x^2+x) %设置曲线名称 结果如图1-5 图1-5 三、弹簧阻尼振荡的simulink仿真 图2-1 如图2-1，若弹簧振子在x轴上受到粘滞阻尼的作用力，则弹簧振子做阻尼的振动。设弹簧振子受到的阻尼力为： （2-1） 式中( 为阻尼系数，与物体的形状以及周围性质有关。 弹簧振子受到的弹力为，则对弹簧振子，有牛顿第二定律有： (2-2) 整理后得： (2-3) 令，，则有： (2-4) 这就是阻尼振动的振动方程。其解即为弹簧振子在时刻t时的振动位移，一阶导数即为弹簧振子在时刻t时振动速度，其二阶导数即为弹簧振子在时刻t时的加速度。 如图2-1所示，我们假设有粘滞阻尼力时，k=43.8N/m, ,，，=1.5513。由此我们可计算得到：，。初始条件设为，连接好个模块后