MATLAB图像处理-张卫平(节选)

第四讲 MATLAB语言(三)应用

【目录】

一、矩  阵运算 2

(1) 概述 2

(2) 矩阵的加减法 2

(3) 矢量积与转置 3

(4) 复数的共轭与转置 3

(5) 矩阵的乘法 4

(6) 矩阵的除法 5

(7) 矩阵的乘幂 6

(8) 解线性方程 6

二、多项式运算 7

(1) 多项式的表示 8

(2) 多项式的值 8

(3) 多项式的根 8

(4) 多项式的系数 9

(5) 多项式的积： 9

(6) 多项式的商： 9

(7) 多项式的导数 10

(8) 多项式的曲线拟合 10

三、字  符运算 11

(1) 字符数组 11

(2) 字符与数值的转换 11

(3) 字符串比较 12

四、符  号运算 12

(1) 符号定义 12

(2) 因式分解 13

(3) 简化 13

(4) 反函数 14

(5) 求和 14

(6) diff(f,v,n)-求导 15

(7) int(f,v,a,b)-积分 15

五、Notebook操作 15

【正文】

一、矩  阵运算

(1) 概述

矩阵与数组在数据形式上是一致的，但在运算上是有区别的。Matlab的线性工具箱在目录matfun中。几个常用的函数：

norm － 矩阵的范数

rank － 矩阵的秩

det  － 方阵的行列式的值

inv  － 逆矩阵

(2) 矩阵的加减法

大小相等的矩阵才能加减：

【例】

A=pascal(3)

B=magic(3)

C=A+B

D=A-B

A =

     1     1     1

     1     2     3

     1     3     6

B =

     8     1     6

     3     5     7

     4     9     2

C =

     9     2     7

     4     7    10

     5    12     8

D =

    -7     0    -5

    -2    -3    -4

    -3    -6     4

(3) 矢量积与转置

【例】

A=[1,2,3]   %行向量

B=[3;5;9]   %列向量

C=A*B       %点积

D=B*A       %叉积

E=D'        %D的转置  

A =

     1     2     3

B =

     3

     5

     9

C =

    40

D =

     3     6     9

     5    10    15

     9    18    27

E =

     3     5     9

     6    10    18

     9    15    27  

(4) 复数的共轭与转置

若A为复数矩阵，则A'表示A的共轭，A.'表示转置。

【例】

A=[3+2i,5+3i;2+8i,9+5i]

B=A'                    % 表示共轭

B=A.'                   % 表示转置   

A =

   3.0000 + 2.0000i   5.0000 + 3.0000i

   2.0000 + 8.0000i   9.0000 + 5.0000i

B =

   3.0000 - 2.0000i   2.0000 - 8.0000i

   5.0000 - 3.0000i   9.0000 - 5.0000i

B =

   3.0000 + 2.0000i   2.0000 + 8.0000i

   5.0000 + 3.0000i   9.0000 + 5.0000i  

(5) 矩阵的乘法

【例】

A=[1,2,3;3,5,6;2,3,1]  % 定义两个3×3 的矩阵

B=[5,2,8;9,6,3;2,5,6]

C=A*B                  % 矩阵的乘积：（C≠D）

D=B*A

E=A.*B                 % 数组的相乘：（E＝F）

F=B.*A  

A =

  Columns 1 through 2 

     1     2

     3     5

     2     3

  Column 3 

     3

     6

     1

B =

  Columns 1 through 2 

     5     2

     9     6

     2     5

  Column 3 

     8

     3

     6

C =

  Columns 1 through 2 

    29    29

    72    66

    39    27

  Column 3 

    32

    75

    31

D =

  Columns 1 through 2 

    27    44

    33    57

    29    47

  Column 3 

    35

    66

    42

E =

  Columns 1 through 2 

     5     4

    27    30

     4    15

  Column 3 

    24

    18

     6

F =

  Columns 1 through 2 

     5     4

    27    30

     4    15

  Column 3 

    24

    18

     6  

(6) 矩阵的除法

【例】

A=[1,2,3;3,5,6;2,3,1]   % 定义两个3×3的矩阵

B=[5,2,8;9,6,3;2,5,6]

C=A/B                   % 矩阵的除法

D=A*B^(-1)

E=A*inv(B)

F=A./B                  % 数组的除法  

A =

     1     2     3

     3     5     6

     2     3     1

B =

     5     2     8

     9     6     3

     2     5     6

C =

    0.0879   -0.0256    0.3956

    0.0769    0.1026    0.8462

   -0.2527    0.2821    0.3626

D =

    0.0879   -0.0256    0.3956

    0.0769    0.1026    0.8462

   -0.2527    0.2821    0.3626

E =

    0.0879   -0.0256    0.3956

    0.0769    0.1026    0.8462

   -0.2527    0.2821    0.3626

F =

    0.2000    1.0000    0.3750

    0.3333    0.8333    2.0000

    1.0000    0.6000    0.1667  

(7) 矩阵的乘幂

A=[1,2,3;3,5,6;2,3,1]  % 定义1个3×3的矩阵

C=A^3                  % 矩阵的乘幂

D=A*A*A

E=A^(-0.1)

F=A.^3                 % 数组的乘幂  

A =

     1     2     3

     3     5     6

     2     3     1

C =

   112   185   183

   267   440   429

   129   211   196

D =

   112   185   183

   267   440   429

   129   211   196

E =

  Columns 1 through 2 

   1.0860 - 0.3139i  -0.1373 + 0.1088i

  -0.1790 + 0.1506i   0.9192 - 0.1466i

   0.0065 + 0.0417i  -0.0352 + 0.0835i

  Column 3 

   0.0368 + 0.0499i

  -0.0883 + 0.1753i

   0.8702 - 0.2132i

F =

     1     8    27

    27   125   216

     8    27     1  

(8) 解线性方程

若：ＡＸ=Ｂ －解为Ｘ＝Ｂ＼Ａ；

若：ＸＡ=Ｂ －解为Ｘ＝Ｂ／Ａ。

【例】

A=[2 3 4;5 4 1;1 3 2]

B=[1;2;3]

X=A/B

A =

     2     3     4

     5     4     1

     1     3     2

B =

     1

     2

     3

X =

   -0.8148

    1.6667

   -0.5926  

【例】

A=[2 3 4;5 4 1;1 3 2]

%B=[4 7 5;1 4 7;5 2 6]

B=2*A

X=A/B