《计算机仿真》试题及答案-工程硕士
计算机仿真试题和答案
控制系统计算机辅助设计综合实验指导实验名称：连续与离散系统校正实验，系统可控性与可观性实验，系统的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－1 仿真遵循的基本原则是什么？相似原理 1 - 2龙格－库塔的思想(省略) 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 值。
（12分）解：⎥⎦⎤⎢⎣⎡=11)0(X ，⎥⎦⎤⎢⎣⎡=2101A ，⎥⎦⎤⎢⎣⎡=01B ，[]01=C 欧拉法：()()⎥⎦⎤⎢⎣⎡=⋅⎪⎪⎭⎫⎝⎛⋅⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡=++⎥⎦⎤⎢⎣⎡=+=3.12.11.010*********)0()0(11))0(,0()0()(h Bu AX hX f X h X 2.1)()(==h CX h y 四阶龙格－库塔法： ()⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡==3201112101)0(,01X f K ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛⎥⎦⎤⎢⎣⎡⋅+⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎪⎭⎫⎝⎛⋅++=4.31.2013221.011210121.0)0(,21.0012K X f K⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛⎥⎦⎤⎢⎣⎡⋅+⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎪⎭⎫⎝⎛⋅++=445.3105.2014.31.221.011210121.0)0(,21.0023K X f K()⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛⎥⎦⎤⎢⎣⎡⋅+⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡=⋅++=8995.32105.201445.3105.21.01121011.0)0(,1.0034K X f K()⎥⎦⎤⎢⎣⎡=++++=3432.12103.1226)0()(4321K K K K hX h X 2103.1)()(==h CX h y5、已知系统传递函数)1(2)()()(+==s s s U s Y s G ，试采用双线性变换法求解系统差分方程，计算步长h ＝0.1s （10分）。
《计算机仿真技术》试题(含完整答案)
一、数值计算，编程完成以下各题〔共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。
2020年秋冬智慧树知道网课《计算机仿真技术》课后章节测试答案
第一章测试1【单选题】(2分)计算机仿真系统的三要素（）A.计算机、模型、系统B.系统、模型、仿真C.计算机、系统、数学公式D.模型、仿真、计算机2【多选题】(2分)系统研究一般有哪些方法（）A.直接描述法B.仿真实验法C.理论解析法D.直接实验法3【单选题】(2分)哪些是以硬件为基础的仿真软件()A.MATHCADB.MAPLEC.MATLABD.Proteus4【判断题】(2分)欠实时仿真，仿真时钟比实际时钟慢。
例如烟火爆炸的仿真分析。
A.错B.对5【多选题】(2分)计算机仿真算法的误差主要有（）A.偶然误差B.截断误差C.舍入误差D.系统误差第二章测试1【单选题】(2分)MATLAB软件中，下面哪个命令是用来显示当前文件夹中文件名字的（）A.whatB.dirC.whoD.whos2【单选题】(2分)下面这些哪个是属于元素群运算（点运算）（）A.x/2B.x+bC.exp(x)D.X*y3【单选题】(2分)哪个是以硬件为基础的仿真软件()A.MAPLEB.MATLABC.ProteusD.MATHCAD4【判断题】(2分)。
计算机仿真技术试题含完整答案
一、数值计算，编程完成以下各题（共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。
《计算机仿真技术》试卷A(含答案)
《计算机仿真技术》考试试卷(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分)输入下面的矩阵，求出特征多项式和特征根，并求出特征多项式的微分。
计算机仿真 习题3 课后答案
习题3
1、已知系统响应函数为)sin(11)(θ+-=-bt e b
t y bt ，其中21l b -=，)1arctan(2
l
l -=θ，要求用不同的线型或颜色，在同一张图上绘制l 取值为0.2,、0.4、0.6、0.8时，系统在[]18,0∈t 区间内的响应曲线，并要求用2.0=l 和8.0=l 对它们相应的两条曲线进行文字标志。
答：刚开始因为没有考虑到相乘或相除矩阵的维度而总出错，后来理解清楚之后就简单多了。
编写程序如下：
得到结果：
2、用plot
3、mesh 、surf 指令，绘制2222)1()1(1y x y x z ++++-=
三维图（x 、y 范围自定）。
答：用plot3时，我定义]100,0[∈x ，]100,0[∈y ；并分别令22)1(y x a +-=，22)1(y x b ++=。
同时需注意将1变换成与a 、b 同维的矩阵。
得到结果：
在利用mesh、surf绘图时，为了图形美观取]
x，x=y.
[-
∈
100
,
100
程序如下：
得到结果：
3、对向量t进行以下运算可以构成3个坐标的值向量：)
y=，
cos(t
x=，)
sin(t
z=。
利用plot3指令，并选用绿色的实线绘制相应的三维曲线。
t
答：取pi
:0(
:1.0
=。
注意绿色是‘g’。
编写程序如下：
t*
)
10
得到结果：。
计算机仿真与建模考试试题
计算机仿真与建模考试试题一、单项选择题1. 在计算机仿真与建模中，以下哪个软件是用于创建和模拟三维计算机图形的？A. MATLABB. AutoCADC. SolidWorksD. 3ds Max2. 计算机仿真与建模常用的数学建模方法有哪些？A. 微分方程建模B. 蒙特卡洛模拟C. 神经网络模型D. 遗传算法优化3. 下面哪个软件是用于分析流体和气体流动的计算机仿真工具？A. ANSYSB. CATIAC. Pro/ENGINEERD. PTC Creo4. 在计算机仿真中，有哪些常见的物理模型？A. 质点模型B. 弹簧-质点模型C. 流体模型D. 电路模型5. 计算机仿真与建模的主要应用领域是什么？A. 机械工程B. 航空航天C. 生物医学D. 所有以上领域二、简答题1. 什么是计算机仿真与建模？计算机仿真与建模是指使用计算机技术和数学方法对现实世界的问题或系统进行建模和模拟实验的过程。
通过对系统各个组成部分的特性进行数学描述和计算机编程实现，可以模拟真实系统的行为，并进行虚拟实验和预测。
2. 计算机仿真与建模的主要步骤是什么？计算机仿真与建模的主要步骤包括：问题建模与分析、数学建模、软件工具选择与开发、模型验证与调试、模拟实验与结果分析等。
3. 计算机仿真与建模的优势有哪些？计算机仿真与建模具有以下优势：- 可以对复杂系统进行全面、高效的分析和研究，避免了大量实际实验和浪费资源；- 提供了研究系统性能、优化设计和决策支持的基础；- 可以预测系统行为并进行虚拟实验，帮助解决实际问题和提高生产效率；- 提供了客观、可重复的实验环境，方便对结果进行验证和比较。
4. 计算机仿真与建模的应用范围有哪些？计算机仿真与建模广泛应用于工程领域，包括机械工程、电气工程、航空航天、建筑工程等。
此外，它还在生物医学、交通运输、环境科学等领域有重要应用。
三、综合题某工程团队需要设计一种新型的风力发电机，用于发电和供电给一个小型城市。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
机械工程学院2011—2012学年第二学期工程硕士《计算机仿真》试题姓名 成绩一、分别用矩阵求逆（inv ）指令和左除（\）指令编程求解如下线性方程组。
（5分）⎪⎩⎪⎨⎧=++=++=++134927753457321321321x x x x x x x x x答：程序如下：a= [7 1 5;3 5 7;2 9 4];b=[4 7 13]';X1=inv(a)*b %求逆指令求解X2=a\b %左除指令求解结果：X =[0.5561 1.4393 -0.2664]’程序运行及结果如下图（程序文件：zyl-1-mfile.m ）：二、画出函数5.05.02cos 2)5(sin )(205.02++-=x x x e x x f x 在区间[-5, 5]的图形，并用fminbnd 命令求出该函数在区间[-5, 5]中的最小值点min x 和函数的最小值min f .（5分）答：程序如下：x=-5:0.01:5; %定义区间f1='sin(5*x).^2.*exp(0.05*x.^2)-2*x.*cos(2*x)+0.5*sqrt((x+0.5).^2)'; %函数表达式 f2= sin(5*x).^2.*exp(0.05*x.^2)-2*x.*cos(2*x)+0.5*sqrt((x+0.5).^2) ; %函数表达式 plot(x,f2,'r-') %用红线画出图形[xmin,fmin]=fminbnd(f1,-5,5) %用fminbnd 命令求极小值结果：min x =-5, min f =-6.0796 程序运行结果及图形如下图（程序文件：zyl-2-mfile.m ）：三、画出函数x e x f 2s in )(=在0～π区间的图形，并用quad 命令编程求⎰π0)(dx x f 的积分值。
（5分）答：程序如下：x=0:pi/1000:pi; %定义区间（1000等分）f=exp((sin(x)).^2); %函数表达式plot(x,f,'r.') %用红线画出图形area=quad('exp((sin(x)).^2)',0,pi) %用quad 命令求积分值结果：积分值=area=5.5084 程序运行结果及图形如下图（程序文件：zyl-3-mfile.m ）：四、画出函数x e x x f x 6.0cos )(2.02-⋅=-在[1，3]区间的图形，并用fzero 命令编程求解该非线性方程0)(=x f 的一个根，设初始点为20=x 。
（5分） 答：程序如下：x=1:0.1:3; %定义自变量区间f=cos(x).^2.*exp(-0.2.*x)-0.6.*sqrt(x.^2); %自变量表达式plot(x,f,'r-') %用红线画出图形z=fzero('cos(x).^2.*exp(-0.2.*x)-0.6.*sqrt(x.^2)',2) %求初始点为2时的0点 结果：z = 0.7570程序运行结果及图形如下图（程序文件：zyl-4-mfile.m ）：五、已知非线性方程如下，编程求其解，设初始点为[2，2]。
（10分）⎪⎩⎪⎨⎧=+-=---2122222121x x x x x x答：程序如下：%程序myfun.mfunction y=myfun(x)y=[2*x(1)-x(2)-2^(-x(1));-x(1)+2*x(2)-2^(-x(2))];%主程序x0=[2;2];[x,f]=fsolve(@myfun,x0)结果：x =0.64120.6412f =1.0e-008 *-0.3702-0.3702程序运行结果如下图（程序文件：zyl-5-mfile.m ，myfun.m ）：六、在曲柄滑块机构中，滑块的位移为：22))sin(()cos(θθa b a s -+=，位移s 是曲柄转动角度θ的函数，a 为曲柄的长度，b 为连杆的长度。
当15=a ，20=b ，用plot 函数绘制出θ-s 的关系曲线图，已知 3600≤≤θ。
（10分）要求：1）对所绘制的曲线设置颜色和线型，并打开网格2）定义X 坐标轴为“曲柄转动角度θ” 定义Y 坐标轴为“滑块的位移s ”。
3）定义图形标题名称为“曲柄滑块θ-s 的关系曲线图”。
4）定义坐标轴范围380~0=X ，30~0=Y 。
5）在（100，15）处标出注释“θ-s 曲线”。
6）用图例注解说明函数曲线为“θ-s 曲线”。
答：程序如下：t=[0:pi/360:360]; %定义变量范围hold onaxis([0,380,0,30]); %定义坐标轴范围s=15*cosd(t)+sqrt(20.^2-(15*sind(t)).^2); %定义函数表达式plot(t,s,'r-') %设置曲线为红色实线xlabel('曲柄转动角度θ'); %定义X 坐标轴ylabel('滑块的位移s'); %定义Y 坐标轴title('曲柄滑块s-θ的关系曲线图'); %定义图形标题名称legend('s-θ曲线') %图例注解说明text(100,15,'s-θ曲线') %在(100,15)处标出注释grid on %打开栅格程序运行结果及图形如下图（程序文件：zyl-6-mfile.m ）：七、在]3,3[,-∈yx区域里,编制程序绘制226yxyez--=的三维网格图、三维曲面图和三维表面图，要求分别绘制成三张图。
（20分）答：程序如下：x=-3:0.2:3;y=x;[X,Y]=meshgrid(x,y); %产生“经纬”矩阵Z=6*Y.*exp(-X.^2-Y.^2);figure(1)mesh(Z) %三维网格图（四边形网格图）命令title('No.1 三维网格图')figure(2)surfl(X,Y,Z) %三维曲面图（四边形曲面图）命令title('No.2 三维曲面图')figure(3)surf(X,Y,Z) %三维表面图（四边形表面图）命令title('No.3 三维表面图')程序运行结果及图形如下图（程序文件：zyl-7-mfile.m）：八、微分方程求解，设1)(8.0)(3.0)(22=++t y dt t dy dtt y d ，0)0(=y ，0)0(=dt dy ，用ode45方法求微分方程的解。
结果用一张图表示，包括位移—时间（t y -）曲线，速度—时间（t y -'）曲线和加速度—时间（t y -''）曲线，仿真时间t=30秒。
（20分） 答：程序如下：%函数表达式文件myfun1.mfunction xdot=myfun1(t,x)xdot=[x(2);1-0.3*x(2)-0.8*x(1)];%主程序zyl-8-mfile.mt=[0 30];x0=[0;0];[tt,yy]=ode45(@myfun1,t,x0);plot(tt,yy(:,1),':b',tt,yy(:,2),'-r')hold onplot(tt,1-0.3*yy(:,2)-0.8*yy(:,1),'-.k')legend('位移曲线','速度曲线','加速度曲线')title('微分方程的解')程序运行结果及图形如下图（程序文件：zyl-8-mfile.m ，myfun1.m ）：九、已知柴油机的特性参数如下表所示，用多项式拟合的方法编程，分别拟合出转矩——转速曲线)M和功率——转速曲线)(nN。
要求（1）用二阶多项式和三阶多(n项式进行转矩——转速曲线)M的拟合，并将特性参数数据点、二阶拟合曲线和三(n阶拟合曲线，在同一张图形中用不同的颜色和不同的线型表达。
（2）用三阶多项式和四阶多项式进行功率——转速曲线)N的拟合，并将特性参数数据点、三阶拟合(n曲线和四阶拟合曲线，在同一张图形中用不同的颜色和不同的线型表达。
（3）分别写出拟合的转矩—转速曲线)M的二阶拟合公式和三阶拟合公式，功率—转速曲线(nN的三阶拟合公式和四阶拟合公式。
（20分）(n)答：（1）转矩-转速曲线拟合程序如下：x=[800:100:2200];y=[550.4,605.9,665.8,705.5,735.5,...745.5,741.5,735.58,728.68,720.8,...708.04,693.26,678.56,653.86,595.26];%绘参数数据点plot(x,y,'.r')hold on%绘二阶拟合曲线p2=polyfit(x,y,2)py2=polyval(p2,x);plot(x,py2,':b')hold on%绘三阶拟合曲线p3=polyfit(x,y,3)py3=polyval(p3,x);plot(x,py3,'g')grid onhold on%图形注释xlabel('转速n');ylabel('转矩M');legend('数据点','二阶拟合曲线','三阶拟合曲线')title('转矩-转速曲线M(n)拟合')结果：当n=2时，p2 =-0.0003 1.0198 -30.8427当n=3时，p3 =0.0000 -0.0011 2.1220 -524.5667 二阶拟合公式：20.0003 1.019830.8427y x x =-+- 三阶拟合公式：20.0011 2.1220524.5667y x x =-+-程序运行结果及图形如下图（程序文件：zyl-9(1)-mfile.m ）：（2）功率-转速曲线拟合程序如下：x=[800:100:2200];y=[48.69,61.28,73.28,86.86,94.81,...104.00,110.98,117.61,124.21,130.09,...135.24,140.38,145.53,150.30,154.35]; %绘参数数据点plot(x,y,'.r')hold on%绘三阶拟合曲线p3=polyfit(x,y,3)py3=polyval(p3,x);plot(x,py3,':b')hold on%绘四阶拟合曲线p4=polyfit(x,y,4)py4=polyval(p4,x);plot(x,py4,'g')grid onhold on%图形注释xlabel('转速n');ylabel('功率N');legend('数据点','三阶拟合曲线','四阶拟合曲线')title('功率-转速曲线N(n)拟合')结果：当n=3时，p3 = 0.0000 -0.0001 0.3141 -129.7814 当n=4时，p4 = -0.0000 0.0000 -0.0002 0.3799 -151.6974 三阶拟合公式：20.00010.3141129.7814y x x =-+- 四阶拟合公式：20.00020.3799151.6974y x x =-+-程序运行结果及图形如下图（程序文件：zyl-9(2)-mfile.m ）：要求：1）答卷中应含程序清单、计算结果和相关结果图。