Функция для расчёта ДН одиночного вибратора (1.44)
, {2.14}
а полная функция для расчёта ДН ЛПВА (1.43)
.{2.15}
В функциях для расчёта ДН аргумент I – вектор-столбец токов в вибраторах [I], остальные аргументы понятны по наименованию.
Функцию для входного сопротивление ЛПВА вычисляем согласно (1.46)
, {2.16}
где Z – матрица ; I –
вектор-столбец токов в вибраторах; N – количество вибраторов.
Расчёт
КСВ производится обычным образом – по заданному входному сопротивлению и реальному входному сопротивлению
.
, {2.17}
где rwx – заданное входное сопротивление; zin – реальное входное сопротивление. В первом операторе функции {2.17} вычисляется комплексный коэффициент отражения, во втором – КСВ.
Для
расчёта КНД, согласно (1.42), необходимо определить модуль максимума ДН и
действительную часть входного сопротивления. Модуль максимума ДН необходим и
для нормирования при выводе графиков. При выбранной нами системе координат
(рис. 1.1) максимум ДН приходится на углы ,
. В результате
. {2.18}
Функцию для расчёта КНД получаем из (1.42)
, {2.19}
где md – модуль максимума ДН; zin – входное сопротивление ЛПВА.
Для трёхчленного закона распределения тока элементы матрицы системы линейных алгебраических уравнений и правой части описываются выражением (1.37б). Составлять программы начнём с подынтегральных выражений. Этих выражений согласно (1.41), как и интегралов, три, поскольку ток описывается тремя гармониками. Выражения состоят из произведения двух функций – гармоник распределения тока (1.41) и ядра (1.31). Функция для расчёта ядра {2.11} уже была рассмотрена ранее. Для гармоник распределения тока получаем
, {2.20}
, {2.21}
. {2.22}
Аргументы этих функций, за исключением
массива , описывались. Массив чисел
следует вычислить заранее, чтобы во
время интегрирования не загружать процессор компьютера бессмысленной работой.
Также нужен и массив чисел
, который будет
необходим при расчёте ДН
, {2.23}
. {2.24}
В аргументах выражений (1.37а) и
(1.37б) имеются числа . Введём функцию
для расчёта этих чисел
. {2.25}
Во втором члене ядра (1.31) в
соответствии с (1.37б) имеется множитель ,
а в (1.37а) и (1.37б) – множитель
. Они
рассчитываются функциями {2.8} и {2.9}, соответственно, с заменой аргумента kh на kh3, а аргумента N – на
3N.
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.