微分で行列式を求める!
ぼーっとしてたら思いついたことを書きます。次のような行列の行列式を求めよ、なんて問題はその辺の線形代数の教科書に載ってたりします。は対角成分が, 非対角成分が, サイズの行列です。
教科書に載ってるぐらいなので別に難しくありません。行・列基本変形を少し行うだけで簡単に求まります。が、気持ち悪い別解を思いついたので紹介しようという魂胆です。
以下のようなについての関数を定義します。は定数だとしましょう。
先のの行列式はのときのとして求まるのでの正体が分かってしまえば勝ちです。
さて、をどうやって求めましょうか...タイトルでネタバレしてるんですよね()。というわけで微分してみまーす。行列式の微分は各列を微分したときの行列式たちの和で書けます。
はい、長い旅でしたがこれでについての関係式が得られました。
とりあえずは順番に積分していけば求まりそうです。ただ積分定数よく分からんし〜〜〜ということでもう少し詰めていきましょう。一旦落ち着こうってやつです。
まずはその定義からについての次式です。また全成分が等しい行列の行列式はなのでです。よって,はたちの線形結合で書けます。さらにのときは簡単に計算できて
となります。以上の話をまとめると、
ここまでくればを出すのは簡単です。実際、
と順番に求まり、一般に
であることがわかります。よって、冒頭に出てきたについて、
であることが分かりました。
BCとCBの固有値の関係
「特異値分解」ってネットで調べると証明っぽいものが出てきたりします。ぽいものが。その中でありがちな議論は「の固有値はの固有値でもある」、また「の固有値はの固有値でもある」ので「との固有値は一致する」みたいなやつです。
え、重複度について何も議論してなくない??????というわけでそのあたりについて適当に勉強したことを適当に紹介します。特異値分解の証明はいろいろありますが日本語でヒットするものは高確率でガバい気がします。線形代数、そんな簡単ではなさそう。
話がそれたので戻ります。および に対して以下が成り立ちます。
ここででの固有多項式を表しました。この式からとの固有値は重複度を含めて一致することがわかります。(サイズ的に足りない分は0で埋める。)
証明は(読むのは)とっても簡単で
と
が実は相似になっていることを示しておしまいです。
この定理を使えば一番上で話した特異値分解の説明も厳密にできます。他にも応用例があって、コーシーの恒等式
を示すことができたり,成分がで与えられる行列の固有値を求めることができます。
ヒルベルト行列の正定値性
任意のでない多項式に対して,
ここまでは当たり前体操です。この左辺を変形すると
よって,以下のように第成分がで与えられるを用意すれば,
任意のについて,
です。これで先のように定義した行列は正定値であることが分かりました。実はには名前がついていてヒルベルト行列と言います。
こっから先は真面目に計算はしていない妄想です。今回はの表現の基底としてをとったため,に対する二次形式が登場しました。同様の議論での表現の基底としてルジャンドル多項式(厳密には区間をアフィン変換したやつ)を持ってくれば対角行列に対する二次形式が登場することがわかります。つまり,それらの基底の変換行列をとおけば
と、対角化っぽいことができます。ルジャンドル多項式はよく調べられているのでは簡単に求まり,上の表式からの逆行列も求まりそうです。詳しくは知りません(無責任)。とりあえず久保くんさんの試合を観ます。
とあるrank 等式の証明
一般化逆行列の説明でが列フルランクだったらは正則みたいな話を見てほんまか?ってなったのでメモです。一般に任意の行列について以下が成り立ちます。
(証明)
より,が従います。次元定理からも分かります。
標準形とかでごちゃごちゃやってもできるかもしれないですがわりと綺麗に示せました。
例えば行階数と列階数が等しいとか、自明っぽい話でもいざ証明しようとなると難しいものです。
恒等式の使い方
中1とかにするとウケる話をします。手品みたいなものです。
突然ですが2桁の整数を3つ思い浮かべてください。思い浮かべたのは ですね。(えー、そういう体でしばらく付き合ってください。これ当たってたらそれはもう手品を超えたなにかです。)
それに対して僕は という4つの整数を返します。「は?」と言うあなた、普通の反応です。「3つの整数に対して4つの整数返してくるのは失礼では?」というあなた、なんかごめんなさい...
冗談はさておき実は下の式が成り立ちます。確認してみてください。
...!?!?!?。すごくないですか?
他の例でも見てみましょう。あなたが と言えば僕は を返します。先と同じように
が成り立ちます。
以下ネタばらしです。考えたい人はまだ見ないでください。上の2つの例を眺めているとなんとなくカラクリは分かるかもしれません。
以上です。この式をじっと眺めていれば何をやっていたか分かるでしょう。恒等式、神々しい......
Q.最後のそれは洒落のつもりですか?
A.はい
log 2 の近似
本日三度目の投稿です。暇なので。
唐突ですが(以下底は自然対数とします)の値が必要になったことありませんか???僕は、ないです。が、中2のときに同級生にしたらほ〜〜んって言われた面白い話を書きます。
を代入して(収束の話はググってください、最近ググれば全部わかるという噂をよく聞くし)
この式自体、あるいはこの式の右辺の収束がまあそれは遅いということはよく知られています。
もっと収束の速い式を作ってみましょう。先のマクローリン展開およびそこでをに置き換えたものを並べると、
辺々引き算して
を代入すると...
「...!!!この右辺の収束普通に早くない?え、ちょっと式いじっただけじゃん!?!?すごい!!!」という反応をしてください、ほ〜〜んとか言わないの。
とにかく簡単な式変形で有意義なことができました。一般に
は全単射になるので,任意のに対して,とすれば
で、
と割と収束が早い感じに計算できます。