原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
2009～2010学年 第 一 学期 研究生课程考试试题 课程名称： MATLAB程序设计与应用 考试形式： 开 卷 试卷：
第 ( 2 ) 页 共 ( 3 ) 页
专业班级： 09研究生 学号： 姓名： 命题教师：
[该项由出卷人填写] 装 订 线
一 二 三 四 五 六 七 八 九 十 总 分 标准分 20 20 20 20 20 100 得 分 1. 给出迭代方程。
先编写求解方程的函数文件，然后调用该函数文件求30000个点上的x，y，最后在素有的（xi，yi）坐标处标记一个点（不要连线）绘出图形。这种图形又称为埃农（Henon）引力线图，它将迭代出来的随机点吸引到一起，最后得出貌似连贯的引力线图。
2.分别利用数值积分法、符号积分法和Simulink仿真求。
3.某气象观测站测得某日6：00-18：00之间每隔2h的室内外温度（℃）如实验表1所示。
实验表1 室内外温度观测结果（℃）
时间/h 6 8 10 12 14 16 18 室内t1 18.0 20.0 22.0 25.0 30.0 28.0 24.0 室外t2 15.0 19.0 24.0 28.0 34.0 32.0 30.0 试用三次样条插值分别求出该日室内外6：30-17：30之间每隔2h各点的近似温度（℃）。
4.有一跷跷板，两板夹角为120°，左边板长尾1.5m，上面的小孩重500N，右边板长为2m，小孩重400N。求当跷跷板平衡时，左边木板与水平方向的夹角α的大小。要求先求解析解，然后给出两种解决方案。
5.某公司投资2000万元建成一条生产线。投产后，在时刻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 on
plot(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/sqrt(2*pi)*exp(-x.^2/2));
I=int(f,0,1);
double(I)
运行结果如下：
ans =
0.3413
③Simulink仿真
3.解：源程序如下：
h=6:2:18;
x=6.5:2:17.5;
t1=[18.0 20.0 22.0 25.0 30.0 28.0 24.0];
t2=[15.0 19.0 24.0 28.0 34.0 32.0 30.0];
T1=interp1(h,t1,x,spline)
T2=interp1(h,t2,x,spline)
运行结果如下：
T1 =
18.5020 20.4986 22.5193 26.3775 30.2051 26.8178
T2 =
15.6553 20.3355 24.9089 29.6383 34.2568 30.9594
4.解：①方案一：迭代法源程序如下
for alpha=0:0.001:pi/3;
if 750*cos(alpha)-800*cos(pi/3-alpha)10e-6
continue
end
break
end
alpha
cos(alpha)
运行结果如下：
alpha =
0.4680
ans =
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
求最佳生产时间的源程序如下：
p=[-1,-3,0,13];
x=roots(p);
t=x.^3
运行结果如下：
t =
3.6768 +21.4316i
3.6768 -2
文档评论（0）