matlab中如何使用小波分析以及小波分析在matlab中得源程序

凌晨24点

matlab中如何使用小波分析以及小波分析在matlab中得源程序！谢谢

edwardli

function [wave,period,scale,coi] = ...
        wavelet(Y,dt,pad,dj,s0,J1,mother,param);

if (nargin < 8), param = -1;, end
if (nargin < 7), mother = -1;, end
if (nargin < 6), J1 = -1;, end
if (nargin < 5), s0 = -1;, end
if (nargin < 4), dj = -1;, end
if (nargin < 3), pad = 0;, end
if (nargin < 2)
        error('Must input a vector Y and sampling time DT')
end

n1 = length(Y);

if (s0 == -1), s0=2*dt;, end
if (dj == -1), dj = 1./4.;, end
if (J1 == -1), J1=fix((log(n1*dt/s0)/log(2))/dj);, end
if (mother == -1), mother = 'MORLET';, end

%....construct time series to analyze, pad if necessary
x(1:n1) = Y - mean(Y);
if (pad == 1)
        base2 = fix(log(n1)/log(2) + 0.4999);   % power of 2 nearest to N
        x = [x,zeros(1,2^(base2+1)-n1)];
end
n = length(x);

%....construct wavenumber array used in transform [Eqn(5)]
k = [1:fix(n/2)];
k = k.*((2.*pi)/(n*dt));
k = [0., k, -k(fix((n-1)/2):-1:1)];

%....compute FFT of the (padded) time series
f = fft(x);    % [Eqn(3)]

%....construct SCALE array & empty PERIOD & WAVE arrays
scale = s0*2.^((0:J1)*dj);
period = scale;
wave = zeros(J1+1,n);  % define the wavelet array
wave = wave + i*wave;  % make it complex

% loop through all scales and compute transform
for a1 = 1:J1+1
        [daughter,fourier_factor,coi,dofmin]=wave_bases(mother,k,scale(a1),param);       
        wave(a1,:) = ifft(f.*daughter);  % wavelet transform[Eqn(4)]
end

period = fourier_factor*scale;
coi = coi*dt*[1E-5,1:((n1+1)/2-1),fliplr((1:(n1/2-1))),1E-5];  % COI [Sec.3g]
wave = wave(:,1:n1);  % get rid of padding before returning

return

mofangbao

fortran区看看去吧 里面的小波的帖子里就有matlab的小波

尽头的尽头

fortran区有小波分析的

edwardli

忘了说了  其实不用什么源程序   我自己分析的时候首先是从小波工具包入手的  这个是可视化的  在matlab窗口中输入wavemenu就出来了。。。闹明白这个之后我就在窗口写了些语句  比如我要分析的数列的处理z，之后的命令是matlab中的函数，例如result=comr（z，1，2,0.5）  具体的函数可以从matlab help里面搜一下

凌晨24点

edwardli 发表于 2011-9-21 15:33

                               
登录/注册后可看大图

忘了说了  其实不用什么源程序   我自己分析的时候首先是从小波工具包入手的  这个是可视化的  在matlab窗口 ...

谢谢，兄台

凌晨24点

edwardli 发表于 2011-9-21 06:49

                               
登录/注册后可看大图

function [wave,period,scale,coi] = ...
        wavelet(Y,dt,pad,dj,s0,J1,mother,param);

凌晨24点

尽头的尽头 发表于 2011-9-20 23:06

                               
登录/注册后可看大图

fortran区有小波分析的

恩，下到了

凌晨24点

mofangbao 发表于 2011-9-20 16:23

                               
登录/注册后可看大图

fortran区看看去吧 里面的小波的帖子里就有matlab的小波

都是PHP文件，打不开啊

尽头的尽头

凌晨24点 发表于 2011-9-21 17:10

                               
登录/注册后可看大图

都是PHP文件，打不开啊

不要用迅雷下载，直接另存为