MATLAB程序设计与地的应用第二版课后地训练题目答案详解

1、第二章3.设矩阵A为A=24 23 9 21 6;65 74 24 11 21;34 5 98 75 21;8 42 42 53 121;43 21 45 6421；(1)B=A(2:5,1:2:5)B =652421349821842121434521(2)A(7)=A =246534843235422192498424521117553 646212112121(3)A+30(4)size(A) ;n dims(A)(5)题目有误(6)reshape(x,3,4)(7)abs(x)(8)char(x)000010000L2 =111110000L3 =000111000L4 =4565.(

2、1)B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*CB =23.000010.0000-0.7780041.0000-45.000065.00005.000032.00005.0000032.0000C =23.000010.000041.0000-45.000032.00005.00006.0000-9.5400D =65.00005.0000032.000054.00003.1400E =1.0e+003 *0.9141-0.22391.20802.71231.1330-0.2103E<Dans =010001E&Dans =110111E|Dans

3、 =111 1DEans =001000fin d(A>=10&A<25)ans =156.all(A)ans =0any(A)ans =1isnan(A)ans =isi nf(A)ans =0011000isfin ite(A)ans =10001117.'成绩1'；'成绩1'；'成绩1 '；'成绩1''成绩1'A(1).x1='学号'A(1).x2='姓名'A(1).x3='专业'；A(1).x4.x4仁A(2).x1='学号'

4、;；A(2).x2='姓名'A(2).x3='专业'；A(2).x4.x4仁A(3).x1='学号'；A(3).x2='姓名'；A(3).x3='专业'；A(3).x4.x4仁A(4) .x1='学号'；A .x2='姓名'；A(4).x3='专业'；A (4) .x4.x 4仁 A(5).x1='学号'；A(5).x2='姓名'；A(5).x3='专业'；A(5).x4.x4仁8.(1)size(B)ans =2 2n

