日志

highspeedlogic仿真-种大气湍流相位屏仿真函数，可准确模拟大气对激光 ...

热度 1已有 1147 次阅读2019-7-13 03:51 |系统分类:其他

%%phase screen 128*128 function alpha=phase_screen(Cn_2,Z,lamda); %%Parameter--phase screen x=linspace(-0.45,0.45,128); nx=length(x); y=x'; ny=length(y); disp('nx');disp(nx); disp('ny');disp(ny); kx=2*pi./x; ky=2*pi./y; K_L=2*pi/1.28; K_O=2*pi/20; S=1.28; V=100; %%phase_screen t=normrnd(0,1,nx,ny); log_t=log(t); t0=normrnd(0,1,1,1); log_t0=log(t0); theta=2*pi*t0; beta=-2*pi*(t0-0.5); %alpha_IF alpha_IF_Re_=((ones(nx,1)*kx).^2+(ky*ones(1,ny)).^2).^(-11/12); alpha_IF_Re_c=2*pi*sqrt(V)/lamda*sqrt(-(2*pi)^3*Cn_2*log_t/30); alpha_IF_Re=alpha_IF_Re_.*alpha_IF_Re_c; alpha_IF_Im=exp(2*pi*j*(t-0.5)); alpha_IF=alpha_IF_Re.*alpha_IF_Im; alpha_IF_conj=conj(alpha_IF); alpha_IF=ifft(alpha_IF_conj); %csvwrite('jp\alpha_IF.dat',alpha_IF); %disp(alpha_IF_f); alpha_IF_abs=abs(alpha_IF); %csvwrite('jp\alpha_IF_abs.dat',alpha_IF_abs); %alpha_tilt a_R=sqrt(-(2*pi)^4*Cn_2*Z*log_t*(K_L^(1/3)-K_O^(1/3))/(10*lamda*lamda)); alpha_tilt_c=ones(nx,1)*x*cos(theta)+y*ones(1,ny)*sin(theta); alpha_tilt=a_R.*alpha_tilt_c; %csvwrite('jp\alpha_tilt.dat',alpha_tilt); %alpha_focus b_R=sqrt(2*(2*pi)^4*Cn_2*Z*(K_L^(7/3)-K_O^(7/3))/(187*lamda*lamda))*sqrt(-log_t); b_R_x=cos(beta)*b_R; b_R_y=sin(beta)*b_R; alpha_focus_1=(cos(theta)*ones(nx,1)*x+sin(theta)*y*ones(1,ny)).^2; alpha_focus_p1=b_R_x.*alpha_focus_1; alpha_focus_2=(-sin(theta)*ones(nx,1)*x+cos(theta)*y*ones(1,ny)).^2; alpha_focus_p2=b_R_y.*alpha_focus_2; alpha_focus=0.5*(alpha_focus_p1+alpha_focus_p2)-S*S*(b_R_x+b_R_y)/24; %csvwrite('jp\alpha_focus.dat',alpha_focus); %alpha alpha=alpha_IF_abs/3+alpha_tilt/3+alpha_focus/3; %csvwrite('jp\alpha.dat',alpha);

phi=exp(j*alpha);

收藏分享邀请举报