MATLAB课后习题.docx

1、MATLAB课后习题第1章 MATLAB操作基础1 MATLAB的主要功能包括数值计算和符号计算功能，绘图功能，编程语言以及应用工具箱。2 MATLAB的优点在于可以在窗口直接输入运算命令，并的到计算结果；有许多现成的函数可以直接调用；可以直接绘制图形对结果更直接的显示。3 MATLAB中help命令在窗口中输入后将会显示当前帮助系统中的所包含的所有项目，即搜索路径中所有的目录项目，所以可以搜索到。4 主要有控制系统工具箱，信号处理工具箱，神经网络工具箱，最优化工具箱，金融工具箱，统计学工具箱，功能性工具箱。5 inv 为求方阵的逆矩阵，plot为对所给函数作图，max为取最大值，round为

2、四舍五入到最临近的整数。第2章 MATLAB矩阵及其运算2-1 合法的变量有global My_exp AbCd AB_C_2-2 直接输入法；利用M文件建立；利用冒号表达式建立一个向量；建立大矩阵。2-3 A*B要求A阵的列数必须与B阵的行数相等，A.*B表示对应元素相乘，矩阵运算左除与右除值一般不同，点除也有类似的规律，但对于标量则都相同。2-4 （1） a=3.5;b=5;c=-9.8; x=(2*pi+(b+c)/(pi+a*b*c)-exp( 2)/tan(b+c)+ax =3.4054 （2） w=sqrt( 2)*(1+0.34245e-6) w =1.4142 （3） t=2

3、1-3i;5 -0.65; z=(exp(2*t)*log(t+sqrt(1+t2)/2 z = 1.0e+004 * 0.0057 - 0.0007i 0.0049 - 0.0027i 1.9884 - 0.3696i 1.7706 - 1.0539i2-5 A=-1,5,-4;0,7,8;3,61,7; B=8,3,-1;2,5,3;-3,2,0; A+6*Bans = 47 23 -10 12 37 26 -15 73 7 I=eye; A2-B+I ans = -18 -216 18 23 533 110 22 868 526 A*Bans = 14 14 16 -10 51 21 1

4、25 328 180 A.*Bans = -8 15 4 0 35 24-9 122 0 B*A ans = -11 0 -15 7 228 53 3 -1 28 A/Bans = 1.2234 -0.9255 2.9787 -0.9468 2.3511 -0.9574 4.6170 3.8723 13.8936 BAans = -0.5106 -8.6170 -1.1277 0.7340 17.5745 1.8085 -0.8830 -21.2128 0.4043 A,Bans = -1 5 -4 8 3 -1 0 7 8 2 5 3 3 61 7 -3 2 02-6 A=23,10,-0.

5、778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14; B=A(1:3,:)B = 23.0000 10.0000 -0.7780 0 41.0000 -45.0000 65.0000 5.0000 32.0000 5.0000 0 32.0000 C=A(:,1:2)C = 23.0000 10.0000 41.0000 -45.0000 32.0000 5.00006.0000 -9.5400 D=A(2:end,3:end)D = 65.0000 5.0000 0 32.0000 54.0000 3.1400 E=B*CE = 1.0e+003 * 0.

6、9141 -0.2239 1.2080 2.7123 1.1330 -0.2103E E&D ans = 1 1 0 1 1 1 E|D,E| ans = 1 1 1 1 1 1 D ans = 0 0 1 0 0 0 find(A=10&A A=34,NaN,Inf,-Inf,-pi,eps,0A = 34.0000 NaN Inf -Inf -3.1416 0.0000 0 all(A)ans = 0 any(A)ans = 1 isnan(A)ans = 0 1 0 0 0 0 0 isinf(A)ans = 0 0 1 1 0 0 0 isfinite(A)ans = 1 0 0 0

7、1 1 12-8 A=1 -1 2 3;5 1 -4 2;3 0 5 2;11 15 0 9; diag(A)ans = 1 1 5 9 triu(A)ans = 1 -1 2 3 0 1 -4 2 0 0 5 2 0 0 0 9 tril(A)ans = 1 0 0 0 5 1 0 0 3 0 5 0 11 15 0 9 inv(A)ans = -0.1758 0.1641 0.2016 -0.0227 -0.1055 -0.1016 -0.0391 0.0664 -0.0508 -0.0859 0.1516 0.0023 0.3906 -0.0313 -0.1813 0.0281 det(

8、A)ans = 1280 rank(A)ans = 4 trace(A)ans = 16 norm(A,2)ans = 21.3005 norm(A,1)ans = 20 norm(A,inf)ans = 35 cond(A,1)ans = 14.4531 cond(A,2)ans = 11.1739 cond(A,inf)ans = 22.0938 B=0.43 43 2;-8.9 4 21; diag(B)ans = 0.4300 4.0000 triu(B)ans = 0.4300 43.0000 2.0000 0 4.0000 21.0000 tril(B)ans = 0.4300 0

9、 0 -8.9000 4.0000 0 pinv(B)ans = 0.0022 -0.0175 0.0234 -0.0017 -0.0035 0.0405 rank(B)ans = 2 trace(B)ans = 4.4300 norm(B,1)ans = 47 norm(B,2)ans = 43.4271 norm(B,inf)ans = 45.43002-9 format rat h=hilb(5)h = 1 1/2 1/3 1/4 1/5 1/2 1/3 1/4 1/5 1/6 1/3 1/4 1/5 1/6 1/7 1/4 1/5 1/6 1/7 1/8 1/5 1/6 1/7 1/8

10、 1/9 Hh=det(h)Hh = 1/266716800000 Th=cond(h)Th = 476607 p=pascal(6)p = 1 1 1 1 1 1 1 2 3 4 5 6 1 3 6 10 15 21 1 4 10 20 35 56 1 5 15 35 70 126 1 6 21 56 126 252 Hp=det(p)Hp = 1 Tp=cond(p)Tp = 332360/3 从矩阵的分析可以看出帕斯卡矩阵的性能要比希尔伯特矩阵的好。2-10 A=-29 6 18;20 5 12;-8 8 5; V,D=eig(A)V = 698/979 395/1409 989/361

11、9 -859/1412 -2430/3089 349/400 538/1543 423/769 578/1427 D = -3595/142 0 0 0 -3755/357 0 0 0 4697/279 A是n阶方阵，若有数和非零向量X，使得AX=X则称数是A的特征值，非零向量X是A对应于特征值的特征向量。2-11 其中6+7i表示一个复数，而6+7*i则表示一个表达式，i是虚数单位，I为一个变量。2-12 （1）A(7)= (2) t(find(t=0)=eps (3)x=reshape(x,3,4) (4)abs(djdfsfksf”) (5)size(A) ones(x,y) （6）B=

12、diag(A)2-12 A=3+randn(500,1)2-13 结构矩阵的建立具体格式为 结构矩阵名.成员名=表达式； 单元矩阵的建立与结构矩阵的建立基本相同只是用大括号把矩阵元素括起来。2-14 稀疏矩阵可以有效的节省存储空间，稀疏矩阵与完全矩阵不同它有独特的调用方式。第三章 MATLAB程序设计3-1 n=input(请输入一个三位数：); a=fix(n/100); b=rem(fix(n/10),10); c=rem(n,10); m=c*100+b*10+a;disp(m)3-2用if语句实现 s=input(请输入成绩:); if s=90&c=8&s=70&s=6&s=0&smax max=A(i+1); end if A(i+1)min min=A(i+1); endenddisp(max)disp(min)3-4 顺序结构 a=input(please input a:)y=0.5*(exp(0.3*a)-exp(-0.3*a)*sin(a+0.3)+log(0.5*(0.3+a)disp(y) 循环结构for a= -3.0:0.1:3.0 y=0.5*(exp(0.3*a)-exp(-0.3*a)