简单分享：matlab小波用法简单翻译

一如既往啊

谢谢大佬，大佬牛逼谢谢大佬，大佬牛逼

高须龙儿

特地来大神这学习小波分析的方法的

豆豆啦啦啦啦噢

学习学习 希望有帮助‘
’

水水的小家

楼主真的太贴心啦

水水的小家

函数wavelet和cwt之间有什么区别呢，有些人用cwt有些用wavelet，分析出来的结果有什么不同吗？小白求大神们指点

Dancing_Piggy

感谢！！！！！！

hdtcl

clear; %清除变量空间
load q.txt  % 装载q.txt文件到MATLAB
sst = q;    % 获取q.txt文件的句柄，返回给变量sst.
%------------------------------------------------------ Computation  计算
variance = std(sst)^2;  % 计算sst的方差
sst = (sst - mean(sst))/sqrt(variance) ;
n = length(sst);      
dt =1;     %每个Y值之间的时间量 即 x轴的间距，影响图形的拉伸
time = [0:length(sst)-1]*dt + 1963.0 ;  %此处时间“1963”可根据需要自行修改
本帖隐藏的内容
xlim = [1962,2011];  %此处时间[1962,2011]可根据需要自行修改，即资料时间尺度
pad = 1;  % 见matlab小波用法简单翻译http://bbs.06climate.com/forum.p ... &extra=page%3D1
dj = 1/12; %同上  
s0 = 2*dt;    %同上
j1 =fix((log(n*dt/s0)/log(2))/dj);   %同上
mother = 'Morlet';  %同上
% Wavelet transform: 小波变换
[wave,period,scale,coi] = wavelet(sst,dt,pad,dj,s0,j1,mother);%见matlab小波用法简单翻译
power = (abs(wave)).^2 ;  % 求小波的功率谱 在此程序里用不上   
realpart=real(wave); %求小波的实部
modulus=abs(wave); %求小波的模 在此程序里用不上
phase=atan2(imag(wave),real(wave)); % 求小波的阶 在此程序里用不上
%------------------------------------------------------ Plotting 画图 （以后有机会再详述）


% Plot wavelet coefficients of time series
subplot('position',[0.10 0.30 0.80 0.60])
Yticks = 2.^(fix(log2(min(period))):fix(log2(max(period))));
[c,h]=contour(time,log2(period),realpart,'-');  
clabel(c,h);
xlabel('Time (year)')
ylabel('Period (years)')
%title('a) Morlet Wavelet Real Part of the Temperature Anomaly Series')
set(gca,'XLim',xlim(:))
set(gca,'YLim',log2([min(period),max(period)]), ...
    'YDir','default', ...
'YTick',log2(Yticks(:)), ...
'YTickLabel',Yticks)

大理段戬

来学习一下

虾米吃小鱼

感谢大神分享，刚开始学习小波

大理段戬

想问一下，小波分析后得到的周期必须是整数么？与采用年平均或者月平均数据有什么区别？