SVD算法的MATLAB代码

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