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

matlab程序设计与应用第二版习题答案

matlab程序设计与应用第二版习题答案Matlab程序设计与应用第二版习题答案Matlab是一种强大的数学软件，广泛应用于科学计算、数据分析和工程设计等领域。

《Matlab程序设计与应用》是一本经典的教材，对于学习和掌握Matlab编程语言具有重要的意义。

本文将为大家提供《Matlab程序设计与应用第二版》中部分习题的答案，帮助读者更好地理解和应用Matlab。

第一章：Matlab基础1.1 基本操作1. a = 3; b = 4; c = sqrt(a^2 + b^2); disp(c);2. x = linspace(-pi, pi, 100); y = sin(x); plot(x, y);3. A = [1 2 3; 4 5 6; 7 8 9]; B = [9 8 7; 6 5 4; 3 2 1]; C = A + B; disp(C);1.2 控制结构1. for i = 1:10disp(i);end2. n = 0; sum = 0; while sum < 100n = n + 1;sum = sum + n;enddisp(n);3. x = 5; if x > 0disp('x is positive');elseif x < 0disp('x is negative');elsedisp('x is zero');end第二章：向量和矩阵运算2.1 向量运算1. A = [1 2 3]; B = [4 5 6]; C = A .* B; disp(C);2. A = [1 2 3]; B = [4 5 6]; C = A ./ B; disp(C);3. A = [1 2 3]; B = [4 5 6]; C = dot(A, B); disp(C);2.2 矩阵运算1. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = A * B; disp(C);2. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = B * A; disp(C);3. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = A .* B; disp(C); 第三章：函数和脚本文件3.1 函数1. function y = myfunc(x)y = x^2 + 3*x + 2;end2. function [y1, y2] = myfunc(x1, x2)y1 = x1^2 + 3*x1 + 2;y2 = x2^2 + 3*x2 + 2;end3. function [y1, y2] = myfunc(x)y1 = x^2 + 3*x + 2;y2 = sin(x);end3.2 脚本文件1. x = linspace(0, 2*pi, 100); y = sin(x); plot(x, y);2. x = linspace(-10, 10, 100); y = x.^2 + 3*x + 2; plot(x, y);3. x = linspace(0, 2*pi, 100); y1 = sin(x); y2 = cos(x); plot(x, y1, x, y2);通过以上习题的答案，读者可以对Matlab程序设计的基本语法和常用函数有一个初步的了解。

matlab课后习题解答第二章

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

