§ 17.5. Математическое моделирование рпу методом комплексной огибающей

Математическая модель смеси на входе РПУ. В рамках метода комплексной огибающей вместо соотношений (17.22), (17.23) используют запись и₂ = Re (Ёе'^') = = Re(£_t.e^/0V) + Re(£„e^/u)n') 4-; ^_е (/?e/^w«'). Отсюда следуют алгоритмы формирования комплексной огибающей смеси:

Математическая модель комплексной огибающей смеси, построенная по алгоритмам (17.33) — (17.35), приведена на рис. 17.10. Здесь, так же как и в схеме рис. 17.5, главным является моделирование низкочастотных процессов А_с, Л_п, А, В.

Математические модели функциональных звеньев РПУ (библиотека моделей). Суть метода комплексной огибающей (51 сводится к замене реальных радиозвеньев с узкополосными колебаниями на входе эквивалентными звеньями (ЭЗ) низкочастотного типа, связывающими комплексные огибающие на входе и выходе. Нахождение математического описания ЭЗ составляет основную задачу метода комплексной огибающей. В [5] приведены правила перехода от реальных звеньев к эквивалентным и составлена библиотека моделей. Напом-

ним основные положения метода комплексной огибающей.

Математические модели линейных радиозвеньев. Для стационарных линейных радиозвеньев (ЛРЗ) (рис. 17.11, а) с заданной импульсной характеристикой h (t) = Re {Н (t) X _xe/[<a_t<-«И*)]} эквивалентное звено (математическая модель) имеет импульсную характеристику (рис. 17.11. б) вида

