请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 668|回复: 4

[源程序] 小波变换图小尺度为空白

[复制链接] |关注本帖

新浪微博达人勋

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
各位大神,请帮忙看一下下面的小波变换图,为什么尺度为1的是空白的?拜托各位大神了,
下面是所用的程序:
clear;
load g.txt
sst =g;
n=length(g);
variance = std(sst)^2;
sst = (sst - mean(sst))/sqrt(variance);
dt = 1;
year = [0:n-1]+ 1912.0 ;
xlim = [1912,2000];   %
pad = 1;      
dj = 1/12;     
s0 = 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 = [-4,-3.5,-3,-2.5,-2.0,-1.5,-1.0,-0.5,0,0.5,1.0,1.5,2.0,2.5,3,3.5,4];
v = [-3.5,-3,-2.5,-2.0,-1.5,-1.0,-0.5,0,0.5,1.0,1.5,2.0,2.5,3,3.5];
Yticks = 0:10:50;
[c,h]=contourf(year,period,real(wave),levels,'k-');
clabel(c,h,v,'fontsize',5);
y=real(wave)
xlabel('年份/year')
ylabel('周期/年 period/year')
set(gca,'XLim',xlim(:))
set(gca,'YLim',[0 50], ...
     'YDir','default', ...
'YTick',Yticks(:), ...
'YTickLabel',Yticks)
hold on
% 画小波方差图
subplot(3,1,2)
plot(period,variance1,'k-')
hold on;
levels= [0,10,20,30,40,50];
title('(b)')
set(gca,'XLim',[1 50], ...        
     'XTick',levels,...
    'XTickLabel',levels)
xlabel('周期/a')
ylabel('方差 variance')
hold on


                               
登录/注册后可看大图





密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-4-11 12:39:19 | 显示全部楼层 |取消关注该作者的回复

求助:小波分析图0-2尺度为空白

各位大神,请帮忙看一下下面的小波变换图,为什么0-2尺度的是空白的?拜托各位大神了,
下面是所用的程序:
clear;
load g.txt
sst =g;
n=length(g);
variance = std(sst)^2;
sst = (sst - mean(sst))/sqrt(variance);
dt = 1;
year = [0:n-1]+ 1912.0 ;
xlim = [1912,2000];   %
pad = 1;      
dj = 1/12;     
s0 = 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 = [-4,-3.5,-3,-2.5,-2.0,-1.5,-1.0,-0.5,0,0.5,1.0,1.5,2.0,2.5,3,3.5,4];
v = [-3.5,-3,-2.5,-2.0,-1.5,-1.0,-0.5,0,0.5,1.0,1.5,2.0,2.5,3,3.5];
Yticks = 0:10:50;
[c,h]=contourf(year,period,real(wave),levels,'k-');
clabel(c,h,v,'fontsize',5);
y=real(wave)
xlabel('年份/year')
ylabel('周期/年 period/year')
set(gca,'XLim',xlim(:))
set(gca,'YLim',[0 50], ...
     'YDir','default', ...
'YTick',Yticks(:), ...
'YTickLabel',Yticks)
hold on
% 画小波方差图
subplot(3,1,2)
plot(period,variance1,'k-')
hold on;
levels= [0,10,20,30,40,50];
title('(b)')
set(gca,'XLim',[1 50], ...        
     'XTick',levels,...
    'XTickLabel',levels)
xlabel('周期/a')
ylabel('方差 variance')
hold on


小波变换图

小波变换图
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-4-17 09:20:23 | 显示全部楼层 |取消关注该作者的回复
逐年数据的话,两年的周期怎么可能分辨得出来。。
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-4-24 10:01:23 | 显示全部楼层 |取消关注该作者的回复
大神,那这个图就是对的吗,
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-4-24 10:02:42 | 显示全部楼层 |取消关注该作者的回复
伽蓝鸟 发表于 2018-4-17 09:20
逐年数据的话,两年的周期怎么可能分辨得出来。。

大神,那这个图就是对的吗,我看了许多论文他们好像都没有这个问题
密码修改失败请联系qq:937062711
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表