〖解答〗c1=3/7+c2=sym(3/7+c3=sym('3/7+')c4=vpa(sym(3/7+)Cs1=class(c1)Cs2=class(c2)Cs3=class(c3)Cs4=class(c4)c1 =c2 =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 = za11a12a 135 求符号矩阵 A a 21 a 22 a 23 的队列式值和逆，所得结果应采纳 “子a31a32a33表达式置换”简短化。

〖目的〗理解 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）经过符号计算求 y(t)sin t 的导数dy。

matlab 部分答案

4. 用三点公式计算 y=f(x)在 x=1.0,1.2 处的导数值， f(x)值由下表给出； x f(x)
5
1.0 0.25
4
1.1
ww w
1.2
2
.k
1.3 0.1890 1.4 0.1736
后
7. 求解下列线性常微分方程的解析解。
答
6．设方程的根为 x = [−3, −5, −8, −9] ，求它们对应的 x 多项式的系数。
y=x(:,1);
20．求解线性常微分方程 3 y + 4 y + 5 y + 6 y = 3u + 0.5u + 4u ，在输入 u(t)为单位脉冲
''' '' ' '' '
并单位阶跃信号时的解析解。
第四章习题
1．若 x(n) = cos(
课
2．求有限长序列 x( n) = 5(0.6) ( 0 ≤ n < 20 )的圆周移位 f (n) = x((n − 10)) 20 R20 ( n) 。
5．创建 3×4 矩阵魔方阵和相应的随机矩阵，将两个矩阵并接起来，然后提取任意两个列向
6．创建一个 4×4 单位阵，提取对角线以上部分。
8．创建一个 5×5 随机阵并求其逆。 9．利用上题的矩阵，计算矩阵的 5 次方。
求 C = A * B ， D = A.* B 。
B=[5,4,3,-2; 6,-2,3,-8; -1,3,-9,7]; C=A*B
部分习题与解答
第二章习题及部分解答
1 1．计算 y = x3 + ( x − 0.98) 2 /( x + 1.35)3 − 5( x + ) ，当 x = 2 和 x = 4 时的值。 x

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）第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编程与工程应用课后答案-第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章)答案培训资料

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章 MATLAB 矩阵运算基础2.1 在MATLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？ >> a=[5 7 3;4 9 1]2.5 计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。

>> a=[5 3 5;3 7 4;7 9 8];>> b=[2 4 2;6 7 9;8 3 6]; >> a+b ans =7 7 7 9 14 13 15 12 142.6 求⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。

>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i];>> x’ ans =4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i2.7 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。

>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8]; >> a.*b ans =12 36 3 8 42 402.9 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。

>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;-0.5118 4.0427 1.33182.10 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。

第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')的值是（）。

MATLAB课后习题集附标准答案

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

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

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

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

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

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没有逆矩阵。

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

第2章
习题
11
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.1
c2=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.5286
c2 =
37/70
c3 =
0.7143
c4 =
0.7143
Cs1 =
double
Cs2 =
sym
Cs3 =
sym
Cs4 =
sym
12
sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)')
w=symsum(subs(h,k,n)*subs(x,k,k-n),n,0,k)%据定义
y1=simple(w)
w =
piecewise([a = b, b^k + b^k*k], [a <> b, (a*a^k - b*b^k)/(a - b)])
y1 =
piecewise([a = b, b^k + b^k*k], [a <> b, (a*a^k - b*b^k)/(a - b)])
〖解答〗
（1）符号积分
syms x clear
syms x
y=exp(-abs(x))*abs(sin(x))
si=vpa(int(y,-10*pi,1.7*pi),64)
y =
abs(sin(x))/exp(abs(x))
si =
1.6
（2）数值计算复验
xx=-10*pi:pi/100:1.7*pi;
15
〖目的〗
符号变量限定性定义的作用。
fourier指令的应用。
〖解答〗
syms A t w
a=sym('a','positive');
f=A*exp(-a*abs(t));
y=fourier(f,t,w)
sn=trapz(exp(-abs(xx)).*abs(sin(xx)))*pi/100
sn =
1.0877
10
〖目的〗
变上限二重积分的符号计算法。
〖解答〗
syms x y
f=x^2+y^2;
r=int(int(f,y,1,x^2),x,1,2)
r =
1006/105
11
〖目的〗
在符号计算中，经常遇到计算结果是特殊经典函数的情况。
如何应用subs获得超过16位有效数字的符号数值结果。
初步尝试ezplot指令的简便。
〖解答〗
（1）符号计算
symstx;
f=sin(t)/t;
y=int(f,t,0,x)%将得到一个特殊经典函数
y5=subs(y,x,sym('4.5'))
ezplot(y,[0,2*pi])
y =
sinint(x)
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)]
（2）变换法（复验）
syms z
X=ztrans(a^k,k,z);
H=ztrans(b^k,k,z);
y2=iztrans(H*X,z,k)%通过Z变换及反变换
y2 =
piecewise([b <> 0, (a*a^k)/(a - b) - (b*b^k)/(a - b)])
〖说明〗
符号计算不同途径产生的结果在形式上有可能不同，而且往往无法依靠符号计算本身的指令是它们一致。此时，必须通过手工解决。
d =
1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)
8
〖目的〗
diff, limit指令的应用。
如何理解运行结果。
〖解答〗
syms t
y=abs(sin(t))
d=diff(y) %求dy/dt
d0_=limit(d,t,0,'left') %求dy/dt|t=0-
dpi_2=limit(d,t,pi/2) %求dy/dt|t=pi/2
y =
abs(sin(t))
d =
sign(sin(t))*cos(t)
d0_ =
-1
dpi_2 =
0
9
〖目的〗
符号积分的解析解和符号数值解。
符号计算和数值计算的相互校验。
y5 =
1.515
（2）数值计算复验
tt=0:0.001:4.5;
tt(1)=eps;
yn=trapz(sin(tt)./tt)*0.001
yn =
1.6541
12
〖目的〗
一般符号解与高精度符号数值解。
〖解答〗
syms x
syms n positive
f=sin(x)^n;
yn=int(f,x,0,pi/2)
y3s=vpa(subs(yn,n,sym('1/3')))
y3d=vpa(subs(yn,n,1/3))
yn =
beta(1/2, n/2 + 1/2)/2
y3s =
1.656
y3d =
1.553
13
〖目的〗
符号离散卷积直接法和变换法。
〖解答〗
（1）直接法
syms a b kn
x=a^k;
h=b^k;
〖目的〗
理解自由符号变量的确认规则。
〖解答〗
symvar(sym('sin(w*t)'),1)
ans =
w
symvar(sym('a*exp(-X)'),1)
ans =
a
symvar(sym('z*exp(j*thLeabharlann '),1)ans =
z
5
〖目的〗
理解subexpr指令。
〖解答〗
A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]')