将原函数的两个端点作为极值怎么画包络线？

李清志

原函数的EMD分解时的包络线程序如下：
t=0:1:500;
s=cos(0.04*pi*t)+0.6*cos(0.08*pi*t)+0.5*sin(0.01*pi*t)
s1=s;
sd=1;
for n=1:3
if(sd>0.3)
   jd=find(diff(sign(diff(s1)))==-2)+1;  %找h1极大值的位置
   jx=find(diff(sign(diff(s1)))==2)+1; %找h1极小值的位置
   ah=length(jd);
   al=length(jx);
   jdz(ah)=0;
   for i=0:(ah+1)
       bh=jd(i);
       jdz(i)=s1(bh);
   end                      %找出h1极大值对应函数值
   jxz(al)=0;
   for i=1:al
       bl=jx(i);
       jxz(i)=s1(bl);
   end                     %找出h1极小值对应函数值
   jsbl=spline(jd,jdz,t);  %极大值拟和的上包络
   jxbl=spline(jx,jxz,t);  %极小值拟和的下包络
   m1=(jsbl+jxbl)/2;       %上下包络均值
   h1=s1-m1;
   sd=sum(((s1-h1)./h1).^2);
end
end
[indmin,indmax,indzer]=extr(s,t)
tt1=indmax
tt2=indmin
max=s(indmax)
min=s(indmin)
tt=indmax(1):indmax(end)
TT=indmin(1):indmin(end)
max1=spline(tt1,max,tt)
min1=spline(tt2,min,TT)
[up,down] = envelope_s(t,s1,'linear')
figure(1)
plot(s1,'c')
hold on
plot(jsbl,'r')
plot(jxbl,'r')
plot(jd,jdz,'o')
plot(tt,max1,'m')
plot(TT,min1,'m')
plot(t,up,'b-.');
plot(t,down,'b-.');
hold off;
要想把原函数的两个端点作为极值点程序该怎么写啊？要想把延拓后的极值点包进去又该怎么写啊？求高手指点。