《电子信息MATLAB系统仿真与设计》由会员分享，可在线阅读，更多相关《电子信息MATLAB系统仿真与设计（16页珍藏版）》请在金锄头文库上搜索。
1、电子信息系统仿真与设计电子信息系统仿真与设计 课程设计报告课程设计报告设计课题: 油价变化系统的模型 姓 名:学 院: 机电与信息工程学院 专 业: 电子信息科学与技术 班 级: 09 级 2 班 学 号:日 期 2010-2011 第三学期 指导教师: 李光明 张军蕊 山东大学威海分校信息工程学院建模：建模：1背景背景设某一星期的油价为p，其中n表示年份，它与上一星期的油价、油价升值速率以及新增资源所能满足的个体数目之间的动力学方程由如下的差分方程所描述：从此差分方程中可以看出，此油价变化系统为一非线性离散系统。如果设油价初始值、油价升值速率、新增资源所能满足的个体数目，要求建立此油价动态变
2、化系统的系统模型，并分析油价在未来100个星期内之间的变化趋势。2 建立油价变化系统的模型建立油价变化系统的模型(1) Discrete模块库Unit Delay模块：其主要功能是将输入信号延迟一个采样时间，它是离散系统的差分方程描述以及离散系统仿真的基础。在仿真时只要设置延迟模块的初始值便可计算系统输出。(2) Discrete模块库Zero-Order Hold模块：其主要功能是对信号进行零阶保持。使用Simulink对离散系统进行仿真时，单位延迟是Discrete模块库中的Unit Delay模块来完成的。对于油价变化系统模型而言，需要将作为Unit Delay模块的输入以得到，然后按照
3、系统的差分方程来建立人口变化系统的模型。z1Unit Delay ScopeProduct -K-Gain11.05Gain1Constant3 系统参数设置系统参数设置系统模型建立之后，首先需要按照系统的要求设置各个模块的参数，如下所述：(1) 增益模块Gain表示油价升值速率，故取值为1.05。(2) 模块Gain1表示新增资源所能满足的个体数目，故取值为1000000。（3）油价初始值设为10$/L(4) Unit Delay模块参数设置。(5) 仿真时间设置：按照系统仿真的要求，设置系统仿真时间范围为0100。(6) 离散求解器与仿真步长设置：对离散系统进行仿真需要使用离散求解器。实验
4、总结及心得体会实验总结及心得体会MATLAB 是一件很强大的工具，在模拟仿真方面有着不可比拟的优势。不仅可以通过语言脚本可以帮助我们解决很多问题，而且 simulink 也是十分强大的。通过十分直观的方式直接按放各模块，很明显地显示出各种逻辑关系，方便快捷，思路清晰。在实际应用中。Simulink 起到了重要作用。通过对 simulink 的学习，我发现我们所学的课本知识是很重要的，只要通过理解变通，就很容易解决实际问题。但是，有个前提就是你要有着扎实的理论知识。所以，我们千万不能忽略了课本知识的重要性，不要浮躁，理解透彻。Simulink对我来说是很陌生的一个东西，通过几天的摸索，我渐渐摸到
5、了他的奇妙之处，其实不如我们想象那么难，只要没仔细分析好，它会是我们工作学习的一个强力助手。当然，由于时间短暂，我还需要更多时间的学习，才能彻底掌握这个仿真软件。附录附录1.利用 simulink 仿真来实现摄氏温度到华氏温度的转换 3259cfTT2.设系统微分方程为，试建立系统模型并仿真 2) 1 (yyxy3.利用 simulink 仿真，取 A=1, )5cos2513cos91(cos8)(2tttAtx2cosTrigonometric Function2cosTrigonometric Function1cosTrigonometric FunctionSum of Elemen
6、tsScope-K-Gain5-K-Gain4-K-Gain3-K-Gain2-K-Gain11/9GainClock2Clock1Clock4.建立如图 1 所示的仿真模型并进行仿真，改变增益，观察 x-y 图形变化，并用浮动的scope 模块观测各点波形。XY GraphSlider Gain1Sine WaveIntegrator1 sFloating Scope图 1. 题目 4改变增益：继续增大增益：5 有初始状态为 0 的二阶微分方程其中 u(t)是单位阶跃函数，)(24 . 05 . 0tuxxx 试建立系统模型并仿真。6 通过构造 SIMULINK 模型求的结果，其中初值分别为
7、 y1(0)=0, y2(0)=1dtty)cos(当 y1(0)=0 时：当 y1(0)=1 时：7. 分析二阶动态电路的零输入响应 图 2 为典型的二阶动态电路，其零输入响应有过阻尼、临界阻尼和欠阻尼三种情况，已知L=0.5H, C=0.02F, R=1, 2, 3, , 13, 初始值求的零输0)0(,1)0(LciVu)()(tituLc和入响应并画出波形。 （1 用 simulink 的方法，2 用脚本文件的方法） LRC图 2 题目 5 二阶动态电路 （1）用 simulink 的方法-uUnary MinusScope1Scope1 sIntegrator11 sIntegrat
8、or-K-Gain350Gain22Gain1Add（2）用脚本文件的方法定义函数文件funcforex123.mfunction xdot=funcforex123(t,x,flag,R,L,C)xdot=zeros(2,1);xdot(1)=-R/L*x(1)-1/L*x(2)+1/L*f(t);xdot(2)=1/C*x(1);function in=f(t)in=0;脚本文件：L=0.5;C=0.02;for R=1 2 3 4 5 6 7 8 9 10 11 12 13t,x=ode45(funcforex123,0 7,0;1,R,L,C);figure(1);plot(t,x(:
9、,1);hold on;xlabel(timesec);text(0.9,0.07,leftarrowi -L(t);grid;figure(2);plot(t,x(:,2);hold on;xlabel(timesec);text(0.5,0.3,leftarrowu -C(t);grid;end电压图：01234567-0.8-0.6-0.4-0.200.20.40.60.81timesecleftarrowu -C(t)leftarrowu -C(t)leftarrowu -C(t)leftarrowu -C(t)leftarrowu -C(t)leftarrowu -C(t)lefta
10、rrowu -C(t)leftarrowu -C(t)leftarrowu -C(t)leftarrowu -C(t)leftarrowu -C(t)leftarrowu -C(t)leftarrowu -C(t)电流图：01234567-0.2-0.15-0.1-0.0500.050.10.15timesecleftarrowi -L(t)leftarrowi -L(t)leftarrowi -L(t)leftarrowi -L(t)leftarrowi -L(t)leftarrowi -L(t)leftarrowi -L(t)leftarrowi -L(t)leftarrowi -L(t)
11、leftarrowi -L(t)leftarrowi -L(t)leftarrowi -L(t)leftarrowi -L(t)8. 一池中有水 2000，含盐 2 kg，以 6/ 分 的速率向池中注入浓度为 0.5 kg / 3m3m3m的盐水，又以 4 / 分的速率从池中流出混合后的盐水，问欲使池中盐水浓度达到 0.2 3mkg / ，需要多长时间？（1 用 simlink 的方法，2 用脚本文件的方法） 【附加：试画出浓3m度 vs 时间的曲线】ScopeProduct1 sIntegrator2Gain34Gain23Gain1-K-Gain2ConstantClockAdd9. 任意
12、选择一个待仿真的实际问题，建立模型并分析仿真结果，或者 MATLAB Simulink demo 里面一个模块进行分析10. 利用 Simulink 画出以下微分方程组的框图：dx/dt=-x2+y,dy/dt=-x-x*y;x(0)=0,y(0)=0运行结果要求传到工作空间中，并画出相位图（横坐标为 x，纵坐标为 y） 。11. 搭建特定的信号源，建立 SIMULINK 仿真模型、显示仿真结果。netsumSwitch2Switch1SwitchScope=Relational OperatorProduct1Product-1Gain12Gain0Constant42*2Constant32Constant22Constant10ConstantClock4Clock3Clock2Clock1Clock
