matlab第二章习题答案

MATLAB程序设计与应用第二版课后题答案2(最新整理)

-0.6863 0.5621 0.4615
-0.0937 -0.6976 0.7103
D=
-0.0166
0
0
0 1.4801
0
0
0 2.5365
第四章
1.a=input('请输入一个 4 位数:');
while (a<1000|a>9999)
a=input('输入错误，请重新输入一个4位数:');
end
9
function f=factor(n,m); y=0; for k=1:n
y=y+k^m; end
10.(1)S=108 (2)x=4 12 20
y=2 4 6
1. (1) x=-10:0.1:10; y=100./(1+x.^2); plot(x,y)
第五章
(2) x=-10:0.1:10; y=1/(2*pi)*exp(-x.^2/2); plot(x,y)
第二章 3.设矩阵 A 为 A=[24 23 9 21 6;65 74 24 11 21;34 5 98 75 21;8 42 42 53 121;43 21 45 64 21]; （1） B=A(2:5,1:2:5)
B=
65 24 21 34 98 21 8 42 121 43 45 21 （2）A(7)=[] A= 24 65 34 8 43 23 5 42 21 9 24 98 42 45 21 11 75 53 64 6 21 21 121 21 （3）A+30 （4）size(A);ndims(A) （5）题目有误（6）reshape(x,3,4) （7）abs(x) （8）char(x) 4. L1 = 000010000 L2 =

MATLAB原理及应用实验报告第二章答案

实验二数组(矩阵)及其运算《MATLAB原理及应用》实验报告1 •课后练习122-1 3 •1、建立/二25-1和B二31-541034-11⑴求矩阵A和B的乘积，矩阵A左除B,以及矩阵A的2次方⑵求数组A和B的乘积，数组A左除B ,及数组A的2次方程序如下:>)A二［1 2 1； 2 5 -1； 4 10 3］；»B=［2 —1 3 ；3 1 -5 ； 4 — 1 1］；}) C1=A*B；>} C2=A\B；>> C3=A~2；/ ? DI = A •卡B：)> D2 = A. \B；>> D3=A「2；)> Cl运行后显示：C1 =12 0 -615 4-2050 3 -35>> C2C2 =9o S000 4. 0 0 0 0 一了. 5 000 19, 3333 8. 6 6 6 7 -15, 666744, 5 000 19o 00 0 0 -3 5.500 0» C3C3 =9 2 2 28 19 -63 6 88 3» D1D1 二2-2 36 5 51 6 -1 0 3} > D2D2 =0. 5000 -2. 0000 0. 3 3330. 6667 5 o 0000 Oo 2 0 0 01。

0 0 00 - 10o 0 0 00 3. 0 00 0} > D3D3 =1 4 14 25 116 1 0 0 916 3213-1一45 10118 79A =9 6812-424 1514 1 -5152.建立矩阵< 1)用两种方法索引出A矩阵第3行第2列的元素.并将其值改为自己的学号加20(2>索引出A矩阵第2行至第4行.第二列至第5列的所有元素程序如下：〉> A二［16 3 2 1 3 —1 —4； 5 10 11 8 7 9； 9 6 8 1 2 -4 2 ：415 14 1 -5 15]A =16 3 2 13 1 -45 10 1 1 8 7 99 6 8 12 —4 24 15 1 4 1 -5 15> > A (3, 2)an s =6>} A(7)ans =6» A(3, 2) = 21A 二1 6 32 13 -1 -45 1 0 11 8 7 g9 2 18 1 2 -4 24 15 1 4 1 -5 15(2) >> B=A(2:4, 2： 5)B =1 011 8 12 1 8 12 -41 5 1 4 1 一53、使用两神方法建立范围为［1°20］的向呈，使得向量中的元素相邻元袤的间隔是2 <1>改变第二个元素的值，并将其赋给一个新的变量〈学号加20〉，并求两个向量的点积(2)从第二个元素开始提取三个元素，并与向量［123］做叉积程序如下:» a=l 0 :2:20 %求 (1) ID 12 16IS20 i nspace(10. 20< 6) 14 IS20b (2) =21 21 20 >> dotb 〉 a ns 1 5 28 b=(3 : 5) %求(2) [12 3] >> cro & s(b,c )4、复数可二彳+及计程序如下:» zl =3+4* i: z2= 1 +2*i： z3=2* e xp z=zl z2/z30o 3349 4 5。

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 的行数.渗釤呛俨匀谔鱉调硯錦。

(完整word版)含答案《MATLAB实用教程》

第二章 MATLAB 语言及应用实验项目实验一 MATLAB 数值计算三、实验内容与步骤1.创建矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a（1（2）用（3）用（42.矩阵的运算（1）利用矩阵除法解线性方程组。

