915230MATLAB第三版实验答案Word下载.docx

1、lsdhKSDLKklsdkl;k=find（A=A&AZA（k）=实验2 MATLAB矩阵处理%第1题E=eye（3）;R=rand（3,2）;O=zeros（2,3）;S=diag（2,3）;A=E,R;O,S;A2B=E,（R+R*S）;O,S2%第2题A=fix（10*rand（5）H=det（A）Trace=trace（A）Rank=rank（A）Norm=norm（A）H=hilb（5）P=pascal（5）Hh=det（H）Hp=det（P）Th=cond（H）Tp=cond（P）A=-29,6,18;20,5,12;-8,8,5V,D=eig（A）%数学意义略%第5题方法一%（

2、1）:A=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6;b=0.95,0.67,0.52x=inv（A）*b%（2）:B=0.95,0.67,0.53x=inv（A）*B%（3）:cond（A）%第5题方法二A=hilb（4）A（:,1）=A（4,:）=B=0.95,0.67,0.52X=inv（A）*BB1=0.95,0.67,0.53X1=inv（A）*B1N=cond（B）N1=cond（B1）Na=cond（A） %矩阵A为病态矩阵实验3 顺序结构程序设计%程序（略）实验4 选择结构程序设计%第1题程序一x=-5.0,-3.0,1.0,2.0,2.5,3.0,

3、5.0;y=; %建立存放所有y值的矩阵for x0=x if x0=0&x05&x0=2&x0=3 y=y,x0*x0-5*x0+6; else y=y,x0*x0-x0-1; endendx %输出所有xy %输出所有y%第1题程序二x=-5,-3,1,2,2.5,3,5;for a=1:7if x（a）x（a）100|xdisp（您输入的成绩不是百分制成绩，请重新输入。if x=90elseif x=89&=80B=79&=70C=69&60DE%第2题程序二s=input（请输入一个成绩（0分到100分之间）： %s用于存放成绩while 1 %判断输入成绩的合理性 if s100 d

4、isp（输入的成绩需在0到100之间，请重新输入：） s=input（ break;switch fix（s/10） %对成绩做出等级判断 case 9,10 case 8 case 7 case 6 otherwisen=input（请输入员工工号：h=input（该员工工作时数是：if h120x=（h-120）*84*（1+0.15）+120*84;elseif h5超出了矩阵的行数，矩阵的最后一行为： a（5,: disp（矩阵的第,num2str（n）,行为：） a（n,:实验5 循环结构程序设计y=0;n=?for i=1:n y=y+1/（2*i-1）/3（2*i-1）;ylog

5、（2）/2i=1:n;y=1./（2*i-1）./3.（2*i-1）;y=sum（y）n=1;while y1e-5 Xn=Xn1; Xn1=a/（b+Xn）; n=n+1; if n=500Xn1r1=（-b+sqrt（b*b+4*a）/2r2=（-b-sqrt（b*b+4*a）/2 if i=1 f（i）=1; elseif i=2 f（i）=0; elseif i=3 f（i）=f（i-1）-2*f（i-2）+f（i-3）;max（f）min（f）sum（f）length（find（f0）length（find（f=0）length（find（fs=0;for i=2:49 b=i*（i

6、+1）-1; m=fix（sqrt（b）; for j=2:m if rem（b,j）=0 break if j=m s=s+b;s实验6 函数文件function y=mat1（x） %建立函数文件mat1.my=exp（x）,log（x）,sin（x）,cos（x）;%在命令窗口调用上述函数文件：y=mat1（1+i）function a,b,N,M=shiyanwu2（m,n,t）A=m*cos（t*pi/180）,-m,-sin（t*pi/180）,0;m*sin（t*pi/180）,0,cos（t*pi/180）,0;0,n,-sin（t*pi/180）,0;0,0,-cos（t*p

7、i/180）,1;B=0,9.8*m,0,9.8*n;C=inv（A）*Ba=C（1）;b=C（2）;N=C（3）;M=C（4）;%在命令窗口调用该函数文件:m1=input（m1=m2=input（m2=theta=input（theta=a1,a2,N1,N2=shiyanwu2（m1,m2,theta）function X=mat2（m1,m2,t）g=9.8;A=m1*cos（t*pi/180）,-m1,-sin（t*pi/180）,0;m1*sin（t*pi/180）,0,cos（t*pi/180）,0;0,m2,-sin（t*pi/180）,0;B=0;m1*g;0;m2*g;X=

8、inv（A）*B;%在命令窗口调用该函数文件：X=mat2（1,1,60）function flag=mat3（x）flag=1;sqrt（x）if rem（x,i）=0 flag=0; break;for i=10:99j=10*rem（i,10）+fix（i/10）;if mat3（i）&mat3（j） disp（i）function y=fx（x）y=1./（x-2）.2+0.1）+1./（x-3）.4+0.01）;y=fx（2）a=1,2;3,4;y=fx（a）function f1=mat5（n）f1=n+10*log（n*n+5）;%在命令窗口中调用该函数文件：y=mat5（40）

9、/（mat5（30）+mat5（20）%（2）方法一function f2=mat6（n）f2=0; f2=f2+i*（i+1）;%在命令窗口中调用该函数文件如：y=mat6（40）/（mat6（30）+mat6（20）%（2）方法二function f2=mat7（n）m=i.*（i+1）;f2=sum（m）;y=mat7（40）/（mat7（30）+mat7（20）实验7 绘图操作x=linspace（0,2*pi,101）;y=（0.5+3*sin（x）./（1+x.2）.*cos（x）;plot（x,y）x=-5:0.01:5; %起始设y为空向量if x0=0 %不能写成x0=0.5

