数学实验二ans
数学实验练习二 参考答案
数学实验二实验内容：学习matlab的m文件编写和函数的编写，体会matlab编程特点，掌握matlab 的编程基本方法。
要求：一．学习ppt教案的例题代码，能正确的输入、运行代码；二．写出如下各段代码的作用，将以下各段循环执行的代码，改为不需要循环的矩阵和数组运行，并使用tic,toc测试不同代码的执行时间：%程序1,文件名：ex2_2_1.mticdx = pi/30;nx = 1 + 2*pi/dx;for i = 1:nxx(i) = (i-1)*dx;y(i) = sin(3*x(i));endtoc以上程序实现将[0,2*pi]间隔pi/30分成60等分，x和y分别为61个元素的数组，y为计算sin(3x)的值。
以上程序可以使用简单的matlab数组计算实现：x2=0:pi/30:2*pi;y2=sin(3*x2);大家可以比较一下，x1和x2完全相同，y和y2也完全相同。
%程序2，文件名：ex2_2_2.mticA=round(2+rand(50,60)*6); 生成一个在[2,8]上均匀分布的50*60随机数组[X,Y]=size(A); 求出其大小；X=50,Y=60minA=A(1,1); 设最小值为矩阵A的第1行1列的元素for i=1:Xfor j=1:Yif A(i,j)<minAminA=A(i,j);minX=i;minY=j;endendend 以上程序按行、列搜索矩阵A的最小值，若当前值A(I,j)小，则将最小值设为当前值；[minA ,minX,minY] 输出矩阵最小值minA及矩阵最小值所在的行minX、列minY。
toc上述程序可以使用find函数及min函数实现；此时只需： minA=min(A(:));[minX,minY]=find(A ’==minA,1);%注意此处需将矩阵A 转置，因为matlab 中是按列优先搜索的，而题目的程序是按行有限搜索。
数学实验(第二版)课后习题答案
贵州师范学院2012级数本一班李刚数学实验课后练习答案习题2.11. syms x y;>> x=-5:0.01:5;>> y=x.^1/2;>> plot(x,y)2. f plot('exp(-x.^2)',[-5,5])3. ezplot('x.^3+y.^3-3*x*y',[-5,5])4 . ezplot('y.^2-x.^3/(1-x)',[-5,5])5.t=0:0.1:2*pi;x=t-sin(t);y=2*(1-cos(t));plot(x,y)6. t=0:0.1:2*pi; x=cos(t).^3; >> y=sin(t).^3;>> plot(t,y)>>7: t=0:0.1:2*pi; x=cos(t); y=2*sin(t); z=3*t; plot3(x,y,z)8: x =0:0.1:2*pi; r=x; polar(x,r)9: x =0:0.1:2*pi; r=exp(x); polar(x,r)10: x=0:0.1:2*pi; r=sqrt(cos(2*x)); polar(x,r)11: x=0:0.1:2*pi; r=sqrt(sin(2*x)); polar(x,r)12: x =0:0.1:2*pi; r=1+cos(x); polar(x,r)练习2.2 1：(1)(2)：syms n; limit('sqrt(n+2)-2*(sqrt(n+1))+sqrt(n)',n,inf)Ans= 0 (3):: (4):(5):(6):2:3:fplot('x.^2*sin(x.^2-x-2)',[-2,2])练习2.3 1：（2）：2：练习2.4 1：（1）（2）：（3）（4）：2：（1）：syms x;int(x^(-x),x,0,1)ans =int(x^(-x),x = 0 .. 1)vpa(ans,10)ans =1.291285997（2）：syms x；int(exp(2*x)*cos(x)^3,x,0,2*pi)ans =-22/65+22/65*exp(4*pi)（3）：syms x; int(exp(x^2/2)/sqrt(2*pi),x,0,1)ans =-1125899906842624/5644425081792261*i*erf(1/2*i*2^(1/2))*pi^(1/2)*2^(1/2) >> vpa(ans,10)ans =.4767191345（4）：syms x;int(x*log(x^4)*asin(1/x^2),x,1,3)ans =int(x*log(x^4)*asin(1/x^2),x = 1 .. 3)>> vpa(ans,10)ans =2.459772128（5）：syms x ;int(exp(x^2/2)/sqrt(2*pi),x,-inf,inf)ans =Inf（6）：syms x ;int(sin(x)/x,x,0,inf)ans =1/2*pi（7）：syms x ;int(tan(x)/sqrt(x),x,0,1)Warning: Explicit integral could not be found. > In sym.int at 58ans =int(tan(x)/x^(1/2),x = 0 .. 1)>> vpa(ans,10)ans =.7968288892（8）：syms x ;int(exp(-x^2/2)/(1+x^4),x,-inf,inf)ans =1/4*pi^(3/2)*2^(1/2)*(AngerJ(1/2,1/2)-2/pi^(1/2)*sin(1/2)+2/pi^(1/2)*cos(1/2)-WeberE(1/2,1/2 ))>> vpa(ans,10)ans =1.696392536（9）：syms x ;int(sin(x)/sqrt(1-x^2),x,0,1)ans =1/2*pi*StruveH(0,1)>> vpa(ans,10)ans =.8932437410练习2.5（1）：syms n;symsum(1/n^2^n,n,1,inf)ans =sum(1/((n^2)^n),n = 1 .. Inf)（2）：s yms n ;symsum(sin(1/n),n,1,inf)ans =sum(sin(1/n),n = 1 .. Inf)（3）：syms n ;symsum(log(n)/n^3,n,1,inf) ans =-zeta(1,3)（4）：syms n ;symsum(1/(log(n))^n,n,3,inf) ans =sum(1/(log(n)^n),n = 3 .. Inf)（5）：syms n;symsum(1/(n*log(n)),n,2,inf) ans =sum(1/n/log(n),n = 2 .. Inf)（6）：yms n;symsum((-1)^n*n/(n^2+1),n,1,inf)ans =-1/4*Psi(1-1/2*i)+1/4*Psi(1/2-1/2*i)-1/4*Psi(1+1/2*i)+1/4*Psi(1/2+1/2*i)第三章练习3.11：（1）：a=-30:1:30;b=-30:1:30;[x,y]=meshgrid(a,b);z=10*sin(sqrt(x.^2+y.^2))./(sqrt(1+x.^2+y.^2)); meshc(x,y,z)（2）：a=-30:1:30;b=-30:1:30;[x,y]=meshgrid(a,b);z=4*x.^2/9+y.^2;meshc(x,y,z)（3）：（4）：a=-30:1:30;b=-30:1:30;[x,y]=meshgrid(a,b); z=x.^2/3-y.^2/3; meshc(x,y,z)（5）：a=-30:1:30;>> b=-30:1:30;>> [x,y]=meshgrid(a,b); >> z=x*y;>> meshc(x,y,z)（6）：（7）：a=-30:1:30;>> b=-30:1:30;>> [x,y]=meshgrid(a,b); >> z=sqrt(x.^2+y.^2); >> meshc(x,y,z)（8）：（9）：a=-30:1:30;>> b=-30:1:30;>> [x,y]=meshgrid(a,b);>> z=atan(x./y);>> meshc(x,y,z)练习3.21;a=-1:0.1:1;>> b=0:0.1:2;>> [x,y]=meshgrid(a,b);>> z=x.*exp(-x.^2-y.^2);>> [px,py]=gradient(z,0.1,0.1);>> contour(a,b,z)>> hold on>> quiver(a,b,px,py)2:a=-2:0.1:1;>> b=-7:0.1:1;>> [x,y]=meshgrid(a,b);>> z=y.^3/9+3*x.^2.*y+9*x.^2+y.^2+x.*y+9; >> plot3(x,y,z)>> grid on3:[x,y]=meshgrid(-2*pi:0.2:2*pi); z=x.^2+2*y.^2;plot3(x,y,z)hold onezplot('x^2+y^2-1',[-2*pi,2*pi]) ; grid on4:t=0:0.03:2*pi;>> s=[0:0.03:2*pi]';>> x=(0*s+1)*cos(t);y=(0*s+1)*sin(t);z=s*(0*t+1); >> mesh(x,y,z)>> hold on>> [x,y]=meshgrid(-1:0.1:1);>> z=1-x+y;>> mesh(x,y,z)5:syms x y z dx dyz=75-x^2-y^2+x*y;zx=diff(z,x),zy=diff(z,y)zx =-2*x+yzy =-2*y+x练习3.31：ezplot('x^2+y^2-2*x',[-2,2]);>> grid onsyms x y ;s=int(int(x+y+1,y,-sqrt(1-(x-1)^2),sqrt(1-(x-1)^2)),x,0,2)s =2*pi2：syms r t ;>> s=int(int(sqrt(1+r^2*sin(t)),r,0,1),t,0,2*pi)s =int(1/2*((1+sin(t))^(1/2)*sin(t)^(1/2)+log(sin(t)^(1/2)+(1+sin(t))^(1/2)))/sin(t)^(1/2),t = 0 .. 2*pi) 3：syms x y z ;>> s=int(int(int(1/(1+x+y+z)^3,z,0,1-x-y),y,0,1-x),x,0,1)s =-5/16+1/2*log(2)4：s=vpa(int(int(x*exp(-x^2-y^2),y,0,2),x,-1,10))s =0.16224980455070416645061789474030练习3.41：（1）：y=dsolve('Dy=x+y','y(0)=1','x')得：y =-1-x+2*exp(x)（2）：y=dsolve('Dy=2*x+y^2','y(0)=0')y =tan(t*x^(1/2)*2^(1/2))*x^(1/2)*2^(1/2)练习4.11：（1）：p=[5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6 8 0 0 0 -5 0 0]; >> x=roots(p)x =0.97680.9388 + 0.2682i0.9388 - 0.2682i0.8554 + 0.5363i0.8554 - 0.5363i0.6615 + 0.8064i0.6615 - 0.8064i0.3516 + 0.9878i0.3516 - 0.9878i-0.0345 + 1.0150i-0.0345 - 1.0150i-0.4609 + 0.9458i-0.4609 - 0.9458i-0.1150 + 0.8340i-0.1150 - 0.8340i-0.7821 + 0.7376i-0.7821 - 0.7376i-0.9859 + 0.4106i-0.9859 - 0.4106i-1.0416-0.7927（2）: p=[8 36 54 23];x=roots(p)x =-1.8969 + 0.6874i-1.8969 - 0.6874i-0.70632:p1=[1 0 -3 -2 -1];p2=[1 -2 5];[q2,r2]=deconv(p1,p2)q2 =1 2 -4r2 =0 0 0 -20 19 3:syms x;f=x^4+3*x^3-x^2-4*x-3;g=3*x^3+10*x^2+2*x-3;p1=factor(f),p2=factor(g)p1 =(x+3)*(x^3-x-1)p2 =(x+3)*(3*x^2+x-1)4:syms x ;f=x^12-1;p=factor(f)p =(-1+x)*(1+x^2+x)*(1+x)*(1-x+x^2)*(1+x^2)*(x^4-x^2+1)5: (1):p=[1 0 1];q=[1 0 0 0 1];[a,b,r]=residue(p,q)a =-0.0000 - 0.3536i-0.0000 + 0.3536i0.0000 - 0.3536i0.0000 + 0.3536ib =0.7071 + 0.7071i0.7071 - 0.7071i-0.7071 + 0.7071i-0.7071 - 0.7071ir =[](2):p=[1];q=[1 0 0 0 1];[a,b,r]=residue(p,q)a =-0.1768 - 0.1768i -0.1768 + 0.1768i0.1768 - 0.1768i0.1768 + 0.1768ib =0.7071 + 0.7071i0.7071 - 0.7071i -0.7071 + 0.7071i -0.7071 - 0.7071ir =[](3):p=[1 0 1];q=[1 1 -1 -1];[a,b,r]=residue(p,q)a =0.5000-1.00000.5000b =-1.0000-1.00001.0000r =[] (4): p=[1 1 0 0 0 -8];[a,b,r]=residue(p,q)a =-4-38b =-11r =1 1 1练习 4.21：（1）：D=[2 1 3 1;3 -1 2 1;1 2 3 2;5 0 6 2];det(D)ans =6（2）：syms a b c dD=[a 1 0 0 ;-1 b 1 0;0 -1 c 1;0 0 -1 d];det(D)ans =a*b*c*d+a*b+a*d+c*d+12：（1）：D=[1 1 1 1; a b c d;a^2 b^2 c^2 d^2;a^3 b^3 c^3 d^3];det(D)ans =b*c^2*d^3-b*d^2*c^3-b^2*c*d^3+b^2*d*c^3+b^3*c*d^2-b^3*d*c^2-a*c^2*d^3+a*d^2*c^3+a *b^2*d^3-a*b^2*c^3-a*b^3*d^2+a*b^3*c^2+a^2*c*d^3-a^2*d*c^3-a^2*b*d^3+a^2*b*c^3+a^ 2*b^3*d-a^2*b^3*c-a^3*c*d^2+a^3*d*c^2+a^3*b*d^2-a^3*b*c^2-a^3*b^2*d+a^3*b^2*c（2）: s yms a b x y zD=[a*x+b*y a*y+b*z a*z+b*x; a*y+b*z a*z+b*x a*x+b*y;a*z+b*x a*x+b*y a*y+b*z];det(D)ans =3*a^3*x*z*y+3*b^3*y*x*z-a^3*x^3-a^3*y^3-b^3*z^3-a^3*z^3-b^3*x^3-b^3*y^33: (1): D=[1 1 1 1;1 2 -1 4;2 -3 -1 -5;3 1 2 11];D1=[5 1 1 1;-2 2 -1 4;-2 -3 -1 -5;0 1 2 11];D2=[1 5 1 1;1 -2 -1 4;2 -2 -1 -5;3 0 2 11];D3=[1 1 5 1;1 2 -2 4;2 -3 -2 -5;3 1 0 11];D4=[1 1 1 5;1 2 -1 -2;2 -3 -1 -2;3 1 2 0];x1=det(D1)/det(D);x2=det(D2)/det(D);x3=det(D3)/det(D);x4=det(D4)/det(D);x1,x2,x3,x4x1 =1x2 =2x3 =3x4 =-1(2):D=[5 6 0 0 0;1 5 6 0 0;0 1 5 6 0;0 0 1 5 6;0 0 0 1 5]; D1=[1 6 0 0 0;0 5 6 0 0;0 1 5 6 0;0 0 1 5 6;1 0 0 1 5]; D2=[5 1 0 0 0;1 0 6 0 0;0 0 5 6 0;0 0 1 5 6;0 1 0 1 5]; D3=[5 6 1 0 0;1 5 0 0 0;0 1 0 6 0;0 0 0 5 6;0 0 1 1 5]; D4=[5 6 0 1 0;1 5 6 0 0;0 1 5 0 0;0 0 1 0 6;0 0 0 1 5]; D5=[5 6 0 0 1;1 5 6 0 0;0 1 5 6 0;0 0 1 5 0;0 0 0 1 1]; x1=det(D1)/det(D);x2=det(D2)/det(D);x3=det(D3)/det(D);x4=det(D4)/det(D);x5=det(D5)/det(D);x1,x2,x3,x4,x5x1 =2.2662x2 =-1.7218x3 =1.0571x4 =-0.5940x5 =0.3188练习 4.3 1：A=[1 2 0;3 4 -1; 1 1 -1];B=[1 2 3;-1 0 1;-2 4 -3];A',2+A,2*A-B,A*B,A^2,A^(-1)ans =1 3 12 4 10 -1 -1ans =3 4 25 6 13 3 1ans =1 2 -37 8 -34 -2 1ans =-1 2 51 2 162 -2 7ans =7 10 -214 21 -33 5 0ans =-3.0000 2.0000 -2.00002.0000 -1.0000 1.0000-1.0000 1.0000 -2.0000 2：（1）：B=[2 4 3];B'ans =243（2）：A=[1 2 3];B=[2 4 3];A.*B,B.*Aans =2 8 9ans =2 8 93：（1）：A=[0 1 0;1 0 0;0 0 1];B=[1 0 0;0 0 1;0 1 0];C=[1 -4 3;2 0 -1;1 -2 0];A^(-1),B^(-1),X=A^(-1)*C*B^(-1) ans =0 1 01 0 00 0 1ans =1 0 00 0 10 1 0X =2 -1 01 3 -41 0 -2（2）：>> A=[1 2 3;2 2 3;3 5 1];B=[1 0 0;2 0 0;3 0 0];A^(-1),x=A^(-1)*Bans =-1.0000 1.0000 0.00000.5385 -0.6154 0.23080.3077 0.0769 -0.1538x =1 0 00 0 00 0 0练习 4.41：（1）：A=[4 2 -1;3 -1 2;11 3 0];b=[2;10;8];B=[A,b];rank(A),rank(B)ans =2ans =3（2）：A=[2 1 -1 1;3 -2 1 -3;1 4 -3 5];b=[1;4;-2];B=[A,b];rank(A),rank(B)ans =2ans =2（3）：A=[ 1 1 1 1; 1 2 -1 4;2 -3 -1 -5;3 1 2 11];b=[5;-2;-2;0];B=[A,b];rank(A),rank(B)ans =4ans =4（4）：A=[ 1 1 2 -1; 2 1 1 -1;2 2 1 2];b=[0;0;0];B=[A,b];rank(A),rank(B)ans =3ans =32：syms a;A=[-2 1 1;1 -2 1;1 1 -2];b=[-2;a;a^2];B=[A,b];rank(A),rank(B)ans =2ans =3练习4.51：（1）：A=[0 1;-1 0];[a,b]=eig(A)a =0.7071 0.70710 + 0.7071i 0 - 0.7071ib =0 + 1.0000i 000 - 1.0000i（2）：A=[0 0 1;0 1 0;1 0 0];[a,b]=eig(A)a =0.7071 0.7071 00 0 -1.0000-0.7071 0.7071 0b =-1 0 00 1 00 0 1（3）：A=[4 1 -1;3 2 -6;1 -5 3];[a,b]=eig(A)a =0.0185 -0.9009 -0.3066-0.7693 -0.1240 -0.7248-0.6386 -0.4158 0.6170b =-3.0527 0 00 3.6760 00 0 8.3766（4）：A=[1 1 1 1;1 1 -1 -1;1 -1 1 -1;1 1 -1 1];[a,b]=eig(A)a =0.5615 0.3366 0.2673 -0.7683-0.5615 -0.3366 0.0000 -0.0000-0.5615 -0.3366 -0.5345 -0.6236-0.2326 0.8125 0.8018 -0.1447b =-1.4142 0 0 00 1.4142 0 00 0 2.0000 00 0 0 2.0000（5）：A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10];[a,b]=eig(A)a =0.8304 0.0933 0.3963 0.3803-0.5016 -0.3017 0.6149 0.5286-0.2086 0.7603 -0.2716 0.55200.1237 -0.5676 -0.6254 0.5209b =0.0102 0 0 00 0.8431 0 00 0 3.8581 00 0 0 30.2887（6）：A=[5 6 0 0 0;1 5 6 0 0 ;0 1 5 6 0 ;0 0 1 5 6; 0 0 0 1 5 ]; [a,b]=eig(A)a =0.7843 -0.7843 -0.9860 -0.9237 -0.92370.5546 0.5546 0.0000 0.3771 -0.37710.2614 -0.2614 0.1643 -0.0000 0.00000.0924 0.0924 0.0000 -0.0628 0.06280.0218 -0.0218 -0.0274 0.0257 0.02579.2426 0 0 0 00 0.7574 0 0 00 0 5.0000 0 00 0 0 2.5505 00 0 0 0 7.4495 2：（1）：A=[0 1;-1 0];[a,b]=eig(A)a =0.7071 0.70710 + 0.7071i 0 - 0.7071ib =0 + 1.0000i 00 0 - 1.0000i>> P=orth(a),B=P'*A*P,P*P'P =-0.7071 -0.70710 - 0.7071i 0 + 0.7071iB =0 + 1.0000i 0 - 0.0000i0 - 0.0000i 0 - 1.0000ians =1.0000 0 + 0.0000i0 - 0.0000i 1.0000>> inv(a)*A*a0 + 1.0000i 000 - 1.0000i3：（1）：A=[2 0 0;0 3 2;0 2 3]; [a,b]=eig(A)a =0 1.0000 0-0.7071 0 0.70710.7071 0 0.7071b =1.0000 0 00 2.0000 00 0 5.0000>> P=orth(a),B=P'*A*P,P*P'P =-1.0000 0 -0.00000.0000 0.7071 0.7071-0.0000 -0.7071 0.7071B =2.0000 0.0000 0.00000.0000 1.0000 00.0000 0 5.0000ans =1.0000 -0.0000 0.0000-0.0000 1.0000 -0.00000.0000 -0.0000 1.0000（2）：A=[1 1 0 -1;1 1 -1 0;0 -1 1 1;-1 0 1 1];[a,b]=eig(A)a =-0.5000 0.7071 0.0000 0.50000.5000 -0.0000 0.7071 0.50000.5000 0.7071 0.0000 -0.5000-0.5000 0 0.7071 -0.5000 b =-1.0000 0 0 00 1.0000 0 00 0 1.0000 00 0 0 3.0000 >> P=orth(a),B=P'*A*P,P*P'P =-0.5000 -0.4998 -0.4783 -0.52100.5000 -0.4822 0.5212 -0.49580.5000 0.4998 -0.4964 -0.5037-0.5000 0.5175 0.5031 -0.4786 B =-1.0000 0.0000 0.0000 0.00000.0000 2.9988 -0.0362 0.03440.0000 -0.0362 1.0007 -0.00060.0000 0.0344 -0.0006 1.0006 ans =1.0000 0.0000 0.0000 -0.00000.0000 1.0000 -0.0000 00.0000 -0.0000 1.0000 0.0000-0.0000 0 0.0000 1.0000练习5.3 1: [m,v]=unifstat(1,11)m =6v =8.33332:[m,v]=normstat(0,16)m =v =256>> s=sqrt(v)s =163:x=randn(200,6);s=std(x)s =0.9094 0.9757 0.9702 0.9393 0.9272 1.09824: x=normrnd(0,16,300,1);hist(x,10)练习 5.61：x=[352 373 411 441 462 490 529 577 641 692 743];y=[166 153 177 201 216 208 227 238 268 268 274];plot(x,y,'*')4：（1）：x=[10 10 10 15 15 15 20 20 20 25 25 25 30 30 30];y=[25.2 27.3 28.7 29.8 31.1 27.8 31.2 32.6 29.7 31.7 30.1 32.3 29.4 30.8 32.8]; plot(x,y,'*')。
数学建模实验二：微分方程模型Matlab求解与分析
实验二： 微分方程模型Matlab 求解与分析一、实验目的[1] 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析； [2] 熟悉MATLAB 软件关于微分方程求解的各种命令；[3] 通过范例学习建立微分方程方面的数学模型以及求解全过程； [4] 熟悉离散 Logistic 模型的求解与混沌的产生过程。
二、实验原理1. 微分方程模型与MATLAB 求解解析解用MATLAB 命令dsolve(‘eqn1’,’eqn2’, ...) 求常微分方程（组）的解析解。
其中‘eqni'表示第i 个微分方程，Dny 表示y 的n 阶导数,默认的自变量为t 。
（1） 微分方程 例1 求解一阶微分方程 21y dxdy+= (1) 求通解 输入：dsolve('Dy=1+y^2')输出：ans =tan(t+C1)（2）求特解 输入：dsolve('Dy=1+y^2','y(0)=1','x')指定初值为1，自变量为x 输出：ans =tan(x+1/4*pi)例2 求解二阶微分方程 221()04(/2)2(/2)2/x y xy x y y y πππ'''++-=='=-原方程两边都除以2x ，得211(1)04y y y x x'''++-= 输入:dsolve('D2y+(1/x)*Dy+(1-1/4/x^2)*y=0','y(pi/2)=2,Dy(pi/2)=-2/pi','x')ans =- (exp(x*i)*(pi/2)^(1/2)*i)/x^(1/2) +(exp(x*i)*exp(-x*2*i)*(pi/2)^(3/2)*2*i)/(pi*x^(1/2))试试能不用用simplify 函数化简 输入: simplify(ans)ans =2^(1/2)*pi^(1/2)/x^(1/2)*sin(x) （2）微分方程组例3 求解 d f /d x =3f +4g ; d g /d x =-4f +3g 。
数学实验 实验一ans
电 子 科 技 大 学实 验 报 告学生姓名： 学 号： 指导教师： 实验时间: 报告评分： 一、实验室名称：应用数学学院数学实验室 二、实验项目名称：抛射曲线实验 三、实验原理：炮弹飞行的简单数学模型为抛射曲线，其参数方程为⎪⎩⎪⎨⎧-⨯=⨯=20021sin cos gt t v y t v x αα 其中v 0为初始速度，α 为发射角，g 是重力加速度（约为 9.8米/秒2）。
四、实验目的：熟悉MATLAB 的集成环境，掌握MATLAB 的基本运算和数学表达式计算方法。
五、实验内容：对已知射程数据，为使炮弹击中目标需计算发射角；反之对已知的发射角数据，也需要计算射程。
设v 0 = 515米/秒，对不同目标计算射程、发射角和炮弹飞行时间。
六、实验器材（设备、元器件）：台式计算机 七、实验步骤及操作：弹道数据包括射程、发射角及炮弹飞行时间。
令参数方程中y 的表达式为零可确定炮弹飞行时间 t (α )，代入 x 的表达式可导出射程计算公式。
（1） 炮弹飞行时间计算公式：gv t ααsin 2)(0=（2） 射程计算公式：ααα2sin sin 2cos 2001gv g v v x =⨯=（3） 发射角计算公式：)/arcsin(21201v gx =α八、实验数据及结果分析：(2分)注意：红色字是填写的数据，黑色字是实验报告原有的数据。
教师参考程序如下%----抛射曲线仿真程序%----2007/01/24g=9.8;v0=515;distance=inline('515*515*sin(2*dalpha)/9.8');alpha=inline('0.5*asin(9.8*s/(515^2))');times=inline('2*515*sin(alpha)/9.8');results1(1,:)=[10,20,27];a=alpha(10000); results1(2,1)=a*180/pi;results1(3,1)=times(a);a=alpha(20000); results1(2,2)=a*180/pi;results1(3,2)=times(a);a=alpha(27000); results1(2,3)=a*180/pi;results1(3,3)=times(a);results1results2(2,:)=[45,60];results2(1,1)=distance(45*pi/180)/1000;results2(1,2)=distance(60*pi/180)/1000;results2(3,1)=times(45*pi/180);results2(3,2)=times(60*pi/180);results2九、实验结论：（1分）十、总结及心得体会：（1分）十一、对本实验过程及方法、手段的改进建议：（1分）。
数学实验2 班级 学号 姓名 一.用MATLAB计算下列极限： (1); (2
>> [x,y,f]=fsolve('eg2_2fun',[0.1,0.1])
回车显示
y =
-0.2484 -0.3676
y =
-0.2484 -0.3676
y =
-0.2484 -0.3676
y =
0.0249 0.0530
y =
0.0249 0.0530
y =
0.0249 0.0530
ans =
exp(m*n)
（6）
clear
>> syms y m n k
>> limit((1+m*y)^(n/y+k),y,0)
ans =
exp(m*n)
（7）
clear
>> syms x m n
>> limit(sin(m*x)/tan(n*x),x,0)
ans =
m/n
（8）
syms x m n
y =
1.0e-003 *
0.3389 0.2246
y =
1.0e-003 *
0.3389 0.2246
y =
1.0e-003 *
0.3389 0.2246
y =
1.0e-007 *
0.4429 0.7265
y =
1.0e-007 *
0.5281 0.7649
y =
1.0e-007 *
0.4575 0.9178
用MATLAB计算极限的命令语句如下：
clear
>> syms x y m n（生成符号变量 ）
limit(f(x),x,a)（求 ）
实验二：定积分的近似计算
0 1 x2
解： a=0, b=1, n=100, f (x) = 1/( 1+x2 )
==> h =1/100=0.01, xi = i*h, yi = f (xi)
==>
1 dx 0 1 x2
h
y0 2
y1
yn1
yn 2
0.78539399673078
相对误差： 0.78539399673078 / 4 5.305 10-6 /4
>> quad('1./(1+x.^2)',0,1,10e-16)
函数表达式一定要用 单引号 括起来！ 涉及的运算一定要用 数组运算！
2021/4/21
22
二重积分的计算
抛物线法计算二重积分： dblquad
dblquad(f,a,b,c,d,tol)
bd
a c f (x, y)dxdy
tol 为计算精度，若不指定，则缺省精度为 10-6
1
中点法：
0
1
dx x
2
n
h
i 1
f
(
xi1 2
xi
)
0.78540024673078
2021/4/21
6
矩形法举例
误差分析
理论值：
1 dx 0 1 x2
arctan x
1 0
π 4
左点法相对误差：0.78789399673078 / 4 0.003178 /4
右点法相对误差：0.78289399673078 / 4 0.003188 /4
i 1
x
2021/4/21
8
梯形法
曲边小梯形的面积可以由直边小梯形的面积来近似
高等数学实验2 微分、积分(含答案)
班级 学号 姓名高等数学实验2 微分、积分一. 用MA TLAB 计算下列导数：diff 函数（1）已知2xy e =，求y '、y ''、(10)y 。
（2）已知nx y e =，求y '''。
（3）已知210x y xe-=,求y '、y ''与(8)y 。
(4)设2sin ()43x f x x x =++，求()f x '、()f x ''及()6f π''。
二．用MA TLAB 解方程。
solve 函数1．一元方程与线性方程（组）（1） 解方程 062=--x x（2）解方程组⎪⎩⎪⎨⎧=-+=-+060622x y y x （3）解方程组⎪⎪⎩⎪⎪⎨⎧=+++=++=++-=++012412324543213214321431x x x x x x x x x x x x x x2．非线性方程（组）（4）解非线性方程组⎩⎨⎧=+-=--0sin 3.0cos 5.00cos 3.0sin 5.0212211x x x x x x 三。
用MA TLAB 计算极值：(1)已知销售额R 是价格P 的函数，且200184R P P ⎛⎫=-⎪+⎝⎭。
当价格P 为何值时， 销售额R 有最大值，且求此最大值。
(2)已知某公司收益函数210xR xe -=，成本函数32(1085)/100C x x =++，其中x 为产（销）量，求最大收益、最低平均成本和最大利润。
四．用MATLAB 计算下列不定积分 int 函数1．ln xdx ⎰； 2。
321x x e dx -⎰； 3． 42(31)sin(21)x x x dx -++⎰； 4．(sin sin cos )ax bx cx dx ⨯⨯⎰； 5．（练习）5(4)ln(32)x x x dx --⎰； 6．（练习）4sin(25)x x e dx +⎰；五．用MATLAB 计算下列定积分 int 函数1．120(1)x xe dx x +⎰ 2。
MATLAB数学实验第二版课后练习题含答案
MATLAB数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下：第一章课后练习题1.编写MATLAB程序，计算并输出下列公式的结果：y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中，x, $\\mu$, $\\sigma$ 分别由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序，求解下列方程的解：4x + y = 11\\\\x + 2y = 7答案如下：A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序，计算下列多项式的值：P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中，x 由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序，绘制下列函数的图像：f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下：x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。
数学实验第二次测验题及参考答案（09级）
数学实验第二次测验题及参考答案（09级）数学实验第二次测验题及参考答案一、写出下列MATLAB指令的运算结果.1. A=[1;2;3]; transpose(A)1 2 31 2 32. A=[1, 2, 3 ; 4, 5, 6 ]; B=A([1 2], [1,3]) , d=size(A)B =1 34 6d =2 33. a=1:3; b=linspace(1,3,3); x=sum(a.*b), y=cross(a, b)x = 14y = 0 0 04. A=[1,2, 3; 4,5,6; 7,8,9]; B=ones(3); C=A-BC =0 1 23 4 56 7 85. v=[1, 2, 3]; A=diag(v); E=eig(A), D=det(A)E=123D =66. x=[1,2,3,4,5]; [mean(x), median(x), range(x),sum(x), prod(x)]3 34 15 1207. x=[2,3,4]; a=cumsum(x) ,b=sort(x)a =2 5 9b =2 3 48．format rat; v=[1, 2, 3]; A=diag(v); inv(A)ans =1 0 00 1/2 00 0 1/39. [m,v]=normstat(1,4) % 求参数为1,4的正态分布的均值与方差m =1, v =16二、写出下列MATLAB指令的实验目的.1. dsolve('x*Dy+y-exp(-x)=0', 'y(1)=2*exp(1)', 'x')求微分方程0=-+'-x e y y x 在初始条件e y x 2|1==下的特解.2. u=[1,2,3],v=[0,3,2], w=[5, 2, 1]; dot(w, cross(u, v))计算向量u, v, w 的混合积.3. A=[1 2 3; 2 2 5; 3 5 1]; b=[1;2;3]; det(A); inv(A)*b利用逆矩阵解线性方程组=++=++=++3532522132321 321321x x x x x x x x x .4. A=[0 0 1; 0 1 1; 1 1 1; 1 0 0]; rref(A)求向量组)1,0,0(1=α,)1,1,0(2=α,)1,1,1(3=α,)0,0,1(4=α的秩.或对矩阵A 做行初等变换。
数学实验2课件市公开课金奖市赛课一等奖课件
其中
Ra=y
r1 ( x1 ) rm ( x1 )
R
,
r1 ( xn ) rm ( xn )
a1
a
,
am
（3）
y1
y
yn
定理: 当RTR可逆时，超定方程组（3）存在最小二乘解， 且即为方程组
RTRa=RTy
解: a=(RTR)-1RTy
第13页
设 y*=a+bxi , 令 δi=yi-y*i=yi-a-bxi, 依据最 小二乘原理, 即使误差平方和达到最小, 也就是令
第22页
用MATLAB作非线性最小二乘拟合 Matlab(优化工具箱)提供了两个求非线性最小二 乘拟合函数: lsqcurvefit和lsqnonlin。 两个命令都要先建立M-文献定义函数, 然后编写 主程序M-文献求解非线性模型未知参数
注意: 二者定义f(x)方式是不同.
第23页
1. lsqcurvefit 已知数据点: xdata=（xdata1，xdata2，…，xdatan），
1.线性最小二乘拟合 2.非线性最小二乘拟合 3.cftool(自学)
ployfit \
lsqcurvefit lsqnonlin
第18页
多项式在x处值y可用下列命令计算: y=polyval（a，x）
用MATLAB作线性最小二乘拟合
1.作多项式f(x)=a1xm+ …+amx+am+1拟合,可利用已有
第28页
课堂小结
x0赋值:注意单位/依据结果重新赋值 ployfit(x,y,m)中x,y赋值: y换行要加… 画图: 非线性模型画出直线
非线性线性化后没有还原
第29页
初中数学实验
初中数学实验数学实验是初中数学教学中一种重要的教学活动，它通过实际操作和观察，培养学生的数学思维和动手能力，提高他们对数学概念的理解和运用能力。
下面是初中数学实验的相关参考内容。
一、实验目的和基本要求1.理解实验目的，明确实验内容，掌握实验步骤和操作方法。
2.能够准确记录实验过程中的数据和观察结果，并能进行数据分析和归纳。
3.能够总结实验的结论，归纳实验中涉及的数学规律和定理。
4.通过实验的操作和观察，发现问题并进行解决，培养数学探究的能力。
二、实验内容和步骤1.实验一：测量线段的长度实验目的：了解测量线段长度的技巧和方法。
实验步骤：准备测量工具（直尺、卷尺等），选取不同长度的线段进行测量。
实验要点：要保持工具的水平、直角，保证测量的准确性；对于过长或过短的线段，应该采取恰当的测量方法。
2.实验二：探究平均数的性质实验目的：通过实验观察，探究平均数的性质。
实验步骤：准备一组数值（可以是学生的身高、体重等），求出其平均数，并观察平均数与原始数值的关系。
实验要点：要确保数值的准确性，使用适当的计算方法，注意观察数值之间的规律。
3.实验三：研究三角形中角的和实验目的：通过实际操作，研究三角形内角和为180°的性质。
实验步骤：利用直角三角板或其他测量工具，制作不同形状的三角形，并测量三角形的内角和。
实验要点：要仔细观察实验数据，总结内角和与三角形形状的关系，归纳出三角形内角和为180°的规律。
三、实验记录和分析1.实验记录：在实验过程中，要详细记录实验的步骤、数据和观察结果，以便于后续的数据分析和结论总结。
2.数据分析：对实验中得到的数据进行统计和分析，找出规律和问题，并进行初步的解释。
3.结果结论：根据实验的数据分析，总结出实验的结论，提炼出涉及的数学规律和定理。
四、实验总结和思考1.实验总结：对整个实验过程进行总结，明确实验的目的和意义，评价实验的优缺点。
2.思考问题：提出实验中的问题或困惑，并尝试进行解答或寻找解决方法。
数学实验2
1.实验内容①数值积分法计算π的近似值；Taylor级数法计算π的近似值；②Monte Carlo法计算π的近似值；2.实验目的利用计算机，结合数学理论、数学软件，计算无理数π的近似值。
3.实验要求①选择函数使其在指定区间上的积分与π相关，利用数值积分法计算积分的近似值，得到π的近似值；研究不同的函数、不同的数值积分方法对π的计算精度的影响。
n20000;y x_:11x2;s14Sum y kn,k,1,n1y0y12s24y0y12Sum y kn,k,1,n14Sum y 2k12n,k,1,n16n;Print N s1,50Print N s2,50Print N,503.14159265358979323846264338326981476415748685972333.1415926535897932384626433832795028841971693993751由以上程序结果比较可知：由辛普森公式计算所得的π值与真值比较精确到小数点后第28位；由梯形公式所得的π值与真值比较精确到小数点后第9位；即由辛普森公式计算所得的π值比由梯形公式所得的π值更精确一点。
②选择适当的函数展成Taylor级数，在函数值与π有关的点计算此级数的有限项和进而得到π的近似值，控制误差并比较不同公式在计算精度上的差异。
Clear x,nT x_,n_:Sum1k x2k12k1,k,0,N4T1,20000,50TimingPrint N4T 12,200T13,200,17Print N16T 15,2004T1239,200Print N,170{0.671Second,3.14164265108988698690008478855370942654501638409 78}3.141592653589793238462643383279502884197169399375105820 97494459230781640628620899862803482534211706798214808651 32823066470942293662986845361440202058819980233921997399 8343.141592653589793238462643383279502884197169399375105820 97494459230781640628620899862803482534211706798214808651 32823066470938446095505822317253594081284811174502841027 0193.141592653589793238462643383279502884197169399375105820 97494459230781640628620899862803482534211706798214808651 32823066470938446095505822317253594081284811174502841027 019Arctan[x]在x=1处的展开值与真值比较精确到小数点后第3位；Arctan[x]在x={1/2,1/3}处的展开值与真值比较精确到小数点后第122位Arctan[x]在x={1/5,1/239}处的展开值与真值比较精确到小数点后第170位以上由此可得，Arctan[x]在x={1/5,1/239}处的展开值最精确。
数学实验报告mathematics
数学实验报告mathematics实验名称: 探究二次函数的特性摘要:本实验主要通过构建和探究二次函数的图像来研究其特性。
实验使用了数学软件进行模拟，并记录了函数的图像和相应的特性。
实验结果表明，二次函数的图像是一个抛物线，其开口的方向和顶点的位置可以通过函数的系数来确定。
引言:二次函数是高中数学中重要的一种函数类型。
了解二次函数的特性对于理解和解决实际问题具有重要意义。
本实验旨在通过构建二次函数的图像，研究其特性，包括顶点、开口方向和对称轴等。
材料与方法:1. 使用数学软件（如Geogebra）创建一个二次函数的图像。
2. 调整二次函数的系数，观察图像的变化。
3. 记录每次调整后的图像特性，如顶点、开口方向和对称轴等。
4. 比较不同系数对图像的影响。
结果与讨论:通过调整二次函数的系数，我们观察到以下结果:1. 系数a的正负决定了二次函数的开口方向。
当a>0时，图像开口向上; 当a<0时，图像开口向下。
2. 顶点的位置可以通过函数的系数b和c来确定。
顶点的横坐标为-x = -b/2a，纵坐标为y = f(-b/2a)。
3. 对称轴是通过顶点且垂直于x轴的直线。
在二次函数的图像中，顶点和对称轴是对称的。
4. 当系数a的绝对值较小时，图像趋于扁平化，开口较宽; 当系数a的绝对值较大时，图像趋于瘦长，开口较窄。
结论:通过本实验，我们深入了解了二次函数的特性。
我们发现二次函数的图像是一个抛物线，其开口方向、顶点位置和对称轴可以通过函数的系数来确定。
这些特性对于解决实际问题和更深入地理解数学概念都具有重要意义。
建议与展望:本实验仅研究了二次函数的基本特性，未涉及其在实际问题中的应用。
进一步的研究可以探讨二次函数在物理学、经济学和工程学等领域的具体应用，并进一步深入研究其特性与实际问题的关联。
大学数学实验教程第二版成丽波课后答案
大学数学实验教程第二版成丽波课后答案1、的值为（）[单选题] *A.-2B. 0C. 1(正确答案)D. 22、15．下列说法中，正确的是（）[单选题] *A．若AP＝PB，则点P是线段AB的中点B．射线比直线短C．连接两点的线段叫做两点间的距离D．过六边形的一个顶点作对角线，可以将这个六边形分成4个三角形(正确答案)3、1.计算-20+19等于（）[单选题] *A.39B.-1(正确答案)C.1D.394、4．﹣3的相反数是（）[单选题] *A.BC -3D 3(正确答案)5、下列是具有相反意义的量是（）[单选题] *A．身高增加1cm和体重减少1kgB．顺时针旋转90°和逆时针旋转45°(正确答案)C．向右走2米和向西走5米D．购买5本图书和借出4本图书6、代数式a3?a2化简后的结果是（）[单选题] *A. aB. a?(正确答案)C. a?D. a?7、若tan（π-α）>0且cosα>0，则角α的终边在（）[单选题] *A.第一象限B.第二象限C.第三象限D.第四象限(正确答案)8、50．式子（2+1）（22+1）（24+1）（28+1）…（21024+1）+1化简的结果为（）[单选题] *A．21024B．21024+1C．22048(正确答案)D．22048+19、60°用弧度制表示为（）[单选题] *π/3(正确答案)π/62π/32π/510、48、如图，△ABC≌△AED，连接BE．若∠ABC＝15°，∠D＝135°，∠EAC＝24°，则∠BEA的度数为（）[单选题] *A．54°B．63°(正确答案)C．64°D．68°11、30°角是（）[单选题] *A、第一象限(正确答案)B、第一象限C、第三象限D、第四象限12、已知二次函数f（x）=2x2-x+2，那么f（2）的值为（）。
高二数学独立重复试验2
保留两个有效数字）；
0.68
（3）最有可能击中目标几次？ 8次
例2 某车间有5台机床，在1小时内每 台机床需要工人照管的概率都是0.25， 求在1小时内这5台机床中至少有2台需要 工人照管的概率.（结果保留两个有效数 字）
0.37
小结作业
1.在独立重复试验中，若每次试验结 果只有事件A发生或不发生两种可能，则 事件A发生的次数服从二项分布；若每次 试验结果有多种可能，则可以根据需要 适当设定事件A，将其转化为二项分布.
；公司起名 起名网 https:/// 公司起名 起名网；
探究（一）：独立重复试验
思考1：在同等条件下，将一枚硬币重复
抛掷100次，记Ai(i＝1，2，…，100)表 示“第i次抛掷硬币正面朝上”，那么事
件A1，A2，…，A100两两之间是否相互独
立？
相互独立
思考2：在同等条件下，某射手连续射击
思考4：投掷一枚图钉，设针尖向上的概
率为p，连续投掷3次，则仅出现1次针尖
向上有哪几种情形？如何计算仅出现1次 针尖向上的概率？
记Ai(i＝1，2，3)表示第i次投掷针尖 向上，则
Hale Waihona Puke 理论迁移例1 某射手每次射击击中目标的概率
都是0.8，若这名射手射击10次，求
（1）恰有8次击中目标的概率；0.3
（2）至少有8次击中目标的概率.（结果
2.2 二项分布及其应用
2.2.3 独立重复试验与二项分布
问题提出
t
p
1 2
5730
1.事件A与事件B相互独立的充要条 件是什么？
事件A与B相互独立 P(AB)＝P(A)P(B)
2.若事件A1，A2，…，An两两之间相 互独立，则P(A1A2…An)等于什么？
数学实验2 班级 学号 姓名 一.用MATLAB计算下列极限： (1); (2
ans =
m/n
（8）
syms x m n
>> limit((tan(m*x)+x)/sin(n*x),x,0)
ans =
(m+1)/n
二．用MATLAB解方程。
1．一元方程与线性方程组
用两个函数solve ( )和linsolve( )来解线性方程，具体格式如下：
>> X=linsolve(A,B)
X =
[ 1]
[ -1]
[ -1]
[ 1]
说明：X=zeros(4,1)表示建立一个4元素列向量X。
1．非线性方程（组）
用命令函数fsolve（）来解非线性方程（组）。具体格式如下：
X=fsolve(fun,x0,options)
参数fun为定义好的待求解的非线性方程（组）的文件名。x0为求解方程的初始向量或矩阵。Options设置命令函数fsolve求解过程的各种参数。通常我们设为optimset（fsolve）。Options的其他参数项可通过帮助文档查询，或直接在MATLAB的命令行输入help optimset查询。
limit(f(x),x,a,'left')（求 ）
例1．用MATLAB计算下列极限：
（1） ；（2） ；
（3） ；（4） ；
（5） ；（6） ；
（7） ；（8） 。
解（1）
clear
>> syms x
>> limit((exp(2*x)-1)/x,x,0)
ans =
2
（2）clear
>> syms x k
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
电子科技大学
实验报告
学生姓名：学号：指导教师：
实验时间: 报告评分：
一、实验室名称：应用数学学院数学实验室
二、实验项目名称：飞行航程计算实验
三、实验原理：
已知球面上点P的经度φ和纬度θ，转换为直角坐标用坐标转化公式
x = R cosθcosφ
y = R cosθsinφ
z = R sinθ
其中，R=6400（km），θ∈[–π /2，π /2]，φ∈[–π，π ]。
过球面上两点P1（x1，y1，z1），P2（x2，y2，z2），的短程线计算公式为：L = R×α。
其中，α是夹角（弧度），满足cos α= (x1，y1，z1)·(x2，y2，z2) / ( ||(x1，y1，z1)|| ||(x2，y2，z2)|| )
四、实验目的：
了解地球上各大城市间的飞行航程的简单计算方法。
五、实验内容：
计算飞行航程：北京——上海——东京——旧金山——纽约
六、实验器材（设备、元器件）：台式计算机
七、实验步骤及操作：
1．利用坐标转化公式将经纬度数据转换为直角坐标；
2．利用球面上过P1和P2两点的短程线长度计算公式顺序计算出各区间的飞行航程数据；3．累加各区间飞行航程数据，写出总航程数据；
4．计算新航线与旧航线所用时间差。
八、实验数据及结果分析：
教师参考程序
%----航程计算程序
%----2008-03-20
city=[40,118;31,122;36,140;37,-123;41,-76];
R=6400+10;
theta=city(:,1)*pi/180;
fai=city(:,2)*pi/180;
x=R*cos(theta).*cos(fai);
y=R*cos(theta).*sin(fai);
z=R*sin(theta);
op1=[x(1:4),y(1:4),z(1:4)];
op2=[x(2:5),y(2:5),z(2:5)];
D=R*acos(op1*op2'/R^2);
format bank
Dnew=D(1,4),Dold=diag(D)'
distance=sum(Dold)
time0=distance/980
time1=Dnew/980
Dtime=time0-time1
format
九、实验结论：
十、总结及心得体会：
十一、对本实验过程及方法、手段的改进建议：。