1、MATLAB习题答案习题二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/A等效于A矩阵的逆右乘B矩阵,即B*inv(
2、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 转换成矩阵。答:reshape(x,3,4);(7) 求一个字符串的ASCII码。答:abs(123); 或double(123);(8) 求一个ASCII码所对应的字符。答:char
3、(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. 已知 完成下列操作:(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=1
4、0&A25)。答: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的值:); 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=
5、,num2str(y);3. 产生20个两位随机整数,输出其中小于平均值的偶数。答: A=fix(10+89*rand(1,20); sum=0; for i=1:20sum=sum+A(i);endB=A(find(A v_max v_max=x;end;if x v_min v_min=x;end; enddisp(最大数为:, num2str(v_max);disp(最小数为:, num2str(v_min);(2) 用max函数、min函数实现:for i=1:5 A(i)=input(请输入第, num2str(i), 数:);enddisp(最大数为:, num2str(max(A
6、);disp(最小数为:, num2str(min(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:100 sum=sum+(-1)(k+1)/k; endsum使用sum函数: x=; for k
7、=1:10000 x=x, (-1)(k+1)/k;end sum(x) (2) 用循环结构实现: sum=0; for k=1:100 sum=sum+(-1)(k+1)/(2*k-1); endsum 使用sum函数:x=;for k=1:100 x=x, (-1)(k+1)/(2*k-1);end sum(x)(3) 用循环结构实现:sum=0; for k=1:100 sum=sum+1/(4k); endsum使用sum函数实现: x=;for k=1:100 x=x, 1/(4k);endsum(x)(4) 用循环结构实现:t=1;for k=1:100 t=t*(2*k)*(2*
8、k)/(2*k-1)*(2*k+1);endt使用prod函数实现:x=;for k=1:100 x=x, (2*k)*(2*k)/(2*k-1)*(2*k+1);endprod(x)7. 编写一个函数文件,求小于任意自然数n的斐波那契(Fibnacci)数列各项。斐波那契数列定义如下: 答: function x=fibnacci(n)for i=1:n if i=2 x(i)=1; else x(i)=x(i-1)+x(i-2); endend8. 编写一个函数文件,用于求两个矩阵的乘积和点乘,然后在命令文件中调用该函数。答: 函数文件myfnc.m: function x, y= myf
9、nc(A, B) tryx=A*B;catch x=;end y=A.*B; 命令文件myexe.m: A=input(请输入矩阵A:);B=input(请输入矩阵B:);x, y=myfnc(A, B);if length(x)=0 display(两矩阵的维数不匹配,无法进行乘积运算!);else disp(矩阵A和矩阵B的乘积为:); xenddisp(矩阵A和矩阵B的点乘为:);y9. 先用函数的递归调用定义一个函数文件求,然后调用该函数文件求。答: 函数文件myfnc.m: function sum=myfnc(n, m) if n=1 sum=1; else sum= myfnc (n-1, m)+nm; end 在命令窗口中调用myfnc.m文件,计算: sum=myfnc(100, 1)+ myfnc(50, 2)+myfnc(10,-1)10. 写出下列程序的输出结果。1 s=0;a=12,13,14;15,16,17;18,19,20;21,22,23;for k=a for j=1:4 if rem(k(j),2)=0 s=s+k(j); endendend s答:执行结果为 s=1082 命令文件exe.m执行后的结果为:x = 4 12 20 y= 2 4 6
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1