Matlab程序设计与应用第二版刘卫国课后实验答案Word下载.doc

1、disp（A+6*B）;disp（A-B+I=disp（A-B+eye（3）;A*B=disp（A*B）;A.*B=disp（A.*B）;A3=disp（A3）;A.3=disp（A.3）;A/B=disp（A/B）;BA=disp（BA）;A,B=disp（A,B）;A（1,3,:B2=disp（A（1,3,:B2）;T3:z=1:25;A=reshape（z,5,5）;B=3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11;C=A*BD=C（3:5,2:3）T4-1:a=100:999;b=find（rem（a,21）=0）;c=length（b）T4-2:a=inp

2、ut（请输入一个字符串：,sb=find（a=A&ab disp（帕萨卡矩阵P性能更好elseif a希尔伯特矩阵H性能更好else两个矩阵性能相同enda=1:A=reshape（a,5,5）行列式的值：disp（det（A）;矩阵的秩：disp（rank（A）;矩阵的迹：disp（trace（A）;矩阵的范数：disp（norm（A）;T4:A=-29,6,18;20,5,12;-8,8,5V,D=eig（A）T5:A=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6B=0.95,0.67,0.52X1=ABB（3）=0.53X2=AB系数矩阵A的条件数：disp（c

3、ond（A）;T6:B1=sqrtm（A）B1=sqrtm（A）B2=sqrt（A）B2=sqrt（A）B1*B1B1*B1B2.*B2B2.*B2实验三：T1:x=-5.0:2:5.0for i=1:length（x）; if （x（i）=0）&（x（i）5）&（x（i）=2）&（x（i）=3）; y（i）=x（i）2-5*x（i）+6; else y（i）=x（i）2-x（i）-1; endyT2-if:s=input（please enter the score:while （s100） disp（the score is not reasonable s=input（if （s=90&

4、s=80&90）;B=70&80）;C=60&70）;DET2-switch:switch fix（s/10）case9,10case8case7case6case0,1,2,3,4,5t=input（请输入工时time=if t120 w=120*84+（t-120）*1.15*84;elseif t60 w=t*84-700; w=84*t;应发工资为：disp（w）;a=10+floor（rand（1）*89）b=10+floor（rand（1）*89）请输入+或-或*或/while（s=+s=-*/）输入的符号错误，请重新输入switch s case c=a+b; c=a-b; c=

5、a*b; c=a/b; cA=rand（5,6）n=input（请输入n:while （n5 B=A（5,: B=A（n,:End实验四：m=0;n; m=m+1/i2;pii=sqrt（6*m）T1-2:m=1./（1:n）.2;s=sum（m）;pii=sqrt（6*s）T2:n=1000;y=0; y=y+1/（2*i-1）; if y=3 disp（最大n值为： i-1对应的y值为： y-1/（2*i-1） break;format long请输入a:b=input（请输入b:x（1）=a/（1.0+b）;i=1;x（2）=a/（x（i）+b）;while （abs（x（i+1）-x（

6、i）0.00001&i500） i=i+1; x（i+1）=a/（x（i）+b）;r1=（-b+sqrt（b2+4*a）/2;r2=（-b-sqrt（b2+4*a）/2;disp（x,num2str（i+1）,num2str（x（i+1）,8）;,num2str（i）,num2str（x（i）,8）;r1=,num2str（r1,8）,r2=,num2str（r2,8）;format shortf（1）=1;f（2）=0;f（3）=1;i=4;while（i0）fn=length（find（ff0=length（find（f=0）j=0;s=0;for i=3:50 m=i*（i-1）-1;

7、if isprime（m）=1 m s=s+m; j=j+1;亲密数的对数：,num2str（j）;所有亲密素数之和：,num2str（s）;实验五：fn1:function f=fn1（n）f=n+10*log（n.2+5）;fn2:function f=fn2（n）f=0;n f=f+i*（i+1）;fx:function f=fx（x）f=1./（x-2）.2+0.1）+1./（x-3）.4+0.01）;x=input（请输入一个复数：fe=exp（x）flg=log（x）fsin=sin（x）fcos=cos（x）m1=input（请输入m1:m2=input（请输入m2:请输入theta:A=m1*cos（t）,-m1,-sin（t）,0;m1*sin（t）,0,. cos（t）,0;0,m2,-sin（t）,0;0,0,-cos（t）,1B=0,m1*9.8,0,m2*9.8a1,a2,N1,N2=disp（ABj=1;n=1;k=1;for i=10:99 % 挑选出1099中所有的素数 if isprime（i）=1 m（j）=i; j=j+1; endfor t=1:length（m）; s（n）=10*rem（m（t）,10）+