%%%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