求各位大神指导【小波分析】程序

24444

有小波分析的程序，不知道对错，请各位大神指正

f=xlsread('year1.xls','sheet1','B1:B60');%读取原始数据

t= xlsread('year1.xls','sheet1','A1:A60');%读取时间

figure(1);title('Stereoscipicgraphy for coef');

subplot(1,1,1);c1=cwt(f,[1:1:60],'morl','plot');%连续的Molet小波

[c,h]=contourf(c1,13);colormap('default');

clabel(c,h)

xlabel('年份');ylabel('周期/年');zlabel('MORLET小波变换系数Coef');title('MORLET时频分布图');

set(gca,'XTickLabel',{1955,1960,1965,1970,1975,1980,1985,1990,1995,2000,2005,2010,2014})

%colorbar('horiz');

print(gcf,'-djpeg','MORLET时频分布图');

figure(2);%绘制小波时频分析图，实际就是coef1/coef2系数等高线图

subplot(1,1,1);

[c,h]=contourf(c1,13);colormap(hsv);

clabel(c,h)

xlabel('年份');ylabel('周期/年');title('MORLET小波实部时频分布等高线');

set(gca,'XTickLabel',{1955,1960,1965,1970,1975,1980,1985,1990,1995,2000,2005,2010,2015})

%colorbar('horiz');

print(gcf,'-djpeg','MORLET时频Coef时频分布分布等高线');

disp(c1);

save('yearMORL.txt','c1','-ASCII');%储存数据

figure(3);

24444

没有人吗