中南大学matlab课后答案-第九章
实验指导
1,
>> figure('Color',[1,0,0],'WindowButtonDownFcn','text(0.5,0.5,''Left Button Pressed'')');axis off 2,
(1)默认属性
>> x=linspace(0,600,100000);y=(log(x+sqrt(1+x.*x)))/2;line(x,y)
句柄操作
>> x=linspace(0,600,100000);
y=(log(x+sqrt(1+x.*x)))/2;
line(x,y,'linewidth',3,'linestyle','-.','color','r');
text(300,3.5,'文字标注')
(2)默认属性
>> t=linspace(0,20,500);x=t.*t;y=5.*t.*t;line(x,y)
句柄操作
>> t=linspace(0,20,500);
x=t.*t;y=5.*t.*t;line(x,y);
line(x,y,'linewidth',3,'linestyle','-.','color','r');
text(250,1600,'文字标注')
3,
(1) >> x=linspace(0,50,1000);
[x,y]=meshgrid(x);
z=x.*x+y.*y-5.*sin(x.*y);
axes('view',[-37.5,30]);
surf(x,y,z);
light('position',[10,20,4000]);
shading interp
(2) >> x=linspace(0,50,1000);
[x,y]=meshgrid(x);
z=y.*y.*y;
axes('view',[-37.5,30]);
surf(x,y,z);
light('position',[30,20,1300000]);
shading interp
4,
>> x=-2*pi:0.01:2*pi;
y1=sin(x);
y2=cos(x);
axes('position',[0.1,0.6,0.2,0.2]);
plot(x,y1);
axes('position',[0.6,0.6,0.2,0.2]);
plot(x,y2);
axes('position',[0.1,0.1,0.2,0.2]);
fplot('tan(x)',[-1.5,1.5]);
axes('position',[0.6,0.1,0.2,0.2]);
fplot('cot(x)',[0.1,3])
5,
>> t=linspace(0,2*pi,6);
k=t(1:length(t)-1);
patch(sin(k),cos(k),1:length(k),'facecolor','interp');
colormap summer;
axis equal;
思考练习
3,
(1)默认属性
>> x=linspace(0,50,1000);y=(1+x.^2)/(1+x.^4);line(x,y)
句柄属性
>> x=linspace(0,50,1000);
y=(1+x.^2)/(1+x.^4);
line(x,y,'linewidth',3,'linestyle','-.','color','b');
text(25,1,'文字标注')
(2)默认属性
>> x=linspace(0,800,5000);y=3.*x+sin(x)-exp(x);axes('XLim',[0,7],'YLim',[-10,5]);line(x,y) 句柄属性
>> x=linspace(0,800,5000);
y=3.*x+sin(x)-exp(x);
axes('XLim',[0,7],'YLim',[-10,5]);line(x,y,'linewidth',3,'linestyle','-.','color','c');
text(3,3,'文字标注')
4,
(1) >> t=linspace(0,2*pi,100);x=cos(t);y=sin(t);z=t;axes('view',[-37.5,30]);plot3(x,y,z);grid
(2) >> x=linspace(0,3,100);
[x,y]=meshgrid(x);
z=x.*y.*exp(-1.*x.*x-y.*y);
axes('view',[-37.5,30]);
surf(x,y,z);
light('position',[0.1,0.1,0.5]);
shading interp
5,
>> [x,y,z]=cylinder(0.5,30);surf(x,y,z);axis equal;light('position',[-0.5,-0.4,0]);shading interp。
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课后习题答案
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. 已知23100.7780414565532503269.5454 3.14A --??=-??完成下列操作：(1) 取出A 的前3行构成矩阵B ，前两列构成矩阵C ，右下角32?子矩阵构成矩阵D ，B 与C 的乘积构成矩阵E 。
中南大学材料院matlab操作题集答案
的491操作题集操作题1 《Matlab与科学计算》P28 例2.1要求计算水在温度为0℃，20℃，40℃，60℃，80℃时的黏度，已知水的黏度随/(1+at+bt2)温度的变化公式为μ=μ为0℃水的黏度，值为1.758*10-3。
其中μ要求分别以长短格式显示计算结果。
操作题2 《Matlab与科学计算》P31 例2.2数值数组【1：5】和字符串转换示例。
操作题3《Matlab与科学计算》P44 例2.9比较用左除和右除法分别求解恰定方程的解。
操作题4 《Matlab与科学计算》 P48 例2.14计算矩阵magic(3)的指数，并比较不同函数的结果。
指数运算expm(a)操作题5 《Matlab与科学计算》 P50 例2.18计算下面矩阵[-149 -50 -154；537 180 546；-27 -9 -25]的特征值条件数。
特征值函数condeig(A)操作题6《Matlab与科学计算》P62 例2.29对矩阵a = 1 1 1 11 2 3 41 3 6 101 4 10 20 进行三角抽取使之分别变成0 0 0 0 1 1 1 01 0 0 0 及 123 41 3 0 0 1 3 6 101 4 10 0 1 4 10 20操作题7《Matlab与科学计算》P62 例2.30用3种方法建立向量v 为[1 2 6 20]多项式的伴随矩阵。
操作题8《Matlab与科学计算》P64 例2.31已知a=[2 1 -3 -1；3 1 0 7；-1 2 4 -2；1 0 -1 5]；求该矩阵的3次方及数组的3次方。
操作题9《Matlab与科学计算》P66 例2.32 数组逻辑运算演示。
已知a=[1:3；4：6；7：9]，b=[0 1 0；1 0 1；0 0 1]， x=5, y=ones （3）*5,求x<=a，ab=a&b，～b操作题10《Matlab与科学计算》P71 例2.41对多项式p =[2 -5 6 -1 9]求根，计算其在x=3的值，进行微分计算并写出其表达式。
matlab第九章答案
第九章5、>> x=[0:1:24];y=[15,14,14,14,14,15,16,18,20,22,23,25,28,31,32,31,29,27,25,24,22,20,18,17,16]; >> p=polyfit(x,y,2)p =-0.0936 2.5943 8.4157拟合为为二次函数形式为：-0.0936*x^2+2.5943*x+8.41576、>> clear;t=[0.5 1 2 3 4 5 7 9];V=[6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63];plot(t,V,'o');fun1=inline('10-(10-a(1))*exp(-t/a(2))','a','t');format short g;a=nlinfit(t,V,fun1,[5 3]);x=min(t):0.1:max(t);y=10-(10-a(1))*exp(-x/a(2));plot(t,V,'o',x,y)aa =5.5577 3.5002012345678966.577.588.599.5108、 >> data=[2 106.42;3 109.20;4 109.58;5 109.50;7 110;8 110.49;11 110.59;14 110.6;15 110.9;15 110.76;18 111;19 111.2];>> h=data(:,1);w=data(:,2);>> subplot(1,3,1);plot(h,w,'o');01020106107108109110111112由散点图可以看出线性关系较明显。
x1=[2 3 4 5 7 8 10 11 14 15 15 18 19]';y=[106.42 109.20 109.58 109.50 110.00 109.93 110.49 110.59 110.60110.90 110.76 111.00111.20]';x=[ones(13,1),x1];[b,bint,r,rint,stats]=regress(y,x,alpha);b,bint,stats,rcoplot(r,rint)b =108.25810.1742bint =107.2794 109.23670.0891 0.2593stats =0.6484 20.2866 0.0009 0.5965线性相关系数较小，线性回归模型在alpha>0.0009成立第一个点为异常点（仅指线性模型下），予以剔除>> alpha=0.05;x1=[3 4 5 7 8 10 11 14 15 15 18 19]';y=[109.20 109.58 109.50 110.00 109.93 110.49 110.59 110.60 110.90 110.76 111.00 111.20]'; x=[ones(12,1),x1];[b,bint,r,rint,stats]=regress(y,x,alpha);b,bint,stats,rcoplot(r,rint)b =109.06680.1159bint =108.8264 109.30720.0958 0.1360stats =0.9428 164.8060 0.0000 0.0267剔除第一个点后线性系数和p值都变得好了很多。
matlab课后习题答案(1-9章)
1 数字1.5e2，1.5e3 中的哪个与1500相同吗？1.5e32 请指出如下5个变量名中，哪些是合法的？abcd-2xyz_33chan 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 观察一个数（在此用@记述）在以下四条不同指令作用下的异同。
第9章MATLAB符号计算习题答案.doc
第9章 MATLAB符号计算习题9一、选择题1．设有a=sym(4)。
则1/a+1/a的值是（）。
BA．0.5 B．1/2 C．1/4+1/4 D．2/a2．函数factor(sym(15))的值是（）。
DA．'15' B．15 C．[ 1, 3, 5] D．[ 3, 5]3．在命令行窗口输入下列命令：>> f=sym(1);>> eval(int(f,1,4))则命令执行后的输出结果是（）。
AA．3 B．4 C．5 D．14．MA TLAB将函数展开为幂级数，所使用的函数是（）。
DA．tailor B．tayler C．diff D．taylor5．MATLAB用于符号常微分方程求解的函数是（）。
CA．solve B．solver C．dsolve D．dsolver二、填空题1．在进行符号运算之前首先要建立，所使用的函数或命令有和。
符号对象，sym，syms2．对于“没有定义”的极限，MATLAB给出的结果为；对于极限值为无穷大的极限，MA TLAB给出的结果为。
NaN，Inf3．在命令行窗口输入下列命令：>> syms n;>> s=symsum(n,1,10)命令执行后s的值是。
554．在MATLAB中，函数solve(s,v)用于代数方程符号求解，其中s代表，v 代表。
符号代数方程，求解变量5．在MA TLAB符号计算中y的二阶导数表示为。
D2y三、应用题1．分解因式。
（1）x9-1 （2）x4+x3+2x2+x+1（3）125x6+75x4+15x2+1 （4）x2+y2+z2+2(xy+yz+zx)(1):2x=sym('x'); A=x^9-1; factor(A) (2):x=sym('x');B=x^4+x^3+2*x^2+x+1; factor(B) 2．求函数的极限。
（1）4586lim 22++x x x x --4→x （2）xx -0→x lim(1):x=sym('x');A=(x^2-6*x+8)/(x^2-5*x+4); limit(A,x,4) (2):x=sym('x'); B=abs(x)/x; limit(B)3．求函数的符号导数。
matlab课后习题及答案
第一章5题已知a=4.96，b=8.11，计算)ln(b a e b a +-的值。
解：clearclca=4.96;b=8.11;exp(a-b)/log(a+b)ans =0.01676题已知三角形的三边a=9.6,b=13.7,c=19.4,求三角形的面积。
提示：利用海伦公式area =))()((c s b s a s s ---计算，其中S=(A+B+C)/2.解：clearclca=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 的值解：clearclcS=0;for i=0:1:63S=S+2^i;endSS =1.8447e+0199题分别用for 和while 循环结构编写程序，计算∑=-1001n 1n 2）（的值。
解：clearclcs=0;for n=1:100s=s+(2*n-1);endss =10000clearclcn=1;s=0;while n<=100s=s+(2*n-1);n=n+1;endss =10000第三章2题在同一坐标下绘制函数x ,,2x -,2x xsin(x)在()∏∈,0x 的曲线。
解：clearclcx=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,'-')00.51 1.52 2.53-10-8-6-4-22468109题用不同的线型和颜色在同一坐标内绘制曲线y1=2e x 5.0 、y2=sin(2∏x )的图形。
解：clearclcx=0:0.01:pi;y1=2*exp(-0.5*x);y2=sin(2*pi*x);plot(x,y1,'b:',x,y2,'r--')00.51 1.52 2.53 3.5-1-0.50.511.52第四章1题对表达式f=3238641+++x x x 进行化简。
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课后习题答案
习题二之阿布丰王创作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的值：');if x>=0.5 & x<1.5y=a*x^2+b*x+c;endif x>=1.5 & x<3.5y=a*((sin(b))^c)+x;endif x>=3.5 & x<5.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*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课后习题答案
习题二之阳早格格创做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.问：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].6.当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]7.用结构体矩阵去保存5名教死的基础情况数据，每名教死的数据包罗教号、姓名、博业战6门课程的结果.问：student(1).id='0001';student(1).name='Tom';student(1).major='computer';student(1).grade=[89,78,67,90,86,85];8.修坐单元矩阵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 的值.问：(1) 用循环结构真止：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)(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/(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('二矩阵的维数没有匹配，无法举止乘积运算！');disp('矩阵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问：真止截止为②下令文献exe.m真止后的截止为：x =41220y=2 4 6。