⎪⎪⎩⎪⎪⎨⎧=+++=-+-=+++=+-12224732258232432143214321421x x x x x x x x x x x x x x x 将方程表示为AX=B ，计算X=A\B 。

（2）利用矩阵的基本运算求解矩阵方程。

已知矩阵A 和B 满足关系式A -1BA=6A+BA ，计算矩阵B 。

其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=7/10004/10003/1A ，Ps: format rata=[1/3 0 0;0 1/4 0;0 0 1/7];b=inv(a)*inv(inv(a)-eye(3))*6*a（3）计算矩阵的特征值和特征向量。

已知矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=1104152021X ，计算其特征值和特征向量。

（4）Page:322利用数学函数进行矩阵运算。

已知传递函数G(s)=1/(2s+1),计算幅频特性Lw=-20lg(1)2(2w )和相频特性Fw=-arctan(2w),w 的范围为[0.01，10]，按对数均匀分布。

3.多项式的运算（1）多项式的运算。

已知表达式G(x)=(x-4)(x+5)(x 2-6x+9),展开多项式形式，并计算当x 在[0，20]内变化时G(x)的值，计算出G(x)=0的根。

Page 324（2）多项式的拟合与插值。

将多项式G(x)=x 4-5x 3-17x 2+129x-180，当x 在[0，20]多项式的值上下加上随机数的偏差构成y1，对y1进行拟合。

对G(x)和y1分别进行插值，计算在5.5处的值。

Page 325 四、思考练习题1.使用logspace 函数创建0～4π的行向量，有20个元素，查看其元素分布情况。

Ps: logspace(log10(0),log10(4*pi),20) (2) sort(c,2) %顺序排列 3.1多项式1）f(x)=2x 2+3x+5x+8用向量表示该多项式，并计算f(10)值. 2）根据多项式的根[-0.5 -3+4i -3-4i]创建多项式。

[整理版]matlab课后习题解答第二章

第2章0000符号运算0000习题2及解答00001 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象？00003/7+0.1; sym(3/7+0.1); sym('3/7+0.1'); vpa(sym(3/7+0.1))0000〖目的〗0000不能从显示形式判断数据类型，而必须依靠class指令。

0000〖解答〗0000c1=3/7+0.10000c2=sym(3/7+0.1)0000c3=sym('3/7+0.1')0000c4=vpa(sym(3/7+0.1))0000Cs1=class(c1)0000Cs2=class(c2)0000Cs3=class(c3)0000Cs4=class(c4) 0000c1 =00000.52860000c2 =000037/700000c3 =00000.528571428571428571428571428571430000c4 =00000.528571428571428571428571428571430000Cs1 =0000double0000Cs2 =0000sym0000Cs3 =0000sym0000Cs4 =0000sym00002 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是自由符号变量.0000sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)')0000〖目的〗0000● 理解自由符号变量的确认规则。

0000〖解答〗0000symvar(sym('sin(w*t)'),1) 0 0 0 0ans = 0 0 0 0 w 0 0 0 0symvar(sym('a*exp(-X)'),1) 0 0 0 0ans = 0 0 0 0 a 0 0 0 0symvar(sym('z*exp(j*th)'),1) 0 0 0 0ans = 0 0 0 0 z 0 0 0 05求符号矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333231232221131211a a a a a a a a a A 的行列式值和逆，所得结果应采用“子表达式置换”简洁化。

matlab第二章习题答案