5、 dims(B)ans =2ans =3x3 doubleB ans =3x3 cellB(3)=B =13x3 double3x3 cellB3=B =13x3 double第三章1.(1)A=eye(3) A=100+100*ra nd(5,6)(3) A=1+sqrt(0.2)*ra ndn( 10,50)(4) B=o nes(size(A)(5) A+30*eye(size(A)(6) B=diag(diag(A)B=rot90(A)C=rot90(A,-1)3.B=inv(A);A的逆矩阵C=det(A);A的行列式的值D=A*BE=B*AD=E 因此A与A-1是互逆的c4.A=4

6、 2 -1;3 -1 2;12 3 0;b=2;10;8;x=i nv(A)*bx =-6.000026.666727.33335.(1)diag(A);主对角线元素ans =11triu(A)上三角阵ans =-1-4tril(A)下三角阵ans =1115ran k(A) ans = no rm(A) ans =21.3005 con d(A) ans =trace(A)；迹ans =16(2) 略6.A=1 1 0.5;1 1 0.25;0.5 0.25 2A =1.00001.00000.50001.00001.00000.25000.50000.25002.0000V,D=eig(A

7、)V =0.72120.44430.53150.68630.56210.46150.0937-0.69760.7103D =-0.0166 0 001.48010002.5365第四章1. a=input('请输入一个4位数:');while (a<1000|a>9999)a=input('输入错误，请重新输入一个4位数:');endb=fix(a/1000);c=rem(fix(a/100),10);d=rem(fix(a/10),10);e=rem(a,10);b=b+7;c=c+7;d=d+7;e=e+7;b=rem(b,10);c=rem(c

8、,10);d=rem(c,10);e=rem(e,10);g=b;b=d;d=g;g=c;c=e;e=g;a=1000*d+100*e+10*b+c;disp('加密后：',num2str(a)2. a=input('请输入 a:');b=input('请输入 b:');c=input('请输入 c:');x=0.5:1:5.5;x1=(x>=0.5 &x<1.5);x2=(x>=1.5 &x<3.5);x3=(x>=3.5 &x<=5.5);y1=a.*(x.A2)+b

9、.*x+c;y2=a*(si n(bFc)+x;y3=log(abs(b+c./x);y=y1.*x1+y1.*x2+y3.*x3;disp(y)3. x=fix(ra nd(1,20)*89)+10;x1=fix(sum(x)/20);disp('平均数是：',num2str(x1)m=(rem(x,2)=0 &x<x1);n=fin d(m);disp('小于平均数的数是：',num2str(x(n);4. A=i nput('请输入20个数的一个行向量：');m=A;a=m;b=m;for m=Aif a>=ma=m;e

10、lseif b<=mb=m;endenddisp('最小数是：',num2str(a)disp('最大数是：',num2str(b)5. s=0;a=0;for b=1:64c=2Aa;a=a+1;s=s+c;enddisp('2的0次方到63次方的和是：',num2str(s)6, sum 1=0;for n=1:100x=(-1)A( n+1)*(1/n);sum仁sum1+x;enddisp('当 n取 100 时：sum=',num2str(sum1) sum2=0;for n=1:1000x=(-1)A( n+1)

11、*(1/n);sum2=sum2+x;enddisp('当 n取 1000 时：sum=',num2str(sum2)sum3=0;for n=1:10000x=(-1)A( n+1)*(1/n);sum3=sum3+x;enddisp('当 n取 10000 时：sum=',num2str(sum3)8. clear alla=input('请输入一个矩阵：');b=input('请再输入一个矩阵：(注意：两矩阵要可以相乘)');f1,f2=juzhenji(a,b);9. clear ally=qiuhe(100,1)+qiu

12、he(50,2)+qiuhe(10,-1);disp(y)10. (1)S=108 x=412 20y=246第五早1.(1) x=-10:0.1:10;y=100./(1+x.A2);Plot(x,y) x=-10:0.1:10;y=1/(2*pi)*exp(-x.A2/2);Plot(x,y) ezplot('xA2+yA2=1')2.t=-10:0.1:10;x=t.A2;y=5*t.A3;Plot(x,y)theta=0:0.01:2*pi; rho=5*cos(theta)+4; polar(theta,rho)theta=0.001:0.1:2*pi;rho=12.

13、/sqrt(theta);polar(theta,rho) theta=0.001:0.1:2*pi; rho=5./cos(theta)-7; polar(theta,rho)theta=0.001:0.1:2*pi;rho=pi/3.*theta42;polar(theta,rho)3.(1)t=0:pi/100:2*pi;x=cos(t);y=si n(t);z=t;plot3(x,y,z)4J 一.-1u=0:pi/100:2*pi;v=0:pi/100:2*pi; x=(1+cos(u).*cos(v); y=(1+cos(u).*si n(v); z=s in( u);plot3(

14、x,y,z)5.plot函数：>> x=li nspace(-10,10,200);>> y=;>> for x0=xif x0>0y=y,x042+(1+x0).A(1/4)+5;elseif x0=0y=y,0;elseif x0<0y=y,x0.A3+sqrt(1-x0)-5;endend>> plot(x,y)fplot函数：fplot('(xv0).*(x.A3+sqrt(1-x)-5)+(x=0).*0+(x>0).*(x.A2+(1+x)A(1/4)+5)',-10,10)第八早1.A=ra ndn

15、 (10,5)（1）mean(A)均值std(A);标准方差（2）max(max(A)；最大元素min( mi n( A);最小元素（3）B=sum(A,2);A每行元素的和sum(B);A全部元素之和（4）sort（A）;A的每列元素按升序排列sort(A,2, ' descend ' ) ;A的每行元素按将序排列2.(1)(2)X=1 4 9 16 25 36 49 64 81 100;X1=1:100;Y1=i nterp1(X,Y,X1,'cubic')3.x=165 123 150 123 141;y=187 126 172 125 148;P=pol

16、yfit(x,y,3)P =1.0e+003 *-0.00000.0013-0.17798.4330所以它的线性拟合曲线为：p(x)=1.3x 2 177.9x+84334.(1)P1=0 3 2;P2=5 -1 2;P3=1 0 -0.5;P=co nv(co nv(P1,P2),P3)P =015.00007.0000-3.50000.5000-2.0000-2.0000所以 P(x)=15x 5+7x 4-3.5x 3+0.5x 2-2x-2roots(P)ans =0.1000 + 0.6245i0.1000 - 0.6245i-0.7071-0.6667i=0:10;xi=0.2*i

17、;polyval(P ,xi)ans =15.0000-2.0000-2.3920-2.6112-1.70242.710442.112094.1408184.9056332.5264560.00005.(1)建立函数文件：fun cti on f=fxy(u)x=u(1);y=u(2);f=3.*x.A2+2*x.*y+y.A2在命令窗口中输入以下命令：U,fmi n=fmi nsearch('fxy',1,1)结果：U =-0.06750.1715fmin =1.9920e-010(2)f=i nlin e('-si n(x)-cos(x42)');fmax=

18、fm inbn d(f,0,pi)fmax =0.7310 6.(1) x=pi/6 pi/4 pi/3;可参见第157页例题6.19f=in li ne('si n(x).A2+cos(x)42'); dx=diff(f(x,5*pi/12)/(pi/12) dx =0 0 0x=pi/2时单独计算：x=pi/2;f=in li ne('si n(x).A2+cos(x)42');diff(f(x,pi)/(pi/2)ans =0(2)f=inlin e('sqrt(x42+1)');dx=diff(f(x,4)结果：dx =0.82190.9

19、2620.96087. (1)f=in li ne('si n(x)45.*si n(5*x)'); quad(f,0,pi)ans =0.0982(2)f=i nlin e('(1+x.A2)./(1+x.A4)'); quad(f,-1,1) ans =2.2214(3)f=in li ne('x.*si n(x)./(1+cos(x).A2)'); quad(f,0,pi) ans =2.4674(4)f=in li ne('abs(cos(x+y)'); dblquad(f,O,pi,O,pi) ans =6.28328.

20、N=64;%采样点数T=5;%采样时间终点t=li nspace(O,T,N);%给出N个采样时间ti (i=1:N )y=exp(-t);%求各米样点样本值ydt=t (2)-t(1);%采样周期f=1/dt;%米样频率Y=fft(y);%计算y的快速傅里叶变换丫F=Y(1:N/2+1);%F(k)=Y(k)f=f*(0:N/2)/N;%使频率轴f从0开始plot(f,abs(F)%绘制振幅-频率图9.(1)矩阵求逆法：A=2 3 5;3 7 4;1 -7 1;b=10；3;5;x=i nv(A)*bx =-1.8060-0.53733.0448矩阵除法法：A=2 3 5;3 7 4;1 -

21、7 1;b=10;3;5;x=Abx =-1.8060-0.53733.0448矩阵分解法：A=2 3 5;3 7 4;1 -7 1;b=10;3;5;L,U=lu(A);x=U(Lb)x =-1.8060-0.53733.0448(2)方法同(1)10.函数文件：lin e_solutio n(A,b)fun cti on x,y=li ne_solutio n( A,b)m, n=size(A);y=;if no rm(b)>0if ran k(A)=ra nk(A,b)if ran k(A)=ndisp('原方程组有唯一解x');x=Ab;elsedisp('原方程组有无穷个解，特解为x，齐次方程组的基础解系为y');x=Ab;y=n ull(A,'r');enddisp('方程组无解');X=；endelsedisp('原方程组有零解x');x=zeros( n,1);if ran k(A) <