SINR自适应波束形成器matlab代码

基于最大信干噪比算法的自适应波束形成器matlab代码

%%%16元线阵波束形成d=namda/2
clear
clc
format long;
j=sqrt（-1）;
v=1;
M=16;%%阵元个数
N=1000;%%%%%%%快拍数
sita0=0*pi/180; %信号方向
sita1=input（‘sita1=‘）*pi/180; %干扰方向
sita2=input（‘sita2=‘）*pi/180;
db1=20;%%%%干扰信号的功率db数
db2=20;
p1=1/（10^（db1*0.1））;%%%干扰信号的归一化功率
p2=1/（10^（db2*0.1））;
%%%期望信号
x1=randn（1N）;
x2=randn（1N）;
adt=x1+j*x2;
%%%干扰信号1
x1=0.5*sqrt（p1）*randn（1N）;
x2=0.5*sqrt（p1）*randn（1N）;
s1=x1+j*x2; 
%%%%干扰信号2
x1=sqrt（0.5*p2）*randn（1N）;
x2=sqrt（0.5*p2）*randn（1N）;
s2=x1+j*x2; 
i=1:1:M;%%%%%%%%%%%%信号的导向矢量：线阵的形式
a0（i1）=exp（-j*（i-1）*pi*sin（sita0））;
a1（i1）=exp（-j*（i-1）*pi*sin（sita1））;
a2（i1）=exp（-j*（i-1）*pi*sin（sita2））;
nr=0.01*randn（MN）;%%%产生M行N列的复噪声矩阵 
ni=0.01*randn（MN）; 
n=nr+j*ni;
R=zeros（MM）;
%%阵列的协方差矩阵
for t=1:N
 x=s1（t）*a1+s2（t）*a2+n（:t）; %阵列对干扰的完整响应
     R=R+x*x‘;%所有时刻干扰的协方差矩阵之和
end
R=R/N;%%%%%%%%%协方差矩阵，所有快拍数的平均
w=pinv（R）*a0;%%%最大SINR算法
%%%%形成波束
for sita=-90:0.5:90;
    for i=1:M
        a（i1）=exp（-1j*（i-1）*pi*sin（sita*pi/180））;
    end
    y（v）=w‘*a;%%%%%%%对信号进行加权，消除干扰
    v=v+1;
end
y_max=max（abs（y））;%%%%%%%%%%%%%%%归一化
y_1=abs（y）/y_max;
y_db=20*log（y_1）;
sita=-90:0.5:90;
plot（sitay_db）;
axis（[-9090-1500]）;
Xlabel（‘sita‘）;
Ylabel（‘天线增益db‘）;

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

     文件       1371  2011-06-28 16:31  SINR.m

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

                 1371                    1