evolMTM(RSEIS)
evolMTM()所属R语言包:RSEIS
Evolutive Multi-taper Spectrum
进化多锥度频谱
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Time varying Multi-taper Spectrum (Gabor Transform)
随时间变化多锥形谱(Gabor变换)
用法----------Usage----------
evolMTM(a, dt = 0, numf = 1024, Ns = 0, Nov = 0, fl = 0, fh = 10)
参数----------Arguments----------
参数:a
Signal
信号
参数:dt
Sample rate interval (s)
采样率间隔(秒)
参数:numf
Number of points in fft
在FFT中的点的数量
参数:Ns
Number of sample in sub-window
在子窗口数目的样
参数:Nov
Number of sample to overlap
的采样次数重叠
参数:fl
low frequency to display
低频显示
参数:fh
high frequency to display
高频显示
Details
详细信息----------Details----------
This is a spectrogram function similar to the Gabor Transform but uses the MTM method for spectrum estimation.
这是一个类似的Gabor变换的频谱功能,但使用的MTM谱估计方法。
值----------Value----------
List
表
参数:sig
input signal
输入信号
参数:dt
deltat
DeltaT的
参数:wpars
input parameters
输入参数
参数:DSPEC
spectrum image
光谱图像
参数:freqs
output frequencies (y axis)
输出频率(y轴)
参数:tims
output times (x-axis)
输出时间(x-轴)
(作者)----------Author(s)----------
Jonathan M. Lees<jonathan.lees@unc.edu>
参考文献----------References----------
analysis: A stand-alone C-subroutine,
参见----------See Also----------
evolfft, MTM.drive
evolfft,MTM.drive
实例----------Examples----------
data(KH)
### swig(KH)[##痛饮(KH)]
Xamp = KH$JSTR[[1]]
dt = KH$dt[1]
plot(seq(from=0, length=length(Xamp), by=dt), Xamp, type='l')
## limit the trace, somewhat[#限制的痕迹,有些]
Xamp = Xamp[12670:22669]
plot(seq(from=0, length=length(Xamp), by=dt), Xamp, type='l')
Nfft=4096 ### fft length[##FFT长度]
Ns=512 ### number of samples in a window[##在一个窗口中的样本数]
Nov=480 ### number of samples of overlap per window[##每个窗口的重叠的样本数]
fl=0 ### low frequency to return[##低频到返回]
fh=12 ### high frequency to return[##高频返回]
EV = evolMTM(Xamp, dt = dt, numf = Nfft, Ns = Ns, Nov = Nov, fl = fl, fh
= fh)
PE = plotevol(EV, log=1, fl=0.01, fh=fh, col=rainbow(100), ygrid=FALSE,
STAMP="", STYLE="ar")
## compare with:[#对比:]
## EVf = evolfft(Xamp, dt = dt, Nfft =Nfft , Ns =Ns , Nov =Nov , fl =fl, fh = fh)[#EVF = evolfft(XAMP,DT = DT,NFFT = NFFT,NS = NS,月月,FL = FL,FH = FH)]
## PE = plotevol(EVf, log=1, fl=fl, fh=fh, col=rainbow(100), ygrid=FALSE,STAMP="", STYLE="fft")[#PE = plotevol(EVF,log= 1,FL = FL,FH = FH,列=彩虹(100),YGRID = FALSE,印花税=“STYLE =”FFT“)]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|