单元均匀线阵自适应波束形成图.doc

单元均匀线阵自适应波束形成图
clear
clc
format long;
v=1;
M=4;
N=1000;
f0=21*10^3;
f1=11*10^3;
f2=15*10^3;
omiga0=2*pi*f0;
omiga1=2*pi*f1;
omiga2=2*pi*f2;
sita0=; %信号方向
sita1=; %干扰方向1
sita2=; %干扰方向2
for t=1:N
adt(t)=sin(omiga0*t/(N*f0));
a1t(t)=sin(omiga1*t/(N*f1));
a2t(t)=sin(omiga2*t/(N*f2));
end
for i=1:M
ad(i,1)=exp(j*(i-1)*pi*sin(sita0));
a1(i,1)=exp(j*(i-1)*pi*sin(sita1));
a2(i,1)=exp(j*(i-1)*pi*sin(sita2));
end
R=zeros(M,M);
for t=1:N
x=adt(t)*ad+a1t(t)*a1+a2t(t)*a2; %阵列对信号的完整响应 R=R+x*x';%信号的协方差矩阵
end
R=R/N;
miu=1/(ad'*inv(R)*ad);
w=miu*inv(R)*ad;
for sita=0:pi/100:pi
for i=1:M
x_(i,1)=exp(j*(i-1)*pi*sin(sita));
end
y(1,v)=w'*x_;
v=v+1;
end
y_max=max(y(:));
y_1=y/y_max;
y_db=20*log(y_1);
sita=0:pi/100:pi;
plot(sita,y)
Xlabel(‘sitaa’)
Ylabel(‘天线增益db’)
4单元均匀线阵自适应波束形成目标
clear
clc
format long;
v=1;
M=4;阵元数
N=100;
f0=21*10^3;
omiga0=2*pi*f0;
sita0=;%信号方向
for t=1:N
adt(t)=sin(omiga0*t/(N*f0)); end
for i=1:M
ad(i,1)=exp(j*(i-1)*pi*sin(sita0)); end
R=zeros(4,4);
r=zeros(4,1);
for t=1:N
x=adt(t)*ad;
R=R+x*x.';
end
R=R/N;
miu=1/(ad.'*inv(R)*ad);
w=miu*inv(R)*ad;
for sita=0:pi/100:pi/2
for i=1:M
a(i,1)=exp(j*(i-1)*pi*sin(sita));
end
y(1,v)=w.'*a;
v=v+1;
end
sita=0:pi/100:pi/2;
plot(sita,y)
xlabel('sita')
ylabel('天线增益’)
这是程序全部,有的比较简单的部分就可以省事的,
如果解释的好,再加分的,大大的。
谢谢了
我有更好的答案
提问者采纳
2010-05-04 18:08
你这里有两个程序,第二个程序与第一个实质上是一样的,区别就是信号与导向矢量的写法有点不同,这里我就不注释了。还有,我下面附了一段我自己的写的程序,里面有SIM算法。G-S正交化算法等。是基于圆阵形式的,你的算法是基于线阵的,他们程序上的区别在于导向矢量的不同。我的算法是某项目中的,保证好使。建议学****波束形成技术,注意把程序分块,例如分成,求导向矢量;最优权值;形成波束等等。
程序如下:
clear
clc
format long;
v=1;
M=4;
N=1000;%%%%%%%快拍数
f0=21*10^3;%%%%%%%%%%%信号与干扰的频率
f1=11*10^3;
f2=15*10^3;
omiga0=2*pi*f0;%%%%%%%信号与干扰的角频率
omiga1=2*pi*f1;
omiga2=2*pi*f2;
sita0=; %信号方向
sita1=; %干扰方向1
sita2=; %干扰方向2
for t=1:N %%%%%%%%%%%%信号
adt