%%%%%%%%%%%%%%%基于線性調(diào)頻信號的ESPRIT算法
clear;
close?all;
j=sqrt（-1）;
c=3e8;
tao=10e-6;%采樣間隔
fs=100e6;%采樣頻率
M=8;
f1_l=29.9e6;
f1_h=30.1e6;%信號的低頻和高頻
f2_l=30.1e6;
f2_h=29.9e6;%此種情況下，中頻相同且不相關(guān)，且為窄帶信號
f1_c=（f1_l+f1_h）/2;%中頻，且相同
B1=f1_h-f1_l;%帶寬
kc1=B1/tao;%調(diào)制斜率
kc2=-B1/tao;
N_fs=fix（tao*fs）;%?采樣點數(shù)
t0=（0:N_fs-1）/fs;
snr=10;
theta=[-10?20];%入射角度
N_source=length（theta）;
lambda=c/f1_c;%波長及陣元間距
d=lambda/2;
dd=（0:M-1）*d;
mecnt=100;

F?=?fs/N_fs;
f_xlabe?=?（（0:N_fs-1）-N_fs/2）*F;

%%%%%%%%%%%流行矩陣
delay=zeros（MN_source）;
for?i=1:N_source
????delay（:i）=dd.‘*sin（theta（i）*pi/180）/c;
end
%?A=zeros（MN_source）;
%?for?i2=1:N_source
%?????A（:i2）=exp（j*2*pi*dd.‘*sin（theta（i）*pi/180）/c）;
%??delay（:i2）=dd.‘*sin（theta（i2）*pi/180）/c
%?end


%%%%%%%%%%%產(chǎn)生線性調(diào)頻信號
X1=zeros（MN_fs）;
X2=zeros（MN_fs）;
for?h1=1:M;
????t1=t0-delay（h11）;
????X1（h1:）=exp（j*2*pi*（f1_l*t1+1/2*kc1*t1.^2））;
end
for?h2=1:M;
????t2=t0-delay（h22）;
????X2（h2:）=exp（j*2*pi*（f2_l*t2+1/2*kc2*t2.^2））;
end

%?plot（f_xlabefftshift（abs（fft（X1（1:）N_fs））））;

X_out=X1+X2;
%?X=awgn（Xsnr‘measured‘）;

%%%%%%%%%%%ESPRIT算法
for?loop=1:mecnt

????X=awgn（X_outsnr‘measured‘）;
????Y1=X（1:M-1:）;
????Y2=X（2:M:）;%子陣1和2
????Y=[Y1;Y2];
????R=（Y*Y‘）/N_fs;%輸出陣列
????[edv?d1]=eig（R）;
????d11=diag（d1）;
????[d_sort?pre]=sort（abs（d11?）‘descend‘）;%按照降序排列
????edv_sort=edv（:pre）;%取出信號子空間
????edv1=edv_sort（:1:N_source）;
????Ex=edv1（1:M-1:）;
????Ey=edv1（M:2*M-2:）;%兩個子陣的信號子空間
????F?=?pinv（Ex）*Ey;??%%%%%%%%%子陣1的逆與子陣2的乘積為不變關(guān)系
????[edv2d2]=eig（F）;????????%?對F進行特征分解得到特征值
????doa（:loop）=?-asin（?lambda*angle（diag（d2））/（2*pi*d））*180/pi;?%角度估計
%?????doa_sort（:loop）?=?sort（doa（:loop）‘a(chǎn)scend‘）;
end
????doa_sort?=?sort（doa‘a(chǎn)scend‘）;
%------蒙特卡洛實驗下對不同角度估計值進行繪圖
%?
figure;
plot（doa_sort（1:）‘-b‘）;
hold?on;
plot（doa_sort（2:）‘-r‘）;
legend（‘-10°‘‘20°‘）;
xlabel（‘實驗次數(shù)‘）;
ylabel（‘doa估計‘）;
title（‘ESPRIT蒙特卡羅實驗‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2144??2018-03-10?10:42??ESPRIT.m
?????文件????????1862??2018-03-12?09:26??MUSIC.m