第一大题：（1）a = 7/3b = sym(7/3)c = sym(7/3,'d')d = sym('7/3')v1=vpa(abs(a-d))v2=vpa(abs(b-d))v3=vpa(abs(c-d))a =2.3333b =7/3c =2.3333333333333334813630699500209d =7/3v1 =0.0v2 =0.0v3 =0.00000000000000014802973661668756666666667788716（2）a = pi/3b = sym(pi/3)c = sym(pi/3,'d')d = sym('pi/3')v1=vpa(abs(a-d))v2=vpa(abs(b-d))v3=vpa(abs(c-d))a =1.0472b =pi/3c =1.047197551196597631317786181171d =pi/3v1 =0.0v2 =0.0v3 =0.00000000000000011483642827992216762806615818554(3)a = pi*3^(1/3)b = sym(pi*3^(1/3))c = sym(pi*3^(1/3),'d')d = sym('pi*3^(1/3)')v1=vpa(abs(a-d))v2=vpa(abs(b-d))v3=vpa(abs(c-d))a =4.5310b =1275352044764433/281474976710656c =4.5309606547207899041040946030989d =pi*3^(1/3)v1 =0.00000000000000026601114166290944374842393221638 v2 =0.00000000000000026601114166290944374842393221638 v3 =0.0000000000000002660111416629094726767991785515第二大题：（1）c1=3/7+0.1c1 =0.5286双精度（2）c2=sym(3/7+0.1)c2 =37/70符号（3）c3=vpa(sym(3/7+0.1))c3 =0.52857142857142857142857142857143完整显示精度第三大题：（1）findsym(sym('sin(w*t)'),1)ans =w(2)findsym(sym('a*exp(-X)' ) ,1)ans =a(3)findsym(sym('z*exp(j*theta)'),1)ans =z第四大题：A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]')A =[ a11, a12, a13][ a21, a22, a23][ a31, a32, a33]DA=det(A)DA =a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31w=inv(A)w =[ (a22*a33 - a23*a32)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 +a12*a23*a31 + a13*a21*a32 - a13*a22*a31), -(a12*a33 -a13*a32)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 +a13*a21*a32 - a13*a22*a31), (a12*a23 - a13*a22)/(a11*a22*a33 -a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)] [ -(a21*a33 - a23*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 +a12*a23*a31 + a13*a21*a32 - a13*a22*a31), (a11*a33 -a13*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 +a13*a21*a32 - a13*a22*a31), -(a11*a23 - a13*a21)/(a11*a22*a33 -a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)] [ (a21*a32 - a22*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 +a12*a23*a31 + a13*a21*a32 - a13*a22*a31), -(a11*a32 -a12*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 +a13*a21*a32 - a13*a22*a31), (a11*a22 - a12*a21)/(a11*a22*a33 -a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)] IAs=subexpr(w,'d')d =1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)IAs =[ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 -a13*a22)][ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 -a13*a21)][ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 -a12*a21)]第六大题：syms ksyms x positives_s=2/(2*k+1)*((x-1)/(x+1))^(2*k+1)s_ss=simple(symsum(s_s,k,0,inf))s_s =(2*((x - 1)/(x + 1))^(2*k + 1))/(2*k + 1)警告: simple will be removed in a future release. Use simplify instead. [> In sym.simple at 41]s_ss =log(x)第八大题：syms x clearsyms xh=exp(-abs(x))*abs(sin(x))si=vpa(int(h,-5*pi,1.7*pi),64)h =abs(sin(x))*exp(-abs(x))si =1.087849417255503701102633764498941389696991336803454392428439159 第九大题：syms x y clearsyms x yr=int(int(x^2+y^2,y,1,x^2),x,1,2)r =1006/105第十大题：syms t x;f=sin(t)/t;y=int(f,t,0,x)y1=subs(y,x,sym('4.5'))ezplot(y,[0,2*pi])y =sinint(x)y1 =syms x clearsyms x ny=sin(x)^nyn=int(y,0,1/2*pi)y31=vpa(subs(yn,n,sym('1/3')))y32=vpa(subs(yn,n,1/3))y =sin(x)^nyn =piecewise([-1 < real(n), beta(1/2, n/2 + 1/2)/2], [real(n) <= -1, int(x^n/(1 - x^2)^(1/2), x, 0, 1)])y31 =1.2935547796148952674767575125656y32 =1.2935547796148952674767575125656第二十题:clearsyms y xy=dsolve('(Dy*y)/5+x/4=0','x')y =2^(1/2)*(C6 - (5*x^2)/8)^(1/2)-2^(1/2)*(C6 - (5*x^2)/8)^(1/2)y1=subs(y,'C6',1)y1 =2^(1/2)*(1 - (5*x^2)/8)^(1/2)-2^(1/2)*(1 - (5*x^2)/8)^(1/2)clfhy1=ezplot(y1(1),[-2,2,-2,2],1)set(hy1,'Color','r')grid onhold onhy2=ezplot(y1(2),[-2,2,-2,2],1)set(hy2,'Color','b')grid onxlabel('Y')ylabel('X')hold offbox onlegend('y(1)','y(2)','Location','Best')hy1 =174.0155hy2 =177.0145。

Matlab习题答案

(1) (3− 5i)(4 + 2i) (2) sin (2 − 8i)
参考答案：（1） >> (3-5*i)*(4+2*i) ans =
22.0000 -14.0000i
（2） >> sin(2-8*i) ans =
1.3553e+003 +6.2026e+002i
5．判断下面语句的运算结果。 (1) 4 < 20 (2) 4 <= 20 (3) 4 == 20 (4) 4 ~= 20 (5) 'b'<'B' 参考答案：（1） >> 4<20 ans =
y_nearest(i)=interp1(x,y,scalar_x(i),'nearest'); y_linear(i) =interp1(x,y,scalar_x(i),'linear'); y_spline(i) =interp1(x,y,scalar_x(i),'spline'); y_cubic(i) =interp1(x,y,scalar_x(i),'cubic'); end subplot(2,2,1),plot(x,y,'*'),hold on,plot(scalar_x,y_nearest),title('method=nearest'); subplot(2,2,2),plot(x,y,'*'),hold on,plot(scalar_x,y_linear),title('method=linear'); subplot(2,2,3),plot(x,y,'*'),hold on,plot(scalar_x,y_spline),title('method=spline'); subplot(2,2,4),plot(x,y,'*'),hold on,plot(scalar_x,y_cubic),title('method=cubic'); 得到结果为：

