怎样调整功率谱分析的参数

mcVS米虫

这是根据家园上下载的程序作出的功率谱图片，想问大家的是，感觉置信区间画出来怪怪的，我自己怀疑是纵坐标周期太大的问题，看文献中画NPO这部分周期最大是32，想问下大家怎样调整？以下是程序中需要用到的参数：
dt = 1.0 %每个数据之间的时间步长。

time = [0:length(sst)-1]*dt + 1900.0 ;  % construct time array
xlim = 1899+[1,n];  % plotting range
pad = 1;      % pad the time series with zeroes (recommended)
dj = 0.05;    % this will do 4 sub-octaves per octave; the spacing between discrete scales. Default is 0.25. a smaller # will give better scale resolution, but be slower to plot.
s0 = 2*dt;    % this says start at a scale of 6 months
j1 = 6/dj;    %%  Default is j1 = (LOG2(n DT/S0))/dj.  this says do 7 powers-of-two with dj sub-octaves each   zl：改变纵轴数值 64 128 256...
lag1 = 0.72;  % lag-1 autocorrelation for red noise background
mother = 'morlet';

······
Cdelta = 0.776;   % this is for the MORLET wavelet

孙玻玻

你的第七行 j1=6/dj，所取得周期是你研究时段总时长的二分之一或三分之一，你用你总时长的二分之一再乘以0.25取得的那个整数就是你要改的数。

mcVS米虫

感谢您的回答