%%%LCMV在多个来波方向约束下波束形成%%%
clc;clear all;close all;ima=sqrt(-1);esp=0.01;
%%天线参数设定%%
N=16; %阵元数
d_lamda=0.5; %阵元间距与波长的比值
theta=-90:0.5:90; %搜索范围确定
theta1=-10; %来波方向1
theta2=0; %来波方向2
theta3=40; %来波方向3
theta_jam=70; %干扰方向
L=512; %采样点数
%%%%%%%%%%%%%%%%%%%%%%%%%
%%信号形成%%
for k=1:L;
a1=10*randn(1);
a2=10*randn(1);
a3=10*randn(1);
ajam=10*randn(1);
an=1;
s(:,k)=a1*exp(ima*2*pi*d_lamda*sin(theta1*pi/180)*[0:N-1]')+...
+a2*exp(ima*2*pi*d_lamda*sin(theta2*pi/180)*[0:N-1]')+...
+a3*exp(ima*2*pi*d_lamda*sin(theta3*pi/180)*[0:N-1]');
jam(:,k)=ajam*exp(ima*2*pi*d_lamda*sin(theta_jam*pi/180)*[0:N-1]');
n(:,k)=an*(randn(N,1)+ima*randn(N,1));
end
%%%%%%%%%%%%%%%%%%%%%%%%%
%最优权矢量产生%
x=s+jam;
Rx=1/L*x*x'; %求信号相关矩阵
R=pinv(Rx); %相关矩阵求逆
a1theta=exp(ima*2*pi*d_lamda*sin(theta1*pi/180)*[0:N-1]');
a2theta=exp(ima*2*pi*d_lamda*sin(theta2*pi/180)*[0:N-1]');
a3theta=exp(ima*2*pi*d_lamda*sin(theta3*pi/180)*[0:N-1]');
C=[a1theta a2theta a3theta];%方向矩阵
F=[1 1 1]';
Wopt=R*C*(inv(C'*R*C))*F;
%%%%%%%%%%%%%%%%%%%%%%%%%
%%最优波束形成%%
for m=1:length(theta);
a=exp(ima*2*pi*d_lamda*sin(theta(m)*pi/180)*[0:N-1]');
y(m)=Wopt'*a;
end
%%%%%%%%%%%%%%%%%%%%%%%%%
Y=20*log10(abs(y)/max(abs(y))+esp);
%%作图%%
plot(theta,Y);hold on;grid on;
axis([-90 90 -50 0]);
plot(theta1,-30:0,'.');
plot(theta2,-30:0,'.');
plot(theta3,-30:0,'.');
plot(theta_jam,-30:0,'.');
xlabel('\theta/o');
ylabel('Amplitude in dB');
title('MVDR波束形成');
MVDR_mvdr_MVDRmatlab_
版权申诉
5星 · 超过95%的资源 137 浏览量
2021-09-30
01:29:05
上传
评论 1
收藏 54KB RAR 举报
周玉坤举重
- 粉丝: 61
- 资源: 4782
最新资源
- Python爬取淘宝热卖商品并可视化分析
- 5152单片机proteus仿真和源码将按键次数写入AT24C02再读出并用1602LCD显示
- SE-SSD复现过程(Det3D的安装教程)
- 基于Python的在线学习与推荐系统设计与实现(论文+源码)-kaic
- 串口通过 YMODEM 协议进行文件传输
- 蓝桥杯2024年第十五届省赛真题-前缀总分
- com.qihoo.appstore_300101305-1.apk
- tensorflow-gpu-2.7.1-cp37-cp37m-manylinux2010-x86-64.whl
- tensorflow-2.7.2-cp37-cp37m-manylinux2010-x86-64.whl
- tensorflow-2.7.1-cp39-cp39-manylinux2010-x86-64.whl
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