MATLAB课后习题

合集下载

Matlab课后习题解答

0.70
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

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文件中.8、如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成.在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上.9、在MA TLAB中有几种获得帮助的途径？答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器.（2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息.（3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数.（4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数.注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息.第3章 MATLAB 数值运算3.1在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？答：在Command Window 窗口输入操作:>> a=[5 7 3;4 9 1]3.2有几种建立矩阵的方法？各有什么优点？答：（1）直接输入法，如a=[1 2 3 4]，优点是输入方法方便简捷；（2）通过M 文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改；（3）由函数建立，如y=sin(x)，可以由MATLAB 的内部函数建立一些特殊矩阵；（4）通过数据文件建立，该方法可以调用由其他软件产生数据.3.3在进行算术运算时，数组运算和矩阵运算各有什么要求？答：进行数组运算的两个数组必须有相同的尺寸.进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数.渗釤呛俨匀谔鱉调硯錦。

MATLAB-实用教程-课后习题标准答案

第二章1.计算复数3+4i与5-6i的乘积。

a=3+4ib=5-6ic=a*b2•构建结构体Students» 属性包含Name, age 和Email,数据包括｛"Zhang", 18,"]｝、「Wang，, 21, []｝和｛'『，[]』｝,构建后读取所有Name属性值,并且修改'Zhang，的Age 属性值为19.Students(l).Age=18Students(l).Email-1/'Students(2).Name=,Wang,Students(2).Age=21Students(2).Email=[]Students(3).Name=,Li,Students(3).Age=[]Students(3).Email=[]Student(l).Age(l)=19Stude nt.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,l,4]/[1,2/4L[1,1,1]A5)4.采用向量构造符得到向量[1,5,941],A=l:4:415.按水平和竖直方向分别合并下述两个矩阵：A=[l 0 0;11 0;0 0 1LB=[2 3 4;5 6 7;8 9 10] A=[10 0;110;0 0 1]B=[2 3 4;5 6 7;8 9 10]C=[A B]D=[A;B]6.分别删除第五题两个结果的第2行。

A=[10 0;110;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列的数值改为12 13]oA=[10 0;110;0 0 1]D=[A;B]C(2/4:6)=[ll 12 13]D(2/:)=[ll 12 13]8.分別查看第5题两个结果的各方向长度A=[10 0;110;0 0 1]B=[2 3 4;5 6 7;8 9 10]C=[A B]D=[A;B]a=size(C)b=size(D)9.分別判断pi是否为字符串和浮点数。

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 )的图形。

(完整word版)matlab经典习题及解答

第1章 MATLAB概论1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点.1.2 MATLAB系统由那些部分组成?MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。

1.4 MATLAB操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立窗口?又如何将脱离出去的窗口重新放置到桌面上?在MATLAB操作桌面上有五个窗口,在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口,在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。

1.5 如何启动M文件编辑/调试器？在操作桌面上选择“建立新文件”或“打开文件"操作时，M文件编辑/调试器将被启动.在命令窗口中键入edit命令时也可以启动M文件编辑/调试器.1.6 存储在工作空间中的数组能编辑吗？如何操作？存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。

1。

7 命令历史窗口除了可以观察前面键入的命令外,还有什么用途?命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。

1。

8 如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件有什么区别？当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的Set Path 菜单项来完成。

在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。

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课后习题与答案

习题二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课后习题

5、利用rand函数产生（0,1）间的均匀分布的10*10随机矩阵A，然后统计A中大于等于的元素的个数。

解：A=rand(10);B=A >= ;C=sum(B);count=sum(C)运行结果（每次运行结果是不同的，仅作参考）：count=326、利用randn函数产生均值为0，方差为1的10*10随机矩阵A，然后统计A中大于且小于的元素的个数。

