我想对一个脑电信号进行时频分析。我发现GSL小波函数可以计算小波系数。如何从这些系数中提取实际的频率带(例如8-12 Hz)? GSL手册上说:
对于正向变换,原始数组的元素被离散小波变换f_i -> w_{j,k}所替换,以紧密的三角形存储布局为基础,其中J是级别的索引j = 0 ... J-1和K是每个级别内系数的索引,k = 0 ... (2^j)-1。总共的级别数量为J = log_2(n)。
输出数据具有以下形式,(s_{-1,0}, d_{0,0}, d_{1,0}, d_{1,1}, d_{2,0}, ..., d_{j,k}, ..., d_{J-1,2^{J-1}-1})
如果我理解正确,输出数组data[]在位置1(例如data[1])包含频率带2^0 = 1Hz的幅度。
对于正向变换,原始数组的元素被离散小波变换f_i -> w_{j,k}所替换,以紧密的三角形存储布局为基础,其中J是级别的索引j = 0 ... J-1和K是每个级别内系数的索引,k = 0 ... (2^j)-1。总共的级别数量为J = log_2(n)。
输出数据具有以下形式,(s_{-1,0}, d_{0,0}, d_{1,0}, d_{1,1}, d_{2,0}, ..., d_{j,k}, ..., d_{J-1,2^{J-1}-1})
如果我理解正确,输出数组data[]在位置1(例如data[1])包含频率带2^0 = 1Hz的幅度。
data[2] = 2^1 Hz
data[3] = 2^1 Hz
data[4] = 2^2 Hz
until
data[7] = 2^2 Hz
data[8] = 2^3 Hz
意思是我只有1Hz、2Hz、4Hz、8Hz、16Hz等频率的振幅。如果我想得到5.3Hz振荡频率的分量的振幅,该怎么办?如果我想要整个频率范围比如8-13Hz的振幅,又该怎么做呢?是否有一些好的建议可以获取良好的时频分布呢?