雷达通信基于matlab粒子群算法优化综合线阵低副瓣方向图含Matlab源码 1962期

Posted 2022-07-12 海神之光

一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【雷达通信】基于matlab粒子群算法优化综合线阵低副瓣方向图【含Matlab源码 1962期】

获取代码方式2：
通过订阅紫极神光博客付费专栏，凭支付凭证，私信博主，可获得此代码。

备注：
订阅紫极神光博客付费专栏，可免费获得1份代码（有效期为订阅日起，三天内有效）；

二、部分源代码

%----------------主程序——------------
clc;
clear;
eps;
c=3e8;                     % 光速
fc=1.3e9;                  % 工作频率（hz）
numda=c/fc;                % 波长 wave length
N=16;                      % 阵列数
d=0.5*numda;               % 阵元间距
L=N*d;                     % 天线长
k=(2*pi)/numda;            % 波数
fs=10;                      % 采样频率
theta=-90:1/fs:90;    % 方位角度范围（采样范围）(rad)
Ns=length(theta);          % 采样点数
B=90;    %主波束宽度为90度
theta_d=0;                % 主波束方向
%---------------------适应值函数------------------------
%------即目标函数--------
fd=[10^(-1.5)*ones(1,800) ones(1,200) 10^(-1.5)*ones(1,801)];
fd=20*log10(fd);
plot(theta,fd,':');hold on;
xlim([-90 90]);
%---------------------粒子群算法------------------------
pop_size = 30;          %   pop_size 种群大小
part_size = N/2;            %   part_size 粒子大小，即阵元的馈电幅值的数目，阵元馈电对称，粒子数目只取一半。
a_max=1;a_min=0.3;  %幅度的范围
c1=1.49445;c2=1.49445;           %学习因子
v_max=1.0;           %最大速度
% w=1;               %惯性权重因子，为可变值，采用不同的权重因子可提高收敛速度。
NN=300;%迭代次数
    %---------------------初始寻优------------------------
    pop_a_ini=rand(pop_size,N/2);
    p_a_b=zeros(pop_size,N/2,NN+1); %粒子局部最优解，即每个粒子在寻优过程中的历史最优值。
    p_a_g=zeros(N/2,NN+1);%粒子全局最优解，所有粒子的最优值。
        %-------初始化粒子-------
        for i=1:pop_size
             f_ini_array(i,:)=present_array(Ns,d,theta,pop_a_ini(i,:));
             erro_ini(i)=array_erro(f_ini_array(i,:),fd);
        end

三、运行结果

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社，2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社，2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社，2018.

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除