对数正态分布杂波仿真

%%%海杂波仿真
clc;
clear;
close all;
%%%log-normal仿真
%%功率谱函数
f_3dB=20;
apha=2*sqrt（log（2））;
BW=5*f_3dB;
N=1024;       %%频域采样点数
f=linspace（-BW/2BW/2-BW/NN）;

S=exp（-（apha*f/f_3dB）.^2）;   %功率谱函数
s=fftshift（ifft（S））;
s=s./max（s）;                 %自相关函数
%%求滤波器
uc=1;
fai=0.6;

Rou=log（1+s*（exp（fai^2）-1））/fai^2;    %非线性变换前的相关系数
G=fftshift（fft（Rou））;                 %非线性变换前的功率谱函数

H_abs=sqrt（G/Rou（1））;
H_ph=exp（1i*2*pi*f）;
H=H_abs.*H_ph;                       %构造滤波器
%%%%%%%%%%%%%
randn（‘state‘sum（100*clock））;   %产生服从U（0，1）的随机序列 
x=randn（1N）;
X=fftshift（fft（x））;

Y1=X.*H;          %滤波
y1=real（fftshift（ifft（Y1）））;
sigmac=std（y1）; 
muc=mean（y1）; 
y1=（y1-muc）/sigmac; 

y2=y1*fai+log（uc）;
y=exp（y2）;                       %产生对数正态分布随机序列
%%%%%%%%%%%%%%
figu