axjack's blog

### axjack is said to be an abbreviation for An eXistent JApanese Cool Klutz ###

線形代数

列空間と左零空間は直交する

準備 行列Aを Aの列空間を Aの左零空間を とする。この時、C(A)の任意のベクトルとN(A')の任意のベクトルは直交する。 確認 x ∈ C(A), y ∈ N(A')を任意に取る。ベクトルが直交することを示すには内積が0となることを確認すれば良い。なお、小文字のoを零ベク…

統計学実践ワークブック 第16章 重回帰分析 pp.125-127について、式を導出したり図示したり行間を埋めたりしました。

試験対策だけであったら、ぶっちゃけ重回帰分析した結果(回帰係数、t値、p値、決定係数)が読み解ければ良いと思われる。がしかし、ワークブックpp.125-127にある線形代数的な解釈・幾何学的な解釈の魅力といったらなんとやら。というわけで、行間を埋めた記…

多変量正規分布の条件付き期待値・条件付き分散

やること 『統計学』東京図書 p.147 練習問題 問5.1の[2]、多変量正規分布の条件付き期待値・条件付き分散を解きます。 問題 が平均ベクトル、分散共分散行列が正値対称行列 の時、X = xおよびZ = zを与えた下でのYの 条件付き期待値 条件付き分散 を求めよ…

AB = E ならば BA = Eをランクを用いて示す

チャート式シリーズ 大学教養 線形代数を買ってから、行列の構造、特に「ランク」の理解が深まりました。ということでタイトルの命題の証明です。 AB = E ならば BA = E AもBもn次正方行列、Eはn次単位行列とします。 ここで、rankA = r ≦ n , rankB = s ≦ n…

行簡約行列をRで

pracmaのrrefを使って行簡約行列を出してみます。これで適当な行列を手計算で簡約化して答え合わせができますね。 ソース # install.packages('pracma') # library('pracma') ## 行ベクトルを4本 a1 <- c(1,2,0) a2 <- c(2,4,0) a3 <- c(0,1,3) a4 <- c(1,3,…

カーネルは線型空間

示したいこと 線型空間V, W に対して線型写像 が与えられているとする。 この時、カーネルは線型空間であることを示す。 証明 およびをそれぞれ取る。 すると、より、 となるので、Ker T の任意の元はベクトル和とスカラー倍に関して閉じていることが分かっ…

対角化

示したいこと 正方行列Aに対してR‪⁻¹‬AR=Λと対角化できるような正則行列Rが存在すると仮定する。この時、ΛはAの固有値を並べた対角行列であることを示そう。即ちΛの対角成分λi はAu=auを満たす固有値aであることを示せばよい。ここでu≠0とする。 証明 Au=au…

axjack is said to be an abbreviation for An eXistent JApanese Cool Klutz.