✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
更多Matlab仿真内容点击👇
⛄ 内容介绍
与相控阵发射波束只有角度依赖性不同,频控阵(FDA)通过在阵元上附加一个频偏增量,使其具有距离-角度二维依赖特性,不仅解决了波束发射时的角度指向问题,还可以解决波束发射时的距离指向问题.通过在FDA中引入MIMO,可以将发射信息整合到接收端,从而实现更加灵活的距离维自由度利用,提升距离-角度二维联合信号处理的能力,这使得FDA-MIMO在抗距离维主瓣欺骗干扰方面具有独特的应用优势.
⛄ 代码
%——————————————-------------------跳频信号(Alpha稳定分布)参数估计,对于alpha>1,无需低阶,直接stft
clear,clc,close all
%------------------------------------产生跳频信号
fs=4*10^3; %采样频率
fk = [1.1 1.3 1.6 1.0 1.7 1.5 1.2 1.4]*10^3; %跳频图案集
th=50*10^(-3); %跳频周期,每跳驻留时间
ze = zeros(1,0); %0表示信号延迟为零
x = []; sig = ze;
t = 0:1/fs:th-1/fs;
for i = fk;
x = cos(2*pi*i*t);
sig = [sig x]; %sig为输出跳频信号
end
N = length(sig);
%----------------------------------加入Alpha稳定分布噪声
U=unifrnd(-pi/2,pi/2,1,1600);
W=exprnd(1,1,1600);
alpha=0.8;
X=(sin(alpha*U)./cos(U).^(1/alpha)).*(cos(U-alpha*U)./W).^((1-alpha)/alpha); %产生Alpha稳定分布
m=2;
sig_temp=sig.*(10.^m/20)./std(sig); %根据信噪比重新调整信号幅度
sig=sig_temp+X;
%------------------------------------------------改进后的低阶矩短时傅里叶变换
M=N/4+1; %汉宁窗窗长
W=hanning(M);
ye=zeros(N+M,1);
ye(M/2+1:M/2+N)=sig;
df=fs/M; %频率分辨率
fm=1900; %最高频率
L=round(fm/df); % the index of the max freq.
f=[0:L-1]*df;
i=1;
while i<=N
p=0.2; %低阶统计量的阶数
ye_p=abs(ye(i:i+M-1)).^p.*sign(ye(i:i+M-1)); %p阶矩运算的一种定义,
s=ye_p.*W;
sa=hilbert(s);
S=fft(sa);
B(:,i)=S(1:L);
i=i+1;
end
SP = (1/N)*abs(B.*B);
SP = SP/max(max(SP)); % 归一化
imagesc(SP)
set(gca,'YDir','normal');
title('FLOSTFT时频图');
xlabel('时域采样点','fontsize',14);
ylabel('频域采样点','fontsize',14);
⛄ 运行结果
⛄ 参考文献
[1] 郑木生. 基于Matlab语言实现电磁场中的FDTD算法编程[J]. 现代电子技术, 2005.
[2] 刘润东, 薛峰涛, 杨赟秀,等. FDA-MIMO雷达实值降维求根MUSIC参数估计算法[J]. 现代雷达, 2022, 44(4):7.
[3] 陈浩, 马建朝, 吕明久,等. 基于压缩感知的FDA-MIMO雷达波束形成算法[J]. 空军预警学院学报, 2020(004):034.
[4] 陈浩、马建朝、吕明久、谢谠、刘亚娜. 基于压缩感知的FDA-MIMO雷达波束形成算法[J]. 空军预警学院学报, 2020, 34(4):5.
⛳️ 代码获取关注我
❤️部分理论引用网络文献,若有侵权联系博主删除
❤️ 关注我领取海量matlab电子书和数学建模资料