10、p=n/30000A=45+51*rand（100,5）;Y,U=max（A）a,b=min（A）m=mean（A）s=std（A）sum（A,2）Y,U=max（ans）a,b=min（ans）%（4）zcj,xsxh=sort（ans）p1=1,2,4,0,5;p2=1,2;p3=1,2,3;p=p1+0,conv（p2,p3） %为使两向量大小相同，所以补0A=roots（p）A=-1,1.2,-1.4;0.75,2,3.5;0,5,2.5;polyval（p,A）polyvalm（p,A）h=6:2:18;x=6.5:17.5;t1=18,20,22,25,30,28,24;t2=15

11、,19,24,28,34,32,30;T1=spline（h,t1,x）T2=spline（h,t2,x）x=1:101;y1=log10（x）;p=polyfit（x,y1,5）y2=polyval（p,x）;:实验9 数值微分与积分for x=1: fx=x,x2,x3;1,2*x,3*x;0,2,6*x; diff（fx）f1=（x） sqrt（cos（x.*x）+4*sin（2*x）.2）+1）;I1=quad（f1,0,2*pi）f2=（x） log（1+x）./（1+x.2）;I2=quadl（f2,0,1）I2=integral（f2,0,1）%第4题（略）%第5题（略）实验10

12、 方程数值求解A=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2;b=-4;13;1;11;x1=Abx2=inv（A）*bL,U=lu（A）;x3=U（Lb）%（1） z=fzero（x） 3*x+sin（x）-exp（x）,1.5）%首先建立函数文件myfun.m。function F=myfun（X）x=X（1）;y=X（2）;z=X（3）;F（1）=sin（x）+y2+log（z）-7;F（2）=3*x+2y-z3+1;F（3）=x+y+z-5;X=fsolve（myfun,1,1,1,optimset（Display, off）f=（x） （x3+cos（x

13、）+x*log（x）/exp（x）;fminbnd（f,0,1）%建立函数文件fx12.m：function f=fx12（u）x1=u（1）;x2=u（2）;f=2*x1.3+4*x1*x2.3-10*x1.*x2+x2.2;U,fmin=fminsearch（fx12,0,0）%建立函数文件ztf.mfunction yy=ztf（x,y）yy=（5*y（1）-y（2）/x;y（1）;x,y=ode23（ztf,0,5,0;0）%建立函数文件ztfun.mfunction yy=ztfun（t,y）yy=y（2）*y（3）;-y（1）*y（3）;-0.51*y（1）*y（2）;t,y=od

14、e23（ztfun,0,5,0;1）plot3（y（:,1）,y（:,2）,y（:,3）实验11 符号计算对象与符号微积分x=sym（6）;y=sym（5）;z=（x+1）/（sqrt（3+x）-sqrt（y）syms x yz=sym（x*x*x*x-y*y*y*y）;factor（z）num=sym（5135）;factor（num）syms beta1 beta2f1=sym（sin（beta1）*sin（beta2）-cos（beta1）*cos（beta2）;simplify（f1）sym xf2=sym（4*x*x+8*x+3）/（2*x+1）;simplify（f2）syms a

15、 b c d e f g h iP1=0,1,0;1,0,0;0,0,1;P2=1,0,0;0,1,0;1,0,1;A=a,b,c;d,e,f;g,h,i;B=P1*P2*Ainv（B）tril（B）det（B）syms x a tf1=（x*（exp（sin（x）+1）-2*（exp（tan（x）-1）/sin（x）3limit（f1,x,0） %求（1）f2=（sqrt（pi）-sqrt（acos（x）/sqrt（x+1）;limit（f2,x,-1,right） %求（2）f3=（1-cos（2*x）/x;diff（f3,x）diff（f3,x,2） %求（3）A=power（a,x）,

16、t3;t*cos（x）,log（x）;diff（A,x）diff（A,t,2）diff（diff（A,x）,t） %求（4）f=（x*x-2*x）*exp（-x*x-y*y-x*y）;%求（5）df51=-diff（f,y）/diff（f,x）df52=diff（diff（f,x）,y）format short eint（1/（1+x4+x8）,x）int（1/（asin（x）2*sqrt（1-x*x）,x）int（x*x+1）/（x4+1）,0,inf）a=int（exp（x）*（1+exp（x）2,0,log（2）eval（a）实验12 级数与方程符号求解syms n x;symsum（1

17、/（2*n-1）,n,1,10）s=symsum（n2*x（n-1）,n,1,inf）x=sym（5）;x=1/x;s=symsum（n2*x（n-1）,n,1,inf）/5x=sym（xf=log（x）;taylor（f,6,1）syms x y;y=solve（log（1+x）-5/（1+sin（x）=2）y=solve（x2+9*sqrt（x+1）-1）y=solve（3*x*exp（x）+5*sin（x）-78.5）x y=solve（sqrt（x2+y2）-100,3*x+5*y-8）y=dsolve（D2y+4*Dy+29*yy（0）=0,Dy（0）=15x,y,z=dsolve（Dx-2*x+3*y-3*zDy-4*x+5*y-3*zDz-4*x+4*y-2*zt实验13 低层绘图操作h=figure（MenuBarfigurecolor