【数理知识】《矩阵论》方保镕老师-第5章-矩阵微积分及其应用
矩陣微積分及其應(yīng)用
- 5.1 向量序列和矩陣序列的極限
- 5.1.1 向量序列的極限
- 5.1.2 矩陣序列的極限
- 5.2 矩陣級數(shù)與矩陣函數(shù)
- 5.2.1 矩陣級數(shù)
- 5.2.2 矩陣函數(shù)
- 5.3 函數(shù)矩陣的微分和積分
- 5.3.1 函數(shù)矩陣對實(shí)變量的導(dǎo)數(shù)
- 定義 5.3.2 函數(shù)矩陣的導(dǎo)數(shù)定義
- 5.3.2 函數(shù)矩陣特殊的導(dǎo)數(shù)
- 1. 數(shù)量函數(shù)對于向量的導(dǎo)數(shù)
- 定義 5.3.3 數(shù)量函數(shù) f(x) 對 x 的導(dǎo)數(shù)
- 例 5.3.4
- 例 5.3.5
- 定義 5.3.4 數(shù)量函數(shù) f(x) 對 A 的導(dǎo)數(shù)
- 例 5.3.6
- 2. 矩陣對于矩陣的導(dǎo)數(shù)
- 5.3.3 矩陣的全微分
- 5.3.4 函數(shù)矩陣的積分
- 5.4 矩陣微分方程
- 5.4.1 常系數(shù)齊次線性微分方程組的解
- 定理 5.4.1
- 定理 5.4.2
- 5.4.2 常系數(shù)非齊次線性微分方程組的解
- 5.4.3 n 階常系數(shù)微分方程的解
5.1 向量序列和矩陣序列的極限
5.1.1 向量序列的極限
5.1.2 矩陣序列的極限
5.2 矩陣級數(shù)與矩陣函數(shù)
5.2.1 矩陣級數(shù)
5.2.2 矩陣函數(shù)
5.3 函數(shù)矩陣的微分和積分
5.3.1 函數(shù)矩陣對實(shí)變量的導(dǎo)數(shù)
定義 5.3.2 函數(shù)矩陣的導(dǎo)數(shù)定義
設(shè) A(t)=(aij(t))m×nA(t) = (a_{ij}(t))_{m\times n}A(t)=(aij?(t))m×n?,若 aij(t)(i=1,2,?,m;j=1,2,?,n)a_{ij}(t)(i=1,2,\cdots,m;j=1,2,\cdots,n)aij?(t)(i=1,2,?,m;j=1,2,?,n) 在 t=t0t=t_0t=t0? 處(或[a, b]上)可導(dǎo),則稱 A(t)A(t)A(t) 在 t=t0t=t_0t=t0? 處(或[a, b]上)可導(dǎo),記為
A′(t0)=dA(t)dt∣t=t0=lim?Δt→0A(t0+Δt)?A(t0)Δt=[a11′(t0)a12′(t0)…a1n′(t0)a21′(t0)a22′(t0)…a2n′(t0)????am1′(t0)am2′(t0)…amn′(t0)]\begin{aligned} A'(t_0) &= \frac{dA(t)}{dt}|_{t=t_0} = \lim_{\Delta t \rightarrow 0} \frac{A(t_0+\Delta t) - A(t_0)}{\Delta t} \\ &= \left[\begin{matrix} a'_{11}(t_0) & a'_{12}(t_0) & \dots & a'_{1n}(t_0)\\ a'_{21}(t_0) & a'_{22}(t_0) & \dots & a'_{2n}(t_0)\\ \vdots & \vdots & \ddots & \vdots \\ a'_{m1}(t_0) & a'_{m2}(t_0) & \dots & a'_{mn}(t_0)\\ \end{matrix}\right] \end{aligned}A′(t0?)?=dtdA(t)?∣t=t0??=Δt→0lim?ΔtA(t0?+Δt)?A(t0?)?=??????a11′?(t0?)a21′?(t0?)?am1′?(t0?)?a12′?(t0?)a22′?(t0?)?am2′?(t0?)?……?…?a1n′?(t0?)a2n′?(t0?)?amn′?(t0?)????????
不難證明函數(shù)矩陣的導(dǎo)數(shù)運(yùn)算有下列性質(zhì):
(1)A(t)A(t)A(t) 為常數(shù)矩陣的充要條件是 A′(t)=0A'(t) = \mathbf{0}A′(t)=0;
(2)設(shè) A(t)=(aij(t))m×nA(t) = (a_{ij}(t))_{m\times n}A(t)=(aij?(t))m×n? 和 B(t)=(bij(t))m×nB(t) = (b_{ij}(t))_{m\times n}B(t)=(bij?(t))m×n? 可導(dǎo),則dd(t)(A(t)±B(t))=A′(t)±B′(t)(5.3.1)\fracze8trgl8bvbq{d(t)}(A(t) \pm B(t)) = A'(t) \pm B'(t) \tag{5.3.1}d(t)d?(A(t)±B(t))=A′(t)±B′(t)(5.3.1)
(3)若 k(t)k(t)k(t) 是可導(dǎo)的實(shí)函數(shù),A(t)A(t)A(t) 可導(dǎo),則
dd(t)(k(t)A(t))=k′(t)A(t)+k(t)A′(t)(5.3.2)\fracze8trgl8bvbq{d(t)}(k(t)A(t)) = k'(t)A(t) + k(t)A'(t) \tag{5.3.2}d(t)d?(k(t)A(t))=k′(t)A(t)+k(t)A′(t)(5.3.2)
(4)設(shè) A(t)A(t)A(t) 和 B(t)B(t)B(t) 都可導(dǎo),則
dd(t)(A(t)B(t))=A′(t)B(t)+A(t)B′(t)(5.3.3)\fracze8trgl8bvbq{d(t)}(A(t)B(t)) = A'(t)B(t) + A(t)B'(t) \tag{5.3.3}d(t)d?(A(t)B(t))=A′(t)B(t)+A(t)B′(t)(5.3.3)
(5)若 A(t)A(t)A(t) 與 A?1(t)A^{-1}(t)A?1(t) 都有導(dǎo)數(shù),則
dA?1(t)dt=?A?1(t)A′(t)A?1(t)(5.3.4)\frac{dA^{-1}(t)}{dt} = -A^{-1}(t)A'(t)A^{-1}(t) \tag{5.3.4}dtdA?1(t)?=?A?1(t)A′(t)A?1(t)(5.3.4)
(6)設(shè)函數(shù)矩陣 A(t)A(t)A(t) 是 ttt 的函數(shù),而 t=f(x)t=f(x)t=f(x) 是 xxx 的實(shí)值函數(shù),且 A(t)A(t)A(t) 與 f(x)f(x)f(x) 均可導(dǎo),則有
dA(t)dx=dA(t)dtf′(x)=f′(x)dA(t)dt(5.3.5)\frac{dA(t)}{dx} = \frac{dA(t)}{dt}f'(x) = f'(x)\frac{dA(t)}{dt} \tag{5.3.5}dxdA(t)?=dtdA(t)?f′(x)=f′(x)dtdA(t)?(5.3.5)
補(bǔ)充性質(zhì):不論 AAA 是任何常量方陣,總有
(1)ddteAt=AeAt=eAtA(5.3.6)\fracze8trgl8bvbq{dt}e^{At} = Ae^{At}=e^{At}A \tag{5.3.6}dtd?eAt=AeAt=eAtA(5.3.6)
此條性質(zhì)可聯(lián)想 (ex)′=ex(e^x)' = e^x(ex)′=ex 記憶
(2)ddtcos?At=?A(sin?At)=?(sin?At)A(5.3.7)\fracze8trgl8bvbq{dt}\cos At = -A(\sin At) = -(\sin At)A \tag{5.3.7}dtd?cosAt=?A(sinAt)=?(sinAt)A(5.3.7)
此條性質(zhì)可聯(lián)想 (cos?x)′=?sin?x(\cos x)' = -\sin x(cosx)′=?sinx 記憶
(3)ddtsin?At=A(cos?At)=(cos?At)A(5.3.8)\fracze8trgl8bvbq{dt}\sin At = A(\cos At) = (\cos At)A \tag{5.3.8}dtd?sinAt=A(cosAt)=(cosAt)A(5.3.8)
此條性質(zhì)可聯(lián)想 (sin?x)′=cos?x(\sin x)' = \cos x(sinx)′=cosx 記憶
5.3.2 函數(shù)矩陣特殊的導(dǎo)數(shù)
1. 數(shù)量函數(shù)對于向量的導(dǎo)數(shù)
在場論中,我們對數(shù)量函數(shù) f(x,y,z)f(x,y,z)f(x,y,z) 定義梯度為
gradf=Δf=(?f?x,?f?y,?f?z)\mathrm{grad}\ f = \Delta f = (\frac{\partial f}{\partial x},\frac{\partial f}{\partial y},\frac{\partial f}{\partial z})grad?f=Δf=(?x?f?,?y?f?,?z?f?)
這可以理解為數(shù)量函數(shù) f(x,y,z)f(x,y,z)f(x,y,z) 對向量 (x,y,z)(x,y,z)(x,y,z) 的導(dǎo)數(shù)。
定義 5.3.3 數(shù)量函數(shù) f(x) 對 x 的導(dǎo)數(shù)
設(shè) x=(x1,x2,?,xn)Tx=(x_1,x_2,\cdots,x_n)^Tx=(x1?,x2?,?,xn?)T,f(x)=f(x1,x2,?,xn)f(x) = f(x_1,x_2,\cdots,x_n)f(x)=f(x1?,x2?,?,xn?) 是以向量 xxx 為自變量的數(shù)量函數(shù),即為 nnn 元函數(shù),則規(guī)定數(shù)量函數(shù) f(x)f(x)f(x) 對于向量 xxx 的導(dǎo)數(shù)為
dfdx=(?f?x1,?f?x2,?,?f?xn)T(5.3.10)\frac{df}{dx} = (\frac{\partial f}{\partial x_1},\frac{\partial f}{\partial x_2},\cdots,\frac{\partial f}{\partial x_n})^T \tag{5.3.10}dxdf?=(?x1??f?,?x2??f?,?,?xn??f?)T(5.3.10)
例 5.3.4
數(shù)量函數(shù) f(x)=xTAxf(x)=x^T A xf(x)=xTAx 對于向量 xxx 的導(dǎo)數(shù)為
dfdt=(A+AT)x=Ax+ATx(5.3.11)\frac{df}{dt} = (A+A^T)x = Ax + A^T x \tag{5.3.11}dtdf?=(A+AT)x=Ax+ATx(5.3.11)
其中, A(t)=(aij)n×nA(t) = (a_{ij})_{n\times n}A(t)=(aij?)n×n? 為常量矩陣,x=(x1,x2,?,xn)x=(x_1,x_2,\cdots,x_n)x=(x1?,x2?,?,xn?)。
有特例如下:
(1)當(dāng) AAA 是實(shí)對稱矩陣時,二次型 xTAxx^TAxxTAx 對 xxx 的導(dǎo)數(shù)為
dxTAxdx=2Ax(5.3.12)\frac{dx^TAx}{dx} = 2Ax \tag{5.3.12}dxdxTAx?=2Ax(5.3.12)
(2)當(dāng) A=IA=IA=I 時,函數(shù) f(x)=x12+x22+?+xn2f(x) = x_1^2 + x_2^2 + \cdots +x_n^2f(x)=x12?+x22?+?+xn2? 對 xxx 的導(dǎo)數(shù)為
df(x)dx=2x(5.3.13)\frac{df(x)}{dx} = 2x \tag{5.3.13}dxdf(x)?=2x(5.3.13)
例 5.3.5
令 x=(ξ1(t),ξ2(t),?,ξn(t))Tx=(\xi_1(t), \xi_2(t), \cdots, \xi_n(t))^Tx=(ξ1?(t),ξ2?(t),?,ξn?(t))T,f(x)=(ξ1,?,ξn)f(x) = (\xi_1, \cdots, \xi_n)f(x)=(ξ1?,?,ξn?)。則有
dfdt=(dfdf)Tdxdt(5.3.14)\frac{df}{dt} = (\frac{df}{df})^T\frac{dx}{dt} \tag{5.3.14}dtdf?=(dfdf?)Tdtdx?(5.3.14)
定義 5.3.4 數(shù)量函數(shù) f(x) 對 A 的導(dǎo)數(shù)
設(shè) A∈Rm×nA\in\mathbb{R}^{m\times n}A∈Rm×n,f(A)f(A)f(A) 為矩陣 AAA 的數(shù)量函數(shù),即看成是 m×nm\times nm×n 元函數(shù),則規(guī)定數(shù)量函數(shù) f(A)f(A)f(A) 對于矩陣 AAA 的導(dǎo)數(shù)為
dfdA=(?f?aij)m×n=[?f?a11??f?a1n????f?an1??f?ann](5.3.15)\frac{df}{dA} = (\frac{\partial f}{\partial a_{ij}})_{m\times n} = \left[\begin{matrix} \frac{\partial f}{\partial a_{11}} & \cdots & \frac{\partial f}{\partial a_{1n}} \\ \vdots & \ddots & \vdots \\ \frac{\partial f}{\partial a_{n1}} & \cdots & \frac{\partial f}{\partial a_{nn}} \\ \end{matrix}\right] \tag{5.3.15}dAdf?=(?aij??f?)m×n?=?????a11??f???an1??f???????a1n??f???ann??f??????(5.3.15)
例 5.3.6
二次型 xTAxx^T A xxTAx 對矩陣 AAA 的導(dǎo)數(shù),其中 A∈Rn×nA\in\mathbb{R}^{n\times n}A∈Rn×n 是實(shí)對稱的,那么有
ddA(xTAx)=(?f?aij∑i=1n∑j=1naijxixj)=(xixj)n×n=xxT\fracze8trgl8bvbq{dA}(x^TAx) = (\frac{\partial f}{\partial a_{ij}} \sum_{i=1}^{n} \sum_{j=1}^{n} a_{ij}x_ix_j) = (x_i x_j)_{n\times n} = xx^TdAd?(xTAx)=(?aij??f?i=1∑n?j=1∑n?aij?xi?xj?)=(xi?xj?)n×n?=xxT
2. 矩陣對于矩陣的導(dǎo)數(shù)
5.3.3 矩陣的全微分
5.3.4 函數(shù)矩陣的積分
5.4 矩陣微分方程
5.4.1 常系數(shù)齊次線性微分方程組的解
定理 5.4.1
一階線性常系數(shù)微分方程組的定解問題
{dxdt=Ax(t)x(0)=(x1(0),x2(0),?,xn(0))T(5.4.3)\left\{\begin{aligned} &\frac{dx}{dt} = Ax(t) \\ &x(0) = (x_1(0),x_2(0),\cdots,x_n(0))^T \end{aligned}\right.\tag{5.4.3}?????dtdx?=Ax(t)x(0)=(x1?(0),x2?(0),?,xn?(0))T?(5.4.3)
有唯一解 x=eAtx(0)x=e^{At}x(0)x=eAtx(0)
同理有
{dxdt=Ax(t)x∣t=t0=x(t0)(5.4.4)\left\{\begin{aligned} &\frac{dx}{dt} = Ax(t) \\ &x|_{t=t_0} = x(t_0) \end{aligned}\right.\tag{5.4.4}?????dtdx?=Ax(t)x∣t=t0??=x(t0?)?(5.4.4)
有唯一解 x=eA(t?t0)x(t0)x=e^{A(t-t_0)}x(t_0)x=eA(t?t0?)x(t0?)
定理 5.4.2
{dXdt=AX(t)X(t)∣t=t0=X(t0)(5.4.5)\left\{\begin{aligned} &\frac{dX}{dt} = AX(t) \\ &X(t)|_{t=t_0} = X(t_0) \end{aligned}\right.\tag{5.4.5}?????dtdX?=AX(t)X(t)∣t=t0??=X(t0?)?(5.4.5)
若未知函數(shù) X(t)X(t)X(t) 不是列向量,而是 n×mn\times mn×m 矩陣,X(t0)X(t_0)X(t0?) 是 n×mn\times mn×m 常數(shù)矩陣,AAA 是給定的 nnn 階常數(shù)方陣,則有唯一解為
X(t)=eA(t?t0)X(t0)X(t) = e^{A(t-t_0)}X(t_0)X(t)=eA(t?t0?)X(t0?)
且解 X(t)X(t)X(t) 的秩與 ttt 的取值無關(guān)。
5.4.2 常系數(shù)非齊次線性微分方程組的解
{dxdt=Ax(t)+f(t)x∣t=t0=x(t0)(5.4.6)\left\{\begin{aligned} &\frac{dx}{dt} = Ax(t) + f(t) \\ &x|_{t=t_0} = x(t_0) \end{aligned}\right.\tag{5.4.6}?????dtdx?=Ax(t)+f(t)x∣t=t0??=x(t0?)?(5.4.6)
有唯一解 x=eA(t?t0)x(t0)+∫t0teA(t?τ)f(τ)dτx=e^{A(t-t_0)}x(t_0) + \int_{t_0}^{t} e^{A(t-\tau)} f(\tau) d\taux=eA(t?t0?)x(t0?)+∫t0?t?eA(t?τ)f(τ)dτ
5.4.3 n 階常系數(shù)微分方程的解
總結(jié)
以上是生活随笔為你收集整理的【数理知识】《矩阵论》方保镕老师-第5章-矩阵微积分及其应用的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 【Paper】2003_Consensu
- 下一篇: 【数理知识】《矩阵论》方保镕老师-第6章