241 words
1 minute
工科高等代数4.7

f(x)=j=1ncjϕj(x)f(x)=\sum_{j=1}^{n}c_j\phi_j(x) {xi,yi}i=1mf{cj}j=1n\{x_i,y_i\}_{i=1}^{m}\to f\to \{c_j\}_{j=1}^{n} f(xi)yif(x_i)\to y_i j=1ncjϕj(xi)yi\sum_{j=1}^{n}c_j\phi_j(x_i)\to y_i

[ϕ1(xi)ϕ2(xi)ϕn(xi)][c1c2cn]=yi\begin{align} \begin{bmatrix} \phi_1(x_i)&\phi_2(x_i)&\dots&\phi_n(x_i)\\ \end{bmatrix} \begin{bmatrix} c_1\\ c_2\\ \dots\\ c_n\\ \end{bmatrix} =y_i \end{align}

y=f(x)+εy=f(x)+\varepsilon i=1myif(xi)2=YAc22\sum_{i=1}^{m}|y_i-f(x_i)|^2=\|Y-Ac\|_2^2 mincRn12YAc22\min_{c\in R^n}\frac{1}{2}\|Y-Ac\|_{2}^{2} f(c)=12YAc22f(c)=\frac{1}{2}\|Y-Ac\|_2^2 f(x)=x12+x1x2f(x)=x_1^2+x_1x_2 fx1=2x1+x2\frac{\partial f}{\partial x_1}=2x_1+x_2 fx2=x1\frac{\partial f}{\partial x_2}=x_1

xf(x)=[2x1+x2x1]\begin{align} \nabla_{x}f(x)= \begin{bmatrix} 2x_1+x_2\\ x_1 \end{bmatrix} \end{align} A=[a11a12a21a22]f(A)=a112+a222+a11a21,AR2×2fa11=2a11+a21,fa12=0fa21=a11,fa22=2a22Af(x)=[2a11+a210a112a22]\begin{align} A= \begin{bmatrix} a_{11}&a_{12}\\ a_{21}&a_{22}\\ \end{bmatrix}\\ f(A)=a_{11}^2+a_{22}^2+a_{11}a_{21},A\in R^{2\times 2}\\ \frac{\partial f}{\partial a_{11}}=2a_{11}+a_{21},& \frac{\partial f}{\partial a_{12}}=0\\ \frac{\partial f}{\partial a_{21}}=a_{11},& \frac{\partial f}{\partial a_{22}}=2a_{22}\\ \nabla_{A}f(x)= \begin{bmatrix} 2a_{11}+a_{21}&0\\ a_{11}&2a_{22}\\ \end{bmatrix} \end{align}

f(X)=tr(AX),A,XRn×nf(X)=tr(AX),A,X\in R^{n\times n} f(X)=i=1nj=1naijxjif(X)=\sum_{i=1}^{n}\sum_{j=1}^{n}a_{ij}x_{ji} fxij=aji\frac{\partial f}{\partial x_{ij}}=a_{ji} Xf(X)=AT\nabla_{X}f(X)=A^{T} f(x)=xTATAx,xRn,ARm×nf(x)=x^TA^TAx,x\in R^n,A\in R^{m\times n}

Ax=[j=1na1jxjj=1na2jxjj=1namjxj]f(x)=i=1m(j=1naijxj)2fxj=\begin{align} Ax= \begin{bmatrix} \sum_{j=1}^{n}a_{1j}x_j\\ \sum_{j=1}^{n}a_{2j}x_j\\ \dots\\ \sum_{j=1}^{n}a_{mj}x_j\\ \end{bmatrix}\\ f(x)=\sum_{i=1}^{m}(\sum_{j=1}^{n}a_{ij}x_j)^2\\ \frac{\partial f}{\partial x_{j}}= \end{align} f(x)=i=1m(j=1naijxj)2=i=1m(xk2aik2+j=1,jkn2aijaikxjxk)fxk=i=1m(2aik2xk+2j=1,jknaijaikxj)=2i=1maikj=1naijxj\begin{align} f(x)&=\sum_{i=1}^{m}(\sum_{j=1}^{n}a_{ij}x_j)^2\\ &=\sum_{i=1}^{m}(x_{k}^2a_{ik}^2+\sum_{j=1,j\neq k}^{n}2a_{ij}a_{ik}x_jx_k)\\ \frac{\partial f}{\partial x_k}&=\sum_{i=1}^{m}(2a_{ik}^2x_k+2\sum_{j=1,j\neq k}^{n}a_{ij}a_{ik}x_j)\\ &=2\sum_{i=1}^{m}a_{ik}\sum_{j=1}^{n}a_{ij}x_j\\ \end{align}

B=ATAB=A^TA

x(xTBx)=(B+BT)x=2ATAx\begin{align} \nabla_{x}(x^TBx)=(B+B^T)x=2A^TAx \end{align} f(c)=12YAx22=12(YAc)T(YAc)=12(YTYYTAccTATY+cTATAc)\begin{align} f(c)&=\frac{1}{2}\|Y-Ax\|_2^2\\ &=\frac{1}{2}(Y-Ac)^T(Y-Ac)\\ &=\frac{1}{2}(Y^TY-Y^TAc-c^TA^TY+c^TA^TAc) \end{align} cf(c)=12(0ATYATY+2ATAc)=ATAcATYATAc=ATY\begin{align} \nabla_cf(c)&=\frac{1}{2}(0-A^TY-A^TY+2A^TAc)\\ &=A^TAc-A^TY\\ A^TAc&=A^TY\\ \end{align}

ARm×n,m>n,rank(A)=rA\in R^{m\times n},m>n,rank(A)=r Ac=YAc=Y 可能没有解,系数多 ATAc=ATYA^TAc=A^TY 好解,规模小 rank(A)=n,c=(ATA)1ATYrank(A)=n,c=(A^TA)^{-1}A^TY A=uΣvTA=u\Sigma v^T

c=(uΣvTvΣTuT)1(vΣTuT)Y=(uΣΣTuT)1vΣTuTY=AY\begin{align} c&=(u\Sigma v^Tv\Sigma^T u^T)^{-1}(v\Sigma^T u^T)Y\\ &=(u\Sigma\Sigma^T u^T)^{-1}v\Sigma^T u^TY\\ &=A^{\dagger}Y\\ \end{align}

κ(A)=σ1(A)σr(A)\kappa(A)=\frac{\sigma_{1}(A)}{\sigma_{r}(A)} κ(ATA)=σ1(ATA)σr(ATA)=σ12(ATA)σr2(ATA)=κ2(A)\kappa(A^TA)=\frac{\sigma_{1}(A^TA)}{\sigma_{r}(A^TA)}=\frac{\sigma_{1}^2(A^TA)}{\sigma_{r}^2(A^TA)}=\kappa^2(A) QRQR分解 A=QRA=QR RTQTQRc=RTQTYR^TQ^TQRc=R^TQ^TY RTRc=RTQTYR^TRc=R^TQ^TY κ(A)=κ(R)\kappa(A)=\kappa(R) 如果rank(A)=n,Rrank(A)=n,R可逆 Rc=QTY,c=R1QTYRc=Q^TY,c=R^{-1}Q^TY 样本PCAPCA 找到数据集的一个低维表示形式 x1,x2,,xNRdx_1,x_2,\dots,x_N\in R^d 考虑线性降维,令z1,z2,,zMz_1,z_2,\dots,z_M是数据集的MM哥线性表示形式 zm=j=1Ncjmxjz_m=\sum_{j=1}^Nc_{jm}x_j 其中c1m,c2m,,cNmc_{1m},c_{2m},\dots,c_{Nm}是常数,若标准化的线性组合中z1z_1具有最大的方差,则称为x1,x2,,xNx_1,x_2,\dots,x_N的第一主成分,对应的单位向量[c11,c21,,cN1]T[c_{11},c_{21},\dots,c_{N1}]^T称为第一主成分的载荷量。