close all; clear;
file_name='f01';% 文件名为***.txt
d=importdata([file_name,'.txt']);%数据形式为N*1
figure(1);plot(d);
% wavelet analysis
d=d-mean(d); % normalized to zero
mother = 'MORLET';% Morlet
param=2*pi;%母小波参数
dj0=0.6;
pad=1;
dt=1/32;%抽样时间
dj=1/32;%离散尺度之间的间隔
s0=12*dt; %小波的最小尺度
J0=12;
smax=200;
J1=fix(log2(smax/s0)/dj);
[wave,period,scale,coi]=wt1(d,dt,pad,dj,s0,J0,J1,mother,param);
fre=fliplr(1./period);amp=abs(fliplr_row_matrix(wave));
ind=find(fre<=2 & fre>=0.0095);fre=fre(ind);amp=amp(ind,:);
% view the matrix
figure(2)
t=0:size(amp,2)-1;
pcolor(t,fre,amp);set(gca,'yscale','log');shading interp;axis tight;colorbar
这里你采用了函数fliplr_row_matrix
fre=fliplr(1./period);amp=abs(fliplr_row_matrix(wave));
如果fliplr_row_matrix不是内置函数的话,说明你的程序差以这个函数名的函数m文件夹,需要编写fliplr_row_matrix函数m文件夹
fliplr_row_matrix你没有这个函数