数字图像处理及应用（MATLAB）第2章习题答案

7.平均值说明f (x ，y )的平均值等于其傅里叶变换F (u ，v )在频率原点的值F (0,0)。

2-3证明离散傅里叶变换的频率位移和空间位移性质。

证明：)(2101),(1),(NvyM ux j M x N y e y x f MN v u F +--=-=∑∑=π),(),(1),(100)(21010)(2)(21010000v v u u F dxdy ey x f MNe ey x f MN y Nv v x M u u j M x N y N yv M x u j Nvy M ux j M x N y --==-+---=-=++--=-=∑∑∑∑πππ因为()()v u F y x f ,,⇔ 所以 ),(),(00)(200v v u u F e y x f N y v M x u j --⇔+π2-4小波变换是如何定义的？小波分析的主要优点是什么？小波之所以小，是因为它有衰减性，即是局部非零的；而称为波，则是因为它有波动性，即其取值呈正负相间的振荡形式，将)(2R L 空间的任意函数f (t )在小波基下展开，称其为函数f (t )的连续小波变换。

小波变换是时间（空间）频率的局部化分析，它通过伸缩平移运算对信号（函数）逐步进行多尺度细化，最终达到高频处时间细分，低频处频率细分，能自动适应时频信号的要求从而可聚焦到信号的任意细节，解决了Fourier 变换的困难问题。

2-5 在图像缩放中，采用最近邻域法进行放大时，如果放大倍数太大，可能会出现马赛克效应，这个问题有没有办法解决，或者有所改善。

可以利用线性插值法，当求出的分数地址与像素点不一致时，求出周围四个像素点的距离比，根据该比率，由四个邻域的像素灰度值进行线性插值。

2-6 复合变换的矩阵等于基本变换的矩阵按顺序依次相乘得到的组合矩阵。

即，T=T N T N-1…T 1。

问矩阵顺序的改变能否影响变换的结果。

矩阵顺序的改变不会影响变换的结果。

matlab 部分答案

⎡ 3sin x ⎤ 11．设 y = cos x ⎢0.5 + ⎥ ，把 x = 0 ~ 2π 间分为 125 点，画出以 x 为横坐标， y 为纵坐标 (1 + x 2 ) ⎦ ⎣
的曲线。 MATLAB：x=0:2*pi/125:2*pi; y=cos(x*(0.5+3*sin(x)/(1+x.*x))); plot(x, y) 12．设 x = z sin 3z ， y = z cos 3 z ，要求画出在 z = −45 ~ 45 区间内画出 x, y, z 三维曲线。 MATLAB: z=-45:0.02:45; x=z.*sin(3*z); y=z.*cos(3*z);
2．在某处测得海洋不同深度处水温如下：
深度(M) 水温 446 7.04 714 4.28 950 3.40 1422 2.54
利用分段线性插值函数，求在深度 500 米、1000 米、1500 米处的水温。 3. 已知四点 P1(1,0) ，P2(0,1)，P1(-1,0)，P4(0,-1)，利用样条插值函数画一通过这四点的圆。
ww w
.k
hd a
w. c
方法二：建立函数文件myfun.m
om
MATLAB：方法一：>> y=quad('exp(-x)+x.^2',0,1)
⎧12 y 2 f ( x, y ) = ⎨ ⎩ 0
求 E ( X ) 、 E (Y ) 、 E ( XY ) 。
0 ≤ y ≤ x ≤1 其它
18．生成一个 4 × 4 的随机矩阵，并对其进行三角分解和正交分解。
10．设有对称实矩阵
⎡2 4 9 ⎤ ⎥ a=⎢ ⎢4 2 4 ⎥ ⎢ ⎣9 4 18⎥ ⎦

matlab编程与工程应用课后答案-第2章

matlab编程与工程应用课后答案-第2章1%计算空心球体积r0 = 3;ri = 2;volume = 4*pi/3*(r0^3-ri^3) —————————————————————————————————————2%计算过氧化氢分子量H = 15.9994;O = 1.0079;H2O2 = 2*H + 2*O —————————————————————————————————————3%计算字符创长度string = input('enter string: ','s');stringlength = length(string) —————————————————————————————————————%4 输入实数，fprintf输出这个变量，格式为两位小数num = input('enter a num: ');fprintf('The num is %.2f\n',num) —————————————————————————————————————6fprintf('不指定宽度: %f\n',12345.6789)不指定宽度: 12345.678900fprintf('10个字符宽度，4个小数位: %10.4f\n',12345.6789)10个字符宽度，4个小数位: 12345.6789fprintf('10个字符宽度，2个小数位: %10.2f\n',12345.6789)10个字符宽度，2个小数位: 12345.68fprintf('6个字符宽度，4个小数位: %6.4f\n',12345.6789)6个字符宽度，4个小数位: 12345.6789fprintf('2个字符宽度，4个小数位: %2.4f\n',12345.6789)2个字符宽度，4个小数位: 12345.6789 —————————————————————————————————————7fprintf('不指定宽度：%int16\n',12345)不指定宽度：12345nt16fprintf('不指定宽度：%d\n',12345)不指定宽度：12345fprintf('5个字符宽度：%5d\n',12345)5个字符宽度：12345fprintf('8个字符宽度：%8d\n',12345)8个字符宽度：12345fprintf('3个字符宽度：%3d\n',12345)3个字符宽度：12345 —————————————————————————————————————8x = 12.34;y = 4.56;fprintf('x is %.3f\n',x)x is 12.340fprintf('x is %2.0f\n',x)x is 12fprintf('y is %.1f\n',y)y is 4.6fprintf('y is %-6.1f!\n',y)y is 4.6 ! —————————————————————————————————————9%计算矩形面积x = input('enter length:');y = input('enter borad:');area = x*y;fprintf('area is %.2f\n',area) —————————————————————————————————————10string = input('what is your name?','s');fprintf('WoW,your name is %s\n',string) —————————————————————————————————————11string = input('enter your string:','s');fprintf('your string was:''%s''\n',string) %输出单引号是应采用''格式—————————————————————————————————————12v = input('enter the flow in m^3/s:');fprintf('a flow rate of %.3f meters per sec\nis equivalent to %.3f feet per sec\n',v,v/0.028) —————————————————————————————————————13incomeyear = input('enter your income per year: ');fprintf('the range of food expenditure per year %f~%f\nthe range of food expenditure per mounth %f~%f\n',incomeyear*0.08,incomeyear*0.1,incomeyear* 0.08/12,incomeyear*0.1/12)—————————————————————————————————————14wight = input('wight of plant: ');area = input('area of wing: ');fprintf('W/A: %f kg/m^2\n',wight/area)—————————————————————————————————————15x= 10;y =22;plot(x,y,'g+')—————————————————————————————————————16x = -2:0.1:2;plot(x,exp(x))xlabel('x')ylabel('y')title('y=e^x') —————————————————————————————————————17x = 1:5:100;y = sqrt(x);figure(1)plot(x,y)figure(2)bar(x,y)—————————————————————————————————————18略—————————————————————————————————————19x1 = linspace(0,pi,10);figure(1)plot(x1,sin(x1))x2 = linspace(0,pi,100);figure(2)plot(x2,sin(x2)) —————————————————————————————————————20mat = [1000 2000 3000 5000 10000;288 281 269 256 223]';x = mat(:,1);y = mat(:,2);plot(x,y)xlabel('high')ylabel('tempture')title('high-tempture') —————————————————————————————————————21mat1 = randi([50,100],3,6)save randfile.dat mat1 -ascii;mat2 = randi([50,100],2,6)save randfile.dat mat2 -ascii -append;load randfile.dat;;randfile —————————————————————————————————————22mat = rand(2,3)*4-1;save testtan.dat mat -ascii;load testtan.dat;mattan = tan(testtan)—————————————————————————————————————23mat = [89 42 49 55 72 63 68 77 82 76 67;90 45 50 56 59 62 68 75 77 75 66;91 44 43 60 60 60 65 69 74 70 70]save hightemp.dat mat -ascii;load hightemp.dathightemp(:,1) = hightemp(:,1)+1900;hightempsave y2ktemp.dat hightemp -ascii—————————————————————————————————————24%24 Calculates y as a function of xfunction y = fn(x)y = x^3-4*x^2+sin(x); —————————————————————————————————————25%25 Converts from MWh to GJfunction gj = mwh_to_gj(mwh)gj = 3.6*mwh;—————————————————————————————————————26%26 converta from inch/h to meter/sfunction meter_sec = converta(inch_hour)meter_sec = inch_hour*5280*0.3048/3600;—————————————————————————————————————27function Tn = fn27(P,i,n)Tn = P*(1+i)*n;—————————————————————————————————————28略—————————————————————————————————————29function V = fn29(Pt,Ps)V = 1.016*sqrt(Pt-Ps);—————————————————————————————————————30function THR = fn30(A)THR = (220-A)*0.6;—————————————————————————————————————31function outdate = fn31(n)outdate = sqrt(2*pi*n)*(n/exp(1))^n;—————————————————————————————————————32%32 脚本n = input('enter the number of units: ');Cn = costn(n);fprintf('the cost for %d units will be $%.2f\n',n,Cn)%32 mygcost函数function Cn = costn(n)Cn = 5*n^2-44*n+11; —————————————————————————————————————33%33 脚本rain = input('enter the rain: ');snow = fn33(rain);fprintf('the snow is %f\n',snow)%33 调用函数function snow = fn33(rain)snow = rain*6.5; —————————————————————————————————————34%34 脚本s = input('enter long_side: ');V = fn34(s);fprintf('volume is %.2f\n',V)%34 函数function V = fn34(s)V = sqrt(2)*s^3/12; —————————————————————————————————————35%35 pickone(x) returns a random element from vector xfunction outdate = pickone(x)outdate = randi([x(1) x(end)],1,1); —————————————————————————————————————36function outdate = vecout(x)outdate = x:1:x+5;—————————————————————————————————————37%37b = input('enter the first side: ');c = input('enter the second side: ');alpha = input('enter the angle between them: ');a = sqrt(b^2+c^2-2*b*c*cosd(alpha)); %切记这里用sindfprintf('the third side is %.3f\n',a)—————————————————————————————————————38略—————————————————————————————————————39%39mat = [90.5792 27.8498 97.0593;12.6987 54.6882 95.7167;91.3376 95.7507 48.5376;63.2359 96.4889 80.02809.7540 15.7613 14.1886];save floatnums.dat mat -ascii;load floatnums.dat;floatnums = round(floatnums')save intnums.dat floatnums -ascii—————————————————————————————————————40%40costssales = [1100 800;1233 650;1111 1001;1222 1300;999 1221];fprintf('there were %d quarters in the file\n',length(costssales))x = costssales(:,1);y = costssales(:,2);plot(x,'ko')hold onplot(y,'b*')mat = fliplr(costssales)';save newfile.dat mat -ascii。

