function?phz=vkolmg（wvlDdzNCNL0l0）
%?D?相位屏尺寸
%N采樣點數(shù)目
%L0?外尺度?l0?內(nèi)尺度?CN常數(shù)?dz傳輸距離
delta=D/N;
x=（-N/2:N/2-1）*delta;
y=x;
[X?Y]=meshgrid（xy）;
del_f=1/（N*delta）;
fx=（-N/2:N/2-1）*del_f;
[kx?ky]=meshgrid（2*pi*fx）;
k=2*pi/wvl;
[th?ka]=cart2pol（kxky）;
km=5.92/l0;
k0=2*pi/L0;
%?r0=0.185*（wvl^2/（dz*CN））^（3/5）;
PSD_phi=0.033*CN*exp（-（ka/km）.^2）./（ka.^2+k0^2）.^（11/6）;
PSD_phi（N/2+1N/2+1）=0;
cn=2*pi*k.^2*dz.*PSD_phi*（2*pi*del_f）.^2;
phz_hi=ift2（（randn（N）+1i*randn（N））.*sqrt（cn）1）;%由于matlab本身的FFT算法的緣故，又cn為傅里葉變換展開的系數(shù)，此時del_f=1;
phz_hi=real（phz_hi）;
%?figure;imagesc（phz_hi）;colorbar;
%%?低頻補償
phz_lo=zeros（size（phz_hi））;
for?p=1:

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1354??2012-05-26?23:35??vkolmg.m

-----------?---------??----------?-----??----

?????????????????1354????????????????????1