(17^

где со_с — частота входного сигнала, а Дсо_с == со,. — <о₀ — частотная расстройка.

Стационарные ЛРЗ можно моделировать спектральным методом. Для этого ЛРЗ с заданной КЧХ К (/со) — = у (co)/jc (со) (рис. 17.12, а) должно иметь в качестве модели КЧХ эквивалентного звена:

(17.37)

где ПРФ — индекс преобразования Фурье.

Комплексные частотные характе -ристики вида (17.37) называются укороченными и для их нахождения разработан ряд методов [51.

Так, для одиночного резонансного контура (усилителя) известна нормированная КЧХ [11]: K_N(M = КШ/К+ Q

где а — 2 — (со — со₀) — обобщенная рас-стройка; Q — io,Lr — добротность; со₀ = = \/~[/LC — резонансная частота. Тогда по формуле (17.37) находим

*эк* </") >«*»Ц «»о + Асос + 0)1 щ

= Во 0'Aco_R) [ 1 + A_t (/Дсо_с) /О]-', (17.38)

где В„ ОДсо_с) = ((o₀/2Q); А_х (/Дш_с) = [1 -|-+ 2 (Q/ci>„) /Дш_с] — комплексные коэффициенты.

Для полосового резонансного усилителя с одинаково настроенными контурами КЧХ

* /■ , — У(М'-Р*>

имеет вид K_N (yw) = _{(] + /Чх)а} _<_ р₂ • где

Р — фактор связи. Тогда формула (17.37) дает следующую укороченную КЧХ:

K_3KN (jQ) = В₀ О'Дсос) [ 1 + (/Дсо_с) /О +

+ Л₂(/Дш_с) (/О)»]-». (17.39)

Здесь В₀ (/Дсо_с), Л;(/Дсо_с), i'=--l, 2—комплексные коэффициенты более сложного вида [5].

Для получения цифровых моделей нереализуемых аналоговых фильтров с КЧХ (17.38), (17.39) тем не менее можно использовать методы оцифровывания, рассмотренные в § 17.3, в частности, алгоритмы (17.8) (17.12), (17.19), (17.20), но с комплексными коэффициентами.

Нестационарные ЛРЗ, описываемые интегралом Дюамеля (17.4), после перехода к комплексным огибающим

(17.40)

приводят к математической модели вида рис. 17.11, б с алгоритмом

оо

Ё_и (г) = J t_x (t — т) h_3K (t, т) dT -о

t

= J Ё_х (т) /i_a„ (г, t — т) dT, где введе-

— во

на комплексная импульсная характеристика ЭЗ Л_зк (/, т) h(t. т) X

X exp (—/СО_ст) Л_эке С- ^т) —

— ih_aK8(t,r), также приводящая к цифровой модели комплексного фильтра рис. 17.3.

Математические модели нелинейных радиозвеньев. В РПУ приходится моделировать безынерционные нелинейные радиозвенья (БНРЗ), состоящие из последовательно соединенного БНЗ у — G (х) и так называемого формирующего звена (ФЗ). пропускающего без искажений (с единичным усилением) спектр сигнала биений на k-u гармонике (£со_с) входного воздействия. Случай k —О соответствует БНРЗ типа детектора, k — 1 — БНРЗ типа нелинейных радиоусилителей или полосовых ограничителей, k > 2 — БНРЗ типа умножи-жителей частоты. В результате можно записать выходное напряжение БНРЗ в виде

С целью нахождения явных решений для у₀ (/), E_yh (t) при входном узкополосном колебании используют методы статистической радиотехники 151.

При подключении к схеме БНРЗ спереди и сзади резонансных усилителей получают структурные схемы типовых звеньев РПУ. Для их моделирования методом комплексной огибающей, очевидно, надо использовать поблочное моделирование на основе приведенного выше математического описания.

При моделировании нефункциональных инерционных нелинейных радиозвеньев пригоден лишь «метод укорочения» исходных дифференциальных уравнений, разработанный С. И. Евтяновым и Р. Л. Стратоновичем '151. Суть метода сводится к следующему.

Пусть нелинейное радиозвено описывается дифференциальным уравнением

где х = х (/, к) = Е (/, X) cos [ш₀г — — гр (t, X)] — входной высокочастотный сигнал, а у (t) —выходной низкочастотный сигнал.

Тогда для получения математической модели проводим решение в несколько этапов.

Этап 1. Представляем входное воздействие х и его производные через комплексные огибающие:

(17.42)

Этап 2. Подставляя (17.42) в (17.41), получаем нелинейное дифференциальное уравнение

Этап 3. Пользуясь «медленностью» изменений комплексной огибающей Ё = Ё (t, X) и ее производных по сравнению с высокочастотным множителем exp /cor = ехр /Ф, проводим усреднение правой части полученного уравнения на интервале Ф — О — 2 я:

Дифференциальное уравнение (17.43) связывает в неявном виде все низкочастотные составляющие: входную комплексную огибающую Ё = = Ее.-№ и выходной эффект у.

Этап 4. В дальнейшем, если это необходимо, проводят укорочение дифференциального уравнения (17.43), полагая

где б — величина высокого порядка малости.

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

Пример 17.3. Модель амплитудного диодного детектора (метод комплексной огибающей). Вернемся вновь к схеме рис. 17.7, а и перепишем дифференциальное уравнение (17.32):

(17.45)

Воспользуемся для решения методом Ван-дер-Поля.

Шаг I. Используя (17.43), усредним выражение (17.45) по фазе Ф:

T,d<yЩШ-гР(£, у) = 0, (17.46)

где

Шаг 2. Будем считать детектор линейным, введем коэффициент передачи по огибающей АГд = у!Е и крутизну вольт-амперной характеристики диода S. Тогда можно записать

(17.47)

где Q (К_я) - V\ - К\ - К_я arccos К_п -нелинейная функция.

Шаг 3. Подставив (17.47) в (17.46) и введя коэффициент нагрузки /С_н ⁼ Srln, получим окончательное решение дифференциального уравнения, связывающего низкочастотные процессы у (/) и Е (I):

t_euy(t)/dt + y(t)-z(y, Е, /С„) = 0. (17.48) Здесь

z(y, Е, Кн) = K_WEQ (<//£)•

Полученное уравнение соответствует модели рис. 17.13, о, где К (р) — передаточная функция rC-цепи нагрузки детектора.

Шаг 4. Проведем укорочение уравнения (17.48), положив т₀(/ = 0. Тогда по-

лучим алгебраическое нелинейное уравнение

Кд-/Сн<Э(Лд) = 0.

(17.49)

Шаг 5. Вычисляя* из уравнения (17.49) коэффициенты Кд (Кн), приходим к простейшей разомкнутой модели (рис. 17.13, б):

У=К_я(К_а)Е.

(17.50)

Примеры моделирования на основе структурных схем. В [5] приведена достаточно полная библиотека моделей различных функциональных звеньев РПУ в рамках метода комплексной огибающей. Особенно просто моделируются схемы демодуляторов (детекторов), которые сводятся к простым математическим операциям над комплексной огибающей Ё --■= а — —]Ь входной смеси:

где. а — a_x(t); Ь ш Ь_х (t); Em Е_х; Ч^{1 388} 4\г-

Важное место занимает математическое моделирование согласованных линейных фильтров (СЛФ) с импульсной характеристикой

согласованной с входным сигналом

"с* (*) = Re {Ё_сх (t) ехр(/со_е/)}.

Тогда приходим к модели рис. 17.11, б, где следует положить

(17.52)

Пример 17.4. Математическая модель супергетеродинного РПУ с произвольной модуляцией. С помощью библиотеки моделей [5] можно «набрать» математическую модель для структурной схемы РПУ вида рис. 17.Н. Полная схема такой модели показана на рис. 17.14.

При реальном моделировании все комплексные ЛФ на рис. 17.14 можно заменять схемами вида рис. 17.3 с расчетом ортогональных проекций а, Ь векторов смесей. Тогда переход к огибающей, фазе, мгновенной частоте должен осуществляться по формулам

Оцифровывание математической дели рис. 17.14 позволяет получить до точную эффективную цифровую мо; универсального РПУ для исследования дач помехозащищенности и помехоуе чивости.

Модель подобного вида требует бс шого разнообразия исходных данных по параметрам смеси, так и по парамет блоков РПУ.