function?[sdsdby]=svddec（x）
N=length（x）;
if?mod（N2）
????m=（N+1）/2;
????n=（N+1）/2;
else?
????m=N/2+1;
????n=N/2;
end
for?i=1:m
????for?j=1:n
????????D（ij）=x（i+j-1）;
????end
end
[usv]=svd（D）;
sd=diag（s）‘;
for?i=1:length（sd）-1
sdb（i）=sd（i）-sd（i+1）;
end
[kl]=max（sdb）