模极大值程序.txt

小波模极大源代码

function md=MaxMode(W,L) 

Wmax=0;%求取第一个小波变换的模极大值 
tfirst=0;%求第一个小波变换模极大值点的时刻 

for i=1:5001 
md(i)=0; 
end 

for(i=1:1:4009) 
if(abs(W(i))>abs(Wmax)) 
Wmax=W(i); 
tfirst=i; 
end 
end 
%求取小波变换模极大值点。 
Wflag=abs(0.02*Wmax);%设定小波变换模极大值的门坎值 
tflag=tfirst+2*L/0.28; 
for(i=1:1:tflag) 
if((W(i+1)>Wflag)&(W(i+1)-W(i)>Wflag)) 
md(i+1)=W(i+1); 
elseif((W(i+1)<-Wflag)&(W(i+1)-W(i)<-Wflag)) 
md(i+1)=W(i+1); 
else 
md(i+1)=0; 
end 

end 

for(i=1:1:tflag) 
if((md(i)>0)&(md(i+1)>md(i))) 
md(i)=0; 
elseif((md(i)<0)&(md(i+1)<md(i))) 
md(i)=0; 
end 
end






% Wflag=abs(0.02*Wmax)是阀值可以自己定义； 
% tflag=tfirst+2*L/0.28是选取信号分析区间，这也是应该个人根据需求而定，不要死搬硬套：）
% 我这个程序里W代表信号，L代表输电线路长度（是用来截取信号分析区间的）