matlab教程第二章课堂练习及答案

Matlab教程第二章符号计算课堂练习1 创建符号变量有几种方法？MATLAB提供了两种创建符号变量和表达式的函数：sym和syms。

sym用于创建一个符号变量或表达式，用法如x=sym(‘x’) 及f=sym(‘x+y+z’)，syms用于创建多个符号变量，用法如syms x y z。

f=sym(‘x+y+z’)相当于syms x y zf= x+y+z2 下面三种表示方法有什么不同的含义？（1）f=3*x^2+5*x+2（2）f='3*x^2+5*x+2'（3）x=sym('x')f=3*x^2+5*x+2（1）f=3*x^2+5*x+2表示在给定x时，将3*x^2+5*x+2的数值运算结果赋值给变量f，如果没有给定x则指示错误信息。

（2）f='3*x^2+5*x+2'表示将字符串'3*x^2+5*x+2'赋值给字符变量f，没有任何计算含义，因此也不对字符串中的内容做任何分析。

（3）x=sym('x')f=3*x^2+5*x+2表示x是一个符号变量，因此算式f=3*x^2+5*x+2就具有了符号函数的意义，f也自然成为符号变量了。

3 用符号函数法求解方程a t2+b*t+c=0。

>> r=solve('a*t^2+b*t+c=0','t')[ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))]4 用符号计算验证三角等式：sin(ϕ1)cos(ϕ2)-cos(ϕ1)sin(ϕ2) =sin(ϕ1-ϕ2) >> syms phi1 phi2;>> y=simple(sin(phi1)*cos(phi2)-cos(phi1)*sin(phi2)) y =sin(phi1-phi2)5 求矩阵⎥⎦⎤⎢⎣⎡=22211211a a a a A 的行列式值、逆和特征根。

