Алгоритм

Идея последовательного приближения

Идея заключается в последовательном удвоении числа элементарных модельных экспериментов. Это продолжиться до тех пор, пока результат текущего моделирования не приблизится к результату моделирования на предыдущей итерации:

$\left | F_i-F_{i-1} \right |<\Delta$ , где

\Delta

- параметр моделирования, задаваемый пользователем.

Однако сравнение соседних итераций не дает достаточного условия на достижение заданной точности (основная причина этого - излучение по направлениям). Поэтому правильнее будет сравнивать итерации через ону, две и т.д.:

$\left | F_i-F_{i-k} \right |<\Delta$ , где

k

также будет задаваться пользователем.

i

здесь - это параметр цикла, стоящего над циклами основной программы,

i=\overline{0:N}

.

При увеличении числа элементарных экспериментов в два раза, только каждый второй будет уникальным. Другая половина будет повторять эксперименты, выполненные на предыдущих итерациях. Поэтому в цикл основной программы введено дополнительное условие для учета этих повторений.

Инициализация геометрической модели

Перерасчет высот с учетом кривизны земли и рефракции радиоволн в тропосфере.

Входной параметр $h$ пересчитывается в соответствии с формулой:

$h'(x,y) = h(x,y)-\frac{-2 R_{eq}+\sqrt{(2 R_{eq})^{2}+4 r^{2}}}{2} \approx h(x,y)-\frac{r^{2}}{2 R_{eq}}$ , где

R_{eq}=\frac{R_0}{1+R_0 \frac{dn}{dh}}

- эквивалентный радиус Земли, где

R_0=6371

- радиус Земли (км),

\frac{dn}{dh}=grad~n

- изменение коэффициента преломления с высотой.

r=\sqrt{x^{2}+y^{2}}

- расстояние до точки с высотой

h(x,y)

.

Инициализация источников

При переотражениях лучевая модель предполагает один входящий луч и два выходящих - отраженный и преломленный. Инициализируем множество источников $SC$ .