数学建模实验小结.doc

例1-1 r=2;V=4/3*pi*r^3 V = 33.5103 例2-1 s=0; for n=1:100 s=s+1/n/n; end s s = 1.6350 例2-5 两个一元函数y=x3-x-1,y=\x\0.2sin(5x)在区间-1x2的复合图 fplot(inline(x^3-x-1),[-1,2]); hold on; x=-1:0.2:2; y=abs(x).^0.2.*sin(5*x); plot(x,y,:ro); hold off; 曲面图 xa=6:8;ya=1:4; [x,y]=meshgrid(xa,ya); z=x.^2+y.^2; mesh(x,y,z) [x,y,z] ans = 6 7 8 1 1 1 37 50 65 6 7 8 2 2 2 40 53 68 6 7 8 3 3 3 45 58 73 6 7 8 4 4 4 52 65 80 例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]);pause mesh(x,y,z); Page40 1. 先在编辑器窗口写下列M函数，保存为ex2_1.m function [xbar,s]=ex2_1(x) n=length(x); xbar=sum(x)/n; s=sqrt((sum(x.^2)-n*xbar^2)/(n-1)); x=[81 70 65 51 76 66 90 87 61 77]; [xbar,s]=ex2_1(x) xbar = 72.4000 s = 12.1124 2. s=log(1);n=0; while s=100 n=n+1; s=s+log(1+n); end m=n 3. F(1)=1;F(2)=1;k=2;x=0; e=1e-8; a=(1+sqrt(5))/2; while abs(x-a)e k=k+1;F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1); end a,x,k m = 37 a = 1.6180 x = 1.6180 k = 21 4. clear;tic;s=0; for i=1:1000000 s=s+sqrt(3)/2^i; end s,toc tic;s=0;i=1; while i=1000000 s=s+sqrt(3)/2^i;i=i+1; end s,toc tic;s=0; i=1:1000000; s=sqrt(3)*sum(1./2.^i); s,toc s = 1.7321 Elapsed time is 2.038973 seconds. s = 1.7321 Elapsed time is 2.948968 seconds. s = 1.7321 Elapsed time is 0.453414 seconds 5. t=0:24; c=[15 14 14 14 14 15 16 18 20 22 23 25 28 ... 31 32 31 29 27 25 24 22 20 18 17 16]; 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,[-2 2]) (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); plo