MATLAB数据分析课后习题答案
.
.
第三章课后习题第二题
(1)做出x,y的散点图，确定所要的回归模型
x1=0.1:0.01:0.18;
x2=[x1,0.2,0.21,0.23]';
y=[42.0,41.5,45.0,45.5,45.5,47.5,49.0,55.0,50.0,55.5,55.5,6
0.5]';
plot(x2,y,'*')
（2）做回归分析
x=[ones(12,1),x2];
[b,bint,r,rint,stats]=regress(y,x)
运行结果
.
.
.
.
结果表明：决定系数R2=0.9229接近1，F1-0.05（1，12-2）
=0.0041<119.7470
P=0.0000<0.05
故回归模型：y=26.8909+141.7415x成立
matlab课后习题答案
第4章数值运算习题 4 及解答1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff或数值梯度gradient指令计算)(ty'曲线绘制y'，然后把)(t y和)(t 在同一张图上，观察数值求导的后果。
（模拟数据从prob_data401.mat获得）〖目的〗●强调：要非常慎用数值导数计算。
●练习mat数据文件中数据的获取。
●实验数据求导的后果●把两条曲线绘制在同一图上的一种方法。
〖解答〗（1）从数据文件获得数据的指令假如prob_data401.mat文件在当前目录或搜索路径上clearload prob_data401.mat（2）用diff求导的指令dt=t(2)-t(1);yc=diff(y)/dt; %注意yc的长度将比y短1plot(t,y,'b',t(2:end),yc,'r')（3）用gradent求导的指令（图形与上相似）dt=t(2)-t(1);yc=gradient(y)/dt;plot(t,y,'b',t,yc,'r') grid on〖说明〗● 不到万不得已，不要进行数值求导。
● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级以上。
● 求导会使数据中原有的噪声放大。
2 采用数值计算方法，画出dt ttx y x⎰=0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。
〖提示〗● 指定区间内的积分函数可用cumtrapz 指令给出。
● )5.4(y 在计算要求不太高的地方可用find 指令算得。
〖目的〗● 指定区间内的积分函数的数值计算法和cumtrapz 指令。
● find 指令的应用。
〖解答〗dt=1e-4; t=0:dt:10;t=t+(t==0)*eps; f=sin(t)./t;s=cumtrapz(f)*dt;plot(t,s,'LineWidth',3) ii=find(t==4.5); s45=s(ii) s45 =3 求函数xex f 3sin )(=的数值积分⎰=π0 )(dx x f s ，并请采用符号计算尝试复算。
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=。
第2章--MATLAB数据及其运算-习题答案教学内容
第2章--M A T L A B数据及其运算-习题答案第2章 MATLAB数据及其运算习题2一、选择题1．下列可作为MATLAB合法变量名的是（）。
D A．合计 B．123 C．@h D．xyz_2a 2．下列数值数据表示中错误的是（）。
CA．+10 B．1.2e-5 C．2e D．2i3．使用语句t=0:7生成的是（）个元素的向量。
A A．8 B．7 C．6 D．54．执行语句A=[1,2,3;4,5,6]后，A(3)的值是（）。
B A．1 B．2 C．3 D．45．已知a为3×3矩阵，则a(:,end)是指（）。
D A．所有元素 B．第一行元素C．第三行元素 D．第三列元素6．已知a为3×3矩阵，则运行a (1)=[]后（）。
A A．a变成行向量 B．a变为2行2列C．a变为3行2列 D．a变为2行3列7．在命令行窗口输入下列命令后，x的值是（）。
B >> clear>> x=i*jA．不确定 B．-1 C．1 D．i*j 8．fix(354/100)+mod(354,10)*10的值是（）。
D A．34 B．354 C．453 D．439．下列语句中错误的是（）。
BA．x==y==3 B．x=y=3C．x=y==3 D．y=3,x=y10．find(1:2:20>15)的结果是（）。
CA．19 20 B．17 19C．9 10 D．8 911．输入字符串时，要用（）将字符括起来。
C A．[ ] B．{ } C．' ' D．" " 12．已知s='显示"hello"'，则s的元素个数是（）。
A A．9 B．11 C．7 D．1813．eval('sqrt(4)+2')的值是（）。
BA．sqrt(4)+2 B．4 C．2 D．2， 214．有3×4的结构矩阵student，每个结构有name（姓名）、scores（分数）两个成员，其中scores是以1×5矩阵表示的5门课的成绩，那么要删除第4个学生的第2门课成绩，应采用的正确命令是（）。
matlab第七章课后题答案
matlab第七章课后题答案第⼀题分解因式syms x y z>> A=x^9-1;>> factor(A)ans =(x-1)*(x^2+x+1)*(x^6+x^3+1) 解（1）>> B=x^4+x^3+2*x^2+x+1;>> factor(B)ans =(x^2+1)*(x^2+x+1) 解（2）> C=125*x^6+75*x^4+15*x^2+1;>> factor(C)ans =(5*x^2+1)^3 解（3）> D=x^2+y^2+z^2+z*(x*y+y*z+z*x);>> factor(D)ans =x^2+y^2+z^2+z*x*y+y*z^2+z^2*x 解（4）第⼆题化简表达式syms x y a b>> s=y/x+x/y;>> simplify(s)ans =(x^2+y^2)/x/y 解（1）s=sqrt(a+sqrt(a^2-b))/2+sqrt(a-sqrt(a^2-b))/2; ans =1/2*(a+(a^2-b)^(1/2))^(1/2)+1/2*(a-(a^2-b)^(1/2))^(1/2) 解（2）s=2*cos(x)^2*x-sin(x)^2*x;>> simplify(s)ans =x*(3*cos(x)^2-1) 解（3）s=sqrt(3+2*(sqrt2))第三题求函数的极限> syms x>> f=(x^2-6*x+8)/(x^2-5*x+4);> limit(f,x,4)ans =2/3 解（1）>> f=abs(x)/x;>> limit(f,x,0)ans =NaN 解（2）f=(sqrt(1+x^2)-1)/x;>> limit(f,x,0)ans =0 解（3）f=(x+1/x)^x;>> limit(f,x,inf,'left')ans =Inf 解（4）第四题求函数的符号导数f=3*(x^2)-5*x+1;>> diff(f)ans =6*x-5 解(1)y’>> diff(f,x,2)ans =6 解（1）> y=sqrt(x+sqrt(x+sqrt(x)));>> diff(y)ans =1/2/(x+(x+x^(1/2))^(1/2))^(1/2)*(1+1/2/(x+x^(1/2))^(1/2)*(1+1/2/x^(1/2))) 解（2）diff(y,x,2)ans =-1/4/(x+(x+x^(1/2))^(1/2))^(3/2)*(1+1/2/(x+x^(1/2))^(1/2)*(1+1/2/x^(1/2)))^2+1/2/(x+(x+x^(1/2) )^(1/2))^(1/2)*(-1/4/(x+x^(1/2))^(3/2)*(1+1/2/x^(1/2))^2-1/8/(x+x^(1/2))^(1/2)/x^(3/2)) 解(2)y=sin(x)-x^2/2;> diff(y)ans =cos(x)-x 解（3）>> diff(y,x,2)ans =-sin(x)-1 解（3）syms x y z>> z=x+y-sqrt(x^2+y^2);>> diff(z,x,y)ans =1-1/(x^2+y^2)^(1/2)*y 解（4）>> diff(y,x)ans =0 解（4）第五题求不定积分x=sym('x');>> f=1/(x+a);>> int(f)ans =log(x+a) 解（1）>> f=(1-3*x)^3;>> int(f)ans =-1/12*(1-3*x)^4 解（2）>> f=(1/(sin(x)^2*cos(x)^2));>> int(f)ans =1/sin(x)/cos(x)-2/sin(x)*cos(x) 解（3）>> f=x^2/(sqrt(a^2+x^2));>> int(f)ans =1/2*x*(a^2+x^2)^(1/2)-1/2*a^2*log(x+(a^2+x^2)^(1/2)) 解（4）第六题求定积分> x=sym('x');> int((x*(2-sin(x)^2))^12,0,1)ans =-13072167041243000966100527033032931/1439431206610157332070400000000000*sin(1)^19 *cos(1)-63988617583073709724938474490679159346183608452999323027852007/820274272 498737105178830959457441284892917760000000000*cos(1)*sin(1)+417844027386435896683 78350956709518241555640967463723429593/1230411408748105657768246439186161927339 376640000000000*sin(1)^3*cos(1)-6287598784304532394386769772554886862775718358540 017487607/1794349971090987417578692723813152810703257600000000000*sin(1)^7*cos(1)-117903417317/3522410053632*sin(1)^23*cos(1)+93129118771020938708526771772323524014639/1529945519744217245425650892800000000000*sin(1)^17*cos(1)-2677966496932891906 2789407028617562008014032541099/6567698857381017597043504106365255680000000000* sin(1)^10-90936661567370530014104030508869215332048771345315401/53375584364747317 6140678428961747763200000000*cos(1)^2+2541573211/146767085568*sin(1)^24+940325057 70279736611460220749/16522396770089041920000000000*sin(1)^20+2704734082846637530 0822998906838403821858600396601/4670363631915390291230936253415292928000000000*sin(1)^8+5542192477209543230894137867604219553503/255755814666761099367961067520 00000000*sin(1)^16+6290548805350754451916704658025155325352197570086424587493/20 18643717477360844776029314289796912041164800000000000*sin(1)^9*cos(1)-74888453896 484988301898479573506809/1756135812378578105401344000000000*sin(1)^18-66275838868 9551809679364987359346538465101954279/842854686697230591620583026983541145600000000*sin(1)^14+4741716006420769418428944378170996543275038000761225969993/153801 4260935132072210308048982702409174220800000000000*sin(1)^5*cos(1)-273698005143037 11474211657412466731397670944577/100266509581957021396215456910540800000000000*sin(1)^15*cos(1)+18044178399358284974551495/22109663333532016620601344*sin(1)^21*cos(1)+159018588498544047612814017616772807534595903672788230889/18504234076875807 7437802687143231383603773440000000000*sin(1)^13*cos(1)-42990929319556261053136947 04999192498272936575705473777243/2220508089225096929253632245718776603245281280 000000000*sin(1)^11*cos(1)+72644795857216400572361363249893595045872001672201/160 1267530942419528422035286885243289600000000*sin(1)^4-179168559345113148705406926 00563709435304263998599/3002376620517036615791316162909831168000000000*sin(1)^6-3 228431702614231399553/6979060395685611307008*sin(1)^22+1809464903223467961506769 2014871302547716716561173/8669362491742943228097425420402137497600000000*sin(1)^12+908212034006674534482628671295497731536658010905190650989336891/1066356554248 3582367324802472946736703607930880000000000解（1）int(x/(x^2+x+1),-1,1)ans =1/2*log(3)-1/6*3^(1/2)*pi 解（2）> int((x*sin(x))^2,0,pi)ans =1/6*pi^3-1/4*pi 解（3）第七题求级数之和n=sym('n');>> s1=symsum((-1)*(2*n+1)/2^n,n,0,inf)s1 =-6 解（1）>> s2=symsum(x^(2*n-1)/2^n-1,n,1,inf)s2 =sum(x^(2*n-1)/(2^n)-1,n = 1 .. Inf) 解（2）s3=symsum(1/(2*n+1)^2,n,0,inf)s3 =1/8*pi^2 解（3）s4=symsum(1/n*(n+1)*(n+1),n,1,inf)s4 =Inf 解（4）第⼋题求泰勒展开式>> x=sym('x');f1=x^4-5*x^3+x^2-3*x+4;f2=(exp(x)+exp(-x))/2;f3=tan(x);f4=sin(x)^2;f5=sqrt(x^3+x^2+5*x+3);taylor(f1,4,4)ans =-140+21*x+37*(x-4)^2+11*(x-4)^3 解（1）taylor(f2,5,0)ans =1+1/2*x^2+1/24*x^4 解（2）taylor(f3,3,2)ans =tan(2)+(1+tan(2)^2)*(x-2)+tan(2)*(1+tan(2)^2)*(x-2)^2 解（3）taylor(f4,8,0)ans =x^2-1/3*x^4+2/45*x^6 解（4）taylor(f5,5,0)ans =3^(1/2)+5/6*3^(1/2)*x-13/72*3^(1/2)*x^2+137/432*3^(1/2)*x^3-2909/10368*3^(1/2)* x^4 解（5）第九题求⾮线性⽅程的解x=solve(‘a*x^2+b*x+c=0’,’x’)x =1/2/a*(-b+(b^2-4*a*c)^(1/2))1/2/a*(-b-(b^2-4*a*c)^(1/2)) 解（1）x=solve(‘2*sin(3*x-pi/4)=1’,’x’)x =5/36*pi 解（2）x=solve(‘sin(x)-sqrt(3)*cos(x)=sqrt(2)’,’x’)x =-atan(2*(1/4*2^(1/2)+1/4*3^(1/2)*2^(1/2))*2^(1/2)/(3^(1/2)-1))+pi-atan(2*(1/4*2^(1/2)-1/4*3^(1/2)*2^(1/2))*2^(1/2)/(1+3^(1/2)))-pi 解（3）x=solve(‘x^2+10*(x-1)*sqrt(x)+14*x+1=0’,’x’)x =(2^(1/2)-1)^2(-4+17^(1/2))^2 解（4）第⼗题求⽅程组的解[x,y]=solve(‘ln(x/y)=9’,’exp(x+y)=3’,’x,y’)x =exp(9)*log(3)/(exp(9)+1)y =log(3)/(exp(9)+1) 解（1）[x,y,z]=solve(‘(4*x^2)/(4*x^2+1)=y’,’(4*y^2)/(4*y^2+1)=z’,’(4*z^2)/(4*z^2+ 1)=x’,’x,y,z’) x =y =0 解（2）z =第⼗⼀题求初值y=dsolve('x*(D2y)+(1-n)*(Dy)+y=0','y(0)=Dy(0)=0','x')第⼗⼆题，求特解[x,y]=dsolve(‘Dx=3*x+4*y’,’Dy=5*x-7*y’,’x(0)=0’,’y(0)=0’,’t’)x =y =解（1）。
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课后习题答案
习题二之阿布丰王创作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。
(完整版)MATLAB)课后实验答案[1]
1 + e2 (2) z = 1 ln( x + 1 + x 2 ) ，其中 x = ⎡⎢ 2⎣-0.45 ⎦2 2 ⎪t 2 - 2t + 1 2 ≤ t <3 ⎨实验一MATLAB 运算基础1. 先求下列表达式的值，然后显示 MATLAB 工作空间的使用情况并保存全部变量。
(1) z = 2sin 8501221 + 2i ⎤5 ⎥(3) z = e 0.3a - e -0.3asin(a + 0.3) + ln 0.3 + a ,a = -3.0, - 2.9, L , 2.9, 3.03⎧t 2 0 ≤ t < 1 (4) z = ⎪t 2 - 11 ≤ t <2 ，其中 t=0:0.5:2.5 4⎩解：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;3=(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)=[]ch=⎣O2⨯3⎥，其中E、R、O、S分别为单位矩阵、随机矩阵、零矩S⎦阵和对角阵，试通过数值计算验证A=⎢⎥。
matlab习题二答案
matlab习题二答案Matlab习题二答案Matlab是一种强大的数值计算和科学计算软件，广泛应用于工程、科学、金融等领域。
在学习和使用Matlab的过程中，习题是一种很好的练习和巩固知识的方式。
本文将为大家提供一些Matlab习题二的答案，希望能对大家的学习有所帮助。
1. 编写一个函数，输入一个矩阵A，输出矩阵A的转置矩阵。
```matlabfunction B = transposeMatrix(A)B = A';end```2. 编写一个函数，输入一个矩阵A，输出矩阵A的每一行的平均值。
```matlabfunction avg = rowAverage(A)avg = mean(A, 2);end```3. 编写一个函数，输入一个矩阵A和一个标量k，输出矩阵A中大于k的元素个数。
```matlabfunction count = countGreaterThanK(A, k)count = sum(A(:) > k);end```4. 编写一个函数，输入一个矩阵A和一个标量k，输出矩阵A中大于k的元素的索引。
```matlabfunction indices = findGreaterThanK(A, k)indices = find(A > k);end```5. 编写一个函数，输入一个矩阵A，输出矩阵A的每一列的方差。
```matlabfunction variances = columnVariance(A)variances = var(A);end```6. 编写一个函数，输入一个矩阵A和一个标量k，将矩阵A中小于k的元素替换为0。
```matlabfunction B = replaceLessThanK(A, k)B = A;B(A < k) = 0;end```7. 编写一个函数，输入一个矩阵A和一个标量k，将矩阵A中小于k的元素替换为k。
```matlabfunction B = replaceLessThanK(A, k)B = A;B(A < k) = k;end```8. 编写一个函数，输入一个矩阵A和一个标量k，将矩阵A中大于k的元素替换为k。
MATLAB)课后实验标准答案[1]
实验一 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)=[]ch =实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。