上机报告0312.doc

佛 山 科 学 技 术 学 院 《 数学建模 》上机报告 专 业 13数学与应用数学 姓 名 叶-- 成 绩 班 级 13 学 号 2013 日 期 项目一 MATLAB程序设计 一、1．掌握M函数文件和M脚本文件的建立； 2．掌握数组和矩阵的表示，能够独立编写Matlab程序；掌握使用Matlab作二维曲线三维曲面图形参数方程作图；．。 、，其中， (2) ，其中，， (3) ，其中， 注: 自变量取值范围自己确定，以能反映指数函数与S型曲线函数增长方式的区别为标准。 2. 用surf,mesh绘制曲面 3. 用ezplot绘制摆线的图形，,a=1 4. 用起泡法对10个数由小到大排序,即将相邻两个数比较,将小的调到前头 5. 有一个矩阵,编程求出其元素最大值及所处位置 6. 有一函数,写一程序,输入自变量的值,输出函数值 三、t=linspace(0,100); x0=3.9; xm=311.9526; r1=0.304; r=0.2798; y1=x0*exp(r1*t); y2=xm./(1+((xm/x0)-1).*exp(-r.*t)); y3=r*(1-(t./xm)).*t; subplot(2,2,1);plot(t,y1),title(y1); subplot(2,2,2);plot(t,y2),title(y2); subplot(2,2,3);plot(t,y3),title(y3) 图像为： 2.(1)代码：x=-3:0.1:3; y=1:0.1:5; [X,Y]=meshgrid(x,y); Z=2.*(X).^2+(Y).^2; surf(X,Y,Z) （2）代码：x=-3:0.1:3; y=1:0.1:5; [X,Y]=meshgrid(x,y); Z=2.*(X).^2+(Y).^2; mesh(X,Y,Z) 3，代码：ezplot(t-sin(t),1-cos(t),[0,2*pi]) 4.代码： a=[3,2,4,1,6,5,9,7,8,0]; for j=9:-1:1 for i=1:j if(a(i)a(i+1)) t=a(i);a(i)=a(i+1);a(i+1)=t;end end end a 答案：a= 0 1 2 3 4 5 6 7 8 9 5.代码： a= 10000*rand(4,5) max=a(1:1,1:1); hang=1; lie=1; for i=1:4; for j=1:5; x=a(i:i,j:j); if xmax; max=x; hang=i; lie=j; end end end max hang lie 结果： a = 1.0e+003 * 8.1472 6.3236 9.5751 9.5717 4.2176 9.0579 0.9754 9.6489 4.8538 9.1574 1.2699 2.7850 1.5761 8.0028 7.9221 9.1338 5.4688 9.7059 1.4189 9.5949 max = 9.7059e+003 hang = 4 lie = 3 a = 1.0e+003 * 6.5574 6.7874 6.5548 2.7692 6.9483 0.3571 7.5774 1.7119 0.4617 3.1710 8.4913 7.4313 7.0605 0.9713 9.5022 9.3399 3.9223 0.3183 8.2346 0.3445 max = 9.5022e+003 hang = 3 lie = 5 6.代码：function f=fun(x,y) x=input(please?input?x=); y=input(please?input?y=); f=x^2+sin(x*y)+2*y 结果：please?input?x=1 please?input?y=1 f = 3.8415 ans =