求解小波变换中模极大值的算法

[小文]

ddw=zeros(size(swd)
pddw=ddw
nddw=ddw
%小波系数的大于0的点
posw=swd.*(swd>0)
%斜率大于0
pdw=(posw( 1points-1)-posw(,2points)<0)
%正极大值点
pddw(,2 points-1)=(pdw(,1points-2)-pdw(,2points-1)>0)
%小波系数小于0的点
negw=swd.*(swd<0)
ndw=(negw(,1:points-1)-negw(,2points)>0)
%负极大值点
nddw(,2 points-1)=(ndw(,1points-2)-ndw(,2:points-1)>0)
%或运算
ddw=pddw|nddw;
ddw(,1)=1;
ddw(,points)=1;
%求出极值点的值,其他点置0
wpeak=ddw.*swd;
wpeak( ,1)=wpeak(,1)+1e-10;
wpeak( ,points)=wpeak(,points)+1e-10;

%画出各尺度下极值点（即模极大值）
figure;
for i=1:level
subplot(level,1,i)
plot(wpeak(i,) axis tight;grid on;
ylabel(strcat('j= ',num2str(i)
end
subplot(4,1,1)
title('ECG信号在j=1,2,3,4尺度下的小波系数的模极大值点')

interva2=zeros(1,points)
intervaqs=zeros(1,points)
Mj1=wpeak(1)
Mj4=wpeak(3)
请问这是如何实现求极大值和极小值后确定模极大值的？朋友们帮忙解释下