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 等效于 A 的逆左乘 B 矩阵，即 inv(A)*B ，而 B/

2、A 等效于 A 矩阵的逆右乘 B 矩阵，即 B*inv(A) 。3. 写出完成下列操作的命令。(1) 将矩阵 A 第25行中第 1, 3, 5 列元素赋给矩阵 B 答： B=A(2:5,1:2: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 转换成 3 4 矩阵。答： reshape(x,3,4);(7) 求一个字符串的

3、ASCII 码。答： abs(123 ); 或 double( 123);(8) 求一个 ASCII 码所对应的字符。答： 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, 0)L2 的值为 (1, 1, 1, 1, 1, 0, 0, 0, 0)L3 的值为 (0, 0, 0, 1, 1, 1, 0, 0, 0)L4 的值为 (4, 5, 6)5. 已知23100.77804145655A32503269.54543.14完成

4、下列操作：(1) 取出A的前3行构成矩阵B，前两列构成矩阵C,右下角 32子矩阵构成矩阵 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) 。01111100答： E=10&A=0.5 & x=1.5 & x=3.5 & x5.5 y=log(abs(b+c/x);enddisp(y=, num2str(y);(2)用switch语句实现计算: a=input(请输入a的值：);b=input(请输入b的值：);c=input(请输入c的值：); x二input(请输入x的值：);switc

5、h fix(x/0.5)case 1,2y=a*xA2+b*x+c;case n um2cell(3:6)y=a*(sin(b)Ac)+x;case num2cell(7:10)y=log(abs(b+c/x);enddisp(y二,num2str(y);3. 产生20个两位随机整数，输岀其中小于平均值的偶数。 答：A=fix(10+89*ra nd(1,20);sum=0;for i=1:20sum=sum+A(i);endB=A(fi nd(A v_maxv_max=x;en d;if x v_minv_mi n=x;en d;enddisp(最大数为：,n um2str(v_max);d

6、isp(最小数为：，nu m2str(v_mi n);(2) 用max函数、min函数实现：for i=1:5A(i)二input(请输入第,num2str(i),数：);enddisp(最大数为：,nu m2str(max(A);disp(最小数为：，nu m2str(mi n(A);5.2已知：s 12223 III632 ，分别用循环结构和调用MATLAB 的sum函数求s的值。答:(1)用循环结构实现:s=0;for i=0:63 s=s+2Ai;ends(2)调用sum函数实现:s=0:63;s=2.As;sum(s)求下列各式的值。6. 当 n 分别取 100、1000、10000

7、 时,（1）III(1)n(ln2)(2)l|(3)1664III4nK 1)(2n) (2n)答:要求分别用循环结构和向量运算（使用sum或prod函数）来实现。(1)用循环结构实现：sum=0;for k=1:100sum=sum+(-1)A(k+1)/k;endsum使用sum函数：x=；for k=1:10000x=x, (-1)A(k+1)/k;endsum(x)(2)用循环结构实现：sum=0;for k=1:100sum=sum+(-1)A(k+1)/(2*k-1);endsum使用sum函数：for k=1:100x=x, (-1)八(k+1)/(2*k-1);endsum(x

8、)(3) 用循环结构实现：sum=0;for k=1:100sum=sum+1/(4Ak);endsum使用 sum 函数实现：x=;for k=1:100x=x, 1/(4Ak);endsum(x)(4) 用循环结构实现：t=1;for k=1:100t=t*(2*k)*(2*k)/(2*k-1)*(2*k+1);end使用 prod 函数实现：x=;for k=1:100x=x, (2*k)*(2*k)/(2*k-1)*(2*k+1);endprod(x)7. 编写一个函数文件，求小于任意自然数n 的斐波那契(Fibnacci )数列各项。斐波那契数列定义如下：答：function x=f

9、ibnacci(n)for i=1:nif i=2x(i)=1;elseend1,1,n1n 2, nx(i)=x(i-1)+x(i-2);end8.编写一个函数文件，用于求两个矩阵的乘积和点乘，然后在命令文件中调用该函数。答：函数文件 myfnc.m :function x, y= myfnc（A, B）tryx=A*B;catchx=；endy=A.*B;命令文件 myexe.m :A=input（请输入矩阵A:）;B=input（请输入矩阵 B :）;x, y=myfnc（A, B）;if le ngth（x）=Odisplay两矩阵的维数不匹配，无法进行乘积运算！）;elsedisp（矩阵A和矩阵B的乘积为：）;xenddisp（矩阵A和矩阵B的点乘为：）;y10050109.先用函数的递归调用定义一个函数文件求然后调用该函数文件求k k2丄k 1 k 1 k 1 k答：函数文件 myfnc.m :fun cti on sum=myf nc(n, m)if n=1sum=1;elsesum= myf nc (n-1, m)+ nm;end10100502 在命令窗口中调用myfnc.m文件，计算kkk 1 k 1sum=myf nc(100, 1)+ myfnc(50, 2)+myfnc(10,-1)