Формирование пространственного спектра (диаграммы направленности) в зоне френеля объектов с помощью линзовых и зеркальных систем. Часть 3. Струков И.Ф. - 27 стр.

                                                            27


    П рим ер 2. Прог ра мма для ра счета одномерног    о ра спределени я и нтенси в ности
спек тра льной плотности поля пря моу г     ольног  о объек та , сформи ров а нной в
фок а льной и ли эк в и в а лентной фок а льной плоск ости ли нзы, с у четом
на пра в ленных св ойств и мпу льсной ха ра к тери сти к и среды cosθ = f                 .
                                                                               x2 + f 2
Ра счет в едется по форму ле (4.30)
                                                                                                               2
                .
                                                                                 ka  x         
               U& ( x, f )                                                2 sin
                                                                                       − sin α 
                                        G& ( ω1 − ω0 )
                                                       2
                                                                                 2f
                                                                  f                             
                                   =                        =           ⋅
           U& ( x, f )                 G (ω1 − ω0 )
                                                      2       x2 + f 2          a x         
                             max                      max                    k  − sin α 
                                                                                  2 f         
                                                                                                               (4.30)

                        −6                              π
     TOL := 10                 λ := 0.4      k := 2 ⋅         f := 100 ⋅λ          imax := 50
                                                        λ       j := 1 .. 4        i := 0 .. 2 ⋅ imax
                                                                                             ( i − imax)
                                                                                   xi :=                 ⋅10
                                                                                                 imax


                                                                     2
                      sink ⋅ x − sin( α )  ⋅ a              
                                                                                   2
       (         )         f                2                             f   
      U x , a , α :=                                                   ⋅        
                      k ⋅ x − sin( α)  ⋅ a                             2   2
                                                                           x +f 
                          f              2                    

      α1 := 0                                                 α2 := 0.17
      a1 := 10 ⋅ λ           a2 := 20 ⋅ λ

                    (
      u1i := U xi , a1 , α1        )                                     (
                                                            u3i := U xi , a1 , α2    )
      u1max := max( u1)                                     u3max := max( u3)
                    u1i                                                   u3i
      u1i :=                                                u3i :=
               u1max                                                  u3max

                    (
      u2i := U xi , a2 , α1        )                                         (
                                                             u4i := U xi , a2 , α2       )
      u2max := max( u2)                                      u4max := max( u4)
                    u2i                                                      u4i
      u2i :=                                                 u4i :=
               u2max                                                     u4max