jRkjukpjpjpjhforj0jNj
fork0kNkRjkRjkujpk矩阵的QR分解
f以下是二次多项式求根doublerootdoublebdoublec
doublemmbb4cretur
m二次多项式求根
以下是求解矩阵的所有特征值voidcharacteristicdoubleANNdoublechaRNdoublechaIN
i
tk0mN1i
tiji
tLdoublestxdoubleMNNBNNi
tf0doubledchdoubleuNwNpNdoubleQNNRNNforL0LMAXL
extifm0chaR0A00chaI00break
iffabsAmm1E
chaRmAmmchaIm0mgoto
extsAm1m1AmmtAm1m1AmmAmm1Am1mifm1xrootstifxExsqrtxchaRms2x2chaRm1s2x2chaIm0chaIm10elsexsqrtfabsxchaRms2chaRm1s2chaImx2chaIm1x2break
fiffabsAm1m2E
xrootstifxExsqrtxchaRms2x2chaRm1s2x2chaIm0chaIm10elsexsqrtfabsxchaRms2chaRm1s2chaImx2chaIm1x2mm2goto
extfori0imiforj0jmj
ifijfork0Mij0kmkMijAikAkjMijMijMijsAijtelsefork0Mij0kmkMijAikAkjMijMijMijsAij以下是M的QR分解fori0imiforj0jmjifijQij1elseQij0fori0imiforj0jmjRijMijfori0imiforji1jmjifRjiEff1iffmico
ti
ueforjid0jmjddRjiRjidsqrtdc1sg
RiidhcccRiiforji1jmjujRjiforj0jijuj0uiRiicforj0jmjfork0wj0kmkwjQjkukwjforj0jmjfork0kmkQjkQjkwjukhforj0jmjforkipj0kmkpjRkjukpjpjpjhforj0jmj
ffork0kmkRjkRjkujpk
forj0jmj
fork0kmkMjkQjk
以上是M的QR分解fori0imi
forj0jmjfork0Bij0kmkBijMkiAkjBijfori0imiforj0jmjfork0Aij0kmkAijBikMkjAij
以下是求矩阵的所有特征值的特征向量voideige
vectordoubleVNNdoubleTN
doubleANNbaozNNguodNdoubleci
tijkmti
tW0fori0iNiforj0jNjbaozijVijfort0t6t
fori0iNiforj0jNjAijbaozij
fori0iNiAiiAiiTt
fori0iN1i

forjijNjiffabsAjiEkjbreak
forjijNjguodjAijAijAkjAkjguodj
forjijNj

cAji
iffabscE
formimNmAjmAjmc

forj0jNj

cAji
fifjiformimNmAjmAjmAimr