spec(seewave)
spec()所属R语言包:seewave
Frequency spectrum of a time wave
频率的频谱的时间波
译者:生物统计家园网 机器人LoveR
描述----------Description----------
This function returns the frequency spectrum (i.e. the relative amplitude of the frequency content) of a time wave. Results can be obtained either as absolute or dB data.
此函数返回的时间波的频谱(即相对振幅的频率内容)。结果可以为绝对或DB数据。
用法----------Usage----------
spec(wave, f, wl = 512, wn = "hanning", fftw = FALSE, PSD = FALSE,
PMF = FALSE, dB = NULL, dBref = NULL,
at = NULL, from = NULL, to = NULL,
identify = FALSE, col = "black", cex = 1,
plot = 1, flab = "Frequency (kHz)",
alab = "Amplitude", flim = NULL,
alim = NULL, type="l",...)
参数----------Arguments----------
参数:wave
an R object.
R对象。
参数:f
sampling frequency of wave (in Hz). Does not need to be specified if embedded in wave.
wave(赫兹)的采样频率。不需要以指定如果嵌入在wave,。
参数:wl
if at is not null, length of the window for the analysis (even number of points, by defaults = 512).
如果at是不空的窗口,长度的分析(偶数点,默认值= 512)。
参数:wn
window name, see ftwindow (by default "hanning").
窗口的名称,请参阅ftwindow(默认情况下,"hanning"“)。
参数:fftw
if TRUE calls the function FFT of the library fftw. See Notes of the spectro. </table>
如果TRUE调用函数FFT的库fftw。见注意的spectro。 </ TABLE>
参数:PSD
if TRUE return Power Spectra Density, i. e. the square of the spectra.
如果TRUE回报功率谱密度。 E。谱的平方。
参数:PMF
if TRUE return Probability Mass Function, i. e. the probability distribution of frequencie.
如果TRUE返回概率密度函数。 E。频率下的概率分布。
参数:dB
a character string specifying the type dB to return: "max0" for a maximum dB value at 0, "A", "B", "C" and "D" for common dB weights.
指定的类型分贝返回一个字符的字符串:“max0”最大的dB值在0,“A”,“B”,“C”和“D”的共同分贝重量。
参数:dBref
a dB reference value when dB is not NULL. NULL by default but should be set to 2*10e-5 for a 20 microPa reference (SPL).
dB参考值,当dB是不是NULL的。 NULL缺省情况下,但应该被设置为2 * 10e的-5 20 microPa参照(SPL)。
参数:at
position where to compute the spectrum (in s).
的位置来计算的频谱(单位:s)。
参数:from
start mark where to compute the spectrum (in s).
开始标记来计算频谱(单位:秒)。
参数:to
end mark where to compute the spectrum (in s).
结束标记来计算频谱(单位:秒)。
参数:identify
to identify frequency and amplitude values on the plot with the help of a cursor.
识别的频率和振幅值的图上的帮助的光标。
参数:col
colour of the spectrum.
色的光谱。
参数:cex
pitch size of the spectrum.
频谱的间距尺寸。
参数:plot
if 1 returns frequency on x-axis, if 2 returns frequency on y-axis, (by default 1).
如果1返回x轴的频率,如果2返回y轴的频率(默认情况下,1“)。
参数:flab
title of the frequency axis.
在频率轴上的标题。
参数:alab
title of the amplitude axis.
标题的振幅轴。
参数:flim
range of frequency axis.
频率轴的范围。
参数:alim
range of amplitude axis.
范围的振幅轴。
参数:type
if plot is TRUE, type of plot that should be drawn. See plot for details (by default "l" for lines).
如果plot是TRUE,应该得出的图。见plot的详细信息(默认情况下,“L”线)。
参数:...
other plot graphical parameters.
其他plot图形参数。
Details
详细信息----------Details----------
If at, from or to are FALSE then spec
如果at,from或to是FALSE然后spec
值----------Value----------
This function returns a two-column matrix, the first column corresponding to the frequency axis, the second column corresponding to the amplitude axis.<br> If identify is TRUE, spec returns a list with two elements:<br>
该函数返回一个两列的矩阵中,第一列相对应的,在频率轴上的第二列对应于振幅轴。<br>物理化学学报如果identifyTRUE,spec返回一个有两个元素的列表:参考
参数:freq
the frequency of the points chosen on the spectrum </table>
的频率选择的点上的频谱</表>
参数:amp
the relative amplitude of the points chosen on the spectrum </table>
所选择的点上的相对幅度的频谱</表>
警告----------Warning----------
The argument peaks is no more available (version > 1.5.6). See the function fpeaks
参数peaks是没有更多的可用(版本1.5.6)。请参阅的功能fpeaks
注意----------Note----------
This function is based on fft.
此函数是基于fft。
(作者)----------Author(s)----------
Jerome Sueur
参见----------See Also----------
meanspec, fpeaks, localpeaks, dynspec,
meanspec,fpeaks,localpeaks,dynspec,
实例----------Examples----------
data(tico)
# spectrum of the whole signal, in absolute or dB amplitude,[整个信号的频谱中,在绝对或3dB振幅,]
# horizontaly or vertically[horizontaly或垂直]
op<-par(mfrow=c(2,2))
spec(tico,f=22050)
spec(tico,f=22050,col="red",plot=2)
spec(tico,f=22050,dB="max0",col="blue")
spec(tico,f=22050,dB="max0",col="green",plot=2)
par(op)
# an indirect way to compare spectra [以间接的方式比较靠谱]
a<-spec(tico,f=22050,wl=512,at=0.2,plot=FALSE)
b<-spec(tico,f=22050,wl=512,at=0.7,plot=FALSE)
c<-spec(tico,f=22050,wl=512,at=1.1,plot=FALSE)
d<-spec(tico,f=22050,wl=512,at=1.6,plot=FALSE)
all<-cbind(a[,2],b[,2],c[,2],d[,2])
matplot(x=a[,1],y=all,yaxt="n",
xlab="Frequency (kHz)",ylab="Amplitude",xaxs="i",type="l")
legend(8,0.8,c("Note A","Note B", "Note C", "Note D"),bty="o",
lty=c(1:4),col=c(1:4))
# spectrum from a particular position to another one[频谱从一个特定的位置,另一个]
op<-par(mfrow=c(2,1))
oscillo(tico,f=22050)
abline(v=c(0.5,0.9),col="red",lty=2)
spec(tico,f=22050,wl=512,from=0.5,to=0.9,col="red")
title("Spectrum of the note B")
par(op)
# spectrum and spectrogram[频谱及频谱图]
data(orni)
orni1<-cutw(orni,f=22050,from=0.32,to=0.39)
layout(matrix(c(1,2),nc=2),widths=c(3,1))
par(mar=c(5,4,3,0.5))
spectro(orni1,f=22050,wl=128,zp=8,ovlp=85,scale=FALSE)
par(mar=c(5,1,3,0.5))
spec(orni1,f=22050,col="red",plot=2,flab="",yaxt="n")
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|