xiaobo.asv

基于hht的大地电磁信号资料处理MATLAB方法实现

function ii=xiaobo(x,t,deta)
%先得到信号X的连续小波分解系数，得到小波尺度图，三维小波系数图，二维小波系数图，频点图。
%load x
%t=0:0.001:1;
%deta=0.001;
wname = 'db10';
scales = [2:0.25:200];
TAB_PF = scal2frq(scales,wname,deta);

figure(1); 
coefs=cwt(x,scales,wname,'plot');
figure(2);
mesh(t,TAB_PF,abs(coefs));
figure(3);
contour(t,TAB_PF,abs(coefs));
figure(4);
plot(TAB_PF);
end