Замена непрерывной передаточной функции дискретной передаточной функцией

Второй путь – замена непрерывного интегрирования дискретным. Для этого в передаточной функции
перейдем от оператора дифференцирования p к оператору интегрирования 1/p, поделив числитель и знаменатель на pn.
Значения интеграла
отстоящие друг от друга на интервал
времени ∆t, связаны соотношением
С учетом обозначений, принятых на рисунке

Простейшие методы численного интегрирования основаны на замене криволинейной трапеции прямоугольником (ABFE – прямоугольник вперед, AGDE – прямоугольник назад) или трапецией ABDE. Получатся следующие приближенные формулы для вычисления интеграла. Для сокращения записи принято y(i∆t) = yi.
Метод прямоугольников вперед yn = yn-1 + SABFE = yn-1 +∆t xn-1.
Метод прямоугольников назад yn = yn-1 + SAGDE = yn-1 +∆t xn.
Метод трапеций
yn = yn-1 + SABDE = yn-1 +∆t (xn+ xn-1)/2
Возьмем Z-преобразование Y(z) = z-1Y(z) + ∆tz-1X(z).
Тогда
Оператор 1/p в непрерывной передаточной функции можно заменить любой из найденных дискретных передаточных функций

Сравним эти методы на примере переходной характеристики интегрирующей цепи.
Как видим, дискретная модель непрерывной системы может быть получена различными путями. Целесообразно выбрать самый точный.
Разностное уравнение интегрирующей цепи. α = ∆t/T, Т – постоянная времени

Как видим, дискретная модель непрерывной системы может быть получена различными путями. Целесообразно выбрать самый точный. Практика моделирования показывает, что наименьшую ошибку обеспечивают методы численного интегрирования дифференциальных уравнений.

Моделирование узкополосных РЭС
Узкополосными называются РЭС, амплитудно-частотная характеристика которых сосредоточена в области высоких частот, причем центральная частота f0 АЧХ много больше полосы пропускания 2∆f .
В узкополосных РЭС сигнал представляется гармоническим колебанием с медленно меняющимися амплитудой Um(t) и фазой φ(t):
u(t) = Um(t)sin(2πf0t +φ(t)).
В зависимости от типа переменной, для которой составляется модель, различают три типа моделей: 1) для несущей, 2) для комплексной огибающей и 3) для информационного параметра.
Моделирование по методу несущей. Переменной моделирования является мгновенное значение сигнала u(t). Математическая модель составляется по D-схеме в соответствии с функциональной или принципиальной электрической схемами РЭС

Составим математическую модель усилителя промежуточной частоты (УПЧ), охваченного системой автоматической регулировки усиления (АРУ).
В математической модели УПЧ представим в виде соединения безинерционного усилительного элемента с коэффициентом усиления K(UАРУ), зависящим от управляющего напряжения UАРУ, и линейного устройства с импульсной характеристикой gУПЧ(t) или передаточной функцией KУПЧ(p).
Детектор АРУ является выпрямителем переменного выходного напряжения с детекторной характеристикой fд(uвых). Фильтр нижних частот описывается импульсной характеристикой gФНЧ(t) или передаточной функцией KФНЧ(p). УПТ считаем широкополосным с постоянным коэффициентом передачи KУПТ.

Главный недостаток модели, составленной по методу несущей, неразумная трата машинного времени. Как уже отмечалось ранее, частота дискретизации при моделировании непрерывных систем берется в 5 – 10 раз выше граничной частоты, что означает для метода несущей fдискр = (5 – 10)f0, то есть 5 – 10 отсчетов за период высокочастотного колебания.
В то же время интересующая разработчика информация заложена не в самом гармоническом колебании на частоте f0, а в отклонении колебания от гармонического, то есть в его медленно изменяющихся амплитуде и фазе.

Модель, построенная по методу комплексной огибающей, как раз и позволяет исследовать амплитудно-фазовые изменения, не моделируя само несущее колебание. Частота дискретизации в этой модели определяется уже не несущей частотой f0, а шириной спектра сигнала, который близок к полосе пропускания АЧХ РЭС : fдискр = (5 – 10)∆f.
Комплексная амплитуда сигнала:
Нелинейные преобразования совершаются над амплитудой сигнала Um(t), а линейные над действительной Uc(t) и мнимой Us(t) составляющими сигнала.
Импульсная характеристика узкополосного УПЧ
Ее комплексная огибающая
При использовании передаточной функции K(p) для моделирования линейного звена ее необходимо заменить низкочастотным эквивалентом для комплексной огибающей.

Математическая модель УПЧ с АРУ, составленная по методу комплексной огибающей:

Таким образом, в линейной модели УПЧ производится квадратурная обработка комплексной амплитуды.
Входное напряжение детектора АРУ определяется как
.
Модель низкочастотной части: ФНЧ и УПТ не изменилась.

Моделирование по методу информационного параметра. Обычно эта модель используется для изучения систем радиоавтоматики, динамические свойства которых практически определяются низкочастотной частью схемы. Поэтому в модели узкополосные высокочастотные устройства считаются безинерционными элементами или вносящими небольшой фазовый сдвиг.
Для системы АРУ информационным параметром является амплитуда выходного сигнала УПЧ.
УПЧ как элемент системы авторегулирования вносит фазовый сдвиг в огибающую сигнала промежуточной частоты. Если АЧХ УПЧ близка к прямоугольной, ФЧХ в полосе пропускания линейна, и УПЧ заменяют звеном идеального запаздывания. В другом крайнем случае АЧХ считают близкой к АЧХ одиночного колебательного контура и УПЧ заменяют инерционным звеном.