微信公众号:EW Frontier
关注可了解更多的雷达、通信、人工智能相关代码。问题或建议,请公众号留言;
如果你觉得EW Frontier对你有帮助,欢迎加入我的知识星球或面包多,更多代码等你来学
知识星球:https://wx.zsxq.com/dweb2/index/group/15552518881412
面包多:https://mbd.pub/o/author-a2mYl2tsbA==/workQQ交流群:729981694
如有侵权请联系删除~
摘要
关键词:多载频MIMO雷达;接收波束形成;虚拟阵列;旁瓣抑制
基于凸优化理论的旁瓣抑制方法
在众多优化方法中,凸优化[45,46]方法近年来备受关注,成为研究的热 点。该优化方法被应用在控制系统、机械工程、信号处理等很多工程领域中。本节将凸优化理论应用于多载频雷达的方向图综合问题中,来降低波束形成 方向图的旁瓣电平。在保持主瓣分辨率不变的条件下,使方向图的旁瓣电平 最低。
MATLAB代码展示
% L波段 (1--2 GHz) 波长 30-15 cm 远程监视(200-500 km)
% S波段 (2--4 GHz) 波长 15-7.5 cm 中程监视(100-200 km) 远程跟踪(50-150km)
% 参考文献 多载频MIMO高频雷达的波束形成方法研究 候宪美
%图3-8
clear
clc
%% 参数 %%
f0 = 3e9;%参考频率
f1 = f0+32e6;
f2 = f1+32e6;
f3 = f2+32e6;
f4 = f3+32e6;
c = 3e8;
Lambda = c/f0;
M = 5;%发射阵元数
N = 13;%接收阵元数
dr = Lambda/2;
dt = 1*dr;
Lambda1 = c/f0;
Lambda2 = c/f1;
Lambda3 = c/f2;
Lambda4 = c/f3;
Lambda5 = c/f4;
theta = -90:90;
theta0 = 0;
%% 方法1 用循环 %%
%% 接收导向矢量 %%
% ar = exp(-j*2*pi*(N-1)*d*sind(theta)/Lambdai)
%接收导向矢量与发射的频率相关
ar1 = zeros(N,1);
katar1 = zeros(N,length(theta));
for jj = 1:length(theta)
for ii = 1:N
ar1(ii,:) = exp(-1i*2*pi*(ii-1)*dr*sind(theta(jj))/Lambda1); % N*1
end
katar1(:,jj) = kron(1,ar1);
end
ar2 = zeros(N,1);
katar2 = zeros(N,length(theta));
for jj = 1:length(theta)
for ii = 1:N
ar2(ii,:) = exp(-1i*2*pi*(ii-1)*dr*sind(theta(jj))/Lambda2); % N*1
end
katar2(:,jj) = kron(exp(1i*2*pi*dt*sind(theta(jj))/Lambda2),ar2);
end