condest(Matrix)
condest()所属R语言包:Matrix
Compute Approximate CONDition number and 1-Norm of (Large) Matrices
计算近似条件的数量和1-范(大)矩阵
译者:生物统计家园网 机器人LoveR
描述----------Description----------
“Estimate”, i.e. compute approximately the CONDition number of a (potentially large, often sparse) matrix A. It works by apply a fast approximation of the 1-norm, norm(A,"1"), through onenormest(.).
“估计”,即计算约(潜在的大,经常稀疏)矩阵A条件。它的工作原理应用1,规范快速逼近,norm(A,"1"),onenormest(.)。
用法----------Usage----------
condest(A, t = min(n, 5), normA = norm(A, "1"),
silent = FALSE, quiet = TRUE)
onenormest(A, t = min(n, 5), A.x, At.x, n,
silent = FALSE, quiet = silent,
iter.max = 10, eps = 4 * .Machine$double.eps)
参数----------Arguments----------
参数:A
a square matrix, optional for onenormest(), where instead of A, A.x and At.x can be specified, see there.
一个方阵,可选的onenormest(),而不是A,A.x和At.x可以指定,看到有。
参数:t
number of columns to use in the iterations.
使用迭代的列数。
参数:normA
number; (an estimate of) the 1-norm of A, by default norm(A, "1"); may be replaced by an estimate.
,默认情况下数(估计)1规范Anorm(A, "1");可取代估计。
参数:silent
logical indicating if warning and (by default) convergence messages should be displayed.
逻辑表明,如果应显示警告(默认)收敛消息。
参数:quiet
logical indicating if convergence messages should be displayed.
逻辑表明,如果收敛消息应显示。
参数:A.x, At.x
when A is missing, these two must be given as functions which compute A %% x, or t(A) %% x, respectively.
时A丢失,这两个必须为计算A %% x或t(A) %% x,分别的功能。
参数:n
== nrow(A), only needed when A is not specified.
== nrow(A)时A未指定,只需要。
参数:iter.max
maximal number of iterations for the 1-norm estimator.
1-范估计的迭代的最大数量。
参数:eps
the relaive change that is deemed irrelevant.
的relaive变化,被认为是无关紧要的。
值----------Value----------
Both functions return a list; onenormest() with components,
这两个函数返回一个listonenormest()组件,
参数:est
a number > 0, the estimated norm(A, "1").
数字> 0,估计norm(A, "1")。
参数:v
the maximal A X column.
最大A X列。
The function condest() returns a list with components,
功能condest()与组件返回一个列表,
参数:est
a number > 0, the estimated condition number k.; when r :=rcond(A), 1/k. ~= r.
数字> 0,估计的条件数k.;r := rcond(A),1/k. ~= r。
参数:v
integer vector length n, with an 1 at the index j with maximal column A[,j] in A.
整数向量的长度n指数1最大列jA[,j]A。
参数:w
numeric vector, the largest A x found.
数字向量,最大A x发现。
参数:iter
the number of iterations used.
使用迭代的数量。
作者(S)----------Author(s)----------
This is based on octave's <code>condest()</code> and
<code>onenormest()</code> implementations with original author
Jason Riedy, U Berkeley; translation to <font face="Courier New,Courier" color="#666666"><b>R</b></font> and
adaption by Martin Maechler.
参考文献----------References----------
A Block Algorithm for Matrix 1-Norm Estimation, with an Application to 1-Norm Pseudospectra. SIAM J. Matrix Anal. Appl. 21, 4, 1185–1201. http://dx.doi.org/10.1137/S0895479899356080
Condition Estimates. SIAM J. Sci. Stat. Comput. 5, 311–316.
参见----------See Also----------
norm, rcond.
norm,rcond。
举例----------Examples----------
data(KNex)
mtm <- with(KNex, crossprod(mm))
system.time(ce <- condest(mtm))
## reciprocal[#倒数]
1 / ce$est
system.time(rc <- rcond(mtm)) # takes ca 3 x longer[需要约3×长]
rc
all.equal(rc, 1/ce$est) # TRUE -- the approxmation was good[TRUE - 良好的approxmation]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|