- 积分
- 397
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-11-13
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2013-9-15 19:58:45
|
显示全部楼层
本帖最后由 nuistcang 于 2013-9-15 20:21 编辑
scaler 发表于 2013-9-10 19:25
真是,以为可以搞定周期置信度检验的,
我试了下,这样应该可以~
数据分析时加:
lag1 = corrcoef(sst(1:(end-1)),sst(2:end)); %lag-1 autocorrelation for red noise background
lag1=lag1(1,2);
[signif,fft_theor] = wave_signif(1.0,dt,scale,0,lag1,-1,-1,'morlet');
sig95 = (signif')*(ones(1,n));
sig95 = power ./ sig95;
画图时加:
[c,h]=contourf(time,log2(period),sig95,[1 1]);
set(h,'linewidth',2)
hold on
plot(time,log2(coi),'k')
hold off结果是这样的:
这其中用到一个wave_signif.m的程序,在附件中给出。
wave_signif.m
(6.58 KB, 下载次数: 37)
|
|