sigconv(RSEIS)
sigconv()所属R语言包:RSEIS
Convolve spikes with wavelets
卷积与小波高峰
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Convolve spikes with wavelets
卷积与小波高峰
用法----------Usage----------
sigconv(wigmat, wavepulse)
参数----------Arguments----------
参数:wigmat
matrix, spikes
矩阵,尖峰
参数:wavepulse
wavelet for convolution
小波的卷积
Details
详细信息----------Details----------
Convolution is done in Frequency domain on each trace
卷积是在频域上的每个跟踪
值----------Value----------
Matrix, waveforms
矩阵,波形
(作者)----------Author(s)----------
Jonathan M. Lees<jonathan.lees@unc.edu>
参见----------See Also----------
wiggleimage, symshot1, genrick
wiggleimage,symshot1,genrick
实例----------Examples----------
S1 = symshot1()
############## S1$THEORY$treflex[#############S1 $理论$ treflex]
d = dim(S1$smograms)
G1 = matrix( rep(0, length=d[1]*d[2]), ncol=d[2], nrow=d[1])
############ set up the spike set for reflexions[###########设立的尖峰为周明华集]
for(i in 1:3){
p = round( S1$THEORY$treflex[i,]/S1$dt );
G1[cbind(p , 1:d[2]) ] = 1
}
#### plot the spikes[###图的尖峰]
wiggleimage(0.1*G1, dt = -S1$dt, dx = S1$x, col = "black")
### make a ricker wavelet[##做的一个Ricker子波]
wavelet = genrick(25,S1$dt,35)
klem = 11
### nwave = -1 + 2 * (wavelet - min(wavelet) )/diff(range(wavelet))[##nwave = -1 + 2 *“(小波 - 分(小波))/差异范围(小波)]
nwave = RESCALE(wavelet, 0, 1, wavelet[1], max(wavelet))
############## convolve the wavelet with the set of spikes[#############卷积小波与组的尖峰]
H1 = sigconv(G1, nwave)
############ plot[###########图]
wiggleimage(0.1*H1, dt = -S1$dt, dx = S1$x, col = "black")
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|