NTtheta1111.zip_matlab_

f0=1/L0; %大气外部尺度的长度L0分之一 (m) dfx=1/L; %相屏大小L分之一 dfy=1/L; Gx=L;%相屏大小L Gy=L; h=zeros(Nx,Ny); %总相位屏幕尺寸像素 Nx，Ny h1=zeros(Nx,Ny); hlf=zeros(6,6); f=zeros(Nx,Ny); p=zeros(Nx,Ny); phf=zeros(Nx,Ny); plf=zeros(Nx,Ny); % %零均值，单位方差的高斯随机数矩阵 % h1=(randn(Nx/2,Ny)+sqrt(-1)*randn(Nx/2,Ny))/sqrt(2); % h(1:Nx/2,1:Ny)=h1; % h1(:,1:Ny/2)=rot90(h1(:,1:Ny/2)); % h1(:,1:Ny/2)=rot90(h1(:,1:Ny/2)); % h1(:,Ny/2+1:Ny)=rot90(h1(:,Ny/2+1:Ny)); % h1(:,Ny/2+1:Ny)=rot90(h1(:,Ny/2+1:Ny)); % % h(Nx/2+1:Nx,Ny/2+1:Ny)=conj(h1(:,1:Ny/2)); % h(Nx/2+1:Nx,1:Ny/2)=conj(h1(:,Ny/2+1:Ny)); Gau=(randn(N,N)+sqrt(-1)*randn(N,N))/sqrt(2);%零均值，单位方差的高斯随机数矩阵 for m1=-Nx/2:Nx/2-1 for n1=-Ny/2:Ny/2-1 f(m1+Nx/2+1,n1+Ny/2+1)=2*pi*sqrt(0.00058*r0^(-5/3)*((m1*dfx)^2+(n1*dfy)^2+f0^2)^(-11/6))/sqrt(Gx*Gy); %if mod(m1-8,16)==0 && mod(n1-8,16)==0 plf(m1+Nx/2+1,n1+Ny/2+1)=PLF(m1+Nx/2+1,n1+Ny/2+1,4,r0,Gx,Gy,Nx,Ny,f0); %end end end % % plf=medfilt2(real(plf))+medfilt2(imag(plf)); phf=ifft2(Gau.*f)*Nx*Ny; p=abs(phf+plf); % xlswrite('G:\研究所任务\毕设结题\自准直仪编程\相位屏模拟结果\surf.xls',surf); a=abs(phf)-mean(abs(phf(:))); ScreamDATA{1,time}=a; figure(1); mesh(abs(phf)-mean(abs(phf(:))));axis tight; colorbar;title('phf'); figure(2); imagesc(abs(phf)-mean(abs(phf(:))));axis tight; colorbar;title('phf'); figure(3); imagesc(abs(plf)-mean(abs(plf(:))));axis tight; colorbar;title('plf'); figure(4); mesh(abs(plf)-mean(abs(plf(:))));axis tight; colorbar;title('plf'); figure(5); imagesc(p-mean(p(:)));axis tight; colorbar;title('p'); figure(6); mesh(p-mean(p(:)));axis tight; colorbar;title('p'); figure(7); mesh(real(phf));axis tight; colorbar;title('phf'); function NT = NTtheta(W,w0) lambda=0.8*10^(-6);%w0=0.01; k=2*pi/lambda; N1=792;N2=792; %W=0.1; r0=w0/W; z=2;Cn2=(r0/0.185)^(5/3)*z/lambda^2; pixel=12.5*10^(-6); for kx=1:N1 for ky=1:N2 Phi_n(kx,ky)=0.033*Cn2*((kx/(pixel*N1))^2+(ky/(pixel*N2))^2)^(-11/6); end end Chi=(randn(N1,N2)+sqrt(-1)*randn(N1,N2))/sqrt(2); Phi=2*pi*k^2*z*Phi_n; MPhi=fftshift(Phi); theta=ifft2(Chi.*MPhi); Theta=abs(theta); Ttheta= Theta/max(max(Theta)); NT=Ttheta(1:600,1:792); %subplot(211),imshow(Ttheta) %subplot(212),imshow(NTtheta)