R中小波的二维图谱

时间:2015-07-30 12:27:36

标签: r wavelet

我是使用小波的新手。我已经做了一些关于如何在R中实现它们的背景研究。为此,我使用了Carmona,Hwang和Torresani的实用时频分析书中的例4.2。代码如下。

library(Rwave)
sing <- numeric(256)
sing[31:60] <- (sqrt(1:30)/sqrt(30))
sing[60:90] <- 1
sing[140] <- 0.5
sing[200] <- 0.5
sing[201] <- -0.5
plot(sing,type="l", main="Test CWT Data")
cwtsing<-cwt(sing,5,12)
cwtsing

我有两个问题: -

  1. 有没有办法使用另一个“基础”小波?这个例子使用Morlet小波但是可以以某种方式使用db4或sym4小波吗?
  2. 该R代码产生的频谱图在y轴上具有log(标度)。有没有一种快速的方法来改变图形的外观,就像在下面的链接中对问题2(参见彩虹色2D图)所做的那样。 https://mathematica.stackexchange.com/questions/34333/extracting-information-from-the-result-of-continuouswavelettransform
  3. 我知道我可以手动编写算法来进行音阶转换,但我希望能有一个快速的解决方案。任何帮助将不胜感激。提前感谢任何人。

0 个答案:

没有答案