经典Esprit空间谱估计算法的代码

经典Esprit空间谱估计算法的代码 适合初学空间谱估计的人

%最小二乘法LS ψ=US1+*US2
clear all;
close all;
dd = 0.35;               %阵元间距
kelm = 8;   
d = 0:dd:（kelm-1）*dd;   
rad=pi/180;             %角度变弧度
snr = 15;               % 信噪比
n = 400;                % 采样数目
iwave = 4;              
theta = [-30 -40 55 68];    
S = randn（iwaven）;     
A = exp（-1i*2*pi*d.‘*sin（theta*rad））;    
X = A*S;
X1 = awgn（Xsnr‘measured‘）; 
R = X1*X1‘/n;                    % 求解接收信号的协方差矩阵
[EVD] = eig（R）;
EVA = （diag（D））.‘;                  % 提取特征值 转置是因为diag得到的是列矢量
% Rs = R-sum（EVA（11:kelm-iwave））/（kelm-iwave）*eye（kelm）;  %去噪
% [EVD] = eig（Rs）;
% EVA = （diag（D））.‘;                  % 提取特征值 转置是因为diag得到的是列矢量
[EVAI] = sort（EVA）;                % 对特征值进行排序
EVA = fliplr（EVA）;                  % 将特征值从小到大排列
EV = fliplr（EV（:I））;                 % 对特征矢量也进行相应的排序

% LS
Us=EV（:1:iwave）;                          %信号子空间
Us1 = Us（1:kelm-1:）;              % Us1取信号子空间前kelm-1行
Us2 = Us（2:kelm:）;                % Us2取信号子空间后kelm-1行
U = pinv（Us1）*Us2;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1238  2015-01-26 14:17  ESPRIT\LS\LS.m

     文件        704  2015-01-07 16:39  ESPRIT\TAM（低噪比时性能不行）\TAM.m

     文件       1754  2015-01-12 10:18  ESPRIT\TLS（低噪比时也不错）\TLS.m

     文件       1090  2015-01-07 20:48  ESPRIT\矩阵束Esprit（已去噪）\array_Esprit.m

     目录          0  2015-02-01 15:04  ESPRIT\LS

     目录          0  2015-02-01 15:04  ESPRIT\TAM（低噪比时性能不行）

     目录          0  2015-02-01 15:04  ESPRIT\TLS（低噪比时也不错）

     目录          0  2015-02-01 15:04  ESPRIT\矩阵束Esprit（已去噪）

     目录          0  2015-02-01 15:05  ESPRIT

----------- ---------  ---------- -----  ----

                 4786                    9