山东科技大学考试复习资料-系统仿真技术与应用(考试大纲+知识点)

1、一共考5个题一、简答题（3个题必考两个共20分）P87步长的选择与应用不能无限制减小步长的原因减慢计算速度增加累积误差应采取的措施选择适当的步长改进近似算法精度采用变步长方法P85描述被积函数的三种方法（给定被积函数为y=1/sqrt(2*pi)*exp(-x.2/2)）匿名函数例如f=(x)1/sqrt(2*pi)*exp(-x.2/2),其中，括号内为自变量列表，后面的表达式求出的是被积函数的值，赋给变量f。这种方法的速度明显高于其他两种方法。M-函数function y=myerrf(x),y=1/sqrt(2*pi)*exp(-x.2/2)该函数的优点在于可以返回多个变量，而另外两种方

2、法则不行，但该方式需要建立一个文件。inline函数f=inline(1/sqrt(2*pi)*exp(-x.2/2),x)其中，函数首先给出了求值表达式，然后给出自变量。该方法的速度是3种方法中最慢的，尽量不要使用。模型转换函数residue:传递函数模型与部分分式模型互换ss2tf:状态空间模型转换为传递函数模型tf2ss: 传递函数模型转换为状态空间模型ss2zp:状态空间模型转换为零极点增益模型zp2ss:零极点增益模型转换为状态空间模型tf2zp:传递函数模型转换为零极点增益模型zp2tf:零极点增益模型转换为传递函数模型二、大题（必考共40分）Lorenz方程的Simulink仿真

3、，其中方程为：仿真模型绘制图如下：（注意对应的常数（箭头指向常数）-8/3、10、-10、28老师出题的时候会改，找出对应的位置修改就行）用MATLAB语言对Van der Pol方程：，即对方程求解(其中=1)（假如老师把改为2，即在下面程序中将mu=1改为mu=2即可）解：法一（隐函数法）：f=(t,x，mu)x(2);-mu*(x(1)2-1)*x(2)-x(1);h_opt=odeset;x0=-0.2;-0.7;t_final=20;mu=1;t,y=ode45(f,0,t_final,x0,h_opt,mu);plot(t,y),figure;plot(y(:,1),y(:,2)附

4、法二（M-函数法）：在子M-文件中写入：function f=van(t,x,mu) f=x(2);-mu*(x(1)2-1)*x(2)-x(1);在主M-文件中写入：clear all clc close all h_opt=odeset;x0=-0.2;-0.7;t_final=20;mu=1;t,y=ode45(van,0,t_final,x0,h_opt,mu);plot(t,y),figure;plot(y(:,1),y(:,2)三、补充题电路图如下，R=1.4欧，L=2亨，C=0.32法，初始状态：电感电流为零，电容电压为0.5V，t=0时刻接入1V的电压，求0t15s时，i(t)

5、，vo(t)的值，并且画出电流与电容电压的关系曲线。解：电路原理公式子M文件中键入主M文件中键入用两种调用子程序的方式求s=1+2+3+m，当1000s3000时的所有m值解：第一种：子M文件中键入主M文件中键入：第二种：子M文件中键入主M文件中键入 在M文件中键入如下程序此处矩阵可以自己写，只要包含a就行any()函数的用法： 用法例解：B=any(A)，如果A是向量，如果向量里有非0的数，则返回1（true），如果A是矩阵，则把矩阵的列当做向量来处理，函数返回每个列向量的逻辑值；all()函数的用法：用法例解：B = all(A)，如果A是一个向量，如果所有的元素都是非零的，则返回1，如果

6、有一个元素为零，则返回0。如果A是一个矩阵，则返回一个行向量，用于检测每一列是否全为非零元素，如果某一列中有一个元素为零，则返回0，如果某一列中全为非零元素，则返回1，由此构成一个行向量。试探性语句结构：try，语句段1，catch，语句段2，end这种结构在实际编程中很实用，例如将一段不保险但速度快的算法放到try段落中，而将一个保险但速度极慢的程序放到catch段落中，这样就能保证原始问题的求解更加可靠，且能使程序高速执行。系统模型分为物理模型和数学模型，系统仿真是根据被研究的真实系统的数学模型研究系统性能的一个学科，现在尤指利用计算机去研究数学模型行为的方法。一维插值问题可用函数interp1（）