function?[?snr?]?=?svd1（?x_n?）
%UNtitleD4?Summary?of?this?function?goes?here
%???Detailed?explanation?goes?here
m=50;
Lk=length（x_n）;

%{
x_n=zeros（1Lk）;
for?i=1:Lk
????x_n（i）=abs（x_n1（i））;
end
%}

ex1=zeros（1m）;
ex2=zeros（1m）;
for?k=1:m
????j=Lk-（k-1）;
????for?i=1:j
????????ex1（k）=ex1（k）+x_n（i）;
????????ex2（k）=ex2（k）+x_n（i+k-1）;
????end
????ex1（k）=ex1（k）/j;
????ex2（k）=ex2（k）/j;
end
%}
r=zeros（1m）;
for?k=1:m
????j=Lk-（k-1）;
????for?i=1:j
????????r（k）=r（k）+（x_n（i）-ex1（k））*conj（x_n（i+k-1）-ex2（k））;
????end
????r（k）=r（k）/j;
????
end


r1=zeros（1m）;
r1（1）=r（1）;
for?k=2:m
????j=Lk-（k-1）;
????for?i=1:j
????????r1（k）=r1（k）+conj（x_n（i）-ex1（k））*（x_n（i+k-1）-ex2（k））;
????end
????r1（k）=r1（k）/j;
end
%}
r_metrix=zeros（mm）;
%r_metrix=toeplitz（r1r）;

r_metrix（1:）=r;
for?h=2:m
???for?f=1:h-1
???????r_metrix（hf）=r1（h+1-f）;
???end
???for?f=h:m
???????r_metrix（hf）=r（f-h+1）;
???end
end
%}
%{
for?i=1:m
????for?