《计算机仿真技术》试题(含完整答案)
2009年计算机仿真技术试卷
中南大学考试试卷2009 -- 2010 学年 1 学期期末考试试题 时间110分钟计算机仿真技术 课程 32 学时 2 学分 考试形式： 开 卷 专业年级： 自动化06级 总分100分，占总评成绩70% 注：请将答案写在答卷上一 填写下列运算结果（本题12分，每空2分）已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+=542031321i A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+=581051301i B ,则 (1) A (1:3,3)＝ 305⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦(2) B (2:3 ) ＝ []11(3) A ’＝112234305i ⎡⎤⎢⎥-⎢⎥⎢⎥⎣⎦(4) B.’＝ 111058305i ⎡⎤⎢⎥⎢⎥⎢⎣⎦+⎥ (5) A.*B ＝ 10931155023225i i +⎡⎤⎢⎥+⎢⎥⎢⎥⎣⎦(6) A*B ＝ 63418415531160312i i i i i +⎡⎤⎢⎥+++⎢⎥⎢⎥+⎣⎦二 判断正误（本题20分，每小题2分）(1) 计算机仿真的三要素为系统、模型和计算机。
（ √ ）(2) 在MATLAB 中，变量名区分字母的大小写，但是MATLAB 提供的标准函数名以及命令名必须用小写字母。
（ × ）(3) MATLAB 中运算符号前带有小点表示执行元素对元素的运算。
（ √ ）(4) 在MATLAB 中，2ⅹ2的矩阵A=[1 2，3 4]是一种正确的输入方式。
（ × ）(5) clc 命令不仅用于清除命令窗口显示的内容同时也删除MATLAB 工作空间中的变量。
（ × ）(6) MATLAB 提供了3种逻辑操作符，其中~表示非。
（ √ ）(7) subplot(211)把屏幕分成两个部分，并把曲线放在左半部分。
（ × ）(8) 数值积分法中，计算步长越小，总误差越小。
（ × ）(9) tf2ss 函数只能将传递函数模型转换成为可控标准型状态方程。
（ √ ）(10) 函数文件中函数名称和文件名可以不同。
计算机仿真试题和答案
控制系统计算机辅助设计综合实验指导实验名称：连续与离散系统校正实验，系统可控性与可观性实验，系统的simulink仿真实验陈茜编实验人：苏建聪学号：200830810122班级：08电气工程及其自动化1班信息工程系实验任务书1. 有一个单位负反馈控制系统，如果控制对象的传递函数为设计要求：① 相角裕度≥45°；② 当系统的输入信号是单位斜坡信号时，稳态误差ess ≤0.04。
③ 要求绘制出校正后系统和未校正系统的Bode 图及其闭环系统的单位阶跃响应曲线，并进行对比。
2. 有一个单位负反馈控制系统，如果控制对象的传递函数为：试设计一个串联滞后校正装置。
设计要求：①相角裕度≥45°；② 当系统的输入信号是单位斜坡信号时，稳态误差ess ≤0.04。
③ 要求绘制出校正后系统和未校正系统的Bode 图及其闭环系统的单位阶跃响应曲线，并进行对比。
3. 有一个单位负反馈控制系统，如果控制对象的传递函数为 ()()4+=s s k s G p试设计一个串联超前滞后校正装置，设计要求： ①相角裕度≥45°；② 当系统的输入信号是单位斜坡信号时，稳态误差ess ≤0.04。
③ 要求绘制出校正后系统和未校正系统的Bode 图及其闭环系统的单位阶跃响应曲线，并进行对比。
4. 系统结构图如图所示，其中，采样周期Ts=0.01s ，被控对象()()110+=s s s G ，()s G h 为零阶保持器。
用W 变换法设计一超前校正装置D(z)，使系统相位裕度γ≥50°，校验设计后系统的性能指标。
5. 系统结构图如图所示，其中，采样周期Ts=0.01s ，被控对象1)s(0.2s k )(+=s G ο，()s G h 为零阶保持器。
用对数频率法设计D(z)，使系统开环增益k ≥30(1/s)，截止频率ωc ≥15(1/s),相位裕度γ≥50 °1使 s 11se -1(s)-Tsh T G +≈=，求出未校正系统的开环系统的开环传递函数(s)(s)G G (s)0s =G ,的传递函数模型参数。
计算机仿真习题及答案
计算机仿真试题1．编写一个函数，使其能够产生如下的分段函数：错误！未找到引用源。
并调用此函数，绘制x=[0,+2]范围内的f(x)*f(x+2) 。
（10分）function y=f(x)if x<=2y=0.5*x;else if x>6y=0.5;else y=1.5-0.25*x;endendx=0:0.05:2;y= f(x)’*f(x+2));plot(x,y)图 1-12．已知4阶龙格-库塔算法如下：试利用该算法求解以下微分方程：（15分）本题可以调用MATLAB函数中龙格-库塔算法函数ode45，首先编写m文件：function dy=func(x,y)dy=-y+1;end再在主窗口调用此文件：[x,y]=ode45('func',[0,5],0)%这里的[0,5]为任取区间，表示方程在此范围的解。
运行结果如下:x =0.00010.00010.00020.00020.00050.00070.00100.00120.00250.00370.00500.00620.01250.01880.02510.0313 0.06270.09410.12550.15690.28190.40690.53190.65690.78190.90691.03191.15691.28191.40691.53191.65691.78191.90692.03192.15692.28192.40692.53192.65692.78192.90693.03193.15693.28193.40693.53193.65693.78193.90694.03194.15694.28194.40694.53194.65694.74274.82854.91425.0000y =0.00010.00010.00020.00020.00050.00070.00100.00120.0025 0.0037 0.0050 0.0062 0.0124 0.0186 0.0248 0.0309 0.0608 0.0898 0.1180 0.1452 0.2457 0.33430.41250.48160.54250.59630.64370.68550.72250.75510.78390.80930.83170.85150.86890.88430.89790.90990.92050.92980.93810.94540.95180.95740.96240.96690.97080.97420.97720.97990.98230.98430.98620.98780.98920.99050.99130.99200.99270.9933为只管起见，我们使用函数命令画出x-y（plot(x,y)）的关系如下图：图1-23．用matlab语言求下列系统的状态方程、传递函数、零极点增益、和部分分式形式的模型参数，并分别写出其相应的数学模型表达式：（15分）（1）G（s）=324327242410355024s s ss s s s+++++++（2）.X=2.25 -5 -1.25 -0.542.25 -4.25 -1.25 -0.2520.25 -0.5 -1.25 -121.25 -1.75 -0.25 -0.75 0X⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥+⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦uy= [0 2 0 2] X解：（1）a)求对应状态方程参数：num=[1 07 24 24]; den=[1 10 35 50 24]; [A,B,C,D]=tf2ss(num,den) 运行结果：A =-10 -35 -50 -241 0 0 00 1 0 00 0 1 0B =1C =1 7 24 24D =故，状态方程为：.X = x+ uY=[1 7 24 24]xb)求对应零极点增益模型参数：num=[1 07 24 24]; den=[1 10 35 50 24]; [Z,P,K]=tf2zp(num,den) 运行结果如下： Z =-2.7306 + 2.8531i -2.7306 - 2.8531i -1.5388P = -4.0000 -3.0000 -2.0000 -1.0000K = 1故变换后的零极点模型为： G(s)=c)求对应部分分式型：num=[1 07 24 24]; den=[1 10 35 50 24]; [R,P,H]=residue(num,den) 运行结果如下： R =4.0000 -6.0000 2.0000 1.0000P =-4.0000 -3.0000 -2.0000 -1.0000H = []故变换后的部分分式模型为：11223644)(+++++-+=s s s s s G（2）由题给条件，知：A=[2.25 -5 -1.25 -0.5; 2.25 -4.25 -1.25 -0.25;0.25 -0.5 -1.25 -1;1.25 -1.75-10 -35 -50 -24 1 0 0 0 0 1 0 0 0 0 1 010 0 0-0.25 -0.75] B=[4;2;2;0] C=[0 2 0 2],D=0 a)求传递函数矩阵： [num,den]=ss2tf(A,B,C,D) 运行结果为： num =0 4.0000 14.0000 22.0000 15.0000 den =1.0000 4.0000 6.2500 5.25002.2500 故，所对应传递函数模型为：25.225.525.641522144)(23423+++++++=s s s s s s s s Gb)求零极点模型：num=[0 4 14 22 15];en=[1 4 6.25 5.25 2.25]; [Z,P,K]=tf2zp(num,den) 运行结果为： Z =-1.0000 + 1.2247i -1.0000 - 1.2247i -1.5000 P =-1.5000 -1.5000 -0.5000 + 0.8660i -0.5000 - 0.8660iK =4.0000故，零极点模型为：)866.05.0()5.1()2247.11)(5.1(4)(2i s s i s s s G ±++±++=c)求对应部分分式模型： [R,P,H]=residue(num,den) 运行结果为： R =4.0000 -0.0000-0.0000 - 2.3094i -0.0000 + 2.3094iP =-1.5000 -1.5000 -0.5000 + 0.8660i -0.5000 - 0.8660iH = []故变换后的部分分式模型为：i s ii s i s s G 866.05.03094.2866.05.03094.25.14)(+++-+-++=4．已知一单位反馈系统开环传递函数为：，试绘制系统Nyquist图，判断闭环系统的稳定性，并求其单位阶跃响应。
计算机仿真试题A 答案
计算机仿真试题A 参考答案四、程序设计与算法题1. (只要是完成题目要求的程序都可以得分)(a)I=imread('flowers.jpg');J=imadd(I,50);K=immultiply(J,0.7);subplot(1,2,1), subimage(I);subplot(1,2,2),subimage(K);(b)I=imread(‘flowers.jpg’);J=imresize(I,1.2);K=imrotate(J,-45,’bilinear’,’crop’);subplot(131),subimage(I);subplot(132),subimage(J);subplot(133),subimage(K);(c)I=imread('flowers.jpg');[x,y,z]=cylinder;warp(x,y,z,I);J=rgb2gray(I);figure,imshow(J);2.主要算法流程如下：程序代码：%Example of PCA analysisclear all, close all;N=1000;fs=500;w=(1:N)*2*pi/fs;t=1:N;x=0.75*sin(w*5);y=sawtooth(w*7,0.5);D(1,:)=.5*y+.5*x+.1*rand(1,N);D(2,:)=.2*y+.7*x+.1*rand(1,N);D(3,:)=.7*y+.2*x+.1*rand(1,N);D(4,:)=-.6*y+-.24*x+.2*rand(1,N);D(5,:)=.6*rand(1,N);plot(t,D(1,:)+0,t,D(2,:)+2,t,D(3,:)+4,t,D(4,:)+6,t,D(5,:)+8);figure;for i=1:5D(i,:)=D(i,:)-mean(D(i,:));end[U,S,pc]=svd(D,0);pc=pc(:,1:5);for i=1:5pc(:,i)=pc(:,i)*sqrt(eigen(i));endfigure;plot(t,pc(:,1)-2,t,pc(:,2)+2);3(a)(注意这里n=4), As x =得：x As 1-= （b ）流程图类似于上题，但注意这里已经有了混合信号，所以无需产生；并且是用ICA 进行分析。
智慧树答案计算机仿真技术知到课后答案章节测试2022年
第一章1.计算机仿真系统的三要素（）答案:计算机、模型、系统2.系统研究一般有哪些方法（）答案:仿真实验法;理论解析法;直接实验法3.哪些是以硬件为基础的仿真软件( )答案:Proteus4.欠实时仿真，仿真时钟比实际时钟慢。
例如烟火爆炸的仿真分析。
答案:对5.计算机仿真算法的误差主要有（）答案:截断误差;舍入误差第二章1.MATLAB软件中，下面哪个命令是用来显示当前文件夹中文件名字的（）答案:what2.下面这些哪个是属于元素群运算（点运算）（）答案:exp(x)3.哪个是以硬件为基础的仿真软件( )答案:Proteus4.MATLAB是（矩阵实验室），其基本变量是矩阵，其数值元素的格式只有双精度（double）这种形式存在于MATLAB的workpace空间中。
答案:错5.在MATLAB中，运行以下代码，b的结果是（）a=[1,0;1,0];b=all(a)答案:[1,0]第三章1.MATLAB中代码文件中边表明文件为函数文件的关键字是（）答案:function2.在MATLAB的命令窗口中，运行以下代码，b的结果是（）a=[1,0,1,0];forn=ax=n(1)&n(2);end答案:出现错误提示3.在MATLAB的fig窗口中，可以直接打开以下（）格式的图形文件答案:fig4.在MATLAB中，对于m函数文件，下面的描述不正确的是（）答案:m函数文件是不能够被simulink来调用的5.全局变量是共有变量，在函数文件中不能够直接访问，需要设置一个输入参数才能访问。
答案:错第四章1.以下说法不正确的是（）答案:不同幂次的两个多项式不能进行乘法运算2.当实验或测试所获得的样本数据有误差时，适合用来估算数据的方法是（）答案:曲线拟合3.若在MATLAB中已有语句a=sym(1); b=1; x=sym('b')，则以下叙述正确的是（）答案:a是符号常量，x是符号变量4.若使用命令taylor(f,x,1,'Order',6)对f进行泰勒展开，则展开式的最高阶为（）答案:55.在命令行窗口输入下列命令，则命令执行后的输出结果是（）>>f=sym(1);>> eval(int(f,1,4))答案:3第五章1.启动Simulink的方法有（）答案:在命令行窗口中输入simulink命令;在“主页”选项卡中单击“文件”命令组中的“新建”命令按钮;在“主页”选项卡中单击SIMULINK命令组中的“Simulink”命令按钮2.在一个模型窗口上按住一个模块并同时按Shift键移动到另一个模型窗口，则（）答案:在两个模型窗口都有这个模块3.为子系统定制参数设置对话框和图标，使子系统本身有一个独立的操作界面，这种操作称为子系统的（）答案:封装4.已知仿真模型如图所示，各模块参数均采用默认设置，则示波器的输出波形是（）答案:5.以下关于S函数的描述中，正确的有（）答案:S函数有现成的模板程序;S函数模块能够被封装;利用S函数可以对Simulink模块库进行扩充第六章1.下列对象中，哪一个不属于MATLAB控制系统工具箱中规定的LTI对象？（）答案:极点留数对象2.tf2zpk函数的作用是？（）答案:零极点增益模型转换为状态空间模型3.在MATLAB命令窗口运行如下代码，返回结果为：（）>>sys = drss(3,4)；isct(sys)答案:14.求系统sys的阶跃相应，返回系统的输出和状态，step函数的调用格式为：（）答案:[y,t,x]=step(sys)5.给定开环系统模型作波特图，频率向量w自动给出，作波特图，并在图上标注幅值裕度Gm和对应的频率wg，相位裕度Pm和对应的频率wp，所用函数及调用格式为：（）答案:margin(sys)第七章1.在MATLAB的m文件中，可以将滞后环节转化为有理数形式的函数是？（）答案:pade2.在simulink中，有多种PID模块，其中下面哪一个是离散的PID模块（）答案:;3.在MATLAB的m文件中，LTI对象中属于滞后特征的属性是( )答案:Outputdelay;IOdelay;Inputdelay4.simulink中，PID模块中的只有Kp，Ki，Kd三个参数可以设置。
《计算机仿真技术》试题(含完整答案)
一、数值计算，编程完成以下各题〔共20分，每题5分〕1、脉冲宽度为d ，周期为T 的矩形脉冲的傅里叶级数如下式描述：)2cos(/)/sin(21[)(1τπππτn T d n T d n T d f n ∑∞=+=当150=n ，41=T d ，2/12/1<<-τ，绘制出函数)(τf 的图形。
解： syms n t;f=((sin(n*pi/4))/(n*pi/4))*cos(2*pi*n*t); s=symsum(f,n,1,150); y=(1+2*s)/4; x=-0.5:0.01:0.5; Y=subs(y,'t',x); plot(x,Y)2、画出函数5505.025.55.15.1cos 5)5(sin )(2x x x x ex x f x +++-=在区间[3,5]的图形，求出该函数在区间[3, 5]中的最小值点m in x 和函数的最小值m in f .解：程序如下 x=3:0.05:5;y=(sin(5*x).^2).*exp(0.05*x.^2)-5*(x.^5).*cos(1.5*x)+1.5*abs(x+5.5)+x.^2.5; mix_where=find(y==min(y)); xmin=x(mix_where); hold on; plot(x,y);plot(xmin,min(y),'go','linewidth',5);str=strcat('(',num2str(xmin),',',num2str(min(y)),')'); text(xmin,min(y),str); Xlabel('x')Ylabel('f(x)')经过运行后得到的图像截图如下：运行后的最小值点m in x =4.6，m in f = -8337.86253、画出函数xe x xf x 5.2cos )(3.02-⋅=-在[1，3]区间的图形，并用编程求解该非线性方程0)(=x f 的一个根，设初始点为20=x .解： x=1:0.02:3; x0=2;y=@(x)(cos(x).^2).*exp(-0.3*x)-2.5*abs(x); fplot(y,[1,3]); Xlabel('x') Ylabel('f(x)')X1=fzero('(cos(x).^2).*exp(-0.3*x)-2.5*abs(x)',x0)运行后求得该方程的一个根为z=0.3256。
计算机仿真技术 试卷及答案 共2套 AB卷
《计算机仿真技术》考试试卷(A卷)班级姓名学号一．填空题。
(每空1分，共20分)1、计算机仿真中系统由_________________、_________________、________________等三个要素组成。
2、常见的系统实验建模法有________________、________________、________________、________________等几种方法。
3、叙述simulink中常见的可变步长的仿真算法有：________________、________________、________________、________________。
(任举四种)4、MA TALB的Comand windows窗口中可以显示workspace中的全部变量的命令是_________________、_________________。
5、计算机仿真算法的算法引起的误差是________________、________________。
6、三维立体绘图中，________________ 是三维曲线绘图命令，________________是三维网格线绘图命令，________________着色表面图绘图命令。
7、符号运算中，泰勒级数计算的函数是________________傅里叶变换的函数________________二．简答题。
(本题5小题，每题8分，共40分)1、计算机仿真的目的和作用2、列举几种计算机仿真软件，并简要介绍一下（要求4种以上）3、matlab主界面主要包括那些窗口，分别有什么样的功能4、局部变量和全局变量之间的区别5、计算机仿真技术中相似性主要包含那几个方面三．编程题 (本题6小题、共40分)根据下面要求，写出程序代码。
1、(5分)设(2E－C-1B)A T=C-1，其中E是4阶单位矩阵，A T是4阶矩阵A的转置。
求矩阵A2、(5分)输入下面的矩阵，求出特征多项式和特征根，并求出特征多项式的微分。
《计算机仿真技术》试卷B(含答案)
《计算机仿真技术》考试试卷(B 卷)班级 姓名 学号一．填空题。
(每空1分，共20分)1、计算机仿真技术三要素是_______________、_________________、________________。
2、计算机仿真中模型能够代表系统来进行研究，模型主要有________________、________________、________________、________________等几类模型。
3、叙述simulink 中常见的固定步长的仿真算法有：________________、________________、________________、________________。
(任举四种)4、MATALB 主要相似性主要包含以下几种相似_________________、________________、___________________。
5、Simulink 中主要包含_________________、________________、_______________等几种模块。
6、符号运算中，泰勒级数计算的函数是________________，拉普拉斯变换的函数________________，傅里叶级数的计算函数是________________。
二．简答题。
(本题5小题，每题8分，共40分)1、简述根据信号类型分类，计算机仿真主要包含哪几种类型。
2、简述蒙特卡洛法建模的步骤。
3、m 函数文件和m 文本文件之间的区别。
4、叙述三维立体绘图常用的几个绘图指令及其功能。
5、简述matlab 软件的基本组成三．编程题 (本题6小题、共40分)根据下面要求，只需写出程序代码，不需计算结果和图形。
1、(5分)写出右面的矩阵D ，求出其大于4的元素有几个，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=059123107D并通过变量列出其在矩阵的位置和序号2、(5分)求解下面矩阵A 的行列式和逆矩阵，并利用命令取出系数矩阵的逆矩阵的第二行的数据，以小数点后面两位数值的形式显示出来。
计算机仿真题库
25. 请比较数值积分法与离散相似法的区别。 数值积分法比较成熟，精度较高，但计算公式复杂，计算量大，适合离线仿真；而离散 相似法计算公式简单，计算量小，速度快，但精度较低，适合在线仿真。 26. 物理仿真 以真实系统的物理性质和几何形状相似为基础而其他性质不变来构造系统的物理模型， 并在物理模型上进行实验的过程称为物理仿真，也称为实体仿真。其优点在于直观、形象； 缺点在于模型改变困难，实验限制较多，投资较大。 27. 虚拟现实 虚拟现实是用高科技手段构造出来的一种人工环境，它具有模仿人的视觉、听觉、触觉 等感知功能的能力，具有使人可以亲身体验沉浸在这种虚拟环境并与之相互作用的能力。 28. 多媒体仿真 将数字、文字、声音、图形、图像和动画等各种媒体有机组合，并与先进的计算机通信 和广播电视技术相结合，形成一个可组织、存储、操纵和控制多媒体信息的集成环境和交互 系统的仿真技术。 29. 根据系统的动态特性，一般可以将系统做哪些分类？ 根据系统的动态特性，可以将系统分为： （1）连续系统和离散系统。时间和各个组成部 分的变量都具有连续形式的系统称为连续系统； 全部或一些组成部分的变量具有离散信号形 式的系统称为离散系统。 （2）线性或非线性系统。如果微分方程的系数是常数或者仅仅是自 变量的函数， 则为线性微分方程， 其符合叠加性原理， 可分为线性定常系统， 线性时变系统； 非线性系统用非线性微分方程表示，其不符合叠加原理。 （3）确定系统和随机系统。确定系 统可以用确定的数学模型进行描述， 输出和输入变量之间有完全确定的函数关系； 随即系统 则指系统内部或环境变量发生不确定变动， 会影响输出变量。 （4） 单变量系统和多变量系统。 这是对输入输出变量个数而言的。 30. 在连续系统的数字仿真中， 仿真算法的选择一般为考虑什么？常用的仿真算法有哪些？ 在连续系统的数字仿真中，选择仿真算法一般会考虑求解精度、速度、数值稳定性、自 启动能力、步长等。常用的仿真方法有：欧拉法、四阶龙格-库塔法、离散相似法等。欧拉 法属于自启动算法，适用于线性系统和非线性系统，方法简单，但误差的积累导致求解精度 降低。四阶龙格库塔法也属于自启动算法，适用于线性和非线性系统，支持单步长和可变步 长，求解速度和精度较高，但仿真时间较长。离散相似法一般用于状态方程，按环节离散相 似法可用于非线性系统，按系统离散相似法仅限于线性系统。 31. 什么是实时仿真？什么情况下需要进行实时仿真？它在算法上有什么要求？ 实时仿真是指仿真模型的时间比例尺等于系统原模型的时间比例尺的一类仿真，对系统 进行仿真试验时，如果仿真系统有实物，由于实物和人是按真实时间变化和运动的，因此就 需要进行实时仿真。实时仿真要求仿真系统接收实时动态输入，并产生实时动态输出，输入 和输出通常为固定采样时间间隔的数列。 32. 什么是数值稳定性问题？如何保证计算中的数值稳定性？ 所谓的数值稳定性问题，是指在扰动影响下，其计算过程中误差的积累是否受到控制的 问题。要保证数值稳定性，应遵循两条原则： （1）保证计算稳定； （2）有一定的计算精度来 选择合理的计算步长。
计算机仿真技术试卷B
题号 一 二 三 四 总分 得分一、填空（共15分，每空1分）１．机械系统的分类是 、 。
２．仿真系统模型的应用分别是 、 、 。
３．空间状态模型主要包括 、 。
４．机械工程领域常见的数字仿真算法是； 、 。
５．功率键合图主要组成部分是 、 、 、 、 。
二、根据系统模型编写C 程序求曲柄滑块机构的位移、速度、加速度（共15分）弹簧质量系统如下图，相关参数为k1=1.2，k2=1.5，k3=1.6。
三、建立空间状态方程（共20分，每题10分）1、 分析质量弹簧系统，列出该系统的微分方程然后根据微分方程写出其状态方程2、将传递函数转化为空间状态方程写出计算步骤得分阅卷人 得分阅卷人 得分阅卷人2四、画出下列系统的功率键合图并写出绘图步骤（共30分，每题10分）１． 绘出右图的键合图并写出步骤２. 绘出右图小车系统键合图并写出绘图步骤3. 一个生长在罐中的细菌的简单模型。
假定细菌的出生率和当前细菌的总数成正比，死亡率和当前的总数的平方成正比，试建立仿真模型，画出键合图。
五、列出系统的状态方程和simulink 模型（共20分）列出右图蹦极系统的状态方程并写出其状态方程。
桥梁距地面的距离是h2=50m ，人距桥梁的距离是h1=30m 。
系统参数为k=20，a1=a2=1，m=70，g=10系统数学模型是2009～2010学年第二学期《计算机仿真技术》期终试卷（B卷）…………………………………装…………订………线…………………………………专业级班姓名学号得分…………………………………密…………封…………线………………………………第3 页共3 页。
计算机仿真思考题答案
沈阳理工大学研究生院 计算机仿真技术 期末复习题思考题1、某个特定的公司生产和销售高尔夫手推车。
每周周末，公司都将一周所生产的手推车转移到仓库（库存）中，卖出的所有手推车都是从库存中提取。
这个过程的一个简单模型为：)()()()1(k S k I k P k I -+=+其中：)(k P =第k 周生产的手推车数量；)(k I =第k 周库存的手推车数量； )(k S =第k 周所卖出的手推车数量；以下为10周计划中的每周销售额：假设每周的产量都基于前一周的销售额，所以有)1()(-=k S k P 。
假设第一周的产量为50辆手推车：即，50)1(=P 。
编写一个MATLAB 程序计算：10周之内每周库存之中的手推车数量或者计算手推车库存数量减少到0为止的时间，并同时绘制图形。
针对以下两种情况运行程序：a.初始库存为50辆手推车，所以50)1(=I ；b.初始库存为30辆手推车，所以30)1(=I 。
P = zeros(1,10);I1 = zeros(1,10); I2 = zeros(1,10); S = [50 55 60 70 70 75 80 80 90 55]; P(1)=50;I1(1)=50; I2(1)=30; for k=2:10 P(k) = S(k-1);I1(k) = P(k-1)+I1(k-1)-S(k-1); I2(k) = P(k-1)+I2(k-1)-S(k-1); endplot(I1,'k*'); hold on ; plot(I2,'ro'); hold off ;legend('出示库存为50辆','出示库存为30辆'); title('每周库存的手推车数');2、用两种方法求解方程组⎪⎪⎩⎪⎪⎨⎧=+++=++=++-=++012412324543213214321431x x x x x x x x x x x x x x 方法1a=[5 0 4 2;1 -1 2 1;4 1 2 0;1 1 1 1]; b=[3;1;1;0];x=inv(a)*b 或者x=a\b 方法2eq1='5*x1+4*x3+2*x4=3'; eq2='x1-x2+2*x3+x4=1'; eq3='4*x1+x2+2*x3=1'; eq4='x1+x2+x3+x4=0';[x1,x2,x3,x4]=solve(eq1,eq2,eq3,eq4) 3、 计算dx x 3101+⎰，精确到小数点下第四位。
计算机仿真技术试题含完整答案
一、数值计算，编程完成以下各题（共20分，每小题5分） 1、脉冲宽度为d ，周期为T 的矩形脉冲的傅里叶级数如下式描述：)2cos(/)/sin(21[)(1τπππτn T d n T d n T d f n ∑∞=+=当150=n ，41=T d ，2/12/1<<-τ，绘制出函数)(τf 的图形。
解： syms n t;f=((sin(n*pi/4))/(n*pi/4))*cos(2*pi*n*t); s=symsum(f,n,1,150); y=(1+2*s)/4; x=-0.5:0.01:0.5; Y=subs(y,'t',x); plot(x,Y)2、画出函数5505.025.55.15.1cos 5)5(sin )(2x x x x ex x f x +++-=在区间[3,5]的图形，求出该函数在区间[3, 5]中的最小值点m in x 和函数的最小值m in f .解：程序如下 x=3:0.05:5;y=(sin(5*x).^2).*exp(0.05*x.^2)-5*(x.^5).*cos(1.5*x)+1.5*abs(x+5.5)+x.^2.5; mix_where=find(y==min(y)); xmin=x(mix_where); hold on; plot(x,y);plot(xmin,min(y),'go','linewidth',5);str=strcat('(',num2str(xmin),',',num2str(min(y)),')'); text(xmin,min(y),str); Xlabel('x')Ylabel('f(x)')经过运行后得到的图像截图如下：运行后的最小值点m in x =4.6，m in f = -8337.86253、画出函数xe x xf x 5.2cos )(3.02-⋅=-在[1，3]区间的图形，并用编程求解该非线性方程0)(=x f 的一个根，设初始点为20=x .解： x=1:0.02:3; x0=2;y=@(x)(cos(x).^2).*exp(-0.3*x)-2.5*abs(x); fplot(y,[1,3]); Xlabel('x') Ylabel('f(x)')X1=fzero('(cos(x).^2).*exp(-0.3*x)-2.5*abs(x)',x0)运行后求得该方程的一个根为z=0.3256。
计算机仿真技术试卷
1、 控制系统的分析与设计方法可以分为 时域法 、 频域法 。
2、 根据仿真时间与系统实际时钟的关系，控制系统仿真可以分为 实时仿真 、 亚实时仿真 、 超实时仿真 。
3、 判断系统稳定性的方法主要有 利用极点判断稳定性 、 利用特征值判断系统稳定性 、利用李雅普洛夫第二法判断系统稳定性 。
4、 常用的离散化方法主要有 简单替换法 、 双线性变换法 、 零阶或一阶保持器法 。
5、 根据采样周期的选择不同，采样控制系统仿真可以分为 同步采样 、 异步采样 。
6、 增量式PID 控制器的表达式为()())2()1(2)()()1()()1()(-+--++--+-=k e k e k e k k e k k e k e k k u k u d i p 。
7、 已知系统脉冲传递函数为2121425)()()(----+-+==zz z z z z U z Y s G ，则该系统的差分方程为)3()2(5)3(4)2(2)(-+-+---=k u k u k y k y k y 。
8、 已知系统的状态方程为X X ⎥⎦⎤⎢⎣⎡--=2221&，则该系统是否稳定 是 （填是或否），其特征值为 2712,1j s ±-=。
1、简述计算机仿真的基本步骤。
答：（1）根据仿真目的确定仿真方案；（2）建立系统的数学模型；（3）选择合适仿真方法，建立仿真模型；（4）编写仿真程序并进行程序调试；（5）进行仿真实验；（6）仿真结果分析。
2、试比较数值积分法与离散相似法的区别。
答：数值积分法：比较成熟，精度较高，计算公式复杂，计算量大，适合离线仿真；离散相似法：计算公式简单，计算量小，速度快，精度较数值积分法低，适合在线仿真。
3、已知系统结构图如下图所示，1、2、3、4、5为典型环节，a 为常数，试写出其连接矩阵W 、W0、Wc （10分）解：由图可知⎪⎪⎪⎩⎪⎪⎪⎨⎧+=+===-=415314231241x x u xx u x u xu ax r u ， 从而：r x x x x x a u u u u u U ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=000010100100101000100000100005432154321 []⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡==54321510000x x x x x x y所以：⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-=01001010*******000010000a W ， ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=000010W ， []10000=c W4、已知系统状态空间表达式为[]⎪⎪⎩⎪⎪⎨⎧=⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡=Xt y t u X X 01)()(012101&，⎥⎦⎤⎢⎣⎡=11)0(X ，计算步长1.0=h ，输入信号1=u （0≥t ）试采用欧拉法，四阶龙格－库塔法计算h t =时对应的y 值。
(完整版)《计算机仿真技术与CAD》习题答案
第0章绪论0-1 什么是仿真？它所遵循的基本原则是什么？答：仿真是建立在控制理论、相似理论、信息处理技术和计算机技术等理论基础之上的，以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验知识、统计数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性科学。
它所遵循的基本原则是相似原理。
0-2 仿真的分类有几种？为什么？答：依据相似原理来分：物理仿真、数学仿真和混合仿真。
物理仿真：就是应用几何相似原理，制作一个与实际系统相似但几何尺寸较小或较大的物理模型（例如飞机模型放在气流场相似的风洞中）进行实验研究。
数学仿真：就是应用数学相似原理，构成数学模型在计算机上进行研究。
它由软硬件仿真环境、动画、图形显示、输出打印设备等组成。
混合仿真又称数学物理仿真，它是为了提高仿真的可信度或者针对一些难以建模的实体，在系统研究中往往把数学仿真、物理仿真和实体结合起来组成一个复杂的仿真系统，这种在仿真环节中有部分实物介入的混合仿真也称为半实物仿真或者半物理仿真。
0-3 比较物理仿真和数学仿真的优缺点。
答：在仿真研究中，数学仿真只要有一台数学仿真设备（如计算机等），就可以对不同的控制系统进行仿真实验和研究，而且，进行一次仿真实验研究的准备工作也比较简单，主要是受控系统的建模、控制方式的确立和计算机编程。
数学仿真实验所需的时间比物理仿真大大缩短，实验数据的处理也比物理仿真简单的多。
与数学仿真相比，物理仿真总是有实物介入，效果直观逼真，精度高，可信度高，具有实时性与在线性的特点；但其需要进行大量的设备制造、安装、接线及调试工作，结构复杂，造价较高，耗时过长，灵活性差，改变参数困难，模型难以重用，通用性不强。
0-4 简述计算机仿真的过程。
答：第一步：根据仿真目的确定仿真方案根据仿真目的确定相应的仿真结构和方法，规定仿真的边界条件与约束条件。
第二步：建立系统的数学模型对于简单的系统，可以通过某些基本定律来建立数学模型。
《计算机仿真技术和CAD》习题答案解析
第0章绪论0-1 什么是仿真？它所遵循的基本原则是什么？答：仿真是建立在控制理论、相似理论、信息处理技术和计算机技术等理论基础之上的，以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验知识、统计数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性科学。
它所遵循的基本原则是相似原理。
0-2 仿真的分类有几种？为什么？答：依据相似原理来分：物理仿真、数学仿真和混合仿真。
物理仿真：就是应用几何相似原理，制作一个与实际系统相似但几何尺寸较小或较大的物理模型（例如飞机模型放在气流场相似的风洞中）进行实验研究。
数学仿真：就是应用数学相似原理，构成数学模型在计算机上进行研究。
它由软硬件仿真环境、动画、图形显示、输出打印设备等组成。
混合仿真又称数学物理仿真，它是为了提高仿真的可信度或者针对一些难以建模的实体，在系统研究中往往把数学仿真、物理仿真和实体结合起来组成一个复杂的仿真系统，这种在仿真环节中有部分实物介入的混合仿真也称为半实物仿真或者半物理仿真。
0-3 比较物理仿真和数学仿真的优缺点。
答：在仿真研究中，数学仿真只要有一台数学仿真设备（如计算机等），就可以对不同的控制系统进行仿真实验和研究，而且，进行一次仿真实验研究的准备工作也比较简单，主要是受控系统的建模、控制方式的确立和计算机编程。
数学仿真实验所需的时间比物理仿真大大缩短，实验数据的处理也比物理仿真简单的多。
与数学仿真相比，物理仿真总是有实物介入，效果直观逼真，精度高，可信度高，具有实时性与在线性的特点；但其需要进行大量的设备制造、安装、接线及调试工作，结构复杂，造价较高，耗时过长，灵活性差，改变参数困难，模型难以重用，通用性不强。
0-4 简述计算机仿真的过程。
答：第一步：根据仿真目的确定仿真方案根据仿真目的确定相应的仿真结构和方法，规定仿真的边界条件与约束条件。
第二步：建立系统的数学模型对于简单的系统，可以通过某些基本定律来建立数学模型。
计算机仿真技术与cad第三版课后练习题含答案
计算机仿真技术与CAD第三版课后练习题含答案前言计算机仿真技术与CAD第三版是近年来广受欢迎的CAD教材。
本书涵盖了CAD 的基础知识、CAD的应用、计算机仿真技术等内容，适合CAD初学者和高级用户使用。
本文主要介绍本书的课后练习题及答案，希望对大家有所帮助。
第一章课后练习题及答案1.1 选择题1.以下哪项不是CAD应用的领域？A. 工程学B. 数学C. 植物学D. 医学答案：C2.CAD是什么的缩写？A. Computer-ded Design（计算机辅助设计）B. Computer-ded Distance（计算机辅助距离）C. Capsule-med Destruction（胶囊定向破坏）D. Capability-Assisted Driving（能力协助驾驶）答案：A1.2 填空题1.CAD软件主要分为三种类型，分别是2D __________ 、3D__________ 、仿真软件。
答案：制图软件、制造软件2.AutoCAD是一款 __________ 系统。
答案：二维制图3.CAD软件通常可以完成______________ 、___________ 、____________ 等操作。
答案：制图、设计、分析1.3 综合应用题设计一件简单的图形，包括设备名称、设备布局图、视图图和详细图。
并在AutoCAD中实现。
答：见图。
简单图形简单图形图中的简单图形包括了设备名称、设备布局图、视图图和详细图。
其中，设备名称为“简单图形”，设备布局图包含了该设备的位置、大小和方向信息，视图图包含了设备的3D视角，详细图则针对设备内部的结构进行了细节描绘。
使用AutoCAD可以方便地实现这件图形的设计和制图。
第二章课后练习题及答案……第三章课后练习题及答案……第四章课后练习题及答案……结语以上是《计算机仿真技术与CAD第三版》的部分课后练习题及答案，希望对大家有所帮助。
本书内容丰富、简明易懂，是一本非常优秀的CAD教材。
计算机仿真题库
仿真是指用模型（物理模型或数学模型）代替实际的系统进行实验和研究，其遵循原理 抽象、相似性原理的原则，相似性原理包括几何相似、性能相似、环境相似等。 40. 数据插值和曲线拟合的区别？
数据插值是研究如何平滑的估算出基准数据之间的其他点的函数值，所以插值所得曲线 必定穿过基准数据；而曲线拟合研究的是如何寻找平滑曲线以最好的表现带噪声的测量数据， 但不要求你和曲线穿过这些测量数据点。 41. 数值运算与符号运算的区别？
21. Matlab 中提供的三种基本逻辑运算为：与(&)、或( | )、非( ~ )。
22. 计算机仿真的三要素为：系统、模型、计算机。
23. 简述龙格-库塔法的基本思想。
用几个点上的函数值的线性组合来代替泰勒展开式中的各阶导数，然后按泰勒级数展开
确定其中的系数，这样既可以避免计算高阶导数，又能提高积分的精度及截断误差阶数。
t=(0:20)/20; r=sin(2*pi*t)+2; [x,y,z]=cylinder(r,40); cx=imread('flowers.tif');
进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运 算规则，如矩阵 a 与 b 相乘（a*b）时必须满足 a 的列数等于 b 的行数。在加、减运算时数 组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即 为数组运算，如 a*b 为矩阵乘，a.*b 为数组乘。 36. 在 Matlab 中什么是图形句柄？图形句柄有什么用途？
在连续系统的数字仿真中，选择仿真算法一般会考虑求解精度、速度、数值稳定性、自 启动能力、步长等。常用的仿真方法有：欧拉法、四阶龙格-库塔法、离散相似法等。欧拉 法属于自启动算法，适用于线性系统和非线性系统，方法简单，但误差的积累导致求解精度 降低。四阶龙格库塔法也属于自启动算法，适用于线性和非线性系统，支持单步长和可变步 长，求解速度和精度较高，但仿真时间较长。离散相似法一般用于状态方程，按环节离散相 似法可用于非线性系统，按系统离散相似法仅限于线性系统。 31. 什么是实时仿真？什么情况下需要进行实时仿真？它在算法上有什么要求？
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、数值计算，编程完成以下各题（共20分，每小题5分） 1、脉冲宽度为d ，周期为T 的矩形脉冲的傅里叶级数如下式描述：)2cos(/)/sin(21[)(1τπππτn T d n T d n T d f n ∑∞=+=当150=n ，41=T d ，2/12/1<<-τ，绘制出函数)(τf 的图形。
解： syms n t;f=((sin(n*pi/4))/(n*pi/4))*cos(2*pi*n*t); s=symsum(f,n,1,150); y=(1+2*s)/4; x=-0.5:0.01:0.5; Y=subs(y,'t',x); plot(x,Y)2、画出函数5505.025.55.15.1cos 5)5(sin )(2x x x x ex x f x +++-=在区间[3,5]的图形，求出该函数在区间[3, 5]中的最小值点m in x 和函数的最小值m in f .解：程序如下 x=3:0.05:5;y=(sin(5*x).^2).*exp(0.05*x.^2)-5*(x.^5).*cos(1.5*x)+1.5*abs(x+5.5)+x.^2.5; mix_where=find(y==min(y)); xmin=x(mix_where); hold on; plot(x,y);plot(xmin,min(y),'go','linewidth',5);str=strcat('(',num2str(xmin),',',num2str(min(y)),')'); text(xmin,min(y),str); Xlabel('x')Ylabel('f(x)')经过运行后得到的图像截图如下：运行后的最小值点m in x =4.6，m in f = -8337.86253、画出函数xe x xf x 5.2cos )(3.02-⋅=-在[1，3]区间的图形，并用编程求解该非线性方程0)(=x f 的一个根，设初始点为20=x .解： x=1:0.02:3; x0=2;y=@(x)(cos(x).^2).*exp(-0.3*x)-2.5*abs(x); fplot(y,[1,3]); Xlabel('x') Ylabel('f(x)')X1=fzero('(cos(x).^2).*exp(-0.3*x)-2.5*abs(x)',x0)运行后求得该方程的一个根为z=0.3256。
4、已知非线性方程组如下，编程求方程组的解，设初始点为[1 0.5 -1].⎪⎩⎪⎨⎧=+=+-=+03352722yz z x x x解：%在新建中建立函数文件fun2_4.mfunction f=fun2_4(x)f=[x(1).^2+x(1)*sqrt(7)+2;x(1)+5*x(3).^2-3;x(2).*x(3)+3]; %非线性方程组求解主程序fxxfcz.mx0=[1 0.5 -1];fsolve(@fun2_4,x0)运行后结果为：ans =-1.3229 3.2264 -0.9298即是 x=-1.3229 y=3.2264 z=-0.9298 .二、控制系统仿真（15分）某控制系统的开环传递函数为：)105.0)(16()112.0)(15.1(6)(++++=s s s s s S G ，要求：编制一个完整的程序完成以下各小题的要求，所绘制的图形分别定义为四张图。
1） 绘制出系统的阶跃信号响应曲线（响应时间为s 30~0） 2） 绘制出系统的脉冲信号响应曲线（响应时间为s 20~0） 3） 绘制出系统的斜坡信号响应曲线（响应时间为s 10~0） 4） 绘制出系统的Bode 图（要求频率范围为2210~10-rad/sec ）解：由传递函数知，该传递函数是将其用零极点描述法描述的，将其化为用传递函数表述的形式为：s s s s s S G ++++=23205.63.0672.908.1)(，所以num=[0 1.08 9.72 6],den=[0.3 6.05 1 0]。
%用传递函数编程求解 num=[0 1.08 9.72 6]; den=[0.3 6.05 1 0]; sys=tf(num,den); t1=0:0.1:30; figure(1)step(sys) %绘制出系统的阶跃信号响应曲线 t2=0:0.1:20; figure(2)impulse(sys) %绘制出系统的脉冲信号响应曲线 t3=0:0.1:10; figure(3) ramp=t3;lsim(sys,ramp,t3);%绘制出系统的斜坡信号响应曲线 figure(4) w=10^(-2):10^2;bode(sys,w);%绘制出系统的Bode 图fig(1)系统的阶跃信号响应曲线fig(2)系统的脉冲信号响应曲线fig(3)系统的斜坡信号响应曲线fig(4)系统的Bode 图三、曲线拟合（15分）已知某型号液力变矩器原始特性参数，要求用多项式拟合的方法编程完成以下各小题： 1）用二阶多项式拟合出)(i K 曲线；用三阶多项式拟合出)(i η曲线；用三阶多项式拟合出)(i B λ曲线。
2）用不同的颜色和不同的线型，将)(i K 的原始特性参数数据点和二阶拟合曲线绘制在同一张图形中；将)(i η的原始特性参数数据点和三阶拟合曲线绘制在同一张图形中；将)(i B λ的原始特性参数数据点和四阶拟合曲线绘制在同一张图形中。
3）运行程序，写出)(i K 曲线的二阶拟合公式、)(i η曲线的三阶拟合公式和)(i B λ曲线的四阶拟合公式。
解：% 曲线拟合（Curve fitting）disp('Input Data--i; Output Data--k(i),\eta(i),\lambdaB(i):')x=[0.065,0.098,0.147,0.187,0.243,0.295,0.344,0.398,0.448,0.499];y1=[2.37,2.32,2.23,2.15,2.05,1.96,1.87,1.78,1.69,1.59];y2=[0.154,0.227,0.327,0.403,0.497,0.576,0.644,0.707,0.757,0.795];y3=[26.775,26.845,27.147,27.549,28.052,28.389,28.645,28.756,28.645,28.243]; figure(1)pf1=polyfit(x,y1,2)px1=polyval(pf1,x)plot(x,px1,'k')gridxlabel('转速比i')ylabel('变矩比K')title('二阶多项式拟合k曲线')%pausefigure(2)pf2=polyfit(x,y2,3)px2=polyval(pf2,x)plot(x,px2,'b')gridxlabel('转速比i')ylabel('效率\eta')title('三阶多项式拟合\eta 曲线')%pausefigure(3)pf3=polyfit(x,y3,4)px3=polyval(pf3,x)plot(x,px3,'-r')gridxlabel('转速比i')ylabel('泵轮转矩系数\lambdaB')title('四阶多项式拟合\lambdaB曲线' )%figure(4)pf1=polyfit(x,y1,2)px1=polyval(pf1,x)plot(x,y1,'or',x,px1,'k')gridxlabel('转速比i')ylabel('变矩比K')title('二阶多项式拟合k曲线')Legend('原始数据','拟合曲线')%将的原始特性参数数据点和二阶拟合曲线绘制在同一张图形中pausefigure(5)pf2=polyfit(x,y2,3)px2=polyval(pf2,x)plot(x,y2,'*m',x,px2,'b') gridxlabel('转速比i') ylabel('效率\eta')title('三阶多项式拟合\eta 曲线') Legend('原始数据','拟合曲线',0)%将的原始特性参数数据点和三阶拟合曲线绘制在同一张图形中 pause figure(6)pf3=polyfit(x,y3,4) px3=polyval(pf3,x) plot(x,y3,'pk',x,px3,'-r') gridxlabel('转速比i')ylabel('泵轮转矩系数\lambdaB') title('四阶多项式拟合\lambdaB 曲线' ) Legend('原始数据','拟合曲线',0)%将的原始特性参数数据点和四阶拟合曲线绘制在同一张图形中 y1=poly2str(pf1,'x') %)(i K 曲线的二阶拟合公式 y2=poly2str(pf2,'x') %)(i η曲线的三阶拟合公式 y3=poly2str(pf3,'x') %)(i B λ曲线的四阶拟合公式 运行后的结果如下：运行后的二阶，三阶，四阶拟合曲线函数为： y1 = 0.01325 x^2 - 1.8035 x + 2.491y2 =-0.12713 x^3 - 1.6598 x^2 + 2.4499 x + 0.0025474y3 =106.7407 x^4 - 199.9852 x^3 + 95.8404 x^2 - 8.7272 x + 26.9754四、微分方程求解。
（25分）自己选择确定一个三阶微分方程，自己设置初始条件，用ode45方法求微分方程的解。
要求:（例如：1)(8)(4)(2)(2233=+++t y dt t dy dtt y d dt t y d ，0)0(=y ，1)0(=dt dy ，0)0(22=dty d ） 1）仿真时间t=30秒2）结果绘制在一张图中，包括t y -曲线，一阶t y -'曲线，二阶t y -''曲线，三阶t y -'''曲线3）用图例命令分别说明四条曲线为“t y -”，“t y -'”，“t y -''” ，“t y -'''” 4）定义横坐标为“时间”，纵坐标为“输出”，图形标题名称为“微分方程的解”解：系统方程为 1)(8)(4)(2)(2233=+++t y dt t dy dtt y d dt t y d , 这是一个单变量三阶常微分方程。