微分方程第一步,吃透基本概念——复数,多项式方程及矩阵理论

康托的天堂 2022-03-17 23:35:46

最近我开启了“量子力学之路”系列,旨在从数理角度从零解释量子力学。正如我在系列的第一篇文章量子力学之路(1)——坚实的数理基础至关重要,没有捷径可走​中提到的那样,学习量子力学有一些先决条件,而一些先决条件并不简单,如很多数学主题,这些主题我在“量子力学之路”系列中一般都会讲到,但不会深入。因此我决定同步开启“微分方程”系列,这是本系列的第一篇文章。

复数

复数就是形如x+iy的数字,其中x和y是实数,i^2=-1。实数x和y分别称为z的实部和虚部,表示为:

如果z=x+iy,则z̅=x-iy称为z的复共轭。很容易得出:

定理(1)

对于复数z和w,有以下7个性质:

笛卡尔和指数形式

复数可以绘制在一个矩形网格上,类似于实数对(x,y)被绘制在一个直角坐标系统上。你可以简单地用一对实数(x,y)来确定复数z=x+iy,并绘制(x,y)。y轴称为虚轴,x轴称为实轴。

也可以用极坐标(θ, r)来表示:

定理(2)(很重要)

对于复数z=x+iy,

多项式方程的根

设p_n(x)和q_n(x)为n次多项式,可以假设p_n(x)的形式为:

系数a_i可能不是实数,但在本系列文章中它们总是实数。在任何情况下,代数基本定理表明,方程p_n(x)=0有n个解。

假设x=r是p_n(x)的根。那么:

如果p_n(x)除以x-r,就得到恒等式:

其中R为常数,q_(n-1)(x)为n-1次多项式。因此,

但上面的式子是关于x的恒等式,通过令x=r,我们可以得到r =0当且仅当p_n(r)=0,也就是说,当且仅当r是p_n(x)的根。

定理(3)

对于每一个多项式

存在n个复数r_1, r_2,…,r_n,称为多项式的根,使得

P_n (r_i)=0,对于所有i=1,2,…,n

p_n(x)=(x-r_1)(x-r_2)…(x-r_n)

而且,对于一些i,如果p_n(r)=0,则r=r_i。

定理(4)

假设r=a+ib, b≠0,是下面多项式的根

那r的共轭( r̅=a-ib)也是多项式的一个根。

矩阵表示

矩阵是一组数字的矩形数组。一般来说,矩阵用黑体字大写字母表示。我们用A来表示p×q矩阵,它的元素是a_ij。也就是

只有一列的矩阵称为向量。我们用黑体小写字母来表示向量,这与我们对矩阵的约定一致。A的列向量就是

对于联立方程组

系数矩阵为

右边常数用向量b表示:

增广矩阵B由A和b表示:

未知量用向量x表示:

行数和列数相同的矩阵称为方阵。方阵的非对角元素是a_ij,其中i≠j。非对角元素都为零的方阵称为对角阵。

I_n矩阵(其中n是正整数),是所有对角元素都是1的对角矩阵。这些矩阵称为单位矩阵。所以:

是单位矩阵。I_n的列被赋予特殊的符号e_1,e_2,…,e_n;也就是:

当上下文明确了I_n的大小时,就可以去掉下标n。那些位于对角线以下的非对角线项为零的方阵称为上三角矩阵(上三角矩阵同理可得)。

我们定义0矩阵O_n为nXn矩阵,它所有的元素都是0。所以O_n也是对角线。

方程组的解

我们可以用消元法来解方程组。注意消去过程很大程度上依赖于每个方程中未知变量的系数。

把方程组的系数和方程组右边的常数放在一起得到一个增广矩阵。化简这个增广矩阵可以得到方程组的解。注意,当系数矩阵简化为单位矩阵时,右边的系数列就是解向量。

从矩阵A到矩阵B的算术步骤叫做初等行运算。这些运算分为三种类型:

交换任意两行。

将一行乘以一个非零标量。

将一行的α倍添加到另一行的β倍。

这里的关键点是初等行运算用另一个方程组替换了一个方程组,后者的解集与前者的解集相同。这种解法称为高斯消去法。

矩阵代数

设m×n矩阵A和B为:

对于所有的i和j,如果a_ij=b_ij,则A=B。因此,两个矩阵相等,意味着矩阵相应项相等。

我们现在定义矩阵A+B和矩阵与任意标量k的乘法:

由上面的定义可以得到下列代数规则:

A+B=B+A

A+(B+C)=(A+B)+C

A+O=A

A+(-1)A=O

0A=O

k(hA)=(kh)A

k(A+B)=kA+kB

(k+h)A=kA+hA

