MATLAB程序设计与应用第二课后题答案

1、作者：Pan Hon glia ng仅供个人学习第二章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 64 21;(1) B=A(2:5,1:2:5)B =652421349821842121434521(2) A(7)=A =246534 84323 54221924984245 211175 53646 21 2112121(3)A+30(4)size(A); ndiims(A)(5)题目有误(6)reshape(x,：3,4)(7)abs(x)(8)char(x)4. L1=0 0001

2、00 00L2 =1 1111 00 00L3 =0 0011 10 00L4 =4565.(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.1

3、330-0.2103E<Dans =010001E&Dans =1 10111E|D ans =111111DE ans =001000fin d(A>=10&A<25) ans =15&all(A)ans =0any(A)ans =1 isnan(A) ans =0 1 00000isi nf(A)ans =0 011000isfi nite(A)ans =1 0001117.A(1).x1='学号'A(1).x2='姓名'A(1).x3='专业'A(1).x4.x41='成绩 1'A(

4、2).x1='学号'A(2).x2='姓名'A(2).x3='专业'A(2).x4.x41='成绩 1 'A(3).x1='学号'A(3).x2='姓名;A(3).x3='专业'A(3).x4.x41='成绩 1'A(4).x 1='学号'A(4).x2='姓名'A.x3='专业'A(4).x4.x41='成绩 1 'A(5).x1='学号'A(5).x2='姓名'A(5).x3=

5、'专业'A(5).x4.x41='成绩 1 ' &(1)size(B)ans =2 2n dims(B)ans =2B(2)ans =3x3 doubleBans =3x3 cellB(3)=B =1 3x3 double 3x3 cellB3=B =13x3 double AVV*第二早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)2.B=rot90(A)

6、C=rot90(A,-1)3.B=inv(A);A地逆矩阵C=det(A) ;A地行列式地值D=A*BE=B*AD=E 因此A与A-1是互逆地4.A=4 2 -1;3 -1 2;12 3 0;b=2;10;8;x=i nv(A)*bx =-6.000026.666727.33335.diag(A);主对角线元素ans =1159triu(A);上三角阵ans =1-12301-4200520009tril(A);下三角阵ans =10 0051003 050111509ran k(A)；秩ans =4norm(A);范数ans =21.3005con d(A);条件数ans =11.1739t

