《MATLAB程序设计应用试卷与答案(共3页).doc》由会员分享，可在线阅读，更多相关《MATLAB程序设计应用试卷与答案(共3页).doc（3页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上 专业班级： 学号： 姓名： 命题教师： 该项由出卷人填写 装 订 线一二三四五总 分标准分2020202020100得 分一、给出迭代方程 先编写求解方程的函数文件，然后调用该函数文件求30 000个点上的x,y,最后在所有的（）坐标处标记一个点（不要连线）会出图形。这种图形又称为埃农（Henon）引力线图，它将迭代出来的随机点吸引到一起，最后得出貌似连贯的引力线图。二、分别利用数值积分法、符号积分法和Simulink仿真求 。三、已知阿波罗（Apollo）卫星的运动轨迹（）满足下列微分方程：试在以上初值下进行数值求解，并绘制出阿波罗卫星位置（）的轨迹。 提示：先选
2、择一组状态变量，写出一阶常微分方程组，并定义相应的函数文件，然后求方程的数值解。四、实验图4所示是一个跷跷板，两板夹角为，左边板长为1.5m，上面的小孩重500N，右边板长为2m，小孩重400N。求当跷跷板平衡时，左边木板与水平方向夹角的大小。要求先求解析解，然后给出两种解决方案。 提示：这是一个力矩平衡问题，可列方程求解析解。可以考虑的两种方案是用迭代法解方程和分别绘制两个小孩所产生力矩随变化的曲线，两曲线的交点即是跷跷板平衡时的。 跷跷板示意图五、某公司投资2000万元建成一条生产线。投产后，在时刻t的追加成本和追加收益分别为（百万元/年），（百万元/年）。试确定该生产线在何时停产可获最大
3、利润？最大利润是多少？提示：利用函数（百万元），由于H（t）-G（t）单调下降，所以H（t）=G（t）时，R（t）取得最大利润。1解： 源程序如下： x(1)=0;y(1)=0;for i=1:30000;x(i+1)=1+y(i)-1.4*x(i)2;y(i+1)=0.3*x(i);hold onplot(x(i),y(i),*b)end埃农（Henon）引力线图如下：2解：数值积分法源程序代码： X=0:0.001:1;Y=(1/sqrt(2*pi)*exp(-X.2/2);trapz(X,Y) 运行结果如下：ans = 0.3413符号积分法源程序代码： x=sym(x);f=(1/sq
4、rt(2*pi)*exp(-x.2/2);I=int(f,0,1);double(I)运行结果如下：ans = 0.3413Simulink仿真3. function dx=appollo(t,x)mu=1/82.45;mustar=1-mu;r1=sqrt(x(1)+mu)2+x(3)2);r2=sqrt(x(1)-mustar)2+x(3)2);dx=x(2)2*x(4)+x(1)-mustar*(x(1)+mu)/r13-mu*(x(1)-mustar)/r23x(4)-2*x(2)+x(3)-mustar*x(3)/r13-mu*x(3)/r23;-x0=1.2;0;0;-1.;%x0
5、(i)对应与xi的初值 options=odeset(reltol,1e-8); tic t,y=ode45(appollo,0,20,x0,options); toc plot(y(:,1),y(:,3) title(Appollo卫星运动轨迹) xlabel(X) ylabel(Y)运行结果如下：Elapsed time is 0. seconds.图如下：4.解：方案一：迭代法源程序如下 for alpha=0:0.001:pi/3;if 750*cos(alpha)-800*cos(pi/3-alpha)10e-6continue专心-专注-专业 专业班级： 学号： 姓名： 该项由出卷
6、人填写装 订 线 专业班级： 学号： 姓名： 该项由出卷人填写装 订 线endbreakendalphacos(alpha)运行结果如下：alpha = 0.4680ans = 0.8925方案二：曲线相交法源程序如下alpha=0:pi/1000:pi/3;y1=750*cos(alpha);y2=800*cos(pi/3-alpha);plot(alpha,y1,g-,alpha,y2,r)grid on运行得到的效果图如下：该图放大后可大致得到alpha和cos(alpha)的值，见下图：5.解：构造函数f(t)=H(t)-G(t)=13-t-3t2/3=0;令t1/3=x,则f(t)=
7、-t3-3t2+13可得矩阵P=-1,-3,0,13求最佳生产时间的源程序如下：p=-1,-3,0,13;x=roots(p);t=x.3运行结果如下：t = 3.6768 +21.4316i 3.6768 -21.4316i 4.6465 再分别将t的三个值带入函数f(t)，比较大小后，得到最大利润与最佳生产时间。求最大利润的程序代码如下： t=3.6768 +21.4316i;x=0:0.01:t;y=13-x-3*x.(2/3);trapz(x,y)运行结果： ans = 25.2583 t=3.6768 -21.4316i;x=0:0.01:t;y=13-x-3*x.(2/3);trapz(x,y)运行结果： ans = 25.2583 t=4.6465;x=0:0.01:t;y=13-x-3*x.(2/3);trapz(x,y)运行结果： ans = 26.3208比较以上三组数据，可知最佳生产时间t=4.6465年，可获得的最大利润 26.3208（百万元/年）。