MATLAB数学实验第二版答案胡良剑.docx
MATLAB数学实验第二版答案胡良剑数学实验答案Chapter1Page20,ex1(5)等于exp(1),exp(2);exp(3),exp(4)(7)3=1*3,8=2*4(8)a为各列最小值，b为最小值所在的行号(10)1>=4,false,2>=3,false,3>=2,ture,4>=1,ture(11)答案表明：编址第2元素知足不等式(30>=20)和编址第4元素知足不等式(40>=10)(12)答案表明：编址第2行第1列元素知足不等式(30>=20)和编址第2行第2列元素知足不等式(40>=10)Page20,ex2(1)a,b,c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b(2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码Page20,ex3>>r=2;p=0.5;n=12;>>T=log(r)/n/log(1+0.01*p)Page20,ex4>>x=-2:0.05:2;f=x.4-2.x;>>fmin,min_index=min(f)最小值最小值点编址>>x(min_index)ans=0.6500最小值点>>f1,x1_index=min(abs(f)求近似根-绝对值最小的点f1=0.0328x1_index=24>>x(x1_index)ans=-0.8500>>x(x1_index)=;f=x.4-2.x;删去绝对值最小的点以求函数绝对值次小的点>>f2,x2_index=min(abs(f)求另一近似根-函数绝对值次小的点f2=0.0630x2_index=65>>x(x2_index)ans=1.2500Page20,ex5>>z=magic(10)z=929918156774515840988071416735557644148188202254566370478587192136062697128869325296168755234172476839042492633652358289914830323966796139597293138457210129496783537444653111810077843643502759>>sum(z)>>sum(diag(z)>>z(:,2)/sqrt(3)>>z(8,:)=z(8,:)+z(3,:)Chapter2Page45ex1先在编辑器窗口写下列M函数，保存为eg2_1.mfunctionxbar,s=ex2_1(x)n=length(x);xbar=sum(x)/n;当前位置：文档视界MATLAB数学实验第二版答案胡良剑MATLAB数学实验第二版答案胡良剑k=k+1;F(k)=F(k-1)+F(k-2);x=F(k)/F(k-1);enda,x,k计算至k=21可知足精度Page45ex4clear;tic;s=0;fori=1:1000000s=s+sqrt(3)/2i;ends,toctic;s=0;i=1;whileii=1:1000000;s=sqrt(3)*sum(1./2.i);s,tocPage45ex5t=0:24;c=15141414141516182022232528.313231292725242220181716;plot(t,c)Page45ex6(1)x=-2:0.1:2;y=x.2.*sin(x.2-x-2);plot(x,y)y=inline('x2*sin(x2-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)(6)theta=linspace(0,2*pi,50);fai=linspace(0,pi/2,20);theta,fai=meshgrid(theta,fai);x=2*sin(fai).*cos(theta);y=2*sin(fai).*sin(theta);z=2*cos(fai);surf(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)page45,ex7x=-1.5:0.05:1.5;y=1.1*(x>1.1)+x.*(x=-1.1)-1.1*(xp=a*exp(-0.75*y.2-3.75*x.2-1.5*x).*(x+y>1);p=p+b*exp(-y.2-6*x.2).*(x+y>-1).*(x+y>A=123;456;780;C=2-5-22;-5-24-56;-22-56-16;>>X=lyap(A,C)X=1.0000-1.0000-0.0000-1.00002.00001.0000-0.00001.00007.0000Chapter3Page65Ex1>>a=1,2,3;b=2,4,3;a./b,a.b,a/b,abans=0.50000.50001.0000ans=221ans=0.6552一元方程组x2,4,3=1,2,3的近似解ans=0000000.66671.33331.0000矩阵方程1,2,3x11,x12,x13;x21,x22,x23;x31,x32,x33=2,4,3的特解Page65Ex2(1)>>A=41-1;32-6;1-53;b=9;-2;1;>>rank(A),rank(A,b)A,b为增广矩阵ans=
