produce_mes.m

来自「雷达目标跟踪matlab源程序 希望对大家能够有所帮助」· M 代码 · 共 40 行

%用于航迹起始测试程序(模拟产生用于航迹起始的多目标量测信息) 

jpda_parameter;

target_total=2;%目标个数
x1=[10 3 0 50 -2 0]';%格式[x vx ax y vy ay]'
x2=[30 -4 0 20 3 0]';
x=[x1 x2];

all_mes=[];
all_mes_total=[];
all_x=[];
for i=1:100
    noise_process=[];
    noise_mea=[];
    for s=1:target_total
       temp=normrnd([0 0]',diag(q));
       temp1=normrnd([0 0]',diag(r));
       noise_process=[noise_process temp];
       noise_mea=[noise_mea temp1];
    end
    x=fai*x+g*noise_process;
    all_x=[all_x x];
    mes=h*x+noise_mea;
    all_mes=[all_mes mes];
    
    mf=poissrnd(arg);%
    if mf>20
       mf=20;%限制杂波数不超过,just for simulation
    end
    all_mes_total(i)=mf+target_total;
    for j=1:mf
       false_alarm=unifrnd(u1,u2); %在监视空域内产生服从均匀分布的虚警
       all_mes=[all_mes false_alarm];%产生虚警   错误量测在跟踪门内服从均匀分布
    end
end
save all_mes all_mes;
save all_mes_total all_mes_total;
save all_x all_x;