wavelet.m

matlab对测井曲线的小波分析变换进行分层处理

x=0:1:38;
y=[3.2,3.1,2.9,3.1,3.6,3.0,3.3,3.2,3.9,7.6,7.9,9.4,8.9,11.5,...
    16.8,25.7,29.6,30.9,30.7,30.4,33.3,29.5,13.0,9.1,8.2,...
    12.5,13.2,14.4,14.0,9.0,4.0,3.6,3.1,3.4,4.2,3.7,2.7,2.4,2.2];
%原始数据
h=1/sqrt(2*pi)*(1-x^2)*exp(-x^2/2);  %小波基函数 (Mexican Hat)
hj=1/sqrt(2*pi)*(1-(x/2^j)^2)*exp(-(x/2^j)^2/2);
%h关于尺度2^j的伸缩变换 hj=h(x/2^j)/2^j
%求f(x)*hj(x)
%j=1
h1=1/sqrt(2*pi)*(1-(x/2)^2)*exp(-(x/2)^2/2)
conv(y,h1)