%%%海雜波仿真
clc;
clear;
close?all;
%%%log-normal仿真
%%功率譜函數(shù)
f_3dB=20;
apha=2*sqrt（log（2））;
BW=5*f_3dB;
N=1024;???????%%頻域采樣點(diǎn)數(shù)
f=linspace（-BW/2BW/2-BW/NN）;

S=exp（-（apha*f/f_3dB）.^2）;???%功率譜函數(shù)
s=fftshift（ifft（S））;
s=s./max（s）;?????????????????%自相關(guān)函數(shù)
%%求濾波器
uc=1;
fai=0.6;

Rou=log（1+s*（exp（fai^2）-1））/fai^2;????%非線性變換前的相關(guān)系數(shù)
G=fftshift（fft（Rou））;?????????????????%非線性變換前的功率譜函數(shù)

H_abs=sqrt（G/Rou（1））;
H_ph=exp（1i*2*pi*f）;
H=H_abs.*H_ph;???????????????????????%構(gòu)造濾波器
%%%%%%%%%%%%%
randn（‘state‘sum（100*clock））;???%產(chǎn)生服從U（0，1）的隨機(jī)序列?
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）;???????????????????????%產(chǎn)生對(duì)數(shù)正態(tài)分布隨機(jī)序列
%%%%%%%%%%%%%%
figu