ESPRIT算法255749

function estimate =  tls_esprit（ddcrLe）
%*******************************************************
% This function calculates TLS-ESPRIT estimator for
% uniform linear array.
%
% Inputs
%    dd       sensor separation in wavelength
%    cr（KK）  array output covariance matrix
%    Le       estimated number of sources ==L=iwave
%   
% Output
%    estimate  estimated angles in degrees
%              estimated powers
%*******************************************************

twpi = 2.0*pi;
derad = pi / 180.0;
radeg = 180.0 / pi;

% eigen decomposition of cr
[KKK] = size（cr）;
[VD]=eig（cr）;
EVA = real（diag（D）‘）;
[EVAI] = sort（EVA）;
%disp（‘Eigenvalues of correlation matrix in TLS-ESPRIT:‘）
EVA=fliplr（EVA）;
EV=fliplr（V（:I））;

%  composition of E_{xy} and E_{