close?all;
????????clear?all;
??????
?????????theat_i=30*pi/180;
?????????N=64;%total?number?of?sample?points%為了作fft，N取2的n次方個值
?????????Die=N^2;
?????????namda=1e-6;
?????????k=2*pi/namda;
?????????L=16*namda;%粗糙面長度
?????????deatax=L/N;
?????????sigma=1*namda;
?????????l=2*namda;
?????????aipcino0=8.8542e-12;
???????????miu0=pi*4e-7;
????????????namda=1e-6;
????????????yita0=sqrt（miu0/aipcino0）;
tic????????
????????XunHuan=20?;%%%%%%%%%%%表面?zhèn)€數(shù)
????????for?Xun=1:XunHuan
????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????????seed=1455*randn（）;
????????[z_incx_inc]=gauss1（NLsigmalseed）;
????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%??
??????Xun
??%??????figureplot（x_incz_inc）;
?????????g=L/4;
?????????w_r=（2.*（x_inc+z_inc.*tan（theat_i））.^2