Using obvious correlations
 |
(14) |
let us present expression (13) as
 |
(15) |
As a result equation (9) starts looking like this
 |
(16) |
According to (16), the currents and voltages of the radiators are excited by EMF 
and MMF 
generators, as well as by the waves of transmission lines 
reflected from the loads of the partial radiators. Let us present the amplitude of the waves reflected from the loads using the current and voltage of falling waves
 |
(17) |
where the coefficients of diagonal matrix 
are defined by expressions (7), (10), 
As a result equation (16) starts looking like this
 |
(18) |
where
 |
(19) |
Equation (18) is solved with the help of the iteration method. At the initial stage of the iteration procedure, the solution we are looking for is represented as a sum of two summands 
that are the solutions of the following equations
(20а) 
(20б)
and equation (20а) is solved. At the K-th stage, the correction for the approximate solution obtained at the previous stage is presented as sum 
, the summands of which are the solution of the following equations
 |
(21а) |
 |
(21б) |
and equation (21а) is solved. The solution of the initial equation is a sum of the following series
 |
(22) |
The first part of equation (21а) is the vector of doubled amplitudes of waves that reflected from loads of partial radiators at the previous stage of the interaction procedure. Thus, algorithm (20)...(22) actually describes the process of successive reflections of waves between radiator apertures and loads.
Iteration procedure (20)...(22) converges because the modules of the reflection coefficients at the apertures of the partial radiators are less than one due to the radiation and losses in constructive elements.
If we implement the iteration procedure numerically, equations (20а), (21а) are solved with the help of the Fourier transform method. After we apply discrete Fourier transform (DFT) to (20а), (21а), we get
 |
(23) |
where the tilde indicates the discrete Fourier transforms (DFTs) of the corresponding values. Switching to solving equations (20а), (21а) is made by using inverse discrete Fourier transform
 |
(24) |
where 
is the operator of the inverse discrete Fourier transform.
