targetechosimulation.m

一个雷达回波模拟的matlab程序包,含有在RCS服从卡方分布,瑞利分布,对数正态分布下目标的的雷达散射截面模拟子程序

function [targetecho]=targetechosimulation(Num_sh,Num_T,samp_T,K,tagetdata,RCS,f0,T,c)
%功能：目标回波仿真
%输入参数：Num_sh船的数量
%         Num_T扫频周期个数
%         samp_T一个扫频周期内所取样本点数(第一次FFT的采样点数，对应距离元)
%         K扫频斜率
%         a目标径向加速度
%         v目标径向速度
%         r目标径向距离
%         RCS目标雷达散射截面
%         f0雷达载频
%         T扫频周期
%         c光速
%输出参数：targetecho目标的雷达回波数据
t_withinT=(0:samp_T-1)*T/samp_T;
targetecho=zeros(Num_T,samp_T);
sr=zeros(Num_T,samp_T,Num_sh);
a=zeros(Num_T,Num_sh);                           %每条船的径向加速度在每个扫频周期之间满足[-1,1]均匀分布(m/s^2)，一个扫频周期内加速度不变%
v=zeros(Num_T,Num_sh);                           %船的径向速度矩阵初始化
r=zeros(Num_T,Num_sh);                           %船的径向距离矩阵初始化 
sita=zeros(Num_T,Num_sh);                        %船得方位角
for L=1:Num_sh
    r(:,L)=tagetdata(:,1,L);
    v(:,L)=tagetdata(:,3,L);
    a(:,L)=tagetdata(:,4,L);
    sita(:,L)=tagetdata(:,2,L);
    for m=1:Num_T
        sr(m,:,L)=RCS(m,L)*exp(2*pi*j*(2*K*r(m,L)*t_withinT/c-2*K*v(m,L)*t_withinT.^2/c-K*a(m,L)*t_withinT.^3/c-2*f0*r(m,L)/c+2*f0*v(m,L)*t_withinT/c+a(m,L)*f0*t_withinT.^2/c-2*K*r(m,L)^2/c^2-2*K*v(m,L)^2*t_withinT.^2/c^2-K*a(m,L)^2*t_withinT.^4/(2*c^2)+4*K*r(m,L)*v(m,L)*t_withinT/c^2+2*K*r(m,L)*a(m,L)*t_withinT.^2/c^2-2*K*a(m,L)*v(m,L)*t_withinT.^3/c^2));
    end
    targetecho=targetecho+sr(:,:,L);
end