vkolmg kolmogorov功率谱反演

kolmogorov功率谱反演，次谐波补偿的方法生成大气湍流相位屏

function phz=vkolmg（wvlDdzNCNL0l0）
% D 相位屏尺寸
%N采样点数目
%L0 外尺度 l0 内尺度 CN常数 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为傅里叶变换展开的系数，此时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