用第一列替换第一行,用第二列替换第二行,以此类推,直到所有的列都变成行。由这个交换得到的矩阵称为原矩阵的转置,[a_ij]^T=[a_ji]。我们用A^T来表示A的转置。

一个向量的转置是一个只有一行的矩阵,有时称为行向量。为了避免混淆,我们用逗号分隔行向量的各个元素。

请注意,转置和加法的定义引出了这样的结论:C=A+B意味着C^T=A^T+B^T。

如果矩阵A等于它自己的转置,即A^T=A,那么它就是对称的;如果A^T=-A,那么它是反对称的。对称矩阵和反对称矩阵必须是方阵。

矩阵乘法

如果A是m×q矩阵,元素为a_ij ;B是q×n矩阵,元素为b_ij,那么乘积C=AB是m×n矩阵,c_ij为:

为了使上面的定义有意义,A的每一行必须有和B的每一列有一样多的元素。这意味着A的列数必须与B的行数相同。

因此,如果A是2×3矩阵, B是3×3矩阵,那么AB是有定义的,而BA没有定义。因此矩阵乘法是不满足交换律的。

以下事实适用于所有大小兼容的A、B和C:

A(BC)=(AB)C

A(A+C)=AB+AC

(B+C)A=BA+CA

能说明矩阵乘法的特性的一个例子:

这表明即使A和B都不为0,AB也可以是0。

两个矩阵乘积的转置,是它们转置的逆序乘积:

这个结果扩展到三个或更多个矩阵的乘积。

矩阵乘法提供了一种将方程组写成紧凑形式的方法。

上面可以写成Ax=b。我们将反复使用这个表达。在不明确A的大小的情况下,我们假设A是m×n,因此x是一个有n个元素的向量,b是一个有m个元素的向量,尽管在大多数应用中m=n。

矩阵的逆

为了达到类似的目的,我们引入了矩阵逆的符号。如果存在一个方阵B,使AB=I=BA,那么方阵A就是非奇异的,或者说有一个逆,或者说是可逆的。很明显,不是所有的矩阵都有逆矩阵。

由于A的逆矩阵只有一个,所以如果AB=I=BA成立,我们称B为A的逆矩阵,并将B写成A^(-1)。用这种表示法,AB=I=BA可以写成AA^(-1)=(A^(-1))(A)=I。

一个不可逆的方阵就是一个A^(-1)不存在的方阵。这样的矩阵称为奇异或不可逆矩阵。

如果A是2x2矩阵

且ad-bc≠0,那么

定理(5)

假设A和B都是可逆的。那么

(A^(-1))^(-1)=A

(AB)^(-1)=(B^(-1))(A^(-1))

(A^T)^(-1)=(A^(-1))^T

定理(6)

假设A是可逆的。那么Ax=b有且只有一个解,x=A^(-1)b。

结论(1)

当且仅当A是奇异阵时,方程组Ax=0有解x≠0。当且仅当A是可逆的,这个方程组只有解x=0。

行列式

定义与基本定理

A的行列式是一个只在方阵中定义的标量,记作det(A)。它有n的阶乘项,每一项是A的元素的正负乘积:

其中第二个下标,由∗表示,是数字{1,2,…,n}之一,其中没有一个被使用两次。指数k是第二个下标的逆序数。因此,

由于数字{1,2,…,n}的每个排列都有一项,所以上面的和包含n的阶乘项。由于这个原因,实际上并不使用上面的求和来计算。如果n=2,定义是容易使用的,

有两种常用的det(A)的求值方法。在本节中,我们将探讨最高效的方法。该方法依赖于两个基本定理:

定理(7)

如果A是上或下三角矩阵,对角元素为a_11,a_22,…,a_nn,那么det(A)=(a_11)(a_22)…(a_nn)。

定理(8)

设A是一个方阵。

如果A的两行元素互换形成B,那么,det(A)=-det(B)

如果A的一行乘以k得到B,那么kdet(A)=det(B)

如果A的一行的倍数加到A的另一行形成B, 那么det(A)=det(B)

这两个定理为计算行列式提供了一种有效的方法。注意,定理8描述了初等行运算对det(A)的影响。

我们可以快速准确地计算初等行变换的结果。所以对于含有已知常数项的矩阵,行化成三角矩阵是计算行列式的首选方法。对于有参数项的矩阵,通常使用其他方法。

由于矩阵乘法和行列式的定义复杂,乘积的行列式和行列式的乘积之间存在着一种简单的关系:

定理(9)

如果A和B是方阵,det(AB)=det(A)det(B)

如果det(A)=0,那么A一定是奇异的。反之亦然。

定理(10)