7、race(A)；迹ans =16(2) 略&A=1 1 0.5;1 1 0.25;0.5 0.25 2 A =1.00001.00000.50001.00001.00000.25000.50000.25002.0000V,D=eig(A)V =0.72120.44430.5315-0.68630.56210.4615-0.0937-0.69760.7103D =-0.01660001.48010002.5365第四章1.a=i nput('请输入一个4位数:');while (a<1000|a>9999)a=input('输入错误，请重新输入一个4位

8、数:'); 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,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=i nput('请输入 a:');b=i nput('请输入 b:');c=i nput('请

9、输入 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.*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('平均数是：,n um2str(x1)m=(rem(x,2)=0 &x<x

10、1);n=fin d(m);disp('小于平均数地数是：,n um2str(x( n);4. A=i nput('请输入20个数地一个行向量：'); m=A;a=m;b=m;for m=Aif a>=ma=m;elseif b<=mb=m;endenddisp('最小数是：',n um2str(a)disp('最大数是：',n um2str(b)5. s=0;a=0;for b=1:64c=2Aa;a=a+1;s=s+c;enddisp('2 地0次方到63次方地和是：',n um2str(s)6. sum

11、1=0;for n=1:100x=(-1)A( n+1)*(1/n);sum1=sum1+x;enddisp('当 n取 100 时：sum=', nu m2str(sum1) sum2=0;for n=1:1000x=(-1)A( n+1)*(1/n);sum2=sum2+x;enddisp('当 n取 1000时：sum=' ,n um2str(sum2) sum3=0;for n=1:10000x=(-1)A( n+1)*(1/n);sum3=sum3+x;enddisp('当 n取 10000时：sum=',num2str(sum3)8.

12、 clear alla=i nput('请输入一个矩阵：');b=input('请再输入一个矩阵：(注意：两矩阵要可以相乘)');f1,f2=juzhenji(a,b);9. clear ally=qiuhe(100,1)+qiuhe(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)(3) ezplot('xA2

13、+yA2=1')(4)t=-10:0.1:10;x=t.A2;y=5*t.A3;plot(x,y)2.(1)theta=0:0.01:2*pi;rho=5*cos(theta)+4;polar(theta,rho)theta=0.001:0.1:2*pi;rho=12./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.*theta.A2;polar(theta,rho)3.(1)t=0:pi/100:2*

14、pi;x=cos(t);y=si n(t);z=t;plot3(x,y,z)u=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(x,y,z)5.plot函数：>> x=li nspace(-10,10,200);>> y=;>> for x0=xif x0>0y=y,x0.A2+(1+x0).A(1/4)+5;elseif x0=0y=y,0;elseif x0<0y=y,x0.A3+sqrt(1-x0)-5;enden

15、d>> plot(x,y)fplot函数：fplot('(xv0).*(x.A3+sqrt(1-x)-5)+(x=0).*0+(x>0).*(x.A2+(1+x)A(l/4)+5)',-10,10)第八早1.A=ra ndn (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'escend);A地每行元素按将序排列2.(1)(2

16、)X=1 4 9 16 25 36 49 64 81 100;Y=1:10;X1=1:100;Y1=i nterp1(X,Y ,X1,'cubic')3.x=165 123 150 123 141;y=187 126 172 125 148;P=polyfit(x,y,3)P =1.0e+003 *-0.00000.0013-0.17798.4330所以它地线性拟合曲线为：p(x)=1.3x2 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.50

17、000.5000-2.0000所以 P(x)=15x5+7x4-3.5x3+0.5x2-2x-2roots(P)ans =0.70710.1000 + 0.6245i0.1000 - 0.6245i-0.7071-0.6667i=0:10;xi=0.2*i;polyval(P,xi)ans =-2.0000-2.3920-2.6112-1.70242.710415.000094.1408 184.9056 332.5264 560.00005.(1)建立函数文件：fun cti on f=fxy(u)x=u(1);y=u(2);f=3.*x.A2+2*x.*y+y.A2在命令窗口中输入以下命令

18、：U,fmi n=fmi nsearch('fxy',1,1)结果：U =1.0e-004 *-0.06750.1715fmin =1.9920e-010(2)f=i nlin e('-si n(x)-cos(x.A2)');fmax=fm inbn d(f,0,pi)fmax =-2.000042.11200.7310(1) x=pi/6 pi/4 pi/3;f=in li ne('si n(x).A2+cos(x)42');dx=diff(f(x,5*pi/12)/(pi/12)可参见第 157 页例题 6.19dx =0 0 0x=pi/2

19、时单独计算：x=pi/2;f=in li ne('si n(x).A2+cos(x)42');diff(f(x,pi)/(pi/2)ans =0(2)x=1:3;f=inlin e('sqrt(x.A2+1)');dx=diff(f(x,4)结果：dx =0.82190.92620.96087. ( 1)f=in li ne('si n(x).A5.*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.2

20、214(3)f=in li ne('x.*si n(x)./(1+cos(x).A2)');quad(f,0,pi)ans =2.4674(4)f=i nli ne('abs(cos(x+y)');dblquad(f,0,pi,0,pi)ans =6.2832&N=64;%采样点数T=5;%采样时间终点t=linspace(0,T,N); %给出 N 个采样时间 ti (i=1:N) y=exp(-t);%求各采样点样本值ydt=t (2)-t(1); f=1/dt; Y=fft(y); F=Y(1:N/2+1); f=f*(0:N/2)/N; plot

21、(f,abs(F) 9.(1) 矩阵求逆法：%采样周期%米样频率%计算y地快速傅里叶变换丫%F(k)=Y(k)%使频率轴f从0开始%绘制振幅-频率图A=2 3 5;3 7 4;1 -7 1; b=10；3;5;x=i nv(A)*b x =-1.8060-0.53733.0448矩阵除法法：A=2 3 5;3 7 4;1 -7 1; b=10;3;5;x=Ab x =-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.函数文件：li

22、n 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)=n disp('原方程组有唯一解x');x=Ab; elsedisp（'原方程组有无穷个解,特解为x,齐次方程组地基础解系为y'）;x=Ab;y=n ull(A,'r');enddisp('方程组无解')；x=;endelsedisp('原方程组有零解x');x=zeros(

23、 n,1);if ran k(A) <ndisp('方程组有无穷个解,基础解系为y');y=n ull(A,'r');endend程序：A=2 1-1 1;4 2 -2 1;2 1-1-1;b=1;2;1;x,y=li ne_solutio n( A,b)结果：一原方程组有无穷个解，特解为x,齐次方程组地基础解系为yWarning: Rank deficient, rank = 2, tol =4.3512e-015.> In lin e_soluti on at 11方程组无解x =y =-0.50000.50001.0000 001.00000

24、011.(1)f=in li ne('x-s in( x)./x');x=fzero(f,0.5)x =0.8767(2)f=in li ne('(si n(x)42).*exp(-0.1.*x)-0.5.*abs(x)'); x=fzero(f,1.5)x =1.673812.函数文件：fun cti on f=fxy(u)x=u(1)y=u(2)f(1)=x-0.6*si n(x)-0.3*cos(y) f(2)=y-0.6*cos(x)+0.3*si n(y)在命令窗口输入以下命令：x=fsolve('fxy',0.5,0.5,optims

25、et('Display','off) 结果：x =0.63540.373415.A=-1 2 0;-1 2 -1;-1 2 -1;-1 2 -1;0 2 -1; d=-1;0;1;B=spdiags(A,d,5,5);b=1 0 0 0 0'x=(i nv(B)*b)'x =0.83330.66670.50000.33330.1667版权申明本文部分内容，包括文字、图片、以及设计等在网上搜集整理。版权为潘宏亮个人所有This article in eludes someparts, in cludi ng text, pictures, and desig n. Copyright is Pan Hon glia ng's pers onal own ership.用户可将本文的内容或服务用于个人学习、研究或欣赏，以及其他非商业性或非盈利性用途，但同时应遵守著作权法及其他相关法律 的规定，不得侵犯本网站及相关权利人的合法权利。除此以外，将本 文任何内容或服务用于其他用途时，须征得本人及相关权利人的书面 许可，并支付报酬。Users may use the contents or services of this articlefor