swerlingiv.m

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

function [x]=swerlingIV(delta,n)
%产生swerlingIV型的雷达回波信号
%delta：实测雷达目标RCS的均值
%
r=ones(1,n);
u=junyun(n);
v=junyun(n); 
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);
y=zeros(1,n);
for k=1:n
    y(k+1)=m(k)+log(1+y(k));
end
x=(delta/2).*y;