%解：A=randn(10);B=(A<&(A>;C=sum(sum(B))运行结果（每次运行结果是不同的，仅作参考）：C=481、解：if and(a<1,b<=语句1;elseif and(a<1,b>`语句2;elseif and(a>=1,b<=语句3;else语句4;2、有一矩阵A，找出矩阵中值等于1的元素，并将它们重新排列成列向量B。

解：A=2*rand(4);k=find(A<=1);￥A(k)=[];%删除下标为k的元素B=A'运行结果（每次运行结果是不同的，仅作参考）B =)3、在一测量矩阵A（100*3）中，存在有奇异值（假设大于100的置认为是奇异值），编程实现删去奇异值所在的行。

解：A=120*randn(10,3);[i,j]=find(A>100);A(i,:)=[] %删去存在奇异值的行【4、在给定的100*100矩阵中，删去整行为0的行，删去整列为0的列。

解：A=diag([1 2 3 4],1)B=any(A)[i,j]=find(B==0)A(:,i)=[] %删除全为0的列B=any(A')[i,j]=find(B==0)A(j,:)=[] %删除全为0的行—运行结果：初始值：A =0 1 0 0 00 0 2 0 00 0 0 3 00 0 0 0 40 0 0 0 0操作后：A =1 0 0 0；0 2 0 00 0 3 00 0 0 41、将窗口分割成四格，分别绘制正弦、余弦、正切和余切函数曲线，并加上适当的标注。

完整word版,Matlab课后习题

习题 11. 执行下列指令，观察其运算结果, 理解其意义： (1) [1 2;3 4]+10-2i(2) [1 2; 3 4].*[0.1 0.2; 0.3 0.4] (3) [1 2; 3 4].\[20 10;9 2] (4) [1 2; 3 4].^2 (5) exp([1 2; 3 4]) (6)log([1 10 100]) (7)prod([1 2;3 4])(8)[a,b]=min([10 20;30 40]) (9)abs([1 2;3 4]-pi)(10) [1 2;3 4]>=[4,3;2 1](11)find([10 20;30 40]>=[40,30;20 10])(12) [a,b]=find([10 20;30 40]>=[40,30;20 10]) (提示：a 为行号，b 为列号) (13) all([1 2;3 4]>1) (14) any([1 2;3 4]>1) (15) linspace(3,4,5) (16) A=[1 2;3 4];A(:,2)2. 执行下列指令，观察其运算结果、变量类型和字节数，理解其意义： (1) clear; a=1,b=num2str(a),c=a>0, a= =b, a= =c, b= =c (2) clear; fun='abs(x)',x=-2,eval(fun),double(fun)3. 本金K 以每年n 次，每次p %的增值率(n 与p 的乘积为每年增值额的百分比)增加，当增加到rK 时所花费的时间为)01.01ln(ln p n rT +=(单位：年)用MA TLAB 表达式写出该公式并用下列数据计算：r =2, p =0.5, n =12.4．已知函数f (x )=x 4-2x 在(-2, 2)内有两个根。

取步长h =0.05, 通过计算函数值求得函数的最小值点和两个根的近似解。

matlab第二版课后习题答案

matlab第二版课后习题答案
《MATLAB第二版课后习题答案》
MATLAB是一种强大的数学软件，被广泛应用于工程、科学和金融等领域。

《MATLAB第二版》是一本经典的教材，为了帮助学生更好地掌握MATLAB的使用，书中提供了大量的课后习题。

下面我们将为大家总结一些MATLAB第二版课后习题的答案，希望能对大家的学习有所帮助。

1. 第一章课后习题答案
第一章主要介绍了MATLAB的基本操作，包括变量的定义、矩阵的运算、函数的使用等。

在课后习题中，有一道题目是要求计算一个矩阵的逆矩阵。

答案是使用MATLAB中的inv函数，将原矩阵作为参数传入即可得到逆矩阵。

2. 第二章课后习题答案
第二章介绍了MATLAB中的绘图功能，包括二维和三维图形的绘制。

有一道课后习题是要求绘制一个正弦曲线和余弦曲线，并在同一张图上显示。

答案是使用MATLAB中的plot函数，分别绘制正弦曲线和余弦曲线，并使用legend函数添加图例。

3. 第三章课后习题答案
第三章介绍了MATLAB中的控制流程，包括if语句、for循环和while循环等。

有一道课后习题是要求编写一个程序，计算1到100之间所有偶数的和。

答案是使用for循环遍历1到100之间的所有数，判断是否为偶数并累加。

通过以上几个例子，我们可以看到MATLAB第二版课后习题的答案涵盖了各种基本和高级的操作，对于学习MATLAB是非常有帮助的。

希望大家在学习MATLAB的过程中能够多加练习，掌握更多的技巧和方法。

MATLAB数学实验第二版课后练习题含答案

MATLAB数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下：第一章课后练习题1.编写MATLAB程序，计算并输出下列公式的结果：y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中，x, $\\mu$, $\\sigma$ 分别由用户输入。

要求输出结果精确至小数点后两位。

答案如下：x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序，求解下列方程的解：4x + y = 11\\\\x + 2y = 7答案如下：A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序，计算下列多项式的值：P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中，x 由用户输入。

要求输出结果精确至小数点后两位。

答案如下：x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序，绘制下列函数的图像：f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下：x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。

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课后习题

第一部分 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.52. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A(5) [A,B]和[A([1,3],:);B^2]3. 设有矩阵A 和B123453166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C 。

(2) 将矩阵C 的右下角3×2子矩阵赋给D 。

(3) 查看MATLAB 工作空间的使用情况。

4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。

(2) 建立一个字符串向量，删除其中的大写字母。

第二部分 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。

2. 产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。

MATLAB课后习题解答

信号与系统MATLAB平时作业第二章（1）代码：t=-4:0.001:4;>> T=2;>> xt=rectpuls(t-1,T);>> plot(t,xt);>> axis([-4,4,-0.5,1.5]) 图形：（2）代码：t=sym('t');>> y=heaviside(t); >> ezplot(y,[-1,1]); >> grid on;>> axis([-1,1,-0.2,1.2]) 图形：(3)代码：> A=10;a=-1;B=5;b=-2;>> t=0:0.001:10;>> xt=A*exp(a*t)-B*exp(b*t); >> plot(t,xt);图形：（4）代码：t=sym('t');y=t*heaviside(t);ezplot(y,[-1,3]);grid on ;axis([-1 3 -0.1 3.1])图形：（5）代码：xt=abs(A*sin(w0*t+phi)); >> plot(t,xt);>> A=2;w0=10*pi;phi=pi/3;图形：（6）代码：> A=1;w0=1;B=1;w1=2*pi;>> t=0:0.001:20;>> xt=A*cos(w0*t)+B*sin(w1*t); >> plot(t,xt);图形：（7）代码：A=4;a=-0.5;w0=2*pi;>> t=0:0.001:10;>> xt=A*exp(a*t).*cos(w0*t); >> plot(t,xt);图形：（8）代码：w0=30;t=-15:0.001:15;xt=sinc(t/pi).*cos(w0*t); plot(t,xt);图形：M2-2（1）xt程序：t=-1:0.001:1;x=4*t.*[t>=0]-4*t.*[t>=1/2]-4*t.*[t>=1/2]+4*[t>=1/2]+4*t.*[t>=1]-4*[t>=1 ]+[t>=-1]-[t>=0];>> plot(t,x);图形：(2)程序：t=-1:0.001:1;x=(4*t.*[t>=0]-4*t.*[t>=1/2]-4*t.*[t>=1/2]+4*[t>=1/2]+4*t.*[t>=1]-4*[t>= 1]+[t>=-1]-[t>=0]).*cos(50*t);>> plot(t,x);图形：M2-3（1）yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5);（2）t=0:0.001:6;yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5);subplot(3,1,1)plot(t,yt)title('x(t)')axis([0,6,-2,2])t=0:0.001:12;yt=((0.5)*t).*(t>=0&t<=4)+2*(t>=4&t<=6)-1*(t>=6&t<=10); subplot(3,1,2)plot(t,yt);axis([0,12,-2,2]);title('(0.5)*t');t=-6:0.001:6;yt=((2-(0.5)*t)).*(t>=0&t<=4)+2*(t>=-2&t<=0)-1*(t>=-6&t<=-2); subplot(3,1,3)plot(t,yt);axis([-6,6,-2,2]);title('2-(0.5)*t');图形：M2-4奇分量程序：t=-2:0.001:2;x1t=(t.*[t>=-1]-t.*[t>=0]+[t>=-1]-[t>=0]+t.*[t<=1]-t.*[t<=0]-[t<=1]+[t<=0]) /2;>> plot(t,x1t);图形：偶分量：程序：t=-2:0.001:2;>>x2t=(t.*[t>=-1]-t.*[t>=0]+[t>=-1]-[t>=0]-t.*[t<=1]+t.*[t<=0]+[t<=1]-[t<=0]) /2;>> plot(t,x2t);图形：M2-5 （1）function [x,k]=impseq(0,-50,50) k=[-50:50];x=[k==0];stem(k,x)图形：（2）k=[-50:50];>> uk=[zeros(1,50),ones(1,51)];>> stem(k,uk)图形：（3）function [f,k]=impseq(0,-50,50)；k=[-50:50];x=10*(1/2).^k.*[k>=0];>> stem(k,x);图形：（4）function [f,k]=impseq(k0,-50,50) k=[-10:10];x1=[(k+2)>=0];x2=[(k-6)>=0];x=x1-x2;stem(k,x);图形：（5）function [f,k]=impseq(k0,-50,50) k=[-50:50];>> x=k.*[k>=0];>> stem(k,x);图形：（6）k=-50:50;xk=(5*0.8.^k).*cos(0.9*pi*k);>> stem(k,xk);图形：M2-6 Ω0=0.1*pi;k=-50:50;>> xk=sin(0.1*pi*k);>> stem(k,xk);图形：Ω0=0.5*pi;k=-50:50;xk=sin(0.5*pi*k);stem(k,xk);图形：Ω0=0.9*pi;k=-50:50;xk=sin(0.9*pi*k); stem(k,xk);图形：Ω0=1.1*pi;k=-50:50;xk=sin(1.1*pi*k); stem(k,xk);图形：Ω0=1.5*pi;k=-50:50;xk=sin(1.5*pi*k); stem(k,xk);图形：Ω0=1.9*pi;k=-50:50;xk=sin(1.9*pi*k); stem(k,xk);图形：结论：随着Ω0的变动，波形也随之而成周期性变换M2-7 （1）程序：t=-2:0.001:2;x1t=cos(6*pi*t);plot(t,x1t);hold onk=-2:2;x1k=cos(0.6*pi*k);stem(k,x1k,'r');hold off图形：（2）程序：t=-1:0.001:1;x1t=cos(14*pi*t); plot(t,x1t);hold onk=-1:1;x1k=cos(1.4*pi*k); stem(k,x1k,'r'); hold off图形：（3）程序：t=-0.5:0.001:0.5;x1t=cos(26*pi*t);plot(t,x1t);hold onk=-0.5:0.5;x1k=cos(2.6*pi*k);stem(k,x1k,'r');hold off图形：M2-8 （1）用square函数：程序：k=0:40;A=1;P=0.6;y=A*square(P*k);stem(k,y)axis([k(1)-1 k(end)+1 -(A+1) (A+1)]); 图形：（2）程序：k=0:40;x2k=sawtooth(0.1*pi*k,0.5);stem(k,x2k);图形：M2-9(1)k=-4:7;xk=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1];stem(k,xk)图形：(2)k=-12:21x=[-3,0,0,-2,0,0,3,0,0,1,0,0,-2,0,0,-3,0,0,-4,0,0,2,0,0,-1,0,0,4,0,0,1,0,0,-1]; subplot(2,1,1)stem(k,x)title('3倍内插)t=-1:2;y=[-2,-2,2,1];subplot(2,1,2)stem(t,y)title('3倍抽取’)axis([-3,4,-4,4])图形：（3）k=-4:7;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; subplot(2,1,1)stem(k+2,x)title('x[k+2]')subplot(2,1,2)stem(k-4,x)title('x[k-4]')图形：（4）k=-4:7;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; stem(-fliplr(k),fliplr(x))title('x[-k]')图形：第三章M3-1(1)代码：ts=0;te=5;dt=0.01;>> sys=tf([2 1],[1 3 2]);>> t=ts:dt:te;>> x=exp(-3*t);>> y=lsim(sys,x,t);>> plot(t,y);>> xlabel('Time(sec)');>> ylabel('y(t)');图形：（2）代码：ts=0;te=5;dt=1; sys=tf([2 1],[1 3 2]); t=ts:dt:te;x=exp(-3*t);y=lsim(sys,x,t); >> plot(t,y); xlabel('Time(sec)'); ylabel('y(t)');图形：M3-2 代码：ts=0;te=5;dt=0.01;>> sys=tf([6],[1 3 6]);>> t=ts:dt:te;>> y=impulse(sys,t);>> plot(t,y);>> xlabel('Time(sec)');>> ylabel('h(t)');图形：（3）代码：ts=0;te=5;dt=0.01; sys=tf([6],[1 3 6]); t=ts:dt:te;y=step(sys,t);plot(t,y);xlabel('Time(sec)'); ylabel('u(t)')图形：M3-3 (1)代码：ts=0;te=5;dt=1;sys=tf([1],[1 0.2 1]);t=ts:dt:te;y=step(sys,t);plot(t,y);xlabel('Time(sec)');ylabel('u(t)')图形：（2）代码：ts=0;te=5;dt=1; sys=tf([1],[1 1 1]); t=ts:dt:te;y=step(sys,t);plot(t,y);xlabel('Time(sec)'); ylabel('u(t)')图形：（3）代码：ts=0;te=5;dt=1;sys=tf([1],[1 2 1]);t=ts:dt:te;y=step(sys,t);plot(t,y);xlabel('Time(sec)');ylabel('u(t)')图形：M3-4 代码：x=[0.85,0.53,0.21,0.67,0.84,0.12];kx=-2:3;h=[0.68,0.37,0.83,0.52,0.71];kh=-1:3;y=conv(x,h);k=kx(1)+kh(1):kx(end)+kh(end);stem(k,y,'r');图形：M3-6h(t)=x(t)x(t)=u(t)*u(t)+u(t-1)*u(t-1)+2u(t)*u(t-1)=r(t)-2r(t-1)+r(t-2)y(t)=x(t) h(t)=1/2*t^2.*u(t)-3/2*(t-1)^2.*u(t-1)-1/2*(t-3)^2.*u(t-3)0, t<0t^2/2, 0<=t<1 Y(t) = -t^2+3*t-3/2 1<=t<2( t-3)^2/2 2<=t<30 t>=3(b)T=0.1;k=-1:T:4;f1=1*((k>=0)&(k<=1));f2=tripuls(k-1,2);y=T*conv(f1,f2);tmin=-2;tmax=8;t1=tmin:0.1:tmax;plot(t1,y)grid on图形：T=0.1：图形：T=0.01T=0.001：M3-7 代码：b=[1 0.42 -0.19];a=[0.31 0.68];k=0:20;x=(0.6).^k.*[k>=0];y=filter(b,a,x);stem(k,y);图形：M3-8 代码：b=[1 0.7 -0.45 -0.6];a=[0.8 -0.44 0.36 0.02];k=0:30;h=impz(b,a,k)stem(k,h);h =1.25001.5625-0.2656-1.6305-0.81630.29140.56840.2019-0.1520-0.1887-0.04040.06650.05950.0038-0.0263-0.01770.00200.00970.0049-0.0017-0.0034-0.00120.00090.00110.0002-0.0004-0.0004-0.00000.00020.0001-0.0000 图形：第四章M4-1周期矩形信号：代码：n=-20:20;X=-j*1/2*sin(n/2*pi).*sinc(n/2);subplot(2,1,1);stem(n,abs(X));title('幅度谱')xlabel('nw');subplot(2,1,2);stem(n,angle(X));title('相位谱')图形：三角波信号：代码：>> n=-20:20;X=sinc(n)-0.5*((sinc(n/2)).^2); subplot(2,1,1);stem(n,abs(X));title('幅度谱')xlabel('nw');subplot(2,1,2);stem(n,angle(X));title('相位谱')图形：M4-6 (4)代码：x=[1,2,3,0,0];X=fft(x,5);subplot(2,1,1);m=0:4;stem(m,real(X)); title('X[m]实部') subplot(2,1,2); >> stem(m,imag(X)); title('X[m]虚部');图形：M4-7 (3)代码：k=0:10;x=0.5.^k;subplot(3,1,1);stem(k,x)title('x[k]')X=fft(x,10);subplot(3,1,2);m=0:9;stem(m,real(X));title('X[m]实部')>> subplot(3,1,3);>> stem(m,imag(X));title('X[m]虚部')图形：M5-2代码：t=0:0.05:2.5;T=1;xt1=rectpuls(t-0.5,T);subplot(2,2,1)plot(t,xt1)title('x(t1)')axis([0,2.5,0,2])xt2=tripuls(t-1,2);subplot(2,2,2)plot(t,xt2)title('x(t2)')axis([0,2.5,0,2])xt=xt1+xt2.*cos(50*t);subplot(2,2,[3,4])plot(t,xt)title('x(t)')figure;b=[10000];a=[1,26.131,341.42,2613.1,10000];[H,w]=freqs(b,a,w);subplot(2,1,1)plot(w,abs(H));set(gca,'xtick',[0 1 2 3 4 5 6 7 8]);set(gca,'ytick',[0.985 0.99 0.995 1 1.005]);grid; title('幅度曲线')subplot(2,1,2)plot(w,angle(H));set(gca,'xtick',[0 1 2 3 4 5 6 7 8])set(gca,'ytick',[-2 -1.5 -1 -0.5 0]);grid;title('相位曲线')figure;sys=tf([10000],[1 26.131 341.42 2613.1 10000]); yt1=lsim(sys,xt,t);subplot(2,1,1);plot(t,yt1);title('y(t1)')yt2=lsim(sys,xt.*cos(50*t),t);subplot(2,1,2);plot(t,yt2);title('y(t2)');图形：第六章 M6-1已知连续时间信号的s 域表示式如下，使用residue 求出X(s)的部分分式展开式，并写出x(t)的实数形式表达式。

matlab仿真课后习题

第一章习题3.请指出以下的变量名（函数名、M文件名）中，哪些是合法的Abc 2004x lil-1 wu_2004 a&b _xyz 解：合法的变量名有：Abc wu_20044．指令窗操作（1）求[12+2×（7-4）]÷32的运算结果解：>> [12+2*(7-4)]/3^2ans =2（2）输入矩阵A=[1，2，3；4，5，6；7，8，9]，观察输出。

解：>> A=[1,2,3;4,5,6;7,8,9]A =1 2 34 5 67 8 9（3）输入以下指令，观察运算结果；clear;x=-8::8;y=x';X=ones(size(y))*x;Y=y*ones(size(x));R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;mesh(X,Y,Z);colormap(hot)xlabel('x'),ylabel('y'),zlabel('z')解：7．指令行编辑（1）依次键入以下字符并运行：y1=2*sin*pi)/(1+sqrt(5))解：>>y1=2*sin*pi)/(1+sqrt(5))y1 =(2)通过反复按键盘的箭头键，实现指令回调和编辑，进行新的计算；y2=2*cos*pi)/(1+sqrt(5))解：>>y2=2*cos*pi)/(1+sqrt(5))y2 =11.编写题4中（3）的M脚本文件，并运行之。

解：第二章习题1.在指令窗中键入x=1::2和y=2::1,观察所生成的数组。

解:>> x=1::2 x =>> y=2::1 y =Empty matrix: 1-by-02．要求在[0，2π]上产生50个等距采样数据的一维数组，试用两种不同的指令实现。

解： y1=0:2*pi/49:2*pi y2=linspace(0,2*pi,50)3.计算e -2t sint ，其中t 为[0，2π]上生成的10个等距采样的数组。

matlab课后习题答案1到6章

欢迎共阅习题二1.如何理解“矩阵是MATLAB 最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。

(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’); 或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⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，11⎡⎤⎢⎥答：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)的值分别是(2)建立5×6随机矩阵A，其元素为[100,200]范围内的随机整数。

Matlab 课后习题集

Matlab程序设计课程实验目录Matlab操作实验（第1章）实验指导思考练习矩阵运算实验（第2章）实验指导思考练习程序设计实验（第3章）实验指导思考练习高级绘图实验（第4章）实验指导思考练习数据分析与多项式计算实验（第5章）实验指导思考练习解方程与最优化问题求解实验（第6章）实验指导思考练习数值运算实验（第7章）实验指导思考练习符号实验（第8章）实验指导思考练习图形句柄实验（第9章）实验指导思考练习图形用户界面设计实验（第10章）实验指导思考练习Simulink动态仿真实验（第11章）实验指导思考练习外部接口技术实验（第12章）实验指导思考练习1、Matlab操作实验一、实验内容1. 第一章实验内容1：path(path,’系统盘名:\matlab文件名\要添加的文件名’)如：path(path,’c:\matlab701\my work’)实验内容-2例1-1：x=[0:0.5:360]*pi/180;plot(x,sin(x),x,cos(x));例1-2p=[3,7,9,0,-23];x=roots(p)例1-3quad(‘x.*log(1+x)’,0,1)例1-4：实验内容3Help 函数名如：help inv help plot等实验内容4（1）x=0:pi/10:2*pi;y=sin(x);二、思考练习思考题1启动：有三种常见方法，（1）、在Windows桌面，单机任务栏上的‚开始‛按钮，选择‚所有程序‛菜单项，然后选择MATLAB程序组中的MATLAB （版本号）程序选项，即可启动。

（2）、在MATLAB的安装路径中找到MATLAB系统启动程序matlab.exe，然后运行它。

（3）、利用建立快捷方式功能，将MATLAB系统启动程序以快捷方式的形式放在Windows桌面上，以后只要在桌面上双击该图标即可启动MATLAB。

关闭：常见的也有3种方法，（1）、在MATLAB主窗口file菜单中选择exit MATLAB命令。

MATLAB课后习题答案

clear allclcsymvar(sym('sin(w*t)'),1)symvar(sym('a*exp(-X)'),1)symvar(sym('z*exp(j*th)'),1)Ch2_3(1)clcclear allreset(symengine)syms x positivesolve(x^3-44.5,x)Ch2_3(2)clcsyms x clearsyms a positivesolve(x^2-a*x+a^2,x)Ch2_4(1)clcx=7/3a=xb=sym(x)c=sym(x,'d')d=sym('7/3')v1=vpa(abs(a-d))v2=vpa(abs(b-d))v3=vpa(abs(c-d))Ch2_5clcA=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33;]') DA=det(A)IA=inv(A);%[IAs,d]=subexpr(IA,d)subexper(IA)Ch2_6clcsyms x kf=x^(k);Z1=symsum(f,k,0,inf)subs(Z1,x,{sym('-1/3'),sym('1/pi'),sym('3^(k)'Ch2_7clcsyms k;x=sym('x','positive');f_k=2/(2*k+1)*((x-1)/(x+1))^(2*k+1);s=simple(symsum(f_k,k,0,inf))Ch2_8clcsyms ty=abs(sin(t))d=diff(y)d0_=limit(d,t,0,'left')dpi_2=limit(d,t,pi/2)Ch2_9(1)clcsyms x clear%syms xy=exp(-abs(x))*abs(sin(x))si=vpa(int(y,-5*pi,1.7*pi))%si=vpa(int(y,-5*pi,1.7*pi),64)Ch2_9(2)xx=-5*pi:pi/100:1.7*pi;sn=trapz(exp(-abs(xx)).*abs(sin(xx)))*pi/100 Ch2_10clcsyms x yf=x^2+y^2;r=int(int(f,y,1,x^2),x,1,2)Ch2_11(1)clcsyms t x;f=sin(t)/t;y=int(f,t,0,x)y5=subs(y,x,sym('4.5'))ezplot(y,[0,2*pi])gridCh2_11(2)clc,clf%tt=0:0.001:4.5;tt=eps:0.001:4.5;yn=trapz(sin(tt)./tt)*0.001Ch2_12clcsyms xsyms n positivef=sin(x)^nyn=int(f,x,0,pi/2)y3s=vpa(subs(yn,n,sym('1/3'))) y3d=vpa(subs(yn,n,1/3))Ch3_1clear,clc方法一：t1=linspace(0,2*pi,10)方法二：t2=0:2*pi/9:2*piCh3_2clear,clcrng('default')A=rand(3,5)[ri,cj]=find(A>0.5);ri=ri'cj=cj'id=find(A>0.5);id=id'Ch3_3方法一：clear,clcrng defaultA=rand(1,1000);a=2*(A>0.5)-1;Na=sum(a==-1)方法二：clcrandn('state',0)B=randn(1,1000);b=2*(B>0)-1;Nb=sum(b==-1)方法三：clcc=randsrc(1,1000,[-1,1]);Nc=sum(c==-1)Ch3_4clear,clcA=[1,2;3,4]B1=A.^0.5B2=A^0.5A1=B1.*B1;A2=B2*B2;norm(A1-A2,'fro')epsCh3_5clear,clf,clct=linspace(0,10,200);N=length(t);y1=zeros(size(t));for k=1:Ny1(k)=1-exp(-0.5*t(k))*cos(2*t(k));endsubplot(1,2,1),plot(t,y1),xlabel('t'),ylabel(' y1')y2=1-exp(-0.5*t).*cos(2*t);subplot(1,2,2),plot(t,y2),xlabel('t'),ylabel(' y2')Ch3_6clear,clcformat longrng('default')A=rand(3,3)B=diag(diag(A))C=A-BCh3_7clear,clc,clfx=-3*pi:pi/15:3*pi;y=x;[X,Y]=meshgrid(x,y);warning offZ=sin(X).*sin(Y)./X./Y;Num_of_NaN=sum(sum(isnan(Z)))subplot(1,2,1),surf(X,Y,Z),shadinginterp,title('fname')Ch3_8clc,clear,for k=4:-1:1A=reshape(1:4*k,k,4);Sa(k,:)=sum(A);if k==1Sa(k,:)=A;endendSaCh4_1clear,clcload prob_data401.matdt=t(2)-t(1);yc=diff(y)/dt;plot(t,y,'b',t(1:end-1),yc,'r')hold onCh4_2clcdt=1e-4;t=0:dt:10;t=t+(t==0)*eps;f=sin(t)./t;s=cumtrapz(f)*dtplot(t,s,'linewidth',3),hold on,plot(t,f) ii=find(t==4.5);s45=s(ii)gridCh4_3clear,clc,format long edx=pi/2000;x=0:dx:pi;s_trapz=trapz(exp(sin(x).^3))*dxs_sum=sum(exp(sin(x).^3))*dxs_quad=quad('exp(sin(x).^3)',0,pi,1e-9) syms xf=exp(sin(x)^3);ss=int(f,x,0,pi)Ch4_4clear,clc,clfy=@(x)exp(-abs(x)).*abs(sin(x));format longs_quadl=quadl(y,-10*pi,1.7*pi,1e-9)x=linspace(-10*pi,1.7*pi,1000);plot(x,y(x)),gridCh4_5clear,clc,clff=@(t)sin(5*t).^2.*exp(0.06*t.*t)+1.8*abs(t+0.5)-1.5*t.*cos(2*t)t=-5:0.001:5;plot(t,f(t))grid on,shgclear,clcx1=-1.5;x2=-1;yx=@(t)sin(5*t).^2.*exp(0.06*t.*t)+1.8*abs(t+0 .5)-1.5*t.*cos(2*t)[xmin,fmin]=fminbnd(yx,x1,x2)Ch4_12clear,clct=-1:0.001:5;y=@(t)-0.5+t-10*exp(-0.2*t).*abs(sin(sin(t))); plot(t,y(t))grid on,shg[tt1,yy1]=ginput(1)[t,yt]=fzero(y,tt1)Ch4_18clear,clc,clfload prob_data418.matP=polyfit(x,y,4);xx=linspace(x(1),x(end),1*length(x));yy=polyyval(P,xx);PP=poly2str(P,'x')plot(x,y,'r.',xx,yy,'ob')。

matlab课后练习习题及答案详解

matlab课后习题及答案详解第1章MATLAB概论与其余计算机语言对比较，MATLAB语言突出的特色是什么？MATLAB拥有功能强盛、使用方便、输入简捷、库函数丰富、开放性强等特色。

MATLAB系统由那些部分构成？MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。

安装MATLAB时，在选择组件窗口中哪些部分一定勾选，没有勾选的部分此后怎样补安装？在安装MATLAB时，安装内容由选择组件窗口中个复选框能否被勾选来决定，能够依据自己的需要选择安装内容，但基本平台（即MATLAB选项）一定安装。

第一次安装没有选择的内容在补安装时只要依据安装的过程进行，不过在选择组件时只勾选要补装的组件或工具箱即可。

MATLAB操作桌面有几个窗口？怎样使某个窗口离开桌面成为独立窗口？又怎样将离开出去的窗口从头搁置到桌面上？在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是封闭窗口的Close按钮，一个是能够使窗口成为独立窗口的Undock 按钮，点击Undock按钮就能够使该窗口离开桌面成为独立窗口，在独立窗口的view菜单中选择Dock,,菜单项就能够将独立的窗口从头防备的桌面上。

怎样启动M文件编写/调试器？在操作桌面上选择“成立新文件”或“翻开文件”操作时，M文件编写/调试器将被启动。

在命令窗口中键入edit命令时也能够启动M文件编写/调试器。

储存在工作空间中的数组能编写吗？怎样操作？储存在工作空间的数组能够经过数组编写器进行编写：在工作空间阅读器中双击要编写的数组名翻开数组编写器，再选中要改正的数据单元，输入改正内容即可。

命令历史窗口除了能够察看前方键入的命令外，还有什么用途？命令历史窗口除了用于查问从前键入的命令外，还能够直接履行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。

怎样设置目前目录和搜寻路径，在目前目录上的文件和在搜寻路径上的文件有什么差别？目前目录能够在目前目录阅读器窗口左上方的输入栏中设置，搜寻路径能够经过选择操作桌面的file菜单中的SetPath菜单项来达成。

MATLAB课后习题集附标准答案

4、MATLAB 操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？聞創沟燴鐺險爱氇谴净。

答：在MATLAB 操作桌面上有五个窗口，在每个窗口的右下角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口称为独立的Undock 按钮，点击Undock 按钮就可以使该窗口脱离桌面称为独立窗口，在独立窗口的view 菜单中选择Dock，菜单项就可以将独立的窗口重新防止的桌面上.残骛楼諍锩瀨濟溆塹籟。

5、如何启动M 文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M 文件编辑/调试器将被启动.在命令窗口中键入edit 命令时也可以启动M 文件编辑/调试器.酽锕极額閉镇桧猪訣锥。

6、存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可.彈贸摄尔霁毙攬砖卤庑。

7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中. 謀荞抟箧飆鐸怼类蒋薔。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

第一部分 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.52. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A(5) [A,B]和[A([1,3],:);B^2]3. 设有矩阵A 和B123453166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C 。

(2) 将矩阵C 的右下角3×2子矩阵赋给D 。

(3) 查看MATLAB 工作空间的使用情况。

4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。

(2) 建立一个字符串向量，删除其中的大写字母。

2. 产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。

为什么？ 3. 建立一个5×5矩阵，求它的行列式值、迹、秩和范数。

4. 已知2961820512885A -⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦求A 的特征值及特征向量，并分析其数学意义。

5. 下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。

(2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。

(3) 计算系数矩阵A 的条件数并分析结论。

6. 建立A 矩阵，试比较sqrtm(A)和sqrt(A)，分析它们的区别。

第三部分选择结构程序设计1. 求分段函数的值。

2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。

2. 输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E 。

其中90分~100分为A ，80分~89分为B ，79分~79分为C ，60分~69分为D ，60分以下为E 。

要求：(1) 分别用if 语句和switch 语句实现。

(2) 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。

3. 硅谷公司员工的工资计算方法如下：(1) 工作时数超过120小时者，超过部分加发15%。

(2) 工作时数低于60小时者，扣发700元。

(3) 其余按每小时84元计发。

试编程按输入的工号和该号员工的工时数，计算应发工资。

4. 设计程序，完成两位数的加、减、乘、除四则运算，即产生两个两位随机整数，再输入一个运算符号，做相应的运算，并显示相应的结果。

5. 建立5×6矩阵，要求输出矩阵第n 行元素。

当n 值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。

第四部分循环结构程序设计1. 根据2222211116123n π=++++ ，求π的近似值。

当n 分别取100、1000、10000时，结果是多少？要求：分别用循环结构和向量运算（使用sum 函数）来实现。

2. 根据11113521y n =++++- ，求： (1) y<3时的最大n 值。

(2) 与(1)的n 值对应的y 值。

3. 考虑以下迭代公式：1n nax b x +=+ 其中a 、b 为正的学数。

(1) 编写程序求迭代的结果，迭代的终止条件为|x n+1-x n |≤10-5，迭代初值x 0=1.0，迭代次数不超过500次。

(2) 如果迭代过程收敛于r ，那么r 的准确值是242b b a-±+，当(a,b)的值取(1,1)、(8,3)、(10,0.1)时，分别对迭代结果和准确值进行比较。

4. 已知12312311021323n n n n f n f n f n f f f f n ---==⎧⎪==⎪⎨==⎪⎪=-+>⎩求f 1~f 100中：(1) 最大值、最小值、各数之和。

(2) 正数、零、负数的个数。

5. 若两个连续自然数的乘积减1是素数，则称这两个边疆自然数是亲密数对，该素数是亲密素数。

例如，2×3-1=5，由于5是素数，所以2和3是亲密数，5是亲密素数。

求[2,50]区间内：(1) 亲密数对的对数。

(2) 与上述亲密数对对应的所有亲密素数之和。

第五部分函数文件1. 定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数文件。

2. 一物理系统可用下列方程组来表示：11121112220cos sin 0sin 0cos 000sin 00cos 1a m m a m g m N m N m g θθθθθθ--⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦从键盘输入m 1、m 2和θ的值，求a 1、a 2、N 1和N 2的值。

其中g 取9.8，输入θ时以角度为单位。

要求：定义一个求解线性方程组AX=B 的函数文件，然后在命令文件中调用该函数文件。

3. 一个自然数是素数，且它的数字位置经过任意对换后仍为素数。

例如13是绝对素数。

试求所有两位绝对素数。

要求：定义一个判断素数的函数文件。

4. 设2411()(2)0.1(3)0.01f x x x =+-+-+，编写一个MATLAB 函数文件fx.m ，使得调用f(x)时，x 可用矩阵代入，得出的f(x)为同阶矩阵。

5. 已知(40)(30)(20)f y f f =+(1) 当f(n)=n+10ln(n 2+5)时，求y 的值。

(2) 当f(n)=1×2+2×3+3×4+...+n ×(n+1)时，求y 的值。

第六部分高层绘图操作1. 设23sin 0.5cos 1x y x x ⎡⎤=+⎢⎥+⎣⎦，在x=0~2π区间取101点，绘制函数的曲线。

2. 已知y 1=x 2，y 2=cos(2x )，y 3=y 1×y 2，完成下列操作：(1) 在同一坐标系下用不同的颜色和线型绘制三条曲线。

(2) 以子图形式绘制三条曲线。

(3) 分别用条形图、阶梯图、杆图和填充图绘制三条曲线。

3. 已知2201ln(1)02x x e y x x x π⎧+≤⎪⎪=⎨⎪++>⎪⎩在-5≤x ≤5区间绘制函数曲线。

4. 绘制极坐标曲线ρ=asin(b+n θ)，并分析参数a 、b 、n 对曲线形状的影响。

5. 绘制函数的曲线图和等高线。

224cos cos x y z x ye+-=其中x 的21个值均匀分布[-5，5]范围，y 的31个值均匀分布在[0，10]，要求使用subplot(2,1,1)和subplot(2,1,2)将产生的曲面图和等高线图画在同一个窗口上。

6. 绘制曲面图形，并进行插值着色处理。

cos cos 3cos sin 0,022sin x s t y s t s t z sππ=⎧⎪⎪=≤≤≤≤⎨⎪=⎪⎩第七部分低层绘图操作1. 建立一个图形窗口，使之背景颜色为红色，并在窗口上保留原有的菜单项，而且在按下鼠标器的左键之后显示出Left Button Pressed 字样。

2. 先利用默认属性绘制曲线y=x 2e 2x ，然后通过图形句柄操作来改变曲线的颜色、线型和线宽，并利用文件对象给曲线添加文字标注。

3. 利用曲面对象绘制曲面v(x,t)=10e -0.01xsin(2000πt-0.2x+π)。

4. 以任意位置子图形式绘制出正弦、余弦、正切和余切函数曲线。

5. 生成一个圆柱体，并进行光照和材质处理。

第八部分数据处理与多项式计算1. 利用MATLAB 提供的rand 函数生成30000个符合均匀分布的随机数，然后检验随机数的性质：(1) 均值和标准方差。

(2) 最大元素和最小元素。

(3) 大于0.5的随机数个数占总数的百分比。

2. 将100个学生5门功课的成绩存入矩阵P 中，进行如下处理：(1) 分别求每门课的最高分、最低分及相应学生序号。

(2) 分别求每门课的平均分和标准方差。

(3) 5门课总分的最高分、最低分及相应学生序号。

(4) 将5门课总分按从大到小顺序存入zcj 中，相应学生序号存入xsxh 。

提示：上机调试时，为避免输入学生成绩的麻烦，可用取值范围在[45,95]之间的随机矩阵来表示学生成绩。

3. 某气象观测得某日6:00~18:00之间每隔2h 的室内外温度（0C ）如实验表1所示。

实验表1 室内外温度观测结果（0C ）时间h68 10 12 14 16 18 室内温度t1 18.0 20.0 22.0 25.0 30.0 28.0 24.0 室外温度t2 15.019.024.028.034.032.030.0试用三次样条插值分别求出该日室内外6:30~18:30之间每隔2h 各点的近似温度（0C ）。

4. 已知lgx 在[1,101]区间10个整数采样点的函数值如实验表2所示。

实验表2 lgx 在10个采样点的函数值 x 1 11 21 31 41 51 61 71 81 91 101 lgx 0 1.0414 1.3222 1.4914 1.6128 1.7076 1.7853 1.8513 1.9085 1.9510 2.0043试求lgx 的5次拟合多项式p(x)，并绘制出lgx 和p(x)在[1,101]区间的函数曲线。

5. 有3个多项式P1(x)=x4+2x3+4x2+5，P2(x)=x+2，P3(x)=x2+2x+3，试进行下列操作：(1) 求P(x)=P 1(x)+P 2(x)P 3(x)。