• 締切済み

行列計算の微分を教えてください。

D=A・{(C^t・C)^-1・C^t・A}^t・[{(C^t・C)^-1・C^t・A}・{(C^t・C)^-1・C^t・A}^t]^-1 ↑の行列式を下のようにそれぞれ偏微分したらどうなるでしょうか。 ∂D/∂A ∂D/∂C 転置と逆行列が入り混じっており、計算が分からなくなってしまいました。よろしくお願いします。

みんなの回答

  • jcpmutura
  • ベストアンサー率84% (311/366)
回答No.1

B=(C^t・C)^-1・C^t・A とすると D=A・B^t・(B・B^t)^-1 B=(C^t・C)^-1・C^t・A ↓(C^t・C)^-1=C^{-1}・(C^t)^{-1} だから B=C^{-1}・(C^t)^{-1}・C^t・A ↓(C^t)^{-1}・C^t=E(単位行列)だから B=C^{-1}・E・A ↓E・A=A だから B=C^{-1}・A ↓ B^{-1}=(C^{-1}・A)^{-1} ↓ B^{-1}=A^{-1}・C D=A・B^t・(B・B^t)^-1 ↓(B・B^t)^-1=(B^t)^{-1}・B^{-1} だから D=A・B^t・(B^t)^{-1}・B^{-1} ↓B^t・(B^t)^{-1}=E だから D=A・E・B^{-1} ↓A・E=A だから D=A・B^{-1} ↓B^{-1}=A^{-1}・C だから D=A・A^{-1}・C ↓A・A^{-1}=E だから D=E・C ↓E・C=C だから D=C Dの行列式は |D|=|C| だから 行列DのAによる偏微分 ∂D/∂A=∂C/∂A=0 は行列式の偏微分ではありません Dの行列式のAによる偏微分は ∂|D|/∂A=∂|C|/∂A=0 です 行列DのCによる偏微分 ∂D/∂C=∂C/∂C は行列式の偏微分ではありません Dの行列式のCによる偏微分は ∂|D|/∂C=∂|C|/∂C です

全文を見る
すると、全ての回答が全文表示されます。