det(A)=0是A是奇异的一个充要条件。

定理(11)

对于每个方阵A,det(A)=det(A^T)

余子式与代数余子式

a_ij的余子式是,去掉A的第i行和第j列形成的矩阵的行列式。

a_ij的代数余子式写成A_ij,等于余子式乘以 (-1)^(i+j)。代数余子式的重要性是由于以下的重要定理:

定理(12)

对于每个i和j,

线性独立(线性无关)

方程组Ax=0可以有无穷多个解。为了描述这种系统的所有解的集合,我们必须首先理解线性无关的概念。

假设已知k个向量a_1, a_2,…,a_k和k个标量c_1 ,c_2…,c_k。考虑到表达式

不是全部为零如果上面的方程对某些标量成立(不是全部为零),那么向量a_1,a_2,…,a_k就是线性相关的,标量c_1,c_2,…,c_k就叫作权值。由上式可知

其中 c_1≠0。上面的方程表明a_1是其他向量的“加权和”。

如果一个给定的向量组不是线性相关的,那么它称为线性无关的。由于线性无关的集合不可能存在依赖关系,

意味着所有的标量系数必须是零。

一般来说,我们不能轻易判定一个向量组是否是线性无关的,也不能轻易求出权值(如果向量组是线性相关的)。

但我们可以把上式重写为矩阵向量的形式,为此,定义一个矩阵A,它的列是向量a_1,a_2,…,a_k,它的项的权值是c_1,c_2,…,a_k。因此:

根据矩阵乘法的定义:

因此,当且仅当Ac=0存在非零解时,矩阵A的列向量才是线性无关的。那么,c的元素就是权值。

如果矩阵A是方阵,那么基于线性相关的行列式的判据是可能和方便的。基本的思想是:如果A是可逆的,那么推论1,系统Ac+0只有平凡解,因为

证明了c一定是0。下面的定理将det(A)与A的行和列的线性无关联系起来。

定理(13)

当且仅当det(A)≠0时,n×n矩阵A的行(列)是线性无关的。

18 阅读:8668
评论列表
  • 2022-04-02 21:14

    线代24年前学的,都快忘光了[笑着哭][笑着哭]

  • onki 4
    2022-05-22 22:47

    非常详细,也比较基础,但是很多已经快忘光了。大部分都是线性代数中的知识。另外,文中有些小错误,希望作者进行校对。比如复数的极坐标表示中,第一个sin应为cos。

  • 2022-05-27 17:13

    豆豆那个样一样,是命,好命位置不一样罢了

  • 2022-05-30 10:48

    有点小错误A(A+C)[笑着哭]

  • 2022-04-22 16:20

    北大博士都去当城管了,微积分能治理小贩?

    火焰山 回复:
    积一下分,也许能算出谁最不遵守规定。
    qw422502 回复:
    哪里北大博士就一二本大学生[呲牙笑]
  • 2022-04-03 12:40

    [点赞]

  • 2022-04-22 12:14

    我大抵是飘了,这样的文章我也敢点进来

  • 2022-05-27 17:12

    那么喜欢听,,人类的命,不是那么回事吗

  • 2022-06-06 17:09

    多讲点微积分现实中的应用。[点赞][点赞][点赞]

  • 2022-04-20 22:13

    线性代数都忘光了[呲牙笑]

  • 2022-05-11 11:08

    从物理或生命科学角度:线性相关与空间可数对等。

  • 2022-05-07 06:29

    昨晚梦中在解数学题看到一堆不认识的数学公式。也有sin,cos

    walksing天一 回复:
    一个月后,从上海到南昌到广州东的路上,成功破解了幻数定理,并完成加密解密可迭代超级商城令牌无登录技术应用。高效率高压缩加密算法为超算领域核心基石。
  • [得瑟]说实话,的确是200300年前的东西,但是大学学习几十个小时,别人是花几百个小时长期从事的事情,所以学不会,一些就忘是正常的。。。恰恰高中之前又不接触,假货们就开始各种卷,骗人我优秀[笑着哭]

  • 2022-05-01 01:21

    太陌生了!

  • 2022-05-08 23:24

    [得瑟]又回忆了一遍

  • 2022-05-25 13:02

    我们的老祖宗花大量的精力是去研究他人,总是琢磨如融入别人构建的体系之中,这是生存之道。对于研究自然宇宙规律,我们的老祖宗淡淡一笑:你整这个玩意有什么用呢?能当饭吃吗?你是不是吃饱饭撑的?![吐舌头眯眼睛笑]

  • 2022-05-12 23:42

    字母和数字都认识,放在一起就不认识了

康托的天堂

简介:科学如此美妙,我想让你知道