数学模型实验打印.doc

实验1(1) 用起泡法对10个数由小到大排序. 即将相邻两个数比较,将小的调到前头.function f=sort1(x)for i=1:10 for j=1:10-i if x(j)x(j+1) m=x(j+1); x(j+1)=x(j); x(j)=m; end end end x运行 x=1 2 1 4 2 7 6 9 4 8; sort1(x)结果x = 1 1 2 2 4 4 6 7 8 9(2) 有一个45矩阵,编程求出其最大值及其所处的位置.function f=maxmax1(A)m=max(max(A)for i=1:4 for j=1:5 if A(i,j)=m x=i,j end endend运行A=2 2 4 8 1;55 2 0 5 7;3 2 7 55 4;9 8 10 12 21;maxmax1(A)结果m = 55x = 2 1x = 3 4(3)编程求m=1;sum=0;for n=1:20 m=m*n; sum=sum+mendSum结果sum = 2.5613e+018(4)一球从100米高度自由落下,每次落地后反跳回原高度的一半,再落下. 求它在第10次落地时,共经过多少米?第10次反弹有多高?n=input(输入落地次数n=); %落地次数x(1)=100; %初始高度100msum=100; for i=2:n*2 %从100m高处落下为第一段路 x（1）-弹起为 x（2），着地一次-直到n次反弹，为 x（2*n）if mod(i,2)=0 %第偶数段路为前一段路的一半 x(i)=x(i-1)/2; else x(i)=x(i-1); %第奇数段路=前一段路 end sum=sum+x(i); %累计路程（包括反弹的那一段）endx(n*2) %输出第n次反弹的高度sum=sum-x(n*2) %输出第n次落地时，总共经历的路程结果输入落地次数n=10ans = 0.0977sum = 299.6094(5)有一函数 ，写一程序,输入自变量的值,输出函数值.function f=hanshu(x)f=x(1)2+sin(x(1)*x(2)+2*x(1);运行 x=1 2; hanshu(x)结果ans =3.9093实验21. 在同一平面中的两个窗口分别画出心形线和马鞍面. 要求 1)在图形上加格栅、图例和标注 2)定制坐标 3)以不同角度观察马鞍面subplot(1,2,1);cta=0:.01*pi:2*pi; rho=2.*(1-cos(cta);% rho=1+cos(cta);polar( cta,rho,r);grid on;xlabel(自变量x);ylabel(自变量y);title(心形线);subplot(1,2,2);x=-4:0.1:4;y=-3:0.1:5;X Y=meshgrid(x,y);Z=X.2-Y.2;surf(X,Y,Z);grid on;xlabel(自变量x);ylabel(自变量y);title(马鞍面);view(-20,20)2. 以不同的视角观察球面和圆柱面所围区域.subplot(1, 2, 1);x,y,z = sphere; % 定义球面等分数，缺省值为20surf(x, y, z); % 绘制图形hidden off; axis equal;hold on;x = 0:0.01:0.5;y = x;r = (x.2+y.2)/x,(x.2+y.2)/x;x,y,z=cylinder(r, 80); % 定义圆柱面surf(x,y,z);xlabel(x轴); ylabel(y轴); zlabel(z轴);grid on; box on;view(-20, 20);subplot(1, 2, 2);x, y, z=sphere; % 定义球面等分数，缺省值为20surf(x, y, z); % 绘制图形hidden off; axis equal;hold on;x = 0:0.01:0.5;y = x;r = (x.2+y.2)/x,(x.2+y.2)/x;x, y, z = cylinder(r