om(1)=24*2*pi/12.4206;
 om(2)=24*2*pi/12.66;
 om(3)=24*2*pi/12.00;
 om(4)=24*2*pi/23.96;
 om(5)=24*2*pi/25.86;
 om(6)=om(1)*2;
 om(7)=om(1)*3;
 om(8)=om(1)*4;

OM=2*pi/32;
NS=71;
[tlow xlow Hhi]=lowpass(t,H,OM,NS);
H=interp1(t,H,tlow);
t=tlow;

[x]=lsq1(om,Hhi,t);
nc=length(om);
X(1:length(t))=x(1);
for i=1:nc;
    X=X+x(2*i)*cos(om(i)*t) + x(2*i+1)*sin(om(i)*t);
end

H=Hhi-X;