趋势分析小波mk_matlab_小波分析_序列数据_MK_

%小波系数实部等值线图 a=xlsread('a.xlsx');%xlsx还是x要根据版本选择,文件放在d盘下面 t=a(:,1); u=a(:,2)'; figure plot(u,'k-') u=[repmat(u(1),1,10),u,repmat(u(end),1,10)] %添加10点，去除端点效应 hold on plot(u,'r-') legend('原始数据','添加端点后'); fs=1; figure cwt1=cwt(u,1:64,'morl');%Morl sym32 小波变换，得出小波系数‘ cwt1(:,1:10)=[];%去除端点效应 cwt1(:,end-10:end)=[]; %去除端点效应 [C,h]=contourf(cwt1);%求等值线 shading flat, colorbar clabel(C,h); %加标值 set(gca,'xtick',10:10:70) set(gca,'XTickLabel',{'1945','1955','1965','1975','1985','1995','2005'}); xlabel('年份') ylabel('周期尺度/a') title('(a)') %%小波系数模 figure [C,h]=contourf(abs(real(cwt1)));%求等值线 shading flat, colorbar clabel(C,h); %加标值 set(gca,'xtick',10:10:70) set(gca,'XTickLabel',{'1945','1955','1965','1975','1985','1995','2005'}); xlabel('年份') ylabel('周期尺度/a') title('(a)') %小波模方 figure [C,h]=contourf(real(cwt1).^2);%求等值线 shading flat, colorbar clabel(C,h); %加标值 set(gca,'xtick',10:10:70) set(gca,'xtick',10:10:70) set(gca,'XTickLabel',{'1945','1955','1965','1975','1985','1995','2005'}); xlabel('年份') ylabel('周期尺度/a') title('(a)') %小波方差 figure var_wa=sum(abs(real(cwt1)).^2,2) plot(var_wa) set(gca,'xtick',10:10:70) set(gca,'XTickLabel',{'10','20','30','40','50','60','70'}); xlabel('周期/年') ylabel('小波方差') title('()')