load initial_track s; % y:initial data,s:data with noiseT=0.1;
% yp denotes the sample value of position% yv denotes the sample value of velocity% Y=[yp(n);yv(n)];% err