%%%海雜波仿真
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