MATLAB运算基础(第2章)答案

实验01讲评、参考答案讲评未交实验报告的同学名单数学：6人（11、12级）信科：12-04, 12-22, 13-47批改情况：问题1：不仔细，式子中出错。

问题2：提交的过程不完整。

问题3：使用语句尾分号(;)不当，提交的过程中不该显示的结果显示。

问题4：截屏窗口没有调整大小。

附参考答案：《MATLAB软件》课内实验王平实验01 MATLAB运算基础（第2章MATLAB数据及其运算）一、实验目的1. 熟悉启动和退出MATLAB 的方法。

2. 熟悉MATLAB 命令窗口的组成。

3. 掌握建立矩阵的方法。

4. 掌握MATLAB 各种表达式的书写规则以及常用函数的使用。

二、实验内容1. 数学表达式计算先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。

1.1 计算三角函数122sin 851z e=+（注意：度要转换成弧度，e 2如何给出）示例：点击Command Window 窗口右上角的，将命令窗口提出来成悬浮窗口，适当调整窗口大小。

命令窗口中的执行过程：1.2 计算自然对数221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦（提示：clc 命令擦除命令窗口，clear 则清除工作空间中的所有变量，使用时注意区别，慎用clear 命令。

应用点乘方）命令窗口中的执行过程：1.3 求数学表达式的一组值0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。

