clearclcclose?all;

ss=importdata（‘6.mat‘）;
s=awgn（ss60‘measured‘）;

%?t=0:1e-3:1;
%?s=2*sawtooth（2*pi*15*t）;%+4*sin（2*pi*10*t）.*sin（2*pi*0.1*t）+sin（2*pi*5*t）;

%?t=0:1e-3:1;
%?s1=2*sin（2*pi*5*t）;
%?s2=0.3*sin（2*pi*100*t）;
%?s2（1:450）=zeros（1450）;
%?s2（550:end）=zeros（1452）;
%?s=s1+s2;

sbk=s;
n=length（s）;x=1:n;
q=0;
while（1）
q=q+1;
qq=0;
while（1）
qq=qq+1;
j=0;k=0;p=0;
for?i=2:（n-1）
????if（s（i）>s（i+1）?&&?s（i）>s（i-1））
????????j=j+1;
????????jd（j）=s（i）;xd（j）=i;
????end
????if（s（i）????????k=k+1;
????????jx（k）=s（i）;xx（k）=i;
????end
????if（s（i-1）*s（i）<=0）
????????p=p+1;
????end
end
%?if（abs（j+k-p）<=1?||?qq>10）
%?????clear?jd?jx?xd?xx?ub?db?mid;
%?????break;
%?end

if?s（1）>s（2）
????dlp=s（1）;xlp=-s（1）;