MATLAB实验指导书（附答案）(3)

solve(f1,f2, …fn) —— 求n个方程的解 12、符号微分方程求解指令：dsolve

命令格式：dsolve(f,g)： f —— 微分方程，可多至12个微分方程的求

解；g为初始条件，默认自变量为 'x',可任意指定自变量't', 'u'等，微分方程的各阶导数项以大写字母D表示

三、实验内容：

1、生成一个3行3列的随机矩阵，并逆时针旋转90°，左右翻转，上下翻转。 >> a=magic(3) a =

8 1 6 3 5 7 4 9 2

>> b=rot90(a) 逆时针旋转90° b =

6 7 2 1 5 9 8 3 4

>> c=fliplr(b) 左右翻转 c =

2 7 6

10

9 5 1 4 3 8

>> d=flipud(c) 上下翻转 d =

4 3 8 9 5 1 2 7 6

2、已知a=[1 2 3]，b=[4 5 6]， 求a.\\b和a./ b a.\\b ans =

4.0000 2.5000 2.0000 a./ b ans =

0.2500 0.4000 0.5000 3、数组和矩阵有何不同？

数组中的元素可以是字符等，矩阵中的只能是数，这是二者最直观的区别。从外观形状和数据结构上看，二维数组和数学中的矩阵没有区别。但是矩阵作为一种变换或映射算子的体现，矩阵运算有着明确而严格的数学规则。而数组运算是Matlab软件所定义的规则，其目的是为了数据管理方便、操作简单、指令形式自然和执行计算的有效。虽然数组运算尚缺乏严谨的数学推理，而且数组运算仍在完善和成熟中，但是它的作用和影响正随着matlab的发展而扩大。

11

4、已知a=[1 2 3;4 5 6;7 8 0]，求其特征多项式并求其根。 求特征多项式： a =

1 2 3 4 5 6 7 8 0 AA=sym(a) AA =

[ 1, 2, 3] [ 4, 5, 6] [ 7, 8, 0]

>> poly(AA) ans =

x^3-6*x^2-72*x-27 b=poly(a) b =

1.0000 -6.0000 -72.0000 -27.0000

>> c=roots(b)

12

c =

12.1229 -5.7345 -0.3884

5、已知多项式a(x)=x2+2x+3，b(x)=4x2+5x+6，求a，b的积并微分。 a=[1,2,3];b=[4,5,6];polyder(a,b) ans =

16 39 56 27

?x1?2x2?86、求解方程 1）?

2x?3x?13 2?1a=[1 2;2 3];b=[8;13];x=a\\b x = 2 3

?x1?2x2?1?2）?2x1?3x2?2

?3x?4x?32?1a=[1 2;2 3;3 4];b=[1;2;3];x=a\\b x =

1.0000 -0.0000

13

3）??x1?2x2?3x3?1

?2x1?3x2?4x3?2 a=[1 2 3;2 3 4];b=[1;2];x=a\\b x = 1 0 0

7、用两种方法创建符号矩阵，A =[ a, 2*b]

[3*a, 0] 方法一：

A=sym('[a,2*b;3*a,0]') A =

[ a, 2*b] [ 3*a, 0] 方法二：

syms a b,A=[a,2*b;3*a,0] A =

[ a, 2*b] [ 3*a, 0]

8、计算二重不定积分

??xe?xydxdysyms x y;f=int(int(x*exp(-x*y),x),y)

14