命令窗口中的执行过程：1.4 求分段函数的一组值2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5提示：用逻辑表达式求分段函数值。

命令窗口中的执行过程：1.5 对工作空间的操作接着显示MATLAB当前工作空间的使用情况并保存全部变量提示：用到命令who, whos, save, clear, load，请参考教材相关内容。

MATLAB运算基础(第2章)答案培训资料

M A T L A B运算基础(第2章)答案实验01讲评、参考答案讲评未交实验报告的同学名单批改情况：问题1：不仔细，式子中出错。

问题2：提交的过程不完整。

问题3：使用语句尾分号(;)不当，提交的过程中不该显示的结果显示。

问题4：截屏窗口没有调整大小。

附参考答案：实验01 MATLAB 运算基础（第2章 MATLAB 数据及其运算）一、实验目的1. 熟悉启动和退出MATLAB 的方法。

2. 熟悉MATLAB 命令窗口的组成。

3. 掌握建立矩阵的方法。

4. 掌握MATLAB 各种表达式的书写规则以及常用函数的使用。

二、实验内容1. 数学表达式计算先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。

1.1 计算三角函数122sin 851z e=+（注意：度要转换成弧度，e 2如何给出）《MATLAB 软件》课内实验王平示例：点击Command Window 窗口右上角的，将命令窗口提出来成悬浮窗口，适当调整窗口大小。

命令窗口中的执行过程：1.2 计算自然对数221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦（提示：clc 命令擦除命令窗口，clear 则清除工作空间中的所有变量，使用时注意区别，慎用clear 命令。

应用点乘方）命令窗口中的执行过程：1.3 求数学表达式的一组值0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。

命令窗口中的执行过程：1.4 求分段函数的一组值2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 提示：用逻辑表达式求分段函数值。

命令窗口中的执行过程：1.5 对工作空间的操作接着显示MATLAB当前工作空间的使用情况并保存全部变量提示：用到命令who, whos, save, clear, load，请参考教材相关内容。

matlab课后习题解答第二章

第2章符号运算习题2及解答1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象？3/7+0.1; sym(3/7+0.1); sym('3/7+0.1'); vpa(sym(3/7+0.1))〖目的〗●不能从显示形式判断数据类型，而必须依靠class指令。

〖解答〗c1=3/7+0.1c2=sym(3/7+0.1)c3=sym('3/7+0.1')c4=vpa(sym(3/7+0.1))Cs1=class(c1)Cs2=class(c2)Cs3=class(c3)Cs4=class(c4)c1 =0.5286c2 =37/70c3 =c4 =Cs1 =doubleCs2 =symCs3 =symCs4 =sym2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是自由符号变量.sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)')〖目的〗●理解自由符号变量的确认规则。

〖解答〗symvar(sym('sin(w*t)'),1)ans =wsymvar(sym('a*exp(-X)'),1)ans =asymvar(sym('z*exp(j*th)'),1) ans = z5求符号矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333231232221131211a a a a a a a a a A 的行列式值和逆，所得结果应采用“子表达式置换”简洁化。

〖目的〗● 理解subexpr 指令。

〖解答〗A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]')DA=det(A) IA=inv(A);[IAs,d]=subexpr(IA,d) A =[ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA =a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d =1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)8（1）通过符号计算求t t y sin )(=的导数dtdy。

matlab答案第二章

第二章1．设A和B是两个同维同大小的矩阵，问：（1）A*B和 A.*B的值是否相等？（2）A./B和B.\A的值是否相等？（3）A/B和B\A的值是否相等？（4）A/B和B\A所代表的数学含义是什么？答：（1）不等（2）相等（3）不等（4）A/B含义： A乘B的逆，即A* inv(B)B\A含义： B的逆乘A，即inv(B)*A2.写出完成下列操作的命令（1）将矩阵A第2~5行中第1，3，5列元素赋给矩阵B。

（2）删除矩阵A的第7号元素。

（3）将矩阵A的每个元素值加30。

（4）求矩阵A的大小和维数。

（5）将向量t的0元素用机器0来代替。

（6）将含有12个元素的向量x转换成3x4矩阵。

（7）求一个字符串的ASCII码。

（8）求一个ASCII码所对应的字符。

