lg.m

这是雷达杂波幅度符合对数正态分布的仿真程序

%---------相参相关对数－正态分布杂波-------------------%
clear all;close all;clc;
miu=0.8;
% sigma=sqrt(0.2);
p=1.6;
sigma=sqrt(2*log(p));
num=8192;
%-----由功率谱和相关系数确定相关高斯序列；
fs=1000;
f=fs/num*(-num/2:num/2-1);%对称的功率谱
sig=1;
amp=1/(sqrt(2*pi)*sig);
Pw=amp*exp(-(f).^2/(2*sig^2));
figure(1);
plot(f,10*log10(abs(Pw)/max(abs(Pw))),'b--');title('期望对数正态分布随机序列的功率谱(高斯谱)');axis([-50 50 -50 0])
hold on


sn=ifft(Pw);
pn=log(1+sn*(exp(sigma^2)-1))/sigma^2;
Xw=fft(pn,num);     %相关高斯随机序列功率谱
abs_Hw=sqrt(abs(Xw)/pn(1));
op=sqrt(-1);
w=rand(1,num);
ww=2*pi*w;
% ww=ww/max(abs(ww));
nois=exp(op*(ww));  %加入的随机相位  相干信号
Hw=abs_Hw.*nois; %成形滤波器频谱

x1=randn(1,num);
y1=randn(1,num);
x2=ifft(fft(x1).*Hw);
y2=ifft(fft(y1).*Hw);

%------求相位角
sita=angle(1./atan(x2./y2));
r=exp(y2*sigma+log(miu));
i=abs(r).*sin(sita);
q=abs(r).*cos(sita);
j=sqrt(-1);
z=i+j*q;
figure(2);
hist(abs(z),60);
% figure(3);
% hist(i,20);
% figure(4);
% hist(q,20);

pw=psd(z,num,'hamming(1024)',512);
pxx=abs(pw);
figure(1);
plot(f,10*log10(fftshift(pxx/max(pxx))));


r=0:0.1:6-0.1;
miu=1;
p=1.2;
pdf=(1./(sqrt(4*pi*log(p))*r)).*exp((-1/4*log(p))*log(r/miu).^2);
figure(3);
plot(r,pdf);