%最小二乘法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