f(x)=∑j=1ncjϕj(x)
{xi,yi}i=1m→f→{cj}j=1n
f(xi)→yi
∑j=1ncjϕj(xi)→yi
[ϕ1(xi)ϕ2(xi)…ϕn(xi)]c1c2…cn=yi
y=f(x)+ε
∑i=1m∣yi−f(xi)∣2=∥Y−Ac∥22
minc∈Rn21∥Y−Ac∥22
f(c)=21∥Y−Ac∥22
f(x)=x12+x1x2
∂x1∂f=2x1+x2
∂x2∂f=x1
∇xf(x)=[2x1+x2x1]
A=[a11a21a12a22]f(A)=a112+a222+a11a21,A∈R2×2∂a11∂f=2a11+a21,∂a21∂f=a11,∇Af(x)=[2a11+a21a1102a22]∂a12∂f=0∂a22∂f=2a22
f(X)=tr(AX),A,X∈Rn×n
f(X)=∑i=1n∑j=1naijxji
∂xij∂f=aji
∇Xf(X)=AT
f(x)=xTATAx,x∈Rn,A∈Rm×n
Ax=∑j=1na1jxj∑j=1na2jxj…∑j=1namjxjf(x)=i=1∑m(j=1∑naijxj)2∂xj∂f=
f(x)∂xk∂f=i=1∑m(j=1∑naijxj)2=i=1∑m(xk2aik2+j=1,j=k∑n2aijaikxjxk)=i=1∑m(2aik2xk+2j=1,j=k∑naijaikxj)=2i=1∑maikj=1∑naijxj
B=ATA
∇x(xTBx)=(B+BT)x=2ATAx
f(c)=21∥Y−Ax∥22=21(Y−Ac)T(Y−Ac)=21(YTY−YTAc−cTATY+cTATAc)
∇cf(c)ATAc=21(0−ATY−ATY+2ATAc)=ATAc−ATY=ATY
A∈Rm×n,m>n,rank(A)=r
Ac=Y
可能没有解,系数多
ATAc=ATY
好解,规模小
rank(A)=n,c=(ATA)−1ATY
A=uΣvT
c=(uΣvTvΣTuT)−1(vΣTuT)Y=(uΣΣTuT)−1vΣTuTY=A†Y
κ(A)=σr(A)σ1(A)
κ(ATA)=σr(ATA)σ1(ATA)=σr2(ATA)σ12(ATA)=κ2(A)
QR分解
A=QR
RTQTQRc=RTQTY
RTRc=RTQTY
κ(A)=κ(R)
如果rank(A)=n,R可逆
Rc=QTY,c=R−1QTY
样本PCA
找到数据集的一个低维表示形式
x1,x2,…,xN∈Rd
考虑线性降维,令z1,z2,…,zM是数据集的M哥线性表示形式
zm=∑j=1Ncjmxj
其中c1m,c2m,…,cNm是常数,若标准化的线性组合中z1具有最大的方差,则称为x1,x2,…,xN的第一主成分,对应的单位向量[c11,c21,…,cN1]T称为第一主成分的载荷量。