初学应用举例（小波分析，实部画图，小波方差画图及小波模画图等）

adnap

心灵之约 发表于 2013-3-21 13:34

                               
登录/注册后可看大图

估计是软件工作路径没有设置好吧

是滴，太粗心了，图已出如下：


多谢大侠指导！

adnap

楼主的程序运行前原数据是否需经标准化或者距平处理？
我发现经过距平处理和没经过距平处理后数据画出来的图都一样......

心灵之约

adnap 发表于 2013-3-21 16:28

                               
登录/注册后可看大图

楼主的程序运行前原数据是否需经标准化或者距平处理？
我发现经过距平处理和没经过距平处理后数据 ...

请问能否将标准化处理的程序语言贴出来，大家共享哈，我也有这个疑问呢。

心灵之约

adnap 发表于 2013-3-21 15:22

                               
登录/注册后可看大图

是滴，太粗心了，图已出如下：

不用客气哈，希望能指点下标准化的陈故乡语言哈。

adnap

心灵之约 发表于 2013-3-21 16:40

                               
登录/注册后可看大图

不用客气哈，希望能指点下标准化的陈故乡语言哈。

目前看到有些文献说小波分析前需要将原数据标准化或距平处理，我目前只做了距平处理，标准化还没做，发现经过距平处理和没经过距平处理后数据画出来的图都一样......所以目前困惑是否需要做距平处理？

星星

赞！超级给力

心灵之约

adnap 发表于 2013-3-22 09:14

                               
登录/注册后可看大图

目前看到有些文献说小波分析前需要将原数据标准化或距平处理，我目前只做了距平处理，标准化还没做，发现 ...

图是一样的啊，那估计就是多余的啦，呵呵

adnap

心灵之约 发表于 2013-3-22 16:19

                               
登录/注册后可看大图

图是一样的啊，那估计就是多余的啦，呵呵

就是想跟楼主确实一下{:2_25:}

╃゛湘τ

谢谢lz分享，正好需要！

学长吉他

clear;
load q.txt
sst =q;
n=length(q);
variance = std(sst)^2;
sst = (sst - mean(sst))/sqrt(variance);
dt = 1;
year = [0:n-1]+ 1951.0 ;
xlim = [1951,2010];   %
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;%计算小波方差
%画小波系数实部等值线图
fig1=figure;
%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:30;
[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 30], ...
    '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] = contour(year,period,real(wave),levels,'r--');
clabel(c,h,v,'fontsize',5);
hold on
% 画小波方差图
fig2=figure;
%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
%画小波系数模
fig3=figure;
%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]=contour(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)