矩量法一维导体散射计算

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 ;%%%%%%%%%%%表面个数
        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