fastkernel(Rwave)
fastkernel()所属R语言包：Rwave
Kernel for Reconstruction from Wavelet Ridges
小波脊重建内核
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Computes the cost from the sample of points on the estimated ridge and the matrix used in the reconstruction of the original signal, using simple trapezoidal rule for integrals.
从样品的估计脊点和重建原始信号中所使用的矩阵计算成本，使用简单的梯形规则积分。
用法----------Usage----------
fastkernel(node, phinode, nvoice, x.inc=1, x.min=node[1],
x.max=node[length(node)], w0=2 * pi, plot=FALSE)
参数----------Arguments----------
参数：node
values of the variable b for the nodes of the ridge.
脊的节点的变量b的值。
参数：phinode
values of the scale variable a for the nodes of the ridge.
规模变量a的脊的节点的值。
参数：nvoice
number of scales within 1 octave.
1个八度内的鳞片数。
参数：x.inc
step unit for the computation of the kernel
步单元，用于计算的内核
参数：x.min
minimal value of x for the computation of Q_2.
最小的x值计算Q_2。
参数：x.max
maximal value of x for the computation of Q_2.
最大的x值计算Q_2。
参数：w0
central frequency of the wavelet
中心频率的小波
参数：plot
if set to TRUE, displays the modulus of the matrix of Q_2. </table>
如果设置为TRUE，则显示的矩阵Q_2的模量。 </ TABLE>
Details
详细信息----------Details----------
Uses trapezoidal rule (instead of Romberg's method) to evaluate the kernel.
使用梯形规则（而不是伯格的方法），以评估的核心。
值----------Value----------
matrix of the Q_2 kernel.
矩阵Q_2内核。
参考文献----------References----------
参见----------See Also----------
kernel, rkernel, gkernel, zerokernel.
kernel，rkernel，gkernel，zerokernel。
实例----------Examples----------
# The function is currently defined as[的功能目前被定义为]
function(node, phinode, nvoice, x.inc = 1, x.min = node[1], x.max = node[length(node)], w0 = 2 * pi, plot = F)
{
#########################################################################[＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃ ＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃]
# fastkernel: [fastkernel：]
# -----------[-----------]
# Same as kernel, except that the kernel is computed[内核相同，不同之处在于计算内核]
# using Riemann sums instead of Romberg integration.[Romberg积分，而不是使用黎曼和。]
#[]
# Input:[输入：]
# ------[------]
# node: values of the variable b for the nodes of the ridge[节点：的变量b的值的节点的脊]
# phinode: values of the scale variable a for the nodes of the ridge[phinode：刻度变量的值的一个节点的脊]
# nvoice: number of scales within 1 octave[nvoice：1个八度的尺度内数]
# x.inc: step unit for the computation of the kernel[x.inc：步骤单元，用于计算的内核]
# x.min: minimal value of x for the computation of Q2[x.min：最小的x值计算Q2]
# x.max: maximal value of x for the computation of Q2[x.max：最大的x值的计算Q2]
# w0: central frequency of the wavelet[w0中：小波的中心频率]
# plot: if set to TRUE, displays the modulus of the matrix of Q2[图：如果设置为TRUE，则显示Q2矩阵的模]
#[]
# Output:[输出：]
# -------[-------]
# ker: matrix of the Q2 kernel[KER：矩阵的Q2内核]
#[]
#########################################################################[＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃ ＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃]
lng <- as.integer((x.max - x.min)/x.inc) + 1
nbnode <- length(node)
b.start <- x.min - 50
b.end <- x.max + 50
ker.r <- matrix(0, lng, lng)
ker.i <- matrix(0, lng, lng)
dim(ker.r) <- c(lng * lng, 1)
dim(ker.i) <- c(lng * lng, 1)
phinode <- 2 * 2^(phinode/nvoice)
z <- .C(fastkernel,
ker.r = as.double(ker.r),
ker.i = as.double(ker.i),
as.integer(x.min),
as.integer(x.max),
as.integer(x.inc),
as.integer(lng),
as.double(node),
as.double(phinode),
as.integer(nbnode),
as.double(w0),
as.double(b.start),
as.double(b.end))
ker.r <- z$ker.r
ker.i <- z$ker.i
dim(ker.r) <- c(lng, lng)
dim(ker.i) <- c(lng, lng)
ker <- matrix(0, lng, lng)
i <- sqrt(as.complex(-1))
ker <- ker.r + i * ker.i
if(plot == T) {
par(mfrow = c(1, 1))
image(Mod(ker))
title("Matrix of the reconstructing kernel (modulus)")
}
ker
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-29 20:44:26
