matlab课后习题答案

1、 习题二1. 如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。因此，矩阵是MATLAB最基本、最重要的数据对象。2. 设A和B是两个同维同大小的矩阵，问：(1) A*B和A.*B的值是否相等？答：不相等。(2) A./B和B.A的值是否相等？答：相等。(3) A/B和BA的值是否相等？答：不相等。(4) A/B和BA所代表的数学含义是什么？答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而BA等效于B矩阵的逆左乘A矩阵，即inv(B)*A。3. 写出完成下

2、列操作的命令。(1) 将矩阵A第25行中第1, 3, 5列元素赋给矩阵B。答：B=A(2:5,1:2:5); 或B=A(2:5,1 3 5)(2) 删除矩阵A的第7号元素。答：A(7)=(3) 将矩阵A的每个元素值加30。答：A=A+30;(4) 求矩阵A的大小和维数。答：size(A);ndims(A);(5) 将向量 t 的0元素用机器零来代替。答：t(find(t=0)=eps;(6) 将含有12个元素的向量 x 转换成矩阵。答：reshape(x,3,4);(7) 求一个字符串的ASCII码。答：abs(123); 或double(123);(8) 求一个ASCII码所对应的字符。答：

3、char(49);4. 下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9;B=10-A;.L1=A=B;L2=A3&A3&A7);答：L1的值为0, 0, 0, 0, 1, 0, 0, 0, 0L2的值为1, 1, 1, 1, 1, 0, 0, 0, 0L3的值为0, 0, 0, 1, 1, 1, 0, 0, 0L4的值为4, 5, 65. 已知完成下列操作：(1) 取出A的前3行构成矩阵B，前两列构成矩阵C，右下角子矩阵构成矩阵D，B与C的乘积构成矩阵E。答：B=A(1:3,:);C=A(:,1:2);D=A(2:4,3:4);E=B*C;(2) 分别求E=10&A25)。答

4、：E=10&A=0.5 & x=1.5 & x=3.5 & x5.5y=log(abs(b+c/x);enddisp(y=,num2str(y);(2) 用switch语句实现计算：a=input(请输入a的值：);b=input(请输入b的值：);c=input(请输入c的值：);x=input(请输入x的值：);switch fix(x/0.5) case 1,2 y=a*x2+b*x+c; case num2cell(3:6) y=a*(sin(b)c)+x; case num2cell(7:10) y=log(abs(b+c/x); enddisp(y=,num2str(y);3. 产

5、生20个两位随机整数，输出其中小于平均值的偶数。答：A=fix(10+89*rand(1,20);sum=0;for i=1:20sum=sum+A(i);endB=A(find(A v_maxv_max=x;end;if x v_minv_min=x;end;enddisp(最大数为：, num2str(v_max);disp(最小数为：, num2str(v_min);(2) 用max函数、min函数实现：for i=1:5A(i)=input(请输入第, num2str(i), 数：);enddisp(最大数为：, num2str(max(A);disp(最小数为：, num2str(m

6、in(A);5. 已知：，分别用循环结构和调用MATLAB的sum函数求s的值。答：(1) 用循环结构实现：s=0;for i=0:63 s=s+2i;ends(2) 调用sum函数实现：s=0:63;s=2.s;sum(s)6. 当n分别取100、1000、10000时，求下列各式的值。(1) (2) (3) (4) 要求分别用循环结构和向量运算（使用sum或prod函数）来实现。答：(1) 用循环结构实现：sum=0;for k=1:100sum=sum+(-1)(k+1)/k;endsum使用sum函数：x=;for k=1:10000x=x, (-1)(k+1)/k;endsum(x)

7、(2) 用循环结构实现：sum=0;for k=1:100sum=sum+(-1)(k+1)/(2*k-1);endsum使用sum函数：x=;for k=1:100x=x, (-1)(k+1)/(2*k-1);endsum(x)(3) 用循环结构实现：sum=0;for k=1:100sum=sum+1/(4k);endsum使用sum函数实现：x=;for k=1:100x=x, 1/(4k);endsum(x)(4) 用循环结构实现：t=1;for k=1:100t=t*(2*k)*(2*k)/(2*k-1)*(2*k+1);endt使用prod函数实现：x=;for k=1:100x=

8、x, (2*k)*(2*k)/(2*k-1)*(2*k+1);endprod(x)7. 编写一个函数文件，求小于任意自然数n的斐波那契（Fibnacci）数列各项。斐波那契数列定义如下：答：function x=fibnacci(n)for i=1:nif i=2x(i)=1;elsex(i)=x(i-1)+x(i-2);endend8. 编写一个函数文件，用于求两个矩阵的乘积和点乘，然后在命令文件中调用该函数。答：函数文件myfnc.m：function x, y= myfnc(A, B)tryx=A*B;catchx=;endy=A.*B;命令文件myexe.m：A=input(请输入矩阵

9、A：);B=input(请输入矩阵B：);x, y=myfnc(A, B);if length(x)=0display(两矩阵的维数不匹配，无法进行乘积运算！);elsedisp(矩阵A和矩阵B的乘积为：);xenddisp(矩阵A和矩阵B的点乘为：);y9. 先用函数的递归调用定义一个函数文件求，然后调用该函数文件求。答：函数文件myfnc.m：function sum=myfnc(n, m)if n=1sum=1;elsesum= myfnc (n-1, m)+nm;end在命令窗口中调用myfnc.m文件，计算：sum=myfnc(100, 1)+ myfnc(50, 2)+myfnc(

10、10,-1)10. 写出下列程序的输出结果。 s=0;a=12,13,14;15,16,17;18,19,20;21,22,23;for k=afor j=1:4if rem(k(j),2)=0s=s+k(j);endendend s答：执行结果为s=108 命令文件exe.m执行后的结果为：x = 41220 y=2 4 61. 若不给自己设限，则人生中就没有限制你发挥的藩篱。2. 若不是心宽似海，哪有人生风平浪静。在纷杂的尘世里，为自己留下一片纯静的心灵空间，不管是潮起潮落，也不管是阴晴圆缺，你都可以免去浮躁，义无反顾，勇往直前，轻松自如地走好人生路上的每一步3. 花一些时间，总会看清一些事。用一些事情，总会看清一些人。有时候觉得自己像个神经病。既纠结了自己，又打扰了别人。努力过