testlog.m

雷达信号杂波仿真中,对数正态杂波的CFAR分析程序

function z=testlog(miu,seta,N)
% clc;
% clear;
% close all;
% seta=0.2;         %方差
% miu=0.5;          %均值
% N=1024;           %采样点数

v=randn(1,N);       %输入的序列，v～N（0,1）,对应功率谱密度为1
Sv=fft(v,N);        %v所对应的频谱 

aifa=1.665;         %定义的常数
w3db=40;            %3db带宽  
w=-40:80/N:40-80/N;
w=w(1:N);
Sij=exp(-((aifa*w/w3db).^2));   %高斯谱
sij=ifft(Sij,N);                %期望得到的相关系数序列

pij=log(1+(exp(seta)-1).*sij)./seta;        %非线性变换前的相关系数序列
Su=fft(pij,N);                              %pij所对应的频谱
Hw=abs(sqrt(Su./pij(1)));                   %归一化的传递函数
u=ifft(Sv.*Hw,N);                           %通过传递函数H(w)后的序列，u～N（0,1）
y=sqrt(seta).*u+log(miu);                   %非线性变换前的序列，y～N（ln(miu),seta^2）
%z=real(exp(y));                             %得到所需序列z
z=exp(y);
z=z/max(abs(z));


%plot(real(z))
% figure(3)
% ee=0.01:0.01:4;
% [gs,zjz]=hist(z,ee);            %z序列统计直方图
% xl=gs./N*90;
% bar(ee,xl);
% hold on
% t=0;
% for x=0.01:0.01:4
%   t=t+1;
% f1=(1/(sqrt(2*pi)*sqrt(seta)*x))*exp(-(log(x/miu))^2/(2*seta));     %对数正态分布概率密度函数表达式
% flog(t)=f1;
% end
% plot(ee,flog,'r')   %对数正态分布pdf理论值
% title('对数正态杂波pdf')
% xlabel('z')
% ylabel('pdf')