ESPRIT算法MATLAB程序

% ESPRIT
function [R_estS_est] = ESPRIT（datanum_tar）
global N_per_frame;
% global n_bomen;
% global MUSIC_start;
% global MUSIC_end;
global fs_s;

c =3e8;
delt_freq = 6.25e6;
% delt_freq = 0;
freq_0 = 0;
z_N = data;

p = 50;
d = num_tar;
N_BP = N_per_frame/2 - p + 1;

R_zz = zeros（pp）;

for k = 1:N_BP
    R_k = z_N（k:k+p-1:）*z_N（k:k+p-1:）‘;
    R_k_x = conj（R_k）;                  
    R_k_x_h =  conj（fliplr（conj（fliplr（R_k_x）‘））‘）;
    R_zz = R_zz + （R_k + R_k_x_h）;               % MSSP 改进型空间平滑处理
%   R_zz = R_zz + R_k ;                        % SSP  空间平滑处理  
end

R_zz = R_zz/（2*N_BP）;
% R_zz = z_N*z_N‘;
[V_R_zz D_R_zz] = eig（R_zz）;

if D_R_zz（11）>D_R_zz（pp）
    E_s = V_R_zz（:1:d）;
else