自适应波束形成的经典算法大全

自适应波束形成的经典算法，有 RLS, LMS以及各算法之间的性能比较和分析

clear all
close all
clc
j=sqrt（-1）;
M=10;
m=0:M-1;
theta=[1040]; %信号源10度，干扰40度
theta=theta*pi/180;
us=0;
A0=exp（j*pi*m‘*sin（theta（1）））;%信号源的方向矢量 
Ai=exp（j*pi*m‘*sin（theta（2）））;%干扰的方向矢量
SNR=20;  % 信噪比（可变）
INR=1; 
P0=10.^（SNR/10）; 
Pi=10.^（INR/10）;
o=pi/180;
xxl=[-90:1:90];
xx=xxl*o;
a=exp（j*pi*m‘*sin（xx））;       % 扫描矢量;
%%%%%%%%%%%最优波束形成%%%%%%%%%%%%%%%%%%%%%%%
L=5000;
Xs=A0*sqrt（P0*2）*wgn（1L1‘complex‘）;%期望
Xi=Ai*sqrt（Pi*2）*wgn（1L1‘complex‘）;%干扰
Xn=sqrt（2）*wgn（ML1‘complex‘）;%噪声
Xx=Xs+Xi+Xn;
 %Xx=Xi+Xw;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure
Rs=P0*A0*A0‘;                                      %期望信号协方差矩阵
Ri=Pi*Ai*Ai‘;                                      %干扰信号协方差矩阵
Rn=eye（M）;                                         %噪声信号协方差矩阵
Rx=Rs+Ri+Rn;    
Pout_min=pinv（A0‘*pinv（Rx）*A0）;
Wopt=Pout_min*pinv（Rx）*A0;
%Wop=pinv（Rx）*A0/（A0‘*pinv（Rx）*A0）;
s1=20*log10（abs（Wopt‘*a）/max（abs（Wopt‘*a）））;
plot（xxls1‘--r‘）
grid on
%%%%%%%%%%%%%%%LMS算法%%%%%%%%%%%%%%%%%%
w=zeros（M1）;
eta=0.01;
miu=0.0001;
e=1;
p=1; 
d=Xs（1:）;
while abs（e）>eta 
 e=d（p）-w‘*Xx（:p）;
 w=w+miu*conj（e）*Xx（:p）;
 e1（p）=e;
 p=p+1;
end
s2=20*log10（abs（w‘*a）/max（abs（w‘*a）））;
hold on
plot（xxls2‘-b‘‘linewidth‘1.5）;
legend（‘最优情况‘‘LMS算法‘）;
xlabel（‘到达角度‘）;
ylabel（‘输出幅度增益/dB‘）;

for k = 1:p-1
    en（k）=（abs（e1（k）））.^2;
end
figure
semilogy（en）;
xlim（[0 500]）
% axis（[0 500 0.001 10]）
xlabel（‘样本数‘）;
ylabel（‘误差e^{2}（n）‘）;
% for k = 1:p-1
%     en（k）=sum（abs（e1（1:k））.^2）/k;
% end
% figure
% semilogy（en）;
% ylim（[0 1000]）
% % axis（[0 500 0.001 10]）
% xlabel（‘样本数‘）;
% ylabel（‘均方误差E[e^{2}（n）]‘）;

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

     文件       1925  2014-05-08 22:30  SMI.m

     文件       1755  2014-05-08 22:30  LMS.m

     文件       1808  2014-05-08 22:30  RLS.m

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

                 5488                    3