计算数学
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$uD <- T$dV <- T$v这里T是list,注意这里的U和V是矩阵,D是向量,想要恢复原矩阵X,需要:Y <- U %*% diag(D) %*...原创 2018-08-30 17:58:35 · 5849 阅读 · 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 · 867 阅读 · 0 评论