小波分析中一些参数的涵义

可乐加盐

dj = 0.25;    % this will do 4 sub-octaves per octave
s0 = 2*dt;    % this says start at a scale of 6 months
j1 = 7/dj;    % this says do 7 powers-of-two with dj sub-octaves each
用的是Torren论文中用的程序， 但其中一些参数的设置从论文中还是找不到。或许是我英文功底太差了吧。希望论坛的大虾揭示下具体含义。我想占整个帖子不单单是说这三个参数后面还有些参数，我也提出来，求大虾指导

可乐加盐

dt = 0.25 ;
time = [0:length(sst)-1]*dt + 1871.0 ;  % construct time array
xlim = [1870,2000];  % plotting range
pad = 1;      % pad the time series with zeroes (recommended)
dj = 0.25;    % this will do 4 sub-octaves per octave
s0 = 2*dt;    % this says start at a scale of 6 months
j1 = 7/dj;    % this says do 7 powers-of-two with dj sub-octaves each
lag1 = 0.72;  % lag-1 autocorrelation for red noise background
比如这几个参数。
lag1=0.72是表示一阶自相关系数为0.72的红噪音背景场
dt是 取样时间间隔
但dj和j1是什么意思呢

可乐加盐

还有
levels = [0.0625,0.125,0.25,0.5,1,2,4,8,16]
contour(time,log2(period),log2(power),log2(levels));
这句
应该是以横坐标为time,纵坐标为Log2(period)，画log2(power)的值大小，是以等值线的方式，并且间隔为log2(levels)
但这个levels的取值是怎么来的呢。在Torrence的论文里面并没有遇见呢

可乐加盐

在信度检验里面
dj的赋值又不一样了。之前的dj是表示在蒙特卡洛检验的情况下赋值为0.25.但在信度检验里面
n1 = length(Y);
J1 = length(scale1) - 1;
scale(1:J1+1) = scale1;
s0 = min(scale);
dj = log(scale(2)/scale(1))/log(2.);
dj不等于0.25

可乐加盐

没人做过这方面的吗？

可乐加盐

@阿华田

黎明的夕阳

不错不错

可乐加盐

黎明的夕阳 发表于 2012-10-30 20:27

                               
登录/注册后可看大图

不错不错

啥不错啊。

xjgkaka

主任的问题没人解决？

可乐加盐

xjgkaka 发表于 2012-11-8 08:39

                               
登录/注册后可看大图

主任的问题没人解决？

老徐！！！