swerlingii.m

来自「雷达swerling 目标的仿真代码」· M 代码 · 共 30 行

function [x]=swerlingII(delta,n) 
%产生目标特性是swelingII分布的雷达回波信号
%调用方式：[x]=swerlingII(n)
%n：序列的长度, delta:实测目标均值
%
r=ones(1,n); 
u=junyun(n); %混合同余法得到随机数
v=junyun(n); %取两个随机变量满足均匀分布，且独立性好
%以下构造u和v的独立性
for i=1:n 
if (u(i)==0) 
u(i)=0.0001; 
else 
continue 
end 
end 
for i=1:n 
if (u(i)==v(i)) 
    u(i)=u(i)+0.0001 
else
    continue
end 
end 

t=-2*log(u); 
h=2*pi*v; 
w=sqrt(t).*cos(h); %w和p满足标准正态分布
p=sqrt(t).*sin(h); 
m=(r/2).*(w.^2+p.^2);%满足自由度为2的卡方分布
x=delta.*m