Matlab课后习题解答
matlab课后习题及答案
第一章 5题已知a=4.96，b=8.11，计算)ln(b a eba +-的值。
解：clear clc a=4.96; b=8.11;exp(a-b)/log(a+b) ans =0.0167 6题已知三角形的三边a=9.6,b=13.7, c=19.4,求三角形的面积。
提示：利用海伦公式area =))()((c s b s a s s ---计算，其中S=(A+B+C)/2. 解：clear clc a=9.6; b=13.7; c=19.4; s=(a+b+c)/2area=sqrt(s*(s-a)*(s-b)*(s-c)) s =21.3500 第二章 8题已知S=1+2+2^2+2^3+……+2^63,求S 的值 解：clear clc S=0;for i=0:1:63 S=S+2^i; end S S =1.8447e+019 9题分别用for 和while 循环结构编写程序，计算∑=-1001n 1n 2）（的值。
解：clear clc s=0;for n=1:100 s=s+(2*n-1); end s s =10000 clear clc n=1; s=0;while n<=100 s=s+(2*n-1); n=n+1; end s s =10000 第三章 2题在同一坐标下绘制函数x ,,2x-,2x xsin(x)在()∏∈,0x 的曲线。
解：clear clcx=0:0.2:pi; y1=x; y2=x.^2; y3=-(x.^2); y4=x.*sin(x);plot(x,y1,'-' ,x,y2,'-' ,x,y3,'-' ,x,y4,'-')0.511.522.53-10-8-6-4-202468109题用不同的线型和颜色在同一坐标内绘制曲线y1=2ex5.0 、y2=sin(2∏x )的图形。
matlab课后习题答案（附图）
matlab课后习题答案（附图）习题2.1画出下列常见曲线的图形y （1）⽴⽅抛物线3x命令：syms x y;ezplot('x.^(1/3)')(2)⾼斯曲线y=e^(-X^2);命令：clearsyms x y;ezplot('exp(-x*x)')(3)笛卡尔曲线命令：>> clear>> syms x y;>> a=1;>> ezplot(x^3+y^3-3*a*x*y)(4)蔓叶线命令：>> clear>> syms x y;>> a=1ezplot(y^2-(x^3)/(a-x))(5)摆线:()()tsin-=,=-by1命令：>> clear>> t=0:0.1:2*pi;>> x=t-sin(t);>>y=2*(1-cos(t)); >> plot(x,y)7螺旋线命令：>> clear >> t=0:0.1:2*pi; >> x=cos(t); >> y=sin(t); >> z=t;>>plot3(x,y,z)(8)阿基⽶德螺线命令：clear>> theta=0:0.1:2*pi;>> rho1=(theta);>> subplot(1,2,1),polar(theta,rho1)(9) 对数螺线命令：cleartheta=0:0.1:2*pi;rho1=exp(theta);subplot(1,2,1),polar(theta,rho1)(12)⼼形线命令：>> clear >> theta=0:0.1:2*pi; >> rho1=1+cos(theta); >> subplot(1,2,1),polar(theta,rho1)练习2.21. 求出下列极限值（1）nnn n3→命令：>>syms n>>limit((n^3+3^n)^(1/n)) ans = 3（2）)121(lim n n n n ++-+∞→命令：>>syms n>>limit((n+2)^(1/2)-2*(n+1)^(1/2)+n^(1/2),n,inf) ans = 0（3）x x x 2cot lim 0→命令：syms x ；>> limit(x*cot(2*x),x,0) ans = 1/2 （4）)(coslimcm xx ∞→命令：syms x m ； limit((cos(m/x))^x,x,inf) ans = 1（5）)111(lim 1--→exx x命令：syms x>> limit(1/x-1/(exp(x)-1),x,1) ans =(exp(1)-2)/(exp(1)-1) （6）)(2lim x x xx -+∞>> limit((x^2+x)^(1/2)-x,x,inf)ans = 1/2练习2.41. 求下列不定积分，并⽤diff 验证：（1）+x dxcos 1>>Clear >> syms x y >> y=1/(1+cos(x)); >> f=int(y,x) f =tan(1/2*x) >> y=tan(1/2*x); >> yx=diff(y ,x); >> y1=simple(yx) y1 =1/2+1/2*tan(1/2*x)^2 (2)+exdx1clear syms x yy=1/(1+exp(x));f=int(y,x) f =-log(1+exp(x))+log(exp(x)) syms x yy=-log(1+exp(x))+log(exp(x)); yx=diff(y,x); y1=simple(yx) y1 = 1/(1+exp(x)) (3)dx x x ?sin 2syms x yy=x*sin(x)^2; >> f=int(y,x) f =x*(-1/2*cos(x)*sin(x)+1/2*x)-1/4*cos(x)^2-1/4*x^2 clearsyms x y y=x*(-1/2*cos(x)*sin(x)+1/2*x)-1/4*cos(x)^2-1/4*x^2; yx=diff(y,x); >> y1=simple(yx) y1 = x*sin(x)^2 (4)xdx ?sec3syms x y y=sec(x)^3;f=int(y,x) f =1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x)) clear syms x yy=1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x)); yx=diff(y,x); y1=simple(yx) y1 =1/cos(x)^32. 求下列积分的数值解 1）dx x-10clearsyms xy=int(x^(-x),x,0,1) y =int(x^(-x),x = 0 .. 1) vpa(y,10) ans =1.291285997 2）xdx e x cos3202?πclearsyms xy=int(exp(2*x)*cos(x)^3,x, clear syms xy=int((1/(2*pi)^(1/2))*exp(-x^2/2),x,0,1) y =7186705221432913/36028797018963968*erf(1/2*2^(1/2))*2^(1/2)*pi^(1/0,2*pi) y =22/65*exp(pi)^4-22/65vpa(ans,10)(3)dx xe21221-π>> clear >> syms x>> y=int(1/(2*pi)^(1/2)*exp(-x^2/2),0,1); >> vpa(y,14) ans =.341344746068552（4）>> clear >> syms x>> y=int(x*log(x^4)*asin(1/x^2),1,3); Warning: Explicit integral could not be found. > In sym.int at 58 >> vpa(y,14) ans = 2.45977212823752（5） >> clear >> syms x1判断下列级数的收敛性，若收敛，求出其收敛值。
matlab课后习题答案
习题二之南宫帮珍创作2.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵, 单个数据（标量）可以看成是仅含一个元素的矩阵, 故向量和单个数据都可以作为矩阵的特例来处置.因此, 矩阵是MATLAB最基本、最重要的数据对象.3.设A和B是两个同维同年夜小的矩阵, 问：(1)A*B和A.*B的值是否相等？答：不相等.(2)A./B和B.\A的值是否相等？答：相等.(3)A/B和B\A的值是否相等？答：不相等.(4)A/B和B\A所代表的数学含义是什么？答：A/B等效于B的逆右乘A矩阵, 即A*inv(B), 而B\A等效于B矩阵的逆左乘A矩阵, 即inv(B)*A.4.写出完成下列把持的命令.(1)将矩阵A第2~5行中第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转换成34 矩阵.答：reshape(x,3,4);(7)求一个字符串的ASCII码.答：abs(‘123’); 或double(‘123’);(8)求一个ASCII码所对应的字符.答：char(49);5.下列命令执行后, L1、L2、L3、L4的值分别是几多？A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：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]6.已知完成下列把持：(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<D、E&D、E|D、~E|~D和find(A>=10&A<25).答：E<D=010001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦, E&D=110111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦, E|D=111111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦, ~E|~D=001000⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦find(A>=10&A<25)=[1; 5].7.当A=[34, NaN, Inf, -Inf, -pi, eps, 0]时, 分析下列函数的执行结果：all(A), any(A), isnan(A), isinf(A), isfinite(A).答：all(A)的值为0any(A)的值为1isnan(A)的值为[0, 1, 0, 0, 0, 0, 0]isinf(A)的值为[0, 0, 1, 1, 0, 0, 0]isfinite(A)的值为[1, 0, 0,0, 1, 1, 1]8.用结构体矩阵来存储5名学生的基本情况数据, 每名学生的数据包括学号、姓名、专业和6门课程的成果.答：student(1).id='0001';student(1).name='Tom';student(1).major='computer';student(1).grade=[89,78,67,90,86,85];9.建立单位矩阵B并回答有关问题.B{1,1}=1;B{1,2}='Brenden';B{2,1}=reshape(1:9,3,3);B{2,2}={12,34,2;54,21,3;4,23,67};(1)size(B)和ndims(B)的值分别是几多？答：size(B)的值为2,2.ndims(B)的值为2.(2)B(2)和B(4)的值分别是几多？答：B(2)=147258369⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦, B(4)=[12][34][2][54][21][3][4][23][67]⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦(3)B(3)=[]和B{3}=[]执行后, B的值分别是几多？答：当执行B(3)=[]后,B={1, [1, 4, 7; 2, 5, 8; 3, 6, 9], {12, 34, 2; 54, 21, 3; 4, 23, 67}}当执行B{3}=[]后,B={1,[]; [1, 4, 7; 2, 5, 8; 3, 6, 9],{12, 34, 2; 54, 21, 3; 4, 23, 67}}习题三1.写出完成下列把持的命令.(1)建立3阶单位矩阵A.答：A=eye(3);(2)建立5×6随机矩阵A, 其元素为[100,200]范围内的随机整数.答：round(100+(200-100)*rand(5,6));(3)发生均值为1, 方差为0.2的500个正态分布的随机数.答：1+sqrt(0.2)*randn(5,100);(4)发生和A同样年夜小的幺矩阵.答：ones(size(A));(5)将矩阵A对角线的元素加30.答：A+eye(size(A))*30;(6)从矩阵A提取主对角线元素, 并以这些元素构成对角阵B.答：B=diag(diag(A));2.使用函数, 实现方阵左旋90o或右旋90o的功能.例如, 原矩阵为A, A左旋后获得B, 右旋后获得C.答：B=rot90(A);C=rot90(A,-1);3.建立一个方阵A, 求A的逆矩阵和A的行列式的值, 并验证A与A-1是互逆的.答：A=rand(3)*10;B=inv(A);C=det(A);先计算B*A, 再计算A*B, 由计算可知B*A=A*B, 即A·A-1=A-1·A是互逆.4.求下面线性方程组的解.答：A=[4,2,-1;3,-1,2;12,3,0];b=[2;10;8];x=inv(A)*b方程组的解为x=6.000026.666727.3333 -⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦5.求下列矩阵的主对角线元素、上三角阵、下三角阵、秩、范数、条件数和迹.(1)112351423052111509A-⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥⎣⎦(2)0.434328.9421B⎡⎤=⎢⎥-⎣⎦答：(1)取主对角线元素：diag(A);上三角阵：triu(A)；下三角阵：tril(A);秩：rank(A);范数：norm(A,1); 或 norm(A);或 norm(A,inf);条件数：cond(A,1); 或 cond(A,2); 或 cond(A,inf)迹：trace(A);(2)【请参考(1)】.6.求矩阵A的特征值和相应的特征向量.答：[V,D]=eig(A);习题四1.从键盘输入一个4位整数, 按如下规则加密后输出.加密规则：每位数字都加上7, 然后用和除以10的余数取代该数字；再把第一位与第三位交换, 第二位与第四位交换.答：a=input('请输入4位整数：');A=[a/1000,a/100,a/10,a];A=fix(rem(A,10));A=rem(A+7,10);b=A(3)*1000+A(4)*100+A(1)*10+A(2);disp(['加密后的值为：',num2str(b)]);2.分别用if语句和switch语句实现以下计算, 其中a、b、c的值从键盘输入.答：(1) 用if语句实现计算：a=input('请输入a的值：');b=input('请输入b的值：');c=input('请输入c的值：');x=input('请输入x的值：');y=a*x^2+b*x+c;endy=a*((sin(b))^c)+x;endy=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*x^2+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.发生20个两位随机整数, 输出其中小于平均值的偶数.答：A=fix(10+89*rand(1,20));sum=0;for i=1:20sum=sum+A(i);endB=A(find(A<(sum/20)));C=B(find(rem(B,2)==0));disp(C);4. 输入20个数, 求其中最年夜数和最小数.要求分别用循环结构和调用MATLAB 的max 函数、min 函数来实现.答：(1) 用循环结构实现：v_max=0;v_min=0;for i=1:20x=input(['请输入第', num2str(i), '数：']);if x>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(min(A))]);5. 已知：236312222s =+++++, 分别用循环结构和调用MATLAB 的sum 函数求s 的值.答：s=0;for i=0:63s=s+2^i;ends(2) 调用sum函数实现：s=0:63;s=2.^s;sum(s)6.当n分别取100、1000、10000时, 求下列各式的值.(1)111111(1)(ln2) 234nn+-+-++-+=(2)1111() 3574π-+-+=(3) 11111() 4166443n+++++=(4)224466(2)(2)133557(21)(21)2n nn nπ⎛⎫⨯⨯⨯⨯⎛⎫⎛⎫⎛⎫⎛⎫=⎪⎪⎪⎪ ⎪⨯⨯⨯-⨯+⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭要求分别用循环结构和向量运算（使用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)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/(4^k);endsum使用sum函数实现：x=[];for k=1:100x=[x, 1/(4^k)];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=[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('请输入矩阵A：');B=input('请输入矩阵B：');[x, y]=myfnc(A, B);if length(x)==0display('两矩阵的维数不匹配, 无法进行乘积运算！');elsedisp('矩阵A 和矩阵B 的乘积为：');xenddisp('矩阵A 和矩阵B 的点乘为：');y9. 先用函数的递归调用界说一个函数文件求1n m i i =∑, 然后调用该函数文件求100501021111k k k k k k ===++∑∑∑. 答：函数文件myfnc.m ：function sum=myfnc(n, m)if n<=1sum=1;elsesum= myfnc (n-1, m)+n^m;end在命令窗口中调用myfnc.m 文件, 计算100501021111k k k k k k ===++∑∑∑： sum=myfnc(100, 1)+myfnc(50, 2)+myfnc(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);endendends答：执行结果为s=108②命令文件exe.m执行后的结果为：x =41220y=。
MATLAB及其应用课后习题解答
5
“MATLAB 及其应用”课程作业
图 4- 2 画出函数图像
第5章
1． 已知椭圆的长、短轴 a 4, b 2 ，用“小红点线”画如下图所示的椭圆
x a cos t 。 y b sin t
（提示：参量 t ；点的大小；axis equal） 答 ： 运 行 :clf; a=4;b=2;t=0:pi/80:2*pi;x=a*cos(t);y=b*sin(t);plot(x,y,’r’,’markersize’,15);axis equal;xlabel(‘x’);ylabel(‘y’);shg;结果如图 5-1 所示：
1
“MATLAB 及其应用”课程作业
a5=sin(sym(pi/4)+exp(sym(0.7+pi/3))); a6=sin(sym(pi/4)+sym(exp(0.7+pi/3))); a7=sin(sym(pi/4+exp(0.7+pi/3))); a8=sym(sin(pi/4+exp(0.7+pi/3))); da2 = vpa(a1-a2,30) da3 = vpa(a1-a3,30) da4 = vpa(a1-a4,30) da5 = vpa(a1-a5,30) da6 = vpa(a1-a6,30) da7 = vpa(a1-a7,30) da8 = vpa(a1-a8,30) 结果如图 2-2 所示： 只有 a2 是精准的，观察可知，这是由 于 sym('Num')和 sym(Num)的区别带来 图 2- 2 vpa 观察误差 的前者以字符串的形式传给符号运算内 核，可以保留完整的精度；而后者经过浮点数运算之后再转换为符号类型，存在精度损 失。 3．在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是独立自由变量。 sym('sin(w*t)') , sym('a*exp(-X)' ) , sym('z*exp(j*th)') 答：分别为 w,a,z，具体图 2-3 所示。
matlab课后习题及答案详解
matlab课后习题及答案详解第1章练习题1.安装matlab时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？在安装matlab时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即matlab选项）必须安装。
第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。
2.matlab操作方式桌面存有几个窗口？如何并使某个窗口瓦解桌面沦为单一制窗口？又如何将瓦解过来的窗口再次置放至桌面上？与其他计算机语言相比较，matlab语言注重的特点就是什么？matlab系统由那些部分共同组成？在matlab操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的close按钮，一个是可以使窗口成为独立窗口的undock按钮，点击undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择dock……菜单项就可以将独立的窗口重新防止的桌面上。
matlab具备功能强大、使用方便、输出简便、库函数多样、开放性弱等特点。
matlab系统主要由开发环境、matlab数学函数库、matlab语言、图形功能和应用程序接口五个部分组成。
3.如何设置当前目录和搜寻路径，在当前目录上的文件和在搜寻路径上的文件存有什么区别？命令历史窗口除了可以观测前面键入的命令外，除了什么用途？当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的setpath菜单项来完成。
在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被matlab运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。
命令历史窗口除了用作查阅以前键入的命令外，还可以轻易执行命令历史窗口中选取的内容、将选取的内容拷贝到剪贴板中、将选取内容轻易拷贝到m文件中。
MATLAB-实用教程-课后习题标准答案
MATLAB-实用教程-课后习题标准答案第二章1.计算复数3+4i与5-6i的乘积。
a=3+4ib=5-6ic=a*b2.构建结构体Students，属性包含Name、age和Email，数据包括｛’Zhang’，18，[‘’，’’]｝、｛’Wang’，21，[]｝和｛’Li’，[],[]｝，构建后读取所有Name属性值，并且修改’Zhang’的Age 属性值为19。
Students(1).Age=18Students(1).Email='',''Students(2).Name='Wang'Students(2).Age=21Students(2).Email=[]Students(3).Name='Li'Students(3).Age=[]Students(3).Email=[]/doc/805835364.html,Student(1).Age(1)=19Student.Age3.用满矩阵和稀疏矩阵存储方式分别构造下属矩阵：A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0]A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0]S=sparse(A)S=sparse([2,1,4],[1,2,4],[1,1,1],4,5)4.采用向量构造符得到向量[1,5,9....,41].A=1:4:415.按水平和竖直方向分别合并下述两个矩阵：A=[1 0 0;1 1 0;0 01],B=[2 3 4;5 6 7;8 9 10]A=[1 0 0;1 1 0;0 0 1]B=[2 3 4;5 6 7;8 9 10]C=[A B]D=[A;B]6.分别删除第五题两个结果的第2行。
A=[1 0 0;1 1 0;0 0 1]B=[2 3 4;5 6 7;8 9 10]C=[A B]D=[A;B]C(2,:)=[]D(2,:)=[]7.分别将第5题两个结果的第2行最后3列的数值改为[11 12 13]。
matlab课后习题答案
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
因此，矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵，问：(1)A*B和A.*B的值是否相等？答：不相等。
(2)A./B和B.\A的值是否相等？答：相等。
(3)A/B和B\A的值是否相等？答：不相等。
(4)A/B和B\A所代表的数学含义是什么？答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而B\A等效于B矩阵的逆左乘A矩阵，即inv(B)*A。
3.写出完成下列操作的命令。
(1)将矩阵A第2~5行中第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转换成34矩阵。
答：reshape(x,3,4);(7)求一个字符串的ASCII码。
答：abs(‘123’); 或double(‘123’);(8)求一个ASCII码所对应的字符。
答：char(49);4.下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：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，右下角32⨯子矩阵构成矩阵D，B与C 的乘积构成矩阵E。
matlab课后习题答案
第2章 MATLAB 矩阵运算基础2.1 在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？ >> a=[5 7 3;4 9 1]2.5 计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。
>> a=[5 3 5;3 7 4;7 9 8]; >> b=[2 4 2;6 7 9;8 3 6];>> a+bans =7 7 7 9 14 13 15 12 142.6 求⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。
>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i]; >> x’ans =4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i2.7 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8]; >> a.*b ans =12 36 3 8 42 402.9 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;-0.5118 4.0427 1.33182.10 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
matlab课后习题解答
Amin=A(1);
Amax=A(1);
fori=1:19
ifA(i+1)>Amax
Amax=A(i+1);
end
ifA(i+1)<Amin
Amin=A(i+1);
end
end
Amax=max(A)
Amin=min(A)
第五题
sum=1;
fori=1:63
sum=sum+2^i;
end
s=sum;
A(5).x1=’学号’;A(5).x2=’姓名’;A(5).x3=’专业’;A(5).x4.x41=’成绩1’;……….
8、建立单元矩阵B并回答有关问题。
>> B={1;'Brenden';reshape(1:9,3,3);[12,34,2;54,21,3;4,23,67]};
>> size(B)
m1=fix(x/1000);
m2=fix((x-m1*1000)/100);
m3=fix((x-m1*1000-m2*100)/10);
m4=rem(x,10);
m1=rem((m1+7),10);
m2=rem((m2+7),10);
m3=rem((m3+7),10);
m4=rem((m4+7),10);
Size(A) ndims(A)
(5)、将向量t的0元素用机器零代替。
t=eps
(6)、将含有12个元素的向量x转换成3*4矩阵。
Reshape（x,3,4）
(7)、求一个字符串的ASCII码。
double( )
(8)、求一个ASCII码所对应的字符。
matlab课后习题答案
第2章 MATLAB 矩阵运算基础2.1 在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？ >> a=[5 7 3;4 9 1]2.5 计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。
>> a=[5 3 5;3 7 4;7 9 8];>> b=[2 4 2;6 7 9;8 3 6]; >> a+b ans =7 7 7 9 14 13 15 12 142.6 求⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。
>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i];>> x’ ans =4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i2.7 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8]; >> a.*b ans =12 36 3 8 42 402.9 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;-0.5118 4.0427 1.33182.10 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
MATLAB课后实验答案
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin851z e =+(2) 221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=-- (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解： M 文件:z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5]; z2=1/2*log(x+sqrt(1+x^2)) a=-3.0:0.1:3.0;z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
解：(1) 结果：m=100:999;n=find(mod(m,21)==0);length(n)ans =43(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：ch='ABC123d4e56Fg9';k=find(ch>='A'&ch<='Z');ch(k)=[]实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
matlab课后习题参考答案
Matlab课后习题部分参考答案习题二2．student=struct('name',{'zhang','wang','li'},'age',{18,21,[] },'email',{['zh ang@','zhang@'],'',''})5、a=[1 0 0;1 1 0;0 0 1];b=[2 3 4;5 6 7;8 9 10];c=[a b]d=[a;b]c =1 0 023 41 1 0 5 6 70 0 1 8 9 10d =1 0 01 1 00 0 12 3 45 6 78 9 1010、reshape(c,2,9)1 0 1 0 1 5 3 9 71 0 0 02 8 6 4 10 reshape(d,2,9)ans =1 0 5 0 0 6 0 1 71 2 8 1 3 9 0 4 10 (注意：重新排列矩阵，是将原来的矩阵按列排序)11、a.’12、a+b a.*ba\b (即inv(a)*b，考试时，要求能自己手算出a的逆)ans =2 3 43 3 38 9 1013、c=4*ones(3);a>=cb>=c14．strcat('The picture is ','very good')a=['The picture is ','very good']18．double('very good')1．norm（A），det（A），rank（A）2．b=[1 1 1 1 1];c=b';a=[17,24,1,8,50;23,5,7,14,49;4,6,13,20,43;10,12,19,21,62;11,18,25,2,56]; x=a\c习题四1、a=14*rand(1,10)-5;for i=1:9max=i;for j=i+1:10if a(j)>a(max)max=j;endendtemp=a(i);a(i)=a(max);a(max)=temp;enda将上述脚本保存为script41.m2、function y=function42(n)3*randn(1,n)+3;for i=1:n-1max=i;for j=i+1:nif a(j)>a(max)max=j;endendtemp=a(i);a(i)=a(max);a(max)=temp;endy=a;将上述函数保存为function42.m 3、r=input(‘please input 1/2: ’);if r= =1script41;elsen=input(‘please input the length: ’);y=funtion42(n);end4．function f=function1(x,y)if y==1f=sin(x);elseif y==2f=cos(x);elsef=sin(x).*cos(x);end在命令窗口调用函数f=function1(1,3) 5．function result=function3(x,n)result=0;for i=1:nresult=result+sin(x.*i)+(-1).^i.*cos(i.*x); endfunction result=function4(n,t)result=0;A=[1 2 3;0 1 2;0 0 1];for i=1:nresult=result+A*i.*exp(A*i.*t);end习题五3．x1=-2:0.01:2;x2=-2:0.01:2;y1=x1.*sin(x2);y2=x2.*cos(x1);plot3(x1,x2,y 1,'-d',x1,x2,y2,'-d')14、x=8*rand(1,100)-2;y=reshape(x,10,10);save mydata.mat yclearclcload mydatak=y>ones(10);totel=0;for i=1:100if k(i)==1totel=totel+1;endend15、x=rand(5); %生成5*5均匀分布的随机矩阵（该处只作举例，未按题目要求）fid=fopen(‘text.txt’,’w’);count=fwrite(fid,x,’int32’); closestatus=fclose(fid);清除内存，关闭所有窗口fid=fopen(‘text.txt’,’r’);x=fread(fid,[5,5],’int32’); closestatus=fclose(fid); inv(x)16．随便生成一个矩阵，x1=-1:0.2:0.8;x2=-1:0.2:0.8;y1=x1.*sin(x2)；y2=x2.*cos(x1)；a=reshape([y1;y2],10,10)；fid=fopen('table.txt','w'); fprintf (fid,‘%f’，a); fclose(fid) ;清除内存，关闭所有窗口a=fscanf('table.txt','%f') ; exp(a) ;fclose(fid) ;习题六1.c=ploy2str(A)；B=[2,0,0,1,3,5]；2.x=1:10;B=ployval(A，x)；。
MATLAB课后习题集附标准答案
第2章MATLAB概论1、与其他计算机语言相比较，MA TLAB 语言突出的特点是什么？答：起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富，可扩展性强.2、MA TLAB 系统由那些部分组成？答：开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB 时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装. 第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可. 矚慫润厲钐瘗睞枥庑赖。
4、MATLAB 操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？聞創沟燴鐺險爱氇谴净。
答：在MATLAB 操作桌面上有五个窗口，在每个窗口的右下角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口称为独立的Undock 按钮，点击Undock 按钮就可以使该窗口脱离桌面称为独立窗口，在独立窗口的view 菜单中选择Dock，菜单项就可以将独立的窗口重新防止的桌面上.残骛楼諍锩瀨濟溆塹籟。
5、如何启动M 文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M 文件编辑/调试器将被启动.在命令窗口中键入edit 命令时也可以启动M 文件编辑/调试器.酽锕极額閉镇桧猪訣锥。
6、存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可.彈贸摄尔霁毙攬砖卤庑。
7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中. 謀荞抟箧飆鐸怼类蒋薔。
matlab第四章课后作业解答
matlab第四章课后作业解答第四章习题解答1、求下列多项式的所有根，并进⾏验算。
（3）267235865x x x x-+-（4）4)32(3-+x 解：>> p=zeros(1,24);>> p(1)=5;p(17)=-6;p(18)=8;p(22)=-5;>> root=roots(p)root =0.97680.9388 + 0.2682i0.9388 - 0.2682i0.8554 + 0.5363i0.8554 - 0.5363i0.6615 + 0.8064i0.6615 - 0.8064i0.3516 + 0.9878i0.3516 - 0.9878i-0.0345 + 1.0150i-0.0345 - 1.0150i-0.4609 + 0.9458i-0.4609 - 0.9458i-0.1150 + 0.8340i-0.1150 - 0.8340i-0.7821 + 0.7376i-0.7821 - 0.7376i-0.9859 + 0.4106i-0.9859 - 0.4106i-1.0416-0.7927>> polyval(p,root)ans =1.0e-012 *-0.07120.0459 - 0.0081i0.0459 + 0.0081i-0.0419 + 0.0444i-0.0419 - 0.0444i0.0509 + 0.0929i0.0509 - 0.0929i-0.2059 + 0.0009i-0.2059 - 0.0009i-0.0340 + 0.0145i-0.0340 - 0.0145i0.1342 + 0.0910i0.1342 - 0.0910i0.0025 + 0.0027i0.0025 - 0.0027i-0.0077 + 0.4643i-0.0077 - 0.4643i-0.3548 - 0.1466i-0.3548 + 0.1466i-0.0251-0.0073(4) >> p1=[2 3];>> p=conv(conv(p1,p1),p1)-[0 0 0 4];>> root=roots(p)root =-1.8969 + 0.6874i-1.8969 - 0.6874i-0.7063>> polyval(p,root)ans =1.0e-014 *-0.7105 - 0.6217i-0.7105 + 0.6217i6、求解下列⽅程组在区域1,0<<βα内的解-=+=.sin 2.0cos 7.0,cos 2.0sin 7.0βαββαα解：以初值)5.0,5.0(),(00=βα进⾏求解>> fun=inline('[0.7*sin(x(1))+0.2*cos(x(2))-x(1),0.7*cos(x(1))-0.2*sin(x(2))-x(2)]');>> [x,f,h]=fsolve(fun,[0.5 0.5])Optimization terminated: first-order optimality is less than options.TolFun.x =0.5265 0.5079f =1.0e-007 *-0.1680 -0.2712h =1因⽽，该⽅程组的近似根为5079.0,5265.0==βα。
matlab课后习题答案(1-9章)要点
1 数字1.5e2，1.5e3 中的哪个与1500相同吗？1.5e32 请指出如下5个变量名中，哪些是合法的？abcd-2xyz_3 3chan a 变量ABCDefgh 2、5是合法的。
3 在MATLAB 环境中，比1大的最小数是多少？ 1+eps4 设 a = -8 , 运行以下三条指令，问运行结果相同吗？为什么？w1=a^(2/3) w2=(a^2)^(1/3) w3=(a^(1/3))^2w1 = -2.0000 + 3.4641i ;w2 = 4.0000 ;w3 =-2.0000 + 3.4641i 5 指令clear, clf, clc 各有什么用处？clear 清除工作空间中所有的变量。
clf 清除当前图形。
clc 清除命令窗口中所有显示。
第二章1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象？3/7+0.1双; sym(3/7+0.1)符; sym('3/7+0.1') 符;; vpa(sym(3/7+0.1)) 符;2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') symvar(sym('sin(w*t)'),1) w a z3 （1）试写出求三阶方程05.443=-x 正实根的程序。
注意：只要正实根，不要出现其他根。
（2）试求二阶方程022=+-a ax x 在0>a 时的根。
（1）reset(symengine)syms x positive solve(x^3-44.5) ans =(2^(2/3)*89^(1/3))/2（2）求五阶方程022=+-a ax x 的实根 syms a positive %注意：关于x 的假设没有去除 solve(x^2-a*x+a^2)Warning: Explicit solution could not be found. > In solve at 83 ans =[ empty sym ]syms x clear syms a positivesolve(x^2-a*x+a^2) ans =a/2 + (3^(1/2)*a*i)/2 a/2 - (3^(1/2)*a*i)/24 观察一个数（在此用@记述）在以下四条不同指令作用下的异同。
matlab课后习题答案到章
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
因此，矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵，问：(1)A*B和A.*B的值是否相等答：不相等。
(2) A./B和B.\A的值是否相等答：相等。
(3)A/B和B\A的值是否相等答：不相等。
(4)A/B和B\A所代表的数学含义是什么答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而B\A等效于B矩阵的逆左乘A矩阵，即inv(B)*A。
3.写出完成下列操作的命令。
(1)将矩阵A第2~5行中第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转换成34⨯矩阵。
答：reshape(x,3,4);(7)求一个字符串的ASCII码。
答：abs(‘123’); 或double(‘123’);(8)求一个ASCII码所对应的字符。
答：char(49);4.下列命令执行后，L1、L2、L3、L4的值分别是多少A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：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，右下角32⨯子矩阵构成矩阵D，B与C的乘积构成矩阵E。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
0.80
0.95
电阻y
15
18
19
21
22.6
23.8
26
>> x=[0.1,0.3,0.4,0.55,0.7,0.8,0.95];
y=[15,18,19,21,22.6,23.8,26];
p1=polyfit(x,y,1);
p3=polyfit(x,y,3);
p5=polyfit(x,y,5);
if x>=90
disp('优秀');
elseif x>=80
disp('良好');
elseif x>=60
disp('及格');
else
disp('不及格');
end
>> x=85
x =
85
良好
Q3：编写函数，计算
>> sum=0;
>> for i=1:50
a=1;
for j=1:i
a=a*j;
end
ans =
0
(4)
>> syms n
>> limit(sqrt(n+2)-2*sqrt(n+1)+sqrt(n),n,inf)
ans =
0
Q2:用MATLAB软件求下列函数极限：
(1)
>> syms x
>> limit((((1+x)^(1/3)-1)/x),x,0)
ans =
1/3
(4)
>> syms x
(2)
>> syms x
>> f=sqrt(x+2)*(3-x)^4/(x+1)^5;
(2)
>> a=[1/sqrt(2),-1/sqrt(2),0,0];
>> b=[-1/sqrt(2),1/sqrt(2),0,0];
>> c=[0,0,1/sqrt(2),-1/sqrt(2)];
>> d=[0,0,-1/sqrt(2),1/sqrt(2)];
>> A=[b;c;d];
>> [R,p]=chol(A)
P16
Q2:计算表达式 在 和 时的函数值。
functiony=jie(x)
y=tan(-x.^2)*acos(x);
>> jie(0.25)
ans =
-0.0825
>> jie(0.78*pi)
ans =
0 + 0.4418i
Q3：编写M命令文件，求 的值。
a=0;b=0;
fori=1:50
a=a+i*i;
>> A=[A;0,0]
A =
1 2
3 4
0 0
>> B=[0;0;0];
>> A=[A,B]
A =
1 2 0
3 4 0
0 0 0
(2)
>> A(3,:)=[1 3 5]
A =
1 2 0
3 4 0
1 3 5
(3)
>> a=A(1,2)
a =
2
>> b=A(3,2)
b =
3
Q10:已知矩阵A= ，B= ，求A+B，A-B，AB，BA， ， 。
0 1 0 0 -2
0 0 1 0 3
0 0 0 1 5
增广矩阵的秩为4，等于系数矩阵的秩，等于未知数的个数。
所以有唯一解。
Q23:通过矩阵LU分解求解矩阵方程AX=b,其中A= ，b=
>> A=[1 0 2 0;0 1 0 1;1 2 4 3;0 1 0 3];
>> [L,U]=lu(A);
>> b=[1;2;-1;5];
ps =
4 x + 5
>> pr=poly2str(pr,'x')
pr =
8 x + 9
以上两个多项式的商为 ，余子式为pr=8 x + 9.
Q8:在钢线碳含量对于电阻的效应的研究中，得到以下数据。分别用一次、三次、五次多项式拟合曲线来拟合这组数据并画出图形。
碳含量x
0.10
0.30
0.40
0.55
0 0 1.7321 -3.4641
0 0 0 2.0000
验证 ：
>> R'*R
ans =
1.0000 -1.0000 2.0000 1.0000
-1.0000 3.0000 0 -3.0000
2.0000 0 9.0000 -6.0000
1.0000 -3.0000 -6.0000 19.0000
>> S=sparse(A)
S =
(1,1) 1
(3,2) 1
(2,3) 2
(1,4) -1
(3,5) 3
>> A1=full(S)
A1 =
1 0 0 -1 0
0 0 2 0 0
0 1 0 0 3
(3)
>> A=[1 0 0 0 2;0 0 0 3 0;0 0 1 0 0;0 3 0 0 0;2 0 0 0 1];
>> x=U\(L\b)
x =
8.5000
0.5000
-3.7500
1.5000
Q25:用QR方法求解下列方程组，然后用其他方法验证解的正确性。
(1)
>> A=[5 4 5;7 8 9;12 3 8];
>> b=[1;2;3];
>> [Q,R]=qr(A)
Q =
-0.3386 -0.2552 -0.9057
end
forj=1:10
b=b+1/j;
end
c=a+b;
>> c
c =
4.2928e+004
P27
Q2：矩阵 , ,计算 , ,并比较两者的区别。
>> A=[1 2 3;4 5 6;7 8 9];
>> B=[4 6 8;5 5 6;3 2 2];
>> A*B
ans =
23 22 26
59 61 74
>> plot(x,y1,'b*',x,y2,'r-.',x,y3,'r-.')
P113
Q8:已知矩阵 ，实现下列操作：
（1）添加零元素使之成为一个 的方阵。
（2）在以上操作的基础上，将第三行元素替换为（1 3 5）。
（3）在以上操作的基础上，提取矩阵中第2个元素以及第3行第2列的元素。
(1)
>> A=[1,2;3,4];
>> S=sparse(A)
S =
(1,1) 1
(5,1) 2
(4,2) 3
(3,3) 1
(2,4) 3
(1,5) 2
(5,5) 1
>> A1=full(S)
A1 =
1 0 0 0 2
0 0 0 3 0
0 0 1 0 0
0 3 0 0 0
2 0 0 0 1
Q20:求解下列方程
(1)
>> A=[1 0 3 10;2 1 4 18;1 -1 2 3];
y3=polyval(p3,x1);
y5=polyval(p5,x1);
plot(x,y,'rp',x1,y1,x1,y3,x1,y5);
legend('拟合点','一次拟合','三次拟合','五次拟合')
一阶拟合函数
f1 =
12.5503 x + 13.9584
三阶拟合函数
f3 =
8.9254 x^3 - 14.6277 x^2 + 19.2834 x + 13.2132
-0.4741 -0.7851 0.3985
-0.8127 0.5643 0.1449
R =
-14.7648 -7.5856 -12.4621
0 -5.6088 -3.8275
0 0 0.2174
>> x=R\(Q\b)
x =
-0.5000
-1.0000
1.5000
检验：
>> [L,U]=lu(A);
-0.8443 0.5338 -0.0473
R =
-14.2127 -7.3174 -12.2426
0 -5.9544 -4.4363
0 0 -0.6617
>> x=R\(Q\b)
x =
-1.8214
-2.8571
4.1786
检验：
>> [L,U]=lu(A);
>> x=U\(L\b)
x =
-1.8214
五阶拟合函数
f5 =
146.1598 x^5 - 386.879 x^4 + 385.5329 x^3 - 178.8558 x^2 + 49.9448 x
+ 11.4481
P165
Q1:用MATLAB软件求下列数列极限：
(1)
>> syms n
>> limit(((-2)^n+3*n)/((-2)^(n+1)+3^(n+1)),n,inf)