答：设矩阵A为A=[24 23 9 21 6;65 74 24 11 21;34 5 98 75 21;8 42 42 53 121;43 21 45 64 21];（1） B=A(2:5,1:2:5) 或者 B=A(2:5,[1 3 5])B =65 24 2134 98 218 42 12143 45 21（2）A(7)=[]A =24 65 34 8 43 23 5 42 21 9 24 98 42 45 21 11 75 53 64 6 21 21 121 21（3）A+30（4）size(A); ndims(A)（5）设向量t=[1,2,0,5,6,0,9];p=find(t==0);t(p)=eps;或者：t(find(t==0))=eps;（6）reshape(x,3,4)（7）abs(x) 或 double(x)（8）char(x)3.下列命令执行后，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 4.已知：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=14.35454.9632053256545410778.01023A 完成下列操作：（1）取出A 的前3行构成矩阵B ，前两列构成矩阵C ，右下角3x2子矩阵构成矩阵D ，B 与C 的乘积构成矩阵E 。

matlab 第二章部分答案

1、答：6+7i中的i为虚数单位，是一个复数常量；而6+7*i中的i为一个表达式，是一个运算量。

i是虚数单位，等于sqrt(-1)；I在MATLAB里没有定义，可根据自己的需要定义为各种变量。

2、答：A*B是矩阵和矩阵相乘，A.*B是指它们的对应元素进行乘法运算，要求两矩阵维参数相同。

A./B表示A矩阵除以B矩阵的对应元素，如果A、B两矩阵具有相同的维数，B.\A 等价于A./B。

如果B矩阵是非奇异方阵，A/B等效于B矩阵的逆右乘A矩阵，也就是A*inv(A)；B\A等效于B的逆左乘A矩阵，也就是inv(B)*A。

如果A和B是两个标量数据，则是它们的对应元素进行相关运算。

3、答：（1）A(7)=[ ];（2）find(t==0)t(ans)=eps（3）a=[0:11]b=reshape(a,3,4)（4）abs('a')（5）A=[0:11]ones[A]（6）A=[1,2,3;4,5,6];diag(A)B=diag([1,5])4、答：a=3+sqrt(1)*randn(50,10)5、答：（1）>> A=[1,-1,2,3;5,1,-4,2;3,0,5,2;11,15,0,9];>> diag(A)ans =1159>> triu(A)ans =1 -12 30 1 -4 20 0 5 20 0 0 9>> tril(A)ans =1 0 0 05 1 0 03 0 5 011 15 0 9>> inv(A)ans =-0.1758 0.1641 0.2016 -0.0227 -0.1055 -0.1016 -0.0391 0.0664 -0.0508 -0.0859 0.1516 0.00230.3906 -0.0313 -0.1813 0.0281 >> det(A)ans =1280>> rank(A)ans =4>> norm(A,1)ans =20>> norm(A)ans =21.3005>> norm(A,inf)ans =35>> cond(A,1)ans =14.4531>> cond(A)ans =11.1739>> cond(A,inf)ans =22.0938>> trace(A)ans =16行列式的值：1280秩：4范数：norm(a,1)=20, norm(a)=21.3005, norm(a,inf)=35条件数：cond(a,1)=14.4531,cond(a)=11.1739,cond(a,inf)=22.0938 迹：16（2）>> B=[0.43,43,2;-8.9,4,21];>> diag(B)ans =0.43004.0000>> triu(B)ans =0.4300 43.0000 2.00000 4.0000 21.0000>> tril(B)ans =0.4300 0 0-8.9000 4.0000 0B有逆矩阵的充要条件是|B|≠0，因为B的行列式等于0，所以B没有逆矩阵。

第2章 MATLAB数据及其运算_习题答案

……………………………………………………………最新资料推荐…………………………………………………第2章 MATLAB数据及其运算习题2一、选择题1．下列可作为MA TLAB合法变量名的是（）。

DA．合计B．123 C．@h D．xyz_2a 2．下列数值数据表示中错误的是（）。

CA．+10 B．1.2e-5 C．2e D．2i3．使用语句t=0:7生成的是（）个元素的向量。

AA．8 B．7 C．6 D．54．执行语句A=[1,2,3;4,5,6]后，A(3)的值是（）。

BA．1 B．2 C．3 D．45．已知a为3×3矩阵，则a(:,end)是指（）。

DA．所有元素B．第一行元素C．第三行元素D．第三列元素6．已知a为3×3矩阵，则运行a (1)=[]后（）。

AA．a变成行向量B．a变为2行2列C．a变为3行2列D．a变为2行3列7．在命令行窗口输入下列命令后，x的值是（）。

B>> clear>> x=i*jA．不确定B．-1 C．1D．i*j8．fix(354/100)+mod(354,10)*10的值是（）。

DA．34 B．354 C．453D．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．输入字符串时，要用（）将字符括起来。

CA．[ ] B．{ } C．' ' D．" "12．已知s='显示"hello"'，则s的元素个数是（）。

AA．9 B．11 C．7 D．181 / 31 / 31 / 31……………………………………………………………最新资料推荐…………………………………………………13．eval('sqrt(4)+2')的值是（）。