整体或局部:手撕矩阵:从零开始

举报
码乐 发表于 2023/11/14 10:28:21 2023/11/14
【摘要】 最珍贵的东西通常是免费的,比如阳光,水,空气,人们一刻也无法离开它们。 彼此彼此,从现实中来,回到现实中去。

承接上节:

最珍贵的东西通常是免费的,比如阳光,水,空气,人们一刻也无法离开它们。
	   ___buddha

彼此彼此,从现实中来,回到现实中去。

2 实例计算

 万剑归宗,万物归一。

2.0 几个例子简单说明

加,减,乘,数乘

加: 具体操作,对应元素相加

A + B 矩阵相加时,必须同类型

A =(-1  2  2 
      2  2  -1 
      2  -1  2 )

B = (1  0  0
     0  1  0
     0  0  1)

A + B = (0  2  2
        2  3  -1
       2  -1  3)

减 与加规则相同

数乘 又称为矩阵广播

A = [-1  2  2
      2  2  -1
      2  -1  2]

2A = [-2  4  4
       4  4  -2
       4  -2  4]

加法和数乘的性质

    设 A,B,C,D是同型矩阵,k,l是数,则
        1), A + B = B + A
        2), (A+B) + C = A + (B+C)
        3), A + O = A
        4), A +(-A) = O
        5), 1A = A
        6), k(lA) = (kl)A
        7), (k+l)A = kA + lA
        8), k(A+B) = kA + kB
        9), kA = O <--> (k=0 或 A = O)
            k !=0 且 A != 0, 则 kA != O, 这与kA=O 矛盾

2.1 矩阵乘

条件: 第一个矩阵的列数 = 第二个矩阵的行数.

当AB 或 BA 不满足矩阵乘的条件时,称AB 或 BA无意义

	 A = (a_ij) m*s 
         B = (b_ij) s*n
         A 与 B的乘积是一个 m * n 的矩阵
         C = (c_ij) m*n,其中
         c_ij = ai1b1j + ai2b2j + ... + a_is*b_sj = 𝞢(k=1,s)a_ik*b_kj
         记为 C = AB
         称AB为 以A左乘B,或 B右乘A

         即 A的第i行 与 B的第j列的 全部元素分别求乘积的和,为新的矩阵的 第i行第j列的 元素值

    AB = O 不能 推断 A 或 B = O

*定义: A 与 B可交换

		AB 和 BA 都意义,且结果相等,称A与B可交换
  • 矩阵乘的性质:

          设k是数,A,B,C是矩阵,如果以下各式,一端有意义,另一端也有意义,并且等式成立
          1,(AB)C = A(BC)
          2,A(B+C) = AB+AC
            (A+B)C = AC+BC
          3,(kA)B = k(AB) = A(kB)
          设A1,A2,...,As 是矩阵,且对于任意 i =1, ..., s-1, A_jA_j+1都有意义,则连乘A1A2...As也有意义		                     
      
      例如:
      A = [-1  2  2
            2  2  -1
           2  -1  2]
    
      B = [0  0  1
           0  1  0
           1  0  0]
    
      各元素值: 
      	ab11 = -1 *0 + 2 * 0 + 2 * 1 = 2   //第一行第一列的值 A的一行各值乘 B 的一列的和
      	ab12 = -1 *0 + 2 * 1 + 2 * 0 = 2   //第一行第二列的值 A的一行各值乘 B 的二列的和
      	ab13 = -1 *1 + 2 * 0 + 2 * 0 = -1   //第一行第三列的值 A的一行各值乘 B 的三列的和
      	
      	如此类推...
    

    最后得到:

      AB = [ 2  2 -1
      	 -1  2  2
      	  2   -1  2
           ]
    

2.2 转换

转置

原矩阵A 的行变列,列变行产生的新矩阵(沿左上到右下对角线翻转) A^T ,称为A的转置

A = (1  2
    3  4)
 
A^t = (1  3
       2  4)

初等变换 矩阵间的三种关系,等价,相似,合同。 形散而神不散

    线性方程组的解决工具,矩阵,行列式,向量

特殊矩阵 比如:方阵 3 x 3

可逆与不可逆

判断矩阵是否可逆,若A可逆,求A^-1

A=(1   2   3
   2   2   1
   3   4   3)

解 初等变换转换为 上三角矩阵,其行列式为

|A| = |1    2   3 |
      |0    -2  -5|
      |0    0   -1|
     = 2 != 0
        
所以矩阵A可逆

2.3 获取矩阵特征值

例如有矩阵如下:

    A = (-1  2  2
          2  2  -1
          2  -1  2)

当且仅当 <=>

	A*->v = λ*->v

这个非零向量 ->v 就是特征向量,λ 就是特征值。
矩阵A乘以向量 v

	单位矩阵再乘以向量v
    
	->O = λI_n*->v - A*->v
	->O = (λI_n - A)*->v

定理:单位矩阵乘以向量v 就等于向量v

那么 -A*->v 就是把矩阵A乘以向量v 从方程两边减掉。

	 (λI_n - A)*->v = ->O    (方程式 I)

	I_n 表示单位矩阵
	->O 表示零向量
	->v 表示向量v

它可以说是,列不是独立线性的,或者行列式为0,或是说 非线性的,或是说矩阵不是可逆的。
λ是矩阵A的特征值,当且仅当以上条件成立。

对于某非零向量,当且仅当λ乘以单位矩阵 I_n 减去矩阵A的行列式等于0,如下:

	det(λI_n - A) = 0      (定义 方程式 I)

对于A 这样的 3 阶矩阵,其单位矩阵为

 I_3 = [1  0  0 
        0  1  0
        0  0  1]

假设A矩阵的特征值为 λ,那么 R3空间的单位矩阵为

λI_3 =[λ  0  0
      0  λ  0
      0  0  λ]

我们根据定义的方程式 I

	λI_3 - A = (λ+1  -2  -2
	            -2   λ-2  λ-2
	            -2   1    λ-2)

通过增广构造3个对角线,需要增广 第一列和第二列

    Right = (1) +(2) +(3) :  (λ+1)(λ-2)(λ-2) +(-2*1*-2) +(-2*-2*1)
Left = (A)+(B)+(C): (-2*-2*(λ-2) + (λ+1) +(-2*-2*(λ-2)))

Right - Left = (λ+1)(λ-2)(λ-2) +4 + 4 - [4(λ-2)+λ+1+4(λ-2)] 
                   
                   = P(λ) = λ^3 -3λ^2 -9 λ + 27

如下图:

矩阵A的特征多项式 为 P(λ) 以上。
当λ真正为矩阵A的特征值时,以上多项式的值为 0 也就是

	P(λ) = 0
	λ^3 -3λ^2 -9 λ + 27 = 0

现在问题变成求3次方程式的根,有两个办法

	1 使用数学定理:如果多项式方程有整数解,那么这个根应该是这个常数的某因子。
	2 数学定理: 变换方程式

我们变化方程式

	λ^2(λ - 3) -9(λ-3) = 0

再次提取因子

	(λ-3)(λ^2-9) =0
	(λ-3)(λ-3)(λ+3)=0

有两个整数解(3是重根),因此特征值 λ 是 3和 -3.

2.4 继上: 获取 特征向量

定义:当特征值矩阵乘以它 等于 0,这个向量就是特征向量。

定理: 不可逆的矩阵才有零空间。
只有行列式为零的矩阵具有零空间。

        A = [-1  2  2
         	  2  2  -1
         	  2  -1  2]

   (λI_n - A)*->v = ->O  <=> A->v = λ->v

矩阵乘以特征向量,对于任何特征值必须为0,代入特征值以取得特征向量。

I 式: 
λI_3 - A = (λ+1  -2  -2
	-2   λ-2  λ-2
	-2   1    λ-2)

当 λ = 3时 I式为

	3I_3 =  [4  -2  -2
		-2  1   1
		-2  1   1]

当 λ = -3时 I 式变成  

  -3I_3 = 	[-2  -2  -2
	-2  -5   1
	-2  1   -5]

矩阵A对于特征值3的特征空间 ->v 或为该矩阵的零空间。 它是λ 乘以单位矩阵减去矩阵A;
所以此矩阵零空间就是特征空间。 因此满足该方程的所有值构成 λ=3 的特征空间的特征向量。

	3I_3*->v = 0

2.5 继上: 获取 特征空间

使用简化行阶矩阵,以此求解零空间,简化3I_3为阶梯矩阵。

	[4  -2  -2      --- 行1除以2
	 -2  1   1      --- *2 + 行1
	 -2  1   1]     --- *2 + 行1

得到

	[2 -1  -1
	 0  0  0
	 0  0  0
	]

设 特征向量 为 ->v  

	[v1
	 v2
	 v3]

期望有:

	[2 -1  -1   * [v1    =   [0
	 0  0  0       v2         0
	 0  0  0       v3         0 
	]             ]           ]

得到一个方程如下:

	2v1 - v2 -v3 =0

	v1 = 1/2(v2+v3) 设 v2 = a, v3 = b

    Eλ=3 = {[v1         [1/2       [1/2
             v2     = a   1    +  b  0       a,b ∈ IR
             v3]          0]         1]
              }

或λ=3的特征空间,写作如下形式:

	Eλ=3 = Spanu([1/2     [1/2
		      1        0
		     0]       1])

同理可计算得出,λ=-3的特征空间,写作如下形式

	Eλ=-3  = [1 0 2     [v1      [0
		   0 1 -1    v2    =   0
		   0 0 0]    v3]       0]

因此有

     v1 = -2v3
     v2 = v3 = t
     v1 = -2t

Eλ=-3 = Spanu([-2 
		 1  
		 1]  )

现在有两个特征值的特征空间:

	Eλ=3 = Spanu([1/2     [1/2
		      1        0
		      0]       1])


	Eλ=-3 = Spanu([-2 
			1  
			1]  )

在已知的任何该平面中的单位特征向量 ->x

当施加变换,比如将其乘以 矩阵 A,由于处于特征空间中,特征值为 3 那就是将 ->x 变成 ->3x

矩阵A乘以向量 ->x 就是 A*->x. R3空间

             ______________
	    / ^           /
	   /  |          /
	  /   /->x      /
	 /   /         /
	/_____________/ 	

乘A变换后的R3空间

     ______________
    / ^          /
   /  |         /
  /   /--->Ax  /
 /   /        /
/____________/ 

本节小结:

本节举例基础计算方式,下一节讨论,计算空间。

不要白嫖,请关注我的持续更新。
【版权声明】本文为华为云社区用户原创内容,未经允许不得转载,如需转载请自行联系原作者进行授权。如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。