- 积分
- 44
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-8-24
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2015-1-3 21:06:17
|
显示全部楼层
请帮我看看这个代码是不是有问题
scales=[1:1:32];
%进行连续小波变换得到小波系数矩阵,选择复morlet小波函数
wf=cwt(qiujijiangyu,scales,'cmor1-1.5','plot');
% 求得系数的实部
shibu=real(wf);
subplot(221);
contourf(shibu,10,'-');
colormap('HSV');
colorbar;
time=1951:6:2013
Xlabel('年份Year');
Ylabel('时间尺度Time scale/a');
set(gca,'XTickLabel', time) %更新XTickLabel
% 小波方差是模的平方的算数平均
mo=abs(wf);
mofang=mo.^2;
fangcha=mean(mofang,2);
subplot(222);
plot(fangcha,'k-','linewidth',0.5);
Xlabel('时间尺度/年');
Ylabel('小波方差')
|
|