今天看書,發(fā)現(xiàn)svd公式,在網(wǎng)上查閱了一下資料,主要用到如下的兩句:
主要是兩條語句:
from numpy import *;
U,Sigma,VT=linalg.svd([[1,1],[7,7]]);
創(chuàng)新互聯(lián)建站長期為上1000家客戶提供的網(wǎng)站建設(shè)服務(wù),團(tuán)隊從業(yè)經(jīng)驗10年,關(guān)注不同地域、不同群體,并針對不同對象提供差異化的產(chǎn)品和服務(wù);打造開放共贏平臺,與合作伙伴共同營造健康的互聯(lián)網(wǎng)生態(tài)環(huán)境。為龍井企業(yè)提供專業(yè)的網(wǎng)站建設(shè)、成都網(wǎng)站建設(shè),龍井網(wǎng)站改版等技術(shù)服務(wù)。擁有十余年豐富建站經(jīng)驗和眾多成功案例,為您定制開發(fā)。
用到的實例(李航《統(tǒng)計學(xué)習(xí)方法》第15章奇異值分解):
from numpy import *
A = array([[1,0,0,0],[0,0,0,4],[0,3,0,0],[0,0,0,0],[2,0,0,0]])
U,sigma,V = linalg.svd(A)
U
array([[ 0. , 0. , -0.4472136 , 0. , -0.89442719],
[-1. , 0. , 0. , 0. , 0. ],
[ 0. , -1. , 0. , 0. , 0. ],
[ 0. , 0. , 0. , 1. , 0. ],
[ 0. , 0. , -0.89442719, 0. , 0.4472136 ]])
sigma
array([ 4. , 3. , 2.23606798, -0. ])
V
array([[-0., -0., -0., -1.],
[-0., -1., -0., -0.],
[-1., -0., -0., -0.],
[-0., -0., -1., -0.]])