関連するQ&A

  • 行列を行列で微分することはできるのか?

    行列を1変数で微分することや1変数を行列で微分することは可能ですが行列を行列で微分することは可能なのでしょうか? たとえばA微分可能な行列として d(A^-1)/dA = -(A^1)^2 というような計算はありえるでしょうか? 以上よろしくお願いします.

  • 行列をベクトルで微分

    以下,大文字を行列,小文字をベクトルとします.また," ' "は転置を示します. ∂(b'a)/∂a=∂(a'b)/∂a=b ∂(a'Ba)/∂=2Ba という公式は知っているのですが,以下の微分ではどうなるかさっぱりわかりません. ∂(Ba'C)/∂a ∂(BaC)/∂a ∂(BaCa'D)/∂a 計算方法を教えていただきたいです.

  • 行列の連立微分方程式

    Aという行列に対して、dy/dt=Ay y(0)=(c0、c1、c2)^T【←転置】という連立微分方程式を解け。 ただしAは α 0 0 0 β 1 0 0 β という3×3の行列であるとする。 という問題なのですが、おそらくy=exp(tA)y(0)という形にして計算するのではないかと思うのですが、ここのexp(tA)を求める過程について質問です。 この問題を解く前に、 e^(αt) 0 0 0 e^(βt) te^(βt) 0 0 e^(βt) という3×3行列Xについて、 X^(-1)dX/dt=A という計算をさせられているので、おそらくこの式を利用して解くのだと思いますが、どのように利用したらいいのか見当がつきません。 exp(tA)=X になるのではないかという予想はあるのですが、どうしたらそうなるのか、過程等わかる方がいましたら、ご教授お願いします。

  • 行列の計算

    A=(6 2)    B=(1 1)    (2 6)      (1 0) があって、(AとBは2行2列の行列です。表記の仕方はそう思ってください) D=A^(-1)+A^(-1)×B×(I -B'×A^(-1)×B)^(-1) ×B'×A^(-1) を求めよ。という問題があったのですが、普通に計算していけば答えは出ます。ただ、もっと上手な解法があるのかが気になります。あるようであれば、教えていただけると幸いです。 (A^(-1) は Aの逆行列、B'は Bの転置行列、Iは単位行列です)

  • 逆行列の行列による微分

    今、あるシステムを構成したく行列方程式を解こうとしているのですが、行き詰ってしまったので質問させていただきます。 YとAを5×1行列、Bを5×3行列、Xをスカラ、Wを5×5行列、Hを3×1行列,B'はBの共役転置行列として、 Y = (A+BH)X・・・(1) H = -inv(B'WB)B'WA・・・(2) とします。 ((2)を(1)に代入した状態で)Y'Y=Jとおいて、∂J/∂W=0となるようなWを求めようとしているのですが、逆行列の行列による微分などが出てきてしまい、まったくわからなくなってしまいました。 これは解けるのでしょうか? 行列をvec関数を使ってベクトルになおしてみたりしたのですがinv()の部分がわからなくて・・・。 逆行列の行列による微分さえわかれば何とかなるような気もするのですが・・・。 解けるとすればどのようになるのでしょうか? 解法のヒントや参考になるページなどありましたら紹介していただきたいです。

  • 微分と行列

    微分と行列の問題なのですが解答がないので チェックをお願いします 次の関係からdy/dx,d^2y/dx^2を求めよ x=acos^3t y=asin^3t (a>0) dy/dt=3asin^2tcost dx/dt=-3acos^2tsint dy/dx=3asin^2tcost/-3acos^2tsint=-tant d^2y/dx^2=(d/t)(dy/dx)(dt/dx) =(d/t)(-tant)(1/-3acos^2tsint) =(-1/cos^2t)(-1/3acos^2tsint) =1/(3acos^4tsint) 余因数展開により次の行列式を求めよ D=|a^2 a 1 bcd|  |b^2 b 1 acd|   |c^2 c 1 abd|   |d^2 d 1 abc|   = |a^2 a 1 bcd | |(b+a)(b-a) b-a 0 -cd(b-a)| |(c+a)(c-a) c-a 0 -bd(c-a)| |(d+a)(d-a) d-a 0 -bc(d-a)| =(b-a)(c-a)(d-a)|(b+a) 1 -cd| |(c+a) 1 -bd| |(d+a) 1 -bc| =(b-a)(c-a)(d-a)|(b+a) 1 -cd | |(c-b) 0 d(c-b)| |(d-b) 0 c(d-b)| =-(b-a)(c-a)(d-a)(c-b)(d-b)|1 d|=(a-b)(c-a)(d-a)(c-b)(d-b)(c-b)      |1 c|   面倒ですがお願いします

  • 行列式の計算

       |C| [A B]| | という行列式はどのように計算すればよいのでしょうか。    |D| また、[A B][C D]はどのように計算すればよいのでしょうか。 すみませんが、お願いしますm(__)m

  • 行列の計算の進め方

    行列P(a,b,c,d)(順に左上,右上,左下,右下の成分です)が P^2-dP=P^(-1)…(1) および ad-bc=1…(2) を満たしている。このときa=-1を示せ。ただしa,b,c,dはいずれも実数である。 という問題なんですが、まず逆行列が嫌なので(1)に両辺Pをかけて逆行列を消し、またハミルトンケーリーで今作った式の次数下げをどんどんやっていきます。そうすると、 {a(a+d)-1}P=(a+1)E という式が出てくると思うのですが、ここで場合分けに入ると、 a(a+d)-1=0 の場合はa=-1と定まるのですが、 a(a+d)-1≠0 の場合に、PはEの実数倍にならないということが言えませんので((2)よりdetP=1)P=kEとしてハミルトンケーリーの式に代入することになると思います。ところ が代入した後に{k^2-(a+d)k+1}E=0となり、kの2次方程式が表れますが、これ以上進めません。 どこか見落としているだけかもしれませんが... 解答では素直に逆行列を求めてから成分を比較して、途中背理法を使うことによってa=-1と示していました。もちろん僕のやり方をするなら後者のほうも調べないといけませんから、このままでは点数になりません。 お聞きしたいことは2点あります。 一つ目は僕が進めなかったところから進めるのか。 二つ目はどういうときに成分計算に入ろうと考えるのか。 特に二つ目は、行列は成分計算が面倒になるからなるべく行列式でもちこみたいと教わっていたので、どうも成分計算に入る決意が出来ません。特に今回のようにもう少しで出来るのでないかと思うようなところまでくるとなおさら「いまさら成分なんか計算してたら時間が無くなってしまう!」と考えてしまいます。それは問題量や経験の差かもしれませんが、何かアドバイスあればよろしくお願いします!

  • 行列の計算がわかりません。

    行列の計算がわかりません。 お世話になっております。転置行列の入っている計算がわからなくて質問しました。 (y-H・x)^t×(y-H・x)=x^t・H^t・H・x-2x^t・H^t・y+y^t・y 答えの二項目がこうなる理由がわかりません。 式からするとy^t・H・x=x^t・H^t・yとなるのですよね?(H・x)^t・y=x^t・H^t・yとなるのはわかるのですが、y^t・H・x=x^t・H^t・yのように転置になる行列が変わっても等式が成り立つなんてことは知りませんでした。これは行列の性質として存在しているものなのでしょうか?調べた限りそのような性質はなかったと思います。それとも自分がなにか勘違いしているのでしょうか? 回答お願いいたします

  • 行列ベクトルの微分

    行列ベクトルの微分 A:(n×n) (1) ω(A)をAで微分   (ちなみにω(A)は( a(11) a(22) ・・・ a(nn) )の転置です) (2) (En(×)A)vec(A)をAで微分   (ちなみに(×)はクロネッカー積) この問題がわかる方、一方だけでもいいので教えてください。お願いします。