计算数学
pure小清新
学习进程记录与总结
展开
-
R语言:SVD分解求解线性方程组AX=b
R语言:SVD分解求解AX=b 函数svd的用法 svd分解X,使用函数svd,返回一个列表T,顺序是d, u, v。 X=UDV′X=UDV′X=UDV' T <- svd(X) U <- T$u D <- T$d V <- T$v 这里T是list,注意这里的U和V是矩阵,D是向量,想要恢复原矩阵X,需要: Y <- U %*% diag(D) %*...原创 2018-08-30 17:58:35 · 5900 阅读 · 0 评论 -
R语言:soft thresholding operator
R语言:矩阵计算 问题:已知矩阵AAA和参数kkk,希望返回矩阵SSS,使得对SSS的每个元素: S(a)=⎧⎩⎨⎪⎪a−k,0,a+k,a>k|a|≤ka<−kS(a)={a−k,a>k0,|a|≤ka+k,a<−kS(a) = \begin{cases} a-k, & a>k \\ 0, &|a|\le k \\ a+k, & aSSS实际上是叫做soft thre...原创 2018-08-30 08:42:33 · 880 阅读 · 0 评论