矩阵函数推着玩

矩阵函数探秘:从幂等到幂单
本文深入探讨矩阵函数,包括幂等矩阵、幂单矩阵和幂零矩阵的特性。通过数学推导揭示了这些矩阵函数的公式,如幂等矩阵的函数公式f(sI+tA)=(I−A)f(s)+Af(s+t),幂单矩阵的函数公式f(sI+tA)=21[(I+A)f(s+t)+(I−A)f(s−t)],以及幂零矩阵的函数公式f(sI+tA)=If(s)+tAf′(s)。这些推导加深了对矩阵函数的理解。

用自己浅薄的数学知识推一下反正也用不上的矩阵函数公式,加深并不存在的理解。
矩阵函数定义:设复变幂级数∑m=0∞cmzm\sum_{m=0}^{\infty}{c_{m}z^{m}}m=0cmzm的收敛半径是R,且在收敛域内f(z)=∑m=0∞cmzmf(z)=\sum_{m=0}^{\infty}{c_{m}z^{m}}f(z)=m=0cmzm,当矩阵A的谱半径ρ(A)<R\rho(A)<Rρ(A)<R,定义f(A)=∑m=0∞cmAmf(A)=\sum_{m=0}^{\infty}{c_{m}A^{m}}f(A)=m=0cmAm,并称f(A)f(A)f(A)为矩阵A的函数 。

幂等矩阵函数:

f(sI+tA)=(I−A)f(s)+Af(s+t)f(sI+tA)=(I-A)f(s)+Af(s+t)f(sI+tA)=(IA)f(s)+Af(s+t)

满足条件:

  1. A(I−A)=(I−A)A=OA(I-A)=(I-A)A=OA(IA)=(IA)A=O
  2. An=A,(n=1,2,3,...)A^{n}=A,(n=1,2,3,...)An=A,(n=1,2,3,...)
  3. I−A为幂等矩阵I-A为幂等矩阵IA

则有
f(sI+tA)=∑m=0∞cm(sI+tA)m=∑m=0∞cm[s(I−A)+(s+t)A]m(拆,使用条件1)=∑m=0∞cm[sm(I−A)m+(s+t)mAm](使用条件2,3)=∑m=0∞cm[sm(I−A)+(s+t)mA]=(I−A)∑m=0∞cmsm+A∑m=0∞cm(s+t)m=(I−A)f(s)+Af(s+t) \begin{aligned} f(sI+tA) &=\sum_{m=0}^{\infty}{c_{m}(sI+tA)^{m}}\\ &=\sum_{m=0}^{\infty}{c_{m}[s(I-A)+(s+t)A]^{m}}(拆,使用条件1)\\ &=\sum_{m=0}^{\infty}{c_{m}[s^{m}(I-A)^{m}+(s+t)^{m}A^{m}]}(使用条件2,3)\\ &=\sum_{m=0}^{\infty}{c_{m}[s^{m}(I-A)+(s+t)^{m}A]}\\ &=(I-A)\sum_{m=0}^{\infty}{c_{m}s^{m}}+A\sum_{m=0}^{\infty}{c_{m}(s+t)^{m}}\\ &=(I-A)f(s)+Af(s+t) \end{aligned} f(sI+tA)=m=0cm(sI+tA)m=m=0cm[s(IA)+(s+t)A]m使1=m=0cm[sm(IA)m+(s+t)mAm]使2,3=m=0cm[sm(IA)+(s+t)mA]=(IA)m=0cmsm+Am=0cm(s+t)m=(IA)f(s)+Af(s+t)

幂单矩阵函数:

f(sI+tA)=12[(I+A)f(s+t)+(I−A)f(s−t)]f(sI+tA)=\frac{1}{2}[(I+A)f(s+t)+(I-A)f(s-t)]f(sI+tA)=21[(I+A)f(s+t)+(IA)f(st)]

满足条件:

  1. A2=AA=IA^2=AA=IA2=AA=I,即I−A2=(I+A)(I−A)=OI-A^2=(I+A)(I-A)=OIA2=(I+A)(IA)=O
  2. 矩阵AAA为幂单矩阵⇔\Leftrightarrow12(A+I)\frac{1}{2}(A+I)21(A+I)为幂等矩阵

则有
f(sI+tA)=∑m=0∞cm(sI+tA)m=∑m=0∞cm[12((I+A)(s+t)+(I−A)(s−t))]m(拆,使用条件1)=∑m=0∞cm[((s+t)m(12(I+A))m+(s−t)m(12(I−A))m)](使用条件2,幂等条件2)=∑m=0∞cm[((s+t)m(12(I+A))+(s−t)m(12(I−A)))]=12(I+A)∑m=0∞cm(s+t)m+12(I−A)∑m=0∞cm(s−t)m=12[(I+A)f(s+t)+(I−A)f(s−t)] \begin{aligned} f(sI+tA) &=\sum_{m=0}^{\infty}{c_{m}(sI+tA)^{m}}\\ &=\sum_{m=0}^{\infty}{c_{m}[\frac{1}{2}((I+A)(s+t)+(I-A)(s-t))]^{m}}(拆,使用条件1)\\ &=\sum_{m=0}^{\infty}{c_{m}[((s+t)^{m}(\frac{1}{2}(I+A))^{m}+(s-t)^{m}(\frac{1}{2}(I-A))^{m})]}(使用条件2,幂等条件2)\\ &=\sum_{m=0}^{\infty}{c_{m}[((s+t)^{m}(\frac{1}{2}(I+A))+(s-t)^{m}(\frac{1}{2}(I-A)))]}\\ &=\frac{1}{2}(I+A)\sum_{m=0}^{\infty}{c_{m}(s+t)^{m}}+\frac{1}{2}(I-A)\sum_{m=0}^{\infty}{c_{m}(s-t)^{m}}\\ &=\frac{1}{2}[(I+A)f(s+t)+(I-A)f(s-t)] \end{aligned} f(sI+tA)=m=0cm(sI+tA)m=m=0cm[21((I+A)(s+t)+(IA)(st))]m使1=m=0cm[((s+t)m(21(I+A))m+(st)m(21(IA))m)]使2,2=m=0cm[((s+t)m(21(I+A))+(st)m(21(IA)))]=21(I+A)m=0cm(s+t)m+21(IA)m=0cm(st)m=21[(I+A)f(s+t)+(IA)f(st)]

幂零矩阵函数:

f(sI+tA)=If(s)+tAf′(s)f(sI+tA)=If(s)+tAf'(s)f(sI+tA)=If(s)+tAf(s)

满足条件:

  1. A2=AA=OA^2=AA=OA2=AA=O

则有
f(sI+tA)=∑m=0∞cm(sI+tA)m(拆,使用条件1)=∑m=0∞cm[smI+msm−1tA]=I∑m=0∞cmsm+tA∑m=0∞cmmsm−1=I∑m=0∞cmsm+tA(∑m=0∞cmzm)′∣z=s=If(s)+tAf′(s) \begin{aligned} f(sI+tA) &=\sum_{m=0}^{\infty}{c_{m}(sI+tA)^{m}}(拆,使用条件1)\\ &=\sum_{m=0}^{\infty}{c_{m}[s^{m}I+ms^{m-1}tA]}\\ &=I\sum_{m=0}^{\infty}{c_{m}s^{m}}+tA\sum_{m=0}^{\infty}{c_{m}ms^{m-1}}\\ &=I\sum_{m=0}^{\infty}{c_{m}s^{m}}+tA(\sum_{m=0}^{\infty}{c_{m}z^{m}})'|_{z=s}\\ &=If(s)+tAf'(s) \end{aligned} f(sI+tA)=m=0cm(sI+tA)m使1=m=0cm[smI+msm1tA]=Im=0cmsm+tAm=0cmmsm1=Im=0cmsm+tA(m=0cmzm)z=s=If(s)+tAf(s)

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值