数学建模实验小结.docx

数学建模实验小结

第一篇：数学建模实验小结例1-1>>r=2;V=4/3*pi*r^3V=33.5103例2-1计算s=...>>s=0;>>forn=1:100s=s+1/n/n;end>>ss=1.6350例2-5两个一元函数y=x3-x-1,y=x.2sin(5x)在区间-1y=abs(x).^0.2.*sin(5*x);plot(x,y,':ro');holdoff;曲面图>>xa=6:8;ya=1:4;>>[x,y]=meshgrid(xa,ya);>>z=x.^2+y.^2;>>mesh(x,y,z)>>[x,y,z]ans=例2-6二元函数图z=xexp(-x2-y2).xa=-2:0.2:2;ya=xa;[x,y]=meshgrid(xa,ya);z=x.*exp(-x.^2-y.^2);mesh(x,y,z);pause;surf(x,y,z);pause;contour(x,y,z,[0.1,0.1]);pausemesh(x,y,z);Page401.先在编辑器窗口写下列M函数，保存为ex2_1.mfunction[xbar,s]=ex2_1(x)n=length(x);xbar=sum(x)/n;s=sqrt((sum(x.^2)-n*xbar^2)/(n-1));>>x=[81706551766690876177];>>[xbar,s]=ex2_1(x)xbar=72.4000s=12.11242.s=log(1);n=0;whiless=s+log(1+n);endm=n3.F(1)=1;F(2)=1;k=2;x=0;e=1e-8;a=(1+sqrt(5))/2;whileabs(x-a)>ek=k+1;F(k)=F(k-1)+F(k-2);x=F(k)/F(k-1);enda,x,km=a=1.6180x=1.6180k=4.clear;tic;s=0;fori=1:1000000s=s+sqrt(3)/2^i;ends,toctic;s=0;i=1;whileis=s+sqrt(3)/2^i;i=i+1;ends,toctic;s=0;i=1:1000000;s=sqrt(3)*sum(1./2.^i);s,tocs=1.7321Elapsedtimeis2.038973seconds.s=1.7321Elapsedtimeis2.948968seconds.s=1.7321Elapsedtimeis0.453414seconds5.t=0:24;c=[15141414141516182022232528...313231292725242220181716];plot(t,c)6.(1)x=-2:0.1:2;y=x.^2.*sin(x.^2-x-2);plot(x,y)y=inline('x^2*sin(x^2-x-2)');fplot(y,[-22])(2)t=linspace(0,2*pi,100);x=2*cos(t);y=3*sin(t);plot(x,y)(3）x=-3:0.1:3;y=x;[x,y]=meshgrid(x,y);z=x.^2+y.^2;surf(x,y,z)（4）x=-3:0.1:3;y=-3:0.1:13;[x,y]=meshgrid(x,y);z=x.^4+3*x.^2+y.^2-2*x-2*y-2*x.^2.*y+6;surf(x,y,z)(5)t=0:0.01:2*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z)7.x=linspace(0,pi,100);y1=sin(x);y2=sin(x).*sin(10*x);y3=-sin(x);plot(x,y1,x,y2,x,y3)%page41,ex7x=-1.5:0.05:1.5;y=1.1*(x>1.1)+x.*(x=-1.1)-1.1*(xPage591.>>a=[123];b=[243];>>a./bans=0.50000.50001.0000>>a.bans=>>a/bans=0.6552>>abans=0000000.66671.33331.00002.(1)>>a=[41-1;32-6;1-53];b=[9;-2;1];>>abans=2.38301.48942.0213(2)>>a=[4-33;32-6;1-53],b=[-1;-2;1]a=-5b=>>abans=-0.4706-0.29410(3)>>a=[41;32;1-5],b=[1;1;1]a=-5b=>>abans=0.3311-0.1219(4)>>a=[21-11;121-1;1121],b=[1;2;3]a=-1b=>>abans=006.(1)>>a=[41-1;32-6;1-53]