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

1、封面作者：PanHongliang仅供个人学习第二章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 =65 242134 98218 42 12143 4521(2)A(7)=A =24 65 34 84323 5 42 21 9 24 98 42 45 21 11 75 53 64 6 21 21 121 213)A+304)size(A)。 ndims(A)5)题目有误6)reshape(x,3,4)7)abs(x)8)ch

2、ar(x)4. L10000 10000L2 =1111 10000L3 =0001 11000L4 =4565.(1)B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*CB =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.0000C =23.0000 10.000041.0000 -45.000032.0000 5.00006.0000 -9.5400D =65.0000 5.00000 32.000054.0000 3.1400E =1.0e+003 *0.914

3、1 -0.22391.2080 2.71231.1330 -0.2103(2) E=10&A25)ans =156.all(A)ans =0any(A)ans =1 isnan(A)ans =0 1 0 0 0 isinf(A)ans =0 0 1isfinite(A) ans =1 0 07.A(1).x1= 学 号 。1. .A(2).x1 = 学 号 。1. .1 0 0 00 1 1 1A(1).x2= 姓 名 A(2).x2= 姓 名 A(1).x3=专业A(2).x3=专业A(1).x4.x4 仁成 绩A(2).x4.x4 仁成绩ans =3x3 cell(3)B(3)=B

4、 =13x3 double 3x3 cellB3=B =13x3 double 1.(1) A=eye(3)(2) A=100+100*rand(5,6)(3) A=1+sqrt(0.2)*randn(10,50)(4) B=ones(size(A)(5) A+30*eye(size(A)(6) B=diag(diag(A)2.B=rot90(A)C=rot90(A,-1)3.B=inv(A) 。 A 地逆矩阵A(3).x1= 学 号 1. .A(4).x1 = 学 号 1. .A(5).x 仁学号1. .8.(1)size(B)ans =2 2 ndims(B)ans =2(2)B(2)an

5、s =3x3 doubleB(4)A(3).x2= 姓名A(4).x2= 姓名A(5).x2= 姓名A(3).x3= 专 业A(4).x3= 专 业A(5).x3= 专 业A(3).x4.x41=A(4).x4.x41=A(5).x4.x41=成绩成绩成绩ATy *第三章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=inv(A)*bx =-6.000026.666727.33335.(1)diag(A) ；主对角线元素ans =1159triu(A) ；上三角阵ans =

6、1-12301-4200520009tril(A) ；下三角阵 ans =10005100305011 15 0 9rank(A) ；秩ans =4norm(A) ；范数 ans =21.3005cond(A) ；条件数 ans =11.1739trace(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)V =0.72120.44430.5315-0.68630.56210.4615-0.0937-0.6976

7、0.7103D =-0.01660001.480100 0 2.5365第四章1.a=input( 请输入一个 4 位数 :) 。 while(a9999) a=input( 输入错误 , 请重新输入一个 4 位数 :) end b=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=10

8、00*d+100*e+10*b+c。 disp( 加密后：,num2str(a)2.a=i nput(请输入 a:)。b=i nput(请输入 b:)。c=in put(请输入 c:)。x=0.5:1:5.5 。x1=(x=0.5&x=1.5&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(rand(1,20)*89)+10 。 x1=fix(sum(x)/20) 。disp( 平均数是： ,nu

9、m2str(x1) m=(rem(x,2)=0&x=m a=m。elseif b x=linspace(-10,10,200)。 y= 。 for x0=xif x00y=y,xO42+(1+xO).A(1/4)+5。elseif x0=0y=y,O。elseif xO plot(x,y)fplot 函数：fplot(xO).*(x.A2+(1+x).A(1/4)+5),-1O,1O)第六章1.A=randn(1O,5)( 1 ) mean(A) 。均值 std(A) 。标准方差( 2) max(max(A)。最大元素 min(min(A)。最小元素(3) B=sum(A,2)。A 每

10、行元素地和 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 1OO。Y=1:1O。X1=1:1OO。Y1=interp1(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.0000 0.0013 -0.1779 8.4330 所以它地线性拟合曲线为：p(x)=1.3x2177.9x+84334.

11、P 仁0 3 2。P2=5 -1 2。P3=1 0 -0.5。P=conv(conv(P1,P2),P3)P =0 15.0000 7.0000 -3.5000 0.5000 -2.0000 -2.0000 所以P(x)=15x5+7x4-3.5x3+0.5x2-2x-2(2)roots(P)ans =0.70710.1000 + 0.6245i0.1000 - 0.6245i-0.7071-0.6667(3)i=0:10。xi=0.2*i 。polyval(P,xi)ans =-2.0000 -2.3920-2.6112 -1.70242.710415.0000 42.1120184.905

12、6 332.5264 560.00005.(1)建立函数文件：function f=fxy(u)x=u(1)。 y=u(2)。f=3.*x.A2+2*x.*y+y.A2在命令窗口中输入以下命令：U,fmin=fminsearch(fxy,1,1)结果：U =1.0e-004 *-0.0675 0.1715fmin =1.9920e-010(2)f=inline(-sin(x)-cos(x42)。94.1408fmax=fminbnd(f,0,pi)fmax =0.73106.(1)x=pi/6 pi/4 pi/3 。f=inline(sin(x).A2+cos(x)42)。dx=diff(f(

13、x,5*pi/12)/(pi/12) 可参见第 157 页例题 6.19 dx =0 0 0 x=pi/2 时单独计算：x=pi/2。f=inline(sin(x).A2+cos(x)42)。 diff(f(x,pi)/(pi/2) ans =0(2)x=1:3。f=inline(sqrt(x.A2+1)0dx=diff(f(x,4)结果：dx =0.8219 0.9262 0.96087.(1)f=inline(sin(x).A5.*sin(5*x) 。quad(f,0,pi)ans =0.0982(2)f=inline(1+x.A2)./(1+x.A4) 。 quad(f,-1,1) an

14、s =2.2214(3)f=inline(x.*sin(x)./(1+cos(x).A2) 。quad(f,0,pi)ans =2.4674(4)f=inline(abs(cos(x+y)。 dblquad(f,0,pi,0,pi)ans =6.28328.N=64。采样点数T=5。 %采样时间终点 t=linspace(0,T,N)。 %给出 N 个采样时间 ti(i=1:N)y=exp(-t) 。 %求各采样点样本值 y dt=t(2)-t(1) 。 %采样周期 f=1/dt 。 % 采样频率Y=fft(y) 。 %计算 y 地快速傅里叶变换 YF=Y(1:N/2+1) 。 %F(k)=Y

15、(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=inv(A)*bx =-1.8060-0.53733.0448 矩阵除法法： A=2 3 5。 3 7 4。 1 -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.函

16、数文件： line_solution(A,b)function x,y=line_solution(A,b)m,n=size(A) 。x =y= if norm(b)0if rank(A)=rank(A,b)if rank(A)=ndisp( 原方程组有唯一解 x) 。 x=Ab 。elsedisp( 原方程组有无穷个解 , 特解为 x, 齐次方程组地基础解系为 x=Ab 。y=null(A,r) 。enddisp( 方程组无解 ) 。x= 。endelsedisp( 原方程组有零解 x) 。x=zeros(n,1) 。if rank(A) In line_solution at 11方程组无

17、解y) 。y =-0.50000.50001.0000 00 1.00000 011.(1)f=inline(x-sin(x)./x)。x=fzero(f,0.5)x =0.8767(2)f=inline(sin(x).A2).*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

18、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。0o1oB=spdiags(A,d,5,5)。b=1 0 0 0 0ox=(i nv(B)*b)x =0.83330.66670.50000.33330.1667版权申明本文部分内容，包括文字、图片、以及设计等在网上搜集整理。版权为潘宏亮个人所有This article in eludes some parts, in cludi ng text, pictures, anddesig n. Copyright is Pan Hon

19、glia ngs pers onal own ership.用户可将本文的内容或服务用于个人学习、研究或欣赏，以及 其他非商业性或非盈利性用途，但同时应遵守著作权法及其他相关 法律的规定，不得侵犯本网站及相关权利人的合法权利。除此以 外，将本文任何内容或服务用于其他用途时，须征得本人及相关权 利人的书面许可，并支付报酬。Users may use the contents or services of this article for personal study, research or appreciati on, and other non-commercial ornon-profit purposes, but at the same time, they shall abide by theprovisions of copyright law and other releva nt laws, and shall notinf