clear
close all
T=1/4;
n=2^10;cn2=1.5e-13;
drtz=15;%%%相位屏的间距
dxy=1/n;
x=n/2;y=n/2;
% nx1=linspace(-n/2-1e-5,n/2-1e-5,n)*dxy;nx1=nx1(1:n);
% ny1=linspace(-n/2-1e-5,n/2-1e-5,n)*dxy;ny1=ny1(1:n);
% [X1,Y1]=meshgrid(nx1,ny1);
nx=linspace(-n/2-1e-5,n/2-1e-5,n)*dxy;nx=nx(1:n);
ny=linspace(-n/2-1e-5,n/2-1e-5,n)*dxy;ny=ny(1:n);
[X,Y]=meshgrid(nx,ny);
[TH,R]=cart2pol(X,-Y);
df=2*pi/dxy/n;
ffx=linspace(-n/2-1e-5,n/2-1e-5,n+1)*df;ffx1=ffx(1:n);
ffy=linspace(-n/2-1e-5,n/2-1e-5,n+1)*df;ffy1=ffy(1:n);
[fx,fy]=meshgrid(ffx1,ffy1);
w=0.01;
A1=2;
B=1;
times=1;
Times=1;%循环次数
Intensity(1:times,1:Times)=0;sci(1:times)=0;
mx=1;
my=7;
fr=(fx.^2+fy.^2).^(0.5);
df1=2*pi/1/n;
ffx2=linspace(-n/2-1e-5,n/2-1-1e-5,n+1)*df1;ffx3=ffx2(1:n);
ffy2=linspace(-n/2-1e-5,n/2-1-1e-5,n+1)*df1;ffy3=ffy2(1:n);
[fx1,fy1]=meshgrid(ffx3,ffy3);
fr1=(fx1.^2+fy1.^2).^(0.5);
s=1e-9;%流体每单位质量湍流动能的耗散率
yita=1e-3;xt=1e-8;%均方温度耗散率
ww=(-2);drtr=8.284*(fr1*yita).^(4/3)+12.978*(fr1*yita).^2;%ww温度和盐度波动的相对强度
at=1.863e-2;as=1.9e-4;ats=9.41e-3;
feye=0.388*10^(-8)*s^(-1/3).*fr1.^(-11/3).*(1+2.35.*(fr1.*yita).^(2/3)).*xt.*(ww^2.*exp(-at.*drtr)+exp(as.*drtr)-2*ww.*exp(-ats.*drtr))./(ww^2);
X1=X(n/2,:);
lamd=0.532e-6;k=2*pi/lamd;
u1=Ai((A1-R)./w).*exp(B(A1-R)./w);
% figure(11)%%%%%彩色光强图
% axis square
% colormap('hot')
% colorbar
% surf(X,Y,Ip1);
% shading interp;
% xlabel('m');
% ylabel('m');
% axis([-0.15 0.15 -0.15 0.15]);
% view(2);
% shading flat;
