小波分析中小波系数的等值线问题

huzengyun

本人刚接触小波分析不久，在论坛上找了一些小波分析的程序。但是对于小波分析中小波系数的等值线图画法有些问题不明白？
第一：如何通过小波系数图（图a）判断周期？
第二：图a和图c中等值线的数值是怎么确定的？就是 levels 里面的数值时固定的吗？
附上代码（论坛中的代码，非本人）
close all
clear all
clc
load 'q.mat'

sst =p;
n=length(p);%时间序列的长度
variance = std(sst)^2;%时间序列的标准差
sst = (sst - mean(sst))/sqrt(variance); % 时间序列标准化
dt = 1;% 采样周期
year = [0:n-1]+ 1962.0 ;
xlim = [1962,2011];   % 画图范围，plotting range.
pad = 1;%     
dj = 1/12;     
s0 = 1/2*dt;   
j1 = 6.5/dj;   
lag1 = 0.72;  
mother = 'Morlet';% 选择的小波函数
[wave,period,scale,coi] = wavelet(sst,dt,pad,dj,s0,j1,mother);
power = (abs(wave)).^2 ;  %计算小波系数的模的平方
modulus=abs(wave);  %计算小波系数的模
variance1=sum(power')/n;%计算小波方差
%画小波系数实部等值线图


subplot(3,1,1)
levels = [0,0.5,1.0,1.5,2.0,2.5];
v = [0,0.5,1.0,1.5];
Yticks = 0:5:45;
[c,h]=contour(year,period,real(wave),levels,'k-');
clabel(c,h,v,'fontsize',5);
xlabel('年份/year')
ylabel('周期/年 period/year')
title('(a)')
set(gca,'XLim',xlim(:))
set(gca,'YLim',[0 50], ...
    'YDir','default', ...
'YTick',Yticks(:), ...
'YTickLabel',Yticks)
hold on
levels = [-0.5,-1.0,-1.5,-2.0,-2.5];
v = [-0.5,-1.0,-1.5];
[c,h] = contourf(year,period,real(wave),levels,'r--');
% contourf中心添加颜色, contour 不添加颜色
clabel(c,h,v,'fontsize',5);
hold on
% 画小波方差图
subplot(3,1,2)
plot(period,variance1,'k-')
hold on;
levels= [1,5,10,15,20,25,30,35,40,45];
title('(b)')
set(gca,'XLim',[1,50], ...        
    'XTick',levels,...
   'XTickLabel',levels)
xlabel('周期/a')
ylabel('方差 variance')
hold on
%画小波系数模
subplot(3,1,3)
levels = [0,0.5,1.0,1.5,2.0,2.5];
v = [0,0.5,1.0,1.5];
Yticks = [0:5:30];
[c,h]=contourf(year,period,abs(wave),levels,'k-');
clabel(c,h,v,'fontsize',5);
title('(c)')
xlabel('年份/year')
ylabel('周期/年 period/year')
set(gca,'XLim',xlim(:))
set(gca,'YLim',[0 30], ...
    'YDir','default', ...
'YTick',Yticks(:), ...
'YTickLabel',Yticks)

kongfeng0824

看看文献                           

huzengyun

看了好多文献还是不明白，请指点下。多谢

fengyu1983

你的图呢？

huzengyun

fengyu1983 发表于 2014-8-26 15:34
你的图呢？

不会上传，不好意思。那个程序可以直接画图. 方便可以加qq请教下：250395497.多谢

zilla

请问你明白了吗？同样求教啊，qq1020459500

tangjian

xiexie!!!!!!!!!

huzengyun

没有细追究。现在用EEMD