Радиотехника и электроника, 2021, T. 66, № 5, стр. 483-489

Дифференциальные уравнения для триггера и мультивибратора

А. Е. Китаев *

Нижегородское научно-производственное объединение им. М.В. Фрунзе
603951 Нижний Новгород, просп. Гагарина, 174, Российская Федерация

* E-mail: kitaev_a_e@mail.ru

Поступила в редакцию 21.10.2020
После доработки 29.10.2020
Принята к публикации 06.11.2020

Полный текст (PDF)

Аннотация

Получены дифференциальные уравнения для симметричного транзисторного триггера и мультивибратора. Показано, что при выводе этих уравнений может использоваться полученное ранее автором выражение для токовой характеристики транзисторного усилителя с общим эмиттером, использующее специальную функцию Ламберта (хотя уравнения могут формулироваться и на основе характеристик другого вида). Исследован характер состояний равновесия для дифференциальных уравнений триггера. Проведено исследование “медленных” и “быстрых” движений для дифференциальных уравнений мультивибратора. Приведен пример численного решения для мультивибратора.

1. ПОСТАНОВКА ЗАДАЧИ

В работах [1, 2] автором были предложены формулы для коллекторного тока усилителя с общим эмиттером (при возможном наличии стабилизирующего резистора в эмиттерной цепи):

(1)
$\begin{gathered} {{i}_{c}} = \alpha {{D}_{{01}}}({{U}_{{{\text{in}}}}}) - {{D}_{{22}}}(\alpha {{D}_{{01}}}({{U}_{{{\text{in}}}}}) \times \\ \times \,\,(R + {{R}_{n}}) - (E - {{U}_{{{\text{in}}}}})). \\ \end{gathered} $

Функции D01 и D22 выражаются через специальную математическую функцию Ламберта и имеют следующий вид:

(2)
$\begin{gathered} {{D}_{{01}}}(x) = - {{I}_{{\text{s}}}} + \frac{\varphi }{{R + {{R}_{{ne}}}}} \times \\ \times \,\,W\left( {{{I}_{{\text{s}}}}\frac{{R + {{R}_{{ne}}}}}{\varphi }\exp \left( {\frac{{x + {{I}_{{\text{s}}}}(R + {{R}_{{ne}}})}}{\varphi }} \right)} \right), \\ {{D}_{{22}}}(x) = - {{I}_{{\text{s}}}} + \frac{\varphi }{{R\left( {1 + \frac{1}{\alpha }} \right) + {{R}_{n}}}} \times \\ \times \,\,W\left( {{{I}_{{\text{s}}}}\frac{{R\left( {1 + \frac{1}{\alpha }} \right) + {{R}_{n}}}}{\varphi }\exp \left( {\frac{{x + {{I}_{{\text{s}}}}\left( {R\left( {1 + \frac{1}{\alpha }} \right) + {{R}_{n}}} \right)}}{\varphi }} \right)} \right). \\ \end{gathered} $

Здесь Uin – это напряжение входного усиливаемого сигнала, R – внутреннее сопротивление p–n-переходов n–p–n-транзистора (предполагается одинаковым для обоих переходов), α – коэффициент передачи тока, Rn – нагрузочное сопротивление в цепи коллектора, Rne – стабилизирующее сопротивление в цепи эмиттера, E – напряжение источника питания (включенного в коллекторную цепь), φ – температурный потенциал, Is – ток насыщения. W(z) – специальная функция Ламберта (ее теорию можно посмотреть в [3]). В работе [1] было показано, что расчетные данные (при расчетах использовалась формула (1)) близки к экспериментальным данным, а также к результатам компьютерного моделирования в приложении Orcad.

Напряжение коллектор-эмиттер Uce в рассмотренном усилителе определяется следующей формулой:

(3)
${{U}_{{ce}}} = E - {{R}_{n}}{{i}_{c}}.$

Если подставить в нее выражение (1), мы получим “усилительную” функцию Uce = F(Uin), которая при нулевом значении Rne имеет вид достаточно резкого перепада от напряжения питания к почти нулевому значению.

Цель статьи – использовать данную функцию (а также другие функции подобного рода) при выводе дифференциальных уравнений для двух симметричных схем, скомпонованных из пары усилителей с общим эмиттером (триггер и мультивибратор). Среди материалов, приведенных в открытых источниках, аналогичных уравнений для транзисторных схем автором не найдено (исследовались лишь уравнения для ламповых приборов, имеющие другую форму – см., например, [4]).

2. МЕТОДЫ И РЕЗУЛЬТАТЫ

Рассмотрим симметричный триггер, содержащий два биполярных транзистора (рис. 1). Пусть величины сопротивлений одинаковы: R1= R2= R. Коллекторные цепи обоих транзисторов аналогичны коллекторной цепи усилителя с общим эмиттером (без стабилизации). Будем считать, что присутствие первого усилителя с высоким входным сопротивлением не влияет на режим работы второго усилителя (и наоборот), поэтому расчеты “выходных” функций, произведенные для усилителя с общим эмиттером (в работах [12]), остаются верными. Если пренебречь током базы, “усилительную” функцию можно выразить через “токовую” функцию для коллекторного тока ic(Ube), определяемую формулой (1):

(4)
$F({{U}_{{be2}}}) = E - R{{i}_{c}}({{U}_{{be2}}}).$
Рис. 1.

Схема триггера.

Если применить закон Кирхгофа к контуру, выделенному на рис. 1 пунктиром (см. также работу [5], в которой рассматривалась статическая задача для триггера), можно получить соотношение

${{U}_{{be1}}} = F({{U}_{{be2}}}).$

Рассмотрев аналогично другую часть симметричной схемы, мы придем к выводу, что

(5)
${{U}_{{be1}}} = F(F({{U}_{{be1}}})).$

По своей форме это есть трансцендентное или алгебраическое (все зависит от вида функции F) уравнение, решив которое, мы найдем несколько значений для напряжения база-эмиттер. Ясно, что уравнение (5) может быть применено не только к триггеру на биполярных транзисторах, но и к другим физическим системам с похожим поведением, и выражение для усилительной функции, имеющей характер нисходящего перепада между двумя стационарными состояниями, может быть аппроксимировано посредством выражений, отличающихся от (4). Для удобства дальнейшей графической иллюстрации решения данного уравнения рассмотрим пример, когда усилительная функция аппроксимируется не формулой (4), а следующим выражением:

(6)
$y = F(x) = {{E}_{0}} + \frac{1}{{\sqrt {1 + B\exp (2Gx)} }}.$

Графики для функции F(x) и для результата ее повторного действия F(F(x)) приведены на рис. 2а, 2б (см. [5], в нашем случае E = 0.02, G = 10, B = = 0.004). Из рис. 2б видно, что имеются три решения уравнения (5). Двум крайним решениям соответствуют устойчивые равновесные значения напряжения триггера. Среднее решение – неустойчивое, система при малейшем отклонении от этой точки “сваливается” в одно из крайних состояний равновесия.

Рис. 2.

График усилительной функции (а) и графическое решение уравнения (б) для точек равновесия.

Чтобы обосновать утверждения об устойчивости состояний равновесия, рассмотрим динамическую задачу. Для этого необходимо учесть емкости p–n-переходов в биполярных транзисторах. Расчетная схема триггера должна принять вид, изображенный на рис. 3. Если пренебречь токами базы, схему с емкостями можно упростить (рис. 4), заменив пару емкостей C1 и C2 на емкость C = C1C2/(C1+ C2) (аналогично – для C3,C4).

Рис. 3.

Схема триггера, учитывающая емкости pn-переходов.

Рис. 4.

Упрощенный учет емкостей.

Для транзистора, изображенного на рис. 4 справа, выполняется соотношение, следующее из закона Кирхгофа:

${{U}_{{ce2}}} = E - R{{i}_{2}}$

(мы помним, что величина сопротивления резистора R2 равна R).

Так как (в пренебрежении током базы первого транзистора)

(7)
${{i}_{2}} = {{i}_{c}} + {{i}_{0}},$

то можно записать

${{U}_{{ce2}}} = E - R{{i}_{c}} - R{{i}_{0}}.$

Ток, проходящий через конденсатор, зависит от производной напряжения коллектор–эмиттер:

${{i}_{0}} = C\frac{{d{{U}_{{ce2}}}}}{{dt}}.$

Поэтому получаем

${{U}_{{ce2}}} = (E - R{{i}_{c}}) - RC\frac{{d{{U}_{{ce2}}}}}{{dt}}.$

Выражение в скобках является введенной нами выше усилительной функцией F(Ube2) (напомню: мы считаем, что подключение “симметричного” усилителя с большим входным сопротивлением оставляет в силе формулы для коллекторного тока исходного усилителя). Вследствие этого последнее уравнение можно записать в виде

${{U}_{{ce2}}} = F({{U}_{{be2}}}) - RC\frac{{d{{U}_{{ce2}}}}}{{dt}}.$

Из закона Кирхгофа следует, что Ube1= Uce2. Аналогично Ube2= Uce1. Поэтому имеем

${{U}_{{ce2}}} = F({{U}_{{ce1}}}) - RC\frac{{d{{U}_{{ce2}}}}}{{dt}}.$

Если рассмотреть первый транзистор (слева на рис. 4), мы получим такое же выражение:

${{U}_{{ce1}}} = F({{U}_{{ce2}}}) - RC\frac{{d{{U}_{{ce1}}}}}{{dt}}.$

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

(8)
$\left\{ {\begin{array}{*{20}{c}} {\frac{{d{{y}_{1}}}}{{dt}} = \frac{1}{{RC}}(F({{y}_{2}}) - {{y}_{1}}),} \\ {\frac{{d{{y}_{2}}}}{{dt}} = \frac{1}{{RC}}(F({{y}_{1}}) - {{y}_{2}}).} \end{array}} \right.$

Стационарные состояния данной системы определяются уже рассмотренным выше уравнением (5):

${{y}_{1}} = F(F({{y}_{1}})).$

Остановимся на значениях емкости C, входящей в уравнения (8). Как известно (см. [6] и [7]), емкости p–n-переходов зависят от напряжения, причем полная емкость перехода равна сумме двух емкостей: барьерной и диффузионной. При обратных и небольших прямых напряжениях доминирует барьерная емкость, а при достаточно больших прямых напряжениях – диффузионная емкость (пропорциональная току через переход, который растет при увеличении прямого напряжения). В данной работе мы имеем дело с составной емкостью C, состоящей из последовательно соединенных емкостей коллекторного и эмиттерного перехода. Предположим, что эта емкость в исследуемом нами триггере является функцией напряжения коллектор–эмиттер и что нам известна соответствующая вольт-фарадная характеристика. Следует ожидать, что при малых значениях Uce, когда оба транзисторных перехода открыты, емкость будет иметь диффузионный характер. При больших же значениях Uce, когда оба перехода закрыты, емкость будет барьерной (и ее величина будет существенно меньшей). Обозначим явно зависимость емкости от напряжения в полученных уравнениях:

(9)
$\left\{ {\begin{array}{*{20}{c}} {\frac{{d{{y}_{1}}}}{{dt}} = \frac{1}{{RC({{y}_{1}})}}(F({{y}_{2}}) - {{y}_{1}}),} \\ {\frac{{d{{y}_{2}}}}{{dt}} = \frac{1}{{RC({{y}_{2}})}}(F({{y}_{1}}) - {{y}_{2}}).} \end{array}} \right.$

Если емкость C не принимает бесконечно больших значений, то стационарные состояния такой усложненной системы будут определяться тем же уравнением (5). Пусть y10 и y20 – одно из стационарных состояний. Для определения его типа используем общую теорию динамических систем, изложенную, например, в [4]. Обозначим правую часть первого уравнения (9) как P(y1, y2), а правую часть второго уравнения – как Q(y1, y2). Из общей теории следует, что типы состояний равновесия определяются корнями λ следующего характеристического уравнения:

$\left| {\begin{array}{*{20}{c}} {a - \lambda }&b \\ c&{d - \lambda } \end{array}} \right| = 0.$

Здесь

$\begin{gathered} a = {{\left. {\frac{{\partial P}}{{\partial {{y}_{1}}}}} \right|}_{{{{y}_{{10}}},{{y}_{{20}}}}}},\,\,\,\,b = {{\left. {\frac{{\partial P}}{{\partial {{y}_{2}}}}} \right|}_{{{{y}_{{10}}},{{y}_{{20}}}}}}, \\ c = {{\left. {\frac{{\partial Q}}{{\partial {{y}_{1}}}}} \right|}_{{{{y}_{{10}}},{{y}_{{20}}}}}},\,\,\,d = {{\left. {\frac{{\partial Q}}{{\partial {{y}_{2}}}}} \right|}_{{{{y}_{{10}}},{{y}_{{20}}}}}}. \\ \end{gathered} $

В нашем случае

$\begin{gathered} a = - \frac{1}{{RC({{y}_{{10}}})}},\,\,\,\,b = \frac{1}{{RC({{y}_{{10}}})}}{{\left. {\frac{{dF({{y}_{2}})}}{{d{{y}_{2}}}}} \right|}_{{{{y}_{{20}}}}}}, \\ c = \frac{1}{{RC({{y}_{{20}}})}}{{\left. {\frac{{dF({{y}_{1}})}}{{d{{y}_{1}}}}} \right|}_{{{{y}_{{10}}}}}},\,\,\,\,d = - \frac{1}{{RC({{y}_{{20}}})}}. \\ \end{gathered} $

Характеристическое уравнение имеет следующий вид:

$\begin{gathered} {{\lambda }^{2}} + \frac{{C({{y}_{{10}}}) + C({{y}_{{20}}})}}{{RC({{y}_{{10}}})C({{y}_{{20}}})}}\lambda + \\ + \,\,\frac{1}{{{{R}^{2}}C({{y}_{{10}}})C({{y}_{{20}}})}}\left( {1 - {{{\left. {\frac{{dF({{y}_{2}})}}{{d{{y}_{2}}}}} \right|}}_{{{{y}_{{20}}}}}}{{{\left. {\frac{{dF({{y}_{1}})}}{{d{{y}_{1}}}}} \right|}}_{{{{y}_{{10}}}}}}} \right) = 0. \\ \end{gathered} $

Его решения:

${{\lambda }_{{1,2}}} = \frac{{C({{y}_{{10}}}) + C({{y}_{{20}}})}}{{2RC({{y}_{{10}}})C({{y}_{{20}}})}}\left\{ { - 1 \pm \sqrt {1 - \frac{{4C({{y}_{{10}}})C({{y}_{{20}}})}}{{{{{(C({{y}_{{10}}}) + C({{y}_{{20}}}))}}^{2}}}}\left( {1 - {{{\left. {\frac{{dF({{y}_{2}})}}{{d{{y}_{2}}}}} \right|}}_{{{{y}_{{20}}}}}}{{{\left. {\frac{{dF({{y}_{1}})}}{{d{{y}_{1}}}}} \right|}}_{{{{y}_{{10}}}}}}} \right){\kern 1pt} } } \right\}.$

Так как функция F имеет “квазиступенчатый” вид (вид перепада от величины, близкой к напряжению питания, к почти нулевой величине), в области крайних стационарных точек производные, стоящие под знаком квадратного корня, близки к нулю (если края функции F достаточно плоские). В центральной же стационарной точке эти производные, напротив, велики (если перепад достаточно крутой). Поэтому, если функция F достаточно близка к идеальной “ступеньке”, двум крайним стационарным состояниям соответствуют действительные отрицательные значения λ, а центральному состоянию – действительные значения, одно из которых отрицательно, а другое – положительно. Из общей теории следует, что крайние стационарные состояния будут устойчивыми узлами, а центральное – седлом (неустойчивым состоянием равновесия).

При численном решении дифференциальных уравнений (9) для начальных оценок разумно использовать одно из постоянных значений емкости в одной из стационарных точек или некоторое среднее постоянное значение. При уточнении расчетов нужно использовать более сложные модели.

Далее рассмотрим схему симметричного мультивибратора (рис. 5). Такое устройство можно рассматривать в качестве “наследника” симметричного триггера – схема рис. 4 дополнена парой конденсаторов и парой сопротивлений с номиналами C0 и R0. Конденсаторы включены в разрыв проводников, соединяющих коллекторы и базы транзисторов. Дополнительные сопротивления соединяют противоположные от коллекторов выводы этих конденсаторов с полюсом питания. Так же как и в триггере, мы учитываем присутствие емкостей p–n-переходов транзистора (пренебрегая токами базы).

Рис. 5.

Схема симметричного мультивибратора.

Важную роль в этой задаче будет играть напряжение на дополнительных конденсаторах C0. Уточним, что под данным напряжением мы будем понимать разность между потенциалом того вывода конденсатора, который ближе к краю схемы (к правому краю для правого конденсатора или к левому краю для левого), и потенциалом другого вывода (того, который ближе к центру). В соответствии с этим на рисунке выбрано направление токов i32 и i31.

Так же как и в случае триггера, мы будем составлять уравнения для правого транзистора (обозначенного индексом “2”). Рассмотрим контур, проходящий через источник питания E, резистор R в правой части схемы, коллектор второго транзистора и его эмиттер (контур направим против часовой стрелки). Тогда имеем

(10)
${{i}_{2}}R + {{U}_{{ce2}}} = E.$

В узле A

${{i}_{2}} = {{i}_{{c2}}} + {{i}_{{02}}} + {{i}_{{32}}}.$

С учетом этого предыдущее соотношение (10) записывается в виде

(11)
${{U}_{{ce2}}} = E - {{i}_{{c2}}}R - RC\frac{{d{{U}_{{ce2}}}}}{{dt}} - {{i}_{{32}}}R.$

Теперь рассмотрим контур, проходящий через базу и эмиттер первого транзистора, источник питания E и правый резистор R0 (контур также направим против часовой стрелки). Из законов Кирхгофа следует:

$ - {{i}_{{42}}}{{R}_{0}} + {{U}_{{be1}}} = E.$

Так как в пренебрежении током базы (для первого транзистора) i42= i32, то соотношение (11) можно записать так:

(12)
${{U}_{{ce2}}} = E - {{i}_{{c2}}}R - RC\frac{{d{{U}_{{ce2}}}}}{{dt}} + \frac{{E - {{U}_{{be1}}}}}{{{{R}_{0}}}}R.$

Первые два слагаемых в правой части соотношения (12) представляют собой уже рассмотренную ранее усилительную функцию (замечания, касающиеся высокого входного сопротивления, здесь сохраняют свою силу):

$E - {{i}_{{c2}}}R = F({{U}_{{be2}}}).$

Используем контур, проходящий через правый конденсатор C0, через коллектор и эмиттер второго транзистора и через базу и эмиттер первого (тоже против часовой стрелки), чтобы выразить Ube1 в выражении (12) через “правые” параметры:

${{U}_{{be1}}} = {{U}_{{ce2}}} - {{V}_{2}},$

где V2 – напряжение на правом конденсаторе C0.

Далее рассмотрим контур, проходящий по часовой стрелке через левый конденсатор C0, через коллектор и эмиттер первого транзистора и через базу и эмиттер второго, чтобы выразить Ube2 через более удобные для нас величины (для этого можно воспользоваться и симметрией схемы):

${{U}_{{be2}}} = {{U}_{{ce1}}} - {{V}_{1}}.$

В итоге получим уравнение

$\begin{gathered} {{U}_{{ce2}}} = - RC\frac{{d{{U}_{{ce2}}}}}{{dt}} + F({{U}_{{ce1}}} - {{V}_{1}}) + \\ + \,\,\frac{R}{{{{R}_{0}}}}{{V}_{2}} - \frac{R}{{{{R}_{0}}}}{{U}_{{ce2}}} + \frac{R}{{{{R}_{0}}}}E. \\ \end{gathered} $

Используем контур, проходящий через источник питания, правый резистор R0, правый конденсатор C0 и через коллектор и эмиттер правого транзистора (против часовой стрелки), чтоб записать еще одно уравнение

${{R}_{0}}{{C}_{0}}\frac{{d{{V}_{2}}}}{{dt}} + {{V}_{2}} = {{U}_{{ce2}}} - E.$

Аналогичные уравнения можно записать и для другой половины схемы (они будут выглядеть точно так же, с точностью до замены индекса “2” на “1”). Приведем все четыре уравнения для величин V1, V2, Uce1, Uce2:

(13)
$\begin{gathered} RC\frac{{d{{U}_{{ce1}}}}}{{dt}} = - {{U}_{{ce1}}} + F({{U}_{{ce2}}} - {{V}_{2}}) - \frac{R}{{{{R}_{0}}}}({{U}_{{ce1}}} - E - {{V}_{1}}), \\ {{R}_{0}}{{C}_{0}}\frac{{d{{V}_{1}}}}{{dt}} + {{V}_{1}} = {{U}_{{ce1}}} - E, \\ RC\frac{{d{{U}_{{ce2}}}}}{{dt}} = - {{U}_{{ce2}}} + F({{U}_{{ce1}}} - {{V}_{1}}) - \frac{R}{{{{R}_{0}}}}({{U}_{{ce2}}} - E - {{V}_{2}}), \\ {{R}_{0}}{{C}_{0}}\frac{{d{{V}_{2}}}}{{dt}} + {{V}_{2}} = {{U}_{{ce2}}} - E. \\ \end{gathered} $

В этих уравнениях, как и в уравнениях для триггера, емкость C зависит от напряжения (в первом уравнении системы (13) C зависит от Uce1, в третьем – от Uce2). Для начальных оценок можно использовать одно из постоянных значений емкости в одной из стационарных точек задачи о триггере, в который мультивибратор превратится, если конденсаторы C0 замкнуть накоротко, а резисторы R0 убрать, или же, как и в задаче о триггере, использовать некоторое среднее постоянное значение емкости.

Рассматриваемая система (13) является системой с двумя степенями свободы (и с четырехмерным фазовым пространством). Отметим, что схожая задача для лампового мультивибратора рассматривалась еще в [4] (но там были не только другие активные элементы, но и другая конфигурация соединительных проводников). В процессе решения авторы [4] перешли к системе с одной степенью свободы (с двумя неизвестными функциями). Здесь мы останемся в рамках модели с двумя степенями свободы и проведем по возможности более простое рассмотрение движений исследуемой нами системы без сокращения числа переменных.

Если R0$ \gg $ R, то уравнения записываются в виде

(14)
$\begin{gathered} RC\frac{{d{{U}_{{ce1}}}}}{{dt}} = - {{U}_{{ce1}}} + F({{U}_{{ce2}}} - {{V}_{2}}), \\ {{R}_{0}}{{C}_{0}}\frac{{d{{V}_{1}}}}{{dt}} + {{V}_{1}} = {{U}_{{ce1}}} - E, \\ ..., \\ .... \\ \end{gathered} $

С целью упрощения вида уравнений рассмотрим возможный предельный случай. Если внутренняя емкость транзистора C (в [4] ее аналогом является внутренняя паразитная емкость лампы) стремится к нулю, то первое уравнение становится не дифференциальным, а конечным трансцендентным (или алгебраическим – в зависимости от вида функции F):

(15)
${{f}_{1}}({{U}_{{ce1}}},{{U}_{{ce2}}},{{V}_{2}}) = - {{U}_{{ce1}}} + F({{U}_{{ce2}}} - {{V}_{2}}) = 0.$

Функция f1 задает фазовую поверхность “вырожденной” системы (как и аналогичная ей функция f2). Именно по этим фазовым поверхностям проходят “медленные” движения системы. При этом уравнения “быстрых” движений можно записать так:

$\begin{gathered} RC\frac{{d{{U}_{{ce1}}}}}{{dt}} = {{f}_{1}}({{U}_{{ce1}}},{{U}_{{ce2}}},{{V}_{2}}), \hfill \\ RC\frac{{d{{U}_{{ce2}}}}}{{dt}} = {{f}_{2}}({{U}_{{ce2}}},{{U}_{{ce1}}},{{V}_{1}}). \hfill \\ \end{gathered} $

Если рассматривать точки фазового пространства с конечными величинами функций f1 и f2, то при стремлении C к нулю производные величин Uce1 и Uce2 должны стремиться к бесконечности. Ясно, что линии движения при этом почти параллельны плоскости V1V2 фазового пространства. Это означает, что V1 и V2 (напряжения на конденсаторах) при таких движениях (скачках) не меняются.

Исследуем “медленные” движения. Выразим Uce1 из уравнения (15) и подставим во второе уравнение (14):

(16)
${{R}_{0}}{{C}_{0}}\frac{{d{{V}_{1}}}}{{dt}} + {{V}_{1}} = F({{U}_{{ce2}}} - {{V}_{2}}) - E.$

Вспомним, что усилительная функция F для усилителя с общим эмиттером имеет вид резкого перепада от напряжения питания E к почти нулевому значению (в несколько десятых долей вольта). Допустим, что медленные движения нашей системы происходят на почти постоянных “хвостах” усилительной функции (а “быстрые” – при перепаде ее значения). Пусть для определенности F в формуле (16) равно E. Тогда

(17)
${{R}_{0}}{{C}_{0}}\frac{{d{{V}_{1}}}}{{dt}} + {{V}_{1}} = 0.$

Решение этого уравнения:

${{V}_{1}} = A\exp \left( { - \frac{t}{{{{R}_{0}}{{C}_{0}}}}} \right).$

Здесь A – произвольная константа. Данное решение соответствует разрядке конденсатора.

Если F = 0, то вместо (17) мы будем иметь уравнение

(18)
${{R}_{0}}{{C}_{0}}\frac{{d{{V}_{1}}}}{{dt}} + {{V}_{1}} = - E.$

Его решение

${{V}_{1}} = - E - A\exp \left( { - \frac{t}{{{{R}_{0}}{{C}_{0}}}}} \right).$

Конденсатор при этом находится в процессе зарядки.

Дифференциальные уравнения (13) имеют решение в виде незатухающих колебаний не для всех значений параметров. В заключение работы мы приведем численное решение, полученное средствами пакета Wolfram Mathematica 9.0 (рис. 6). Усилительная функция F при этом определяется формулой (6), параметры задачи следующие: E0= 0.02, G = 5, B = 0.0004, R0= 1.1, C0= = 100, R = 0.01, C = 0.0033, E = 1 + E0.

Рис. 6.

Численные решения для напряжений Uce1(t) (кривая 1), V1(t) (кривая 2), V2(t) (кривая 3), Uce2(t) (кривая 4).

Из графиков видно, что при данных значениях параметров предположение о постоянстве Uce1(t) (и Uce2(t)) в пределах полупериода колебаний (фактически сделанное ранее) является чересчур грубым приближением. Но движения все равно подразделяются на “быстрые” и “медленные”.

Численное решение уравнений с усилительной функцией, которая записывается через формулу (1), имеет похожий вид.

ВЫВОДЫ

В работе получены дифференциальные уравнения для симметричного транзисторного триггера и мультивибратора. При этом использовалось полученное ранее автором выражение для токовой характеристики транзисторного усилителя с общим эмиттером (хотя можно использовать и выражения другого вида – в работе приведен пример подобной аппроксимации). Уравнения триггера в общем виде могут быть записаны следующим образом (для неизвестных функций y1, y2):

$\left\{ {\begin{array}{*{20}{c}} {\frac{{d{{y}_{1}}}}{{dt}} = g(F({{y}_{2}}) - {{y}_{1}}),} \\ {\frac{{d{{y}_{2}}}}{{dt}} = g(F({{y}_{1}}) - {{y}_{2}}).} \end{array}} \right.$

Здесь g – числовой параметр (равный 1/RC) или же функция, зависящая от y1 в первом уравнении и от y2 во втором уравнении, а F(y) – “квазиступенчатая” функция, график которой выглядит как достаточно резкий перепад от большего значения к меньшему.

Уравнения мультивибратора (для неизвестных функций y1, y2, z1, z2) имеют следующий вид:

$\left\{ \begin{gathered} \frac{{d{{y}_{1}}}}{{dt}} = g\{ (F({{y}_{2}} - {{z}_{2}}) - {{y}_{1}}) - \frac{R}{{{{R}_{0}}}}({{y}_{1}} - {{z}_{1}} - E)\} , \hfill \\ \frac{{d{{y}_{2}}}}{{dt}} = g\{ (F({{y}_{1}} - {{z}_{1}}) - {{y}_{2}}) - \frac{R}{{{{R}_{0}}}}({{y}_{2}} - {{z}_{2}} - E)\} , \hfill \\ \frac{{d{{z}_{1}}}}{{dt}} = - {{g}_{0}}({{z}_{1}} - {{y}_{1}} + E), \hfill \\ \frac{{d{{z}_{2}}}}{{dt}} = - {{g}_{0}}({{z}_{2}} - {{y}_{2}} + E). \hfill \\ \end{gathered} \right.$

Здесь g0= 1/R0C0, R, R0 и E – числовые параметры.

Для дифференциальных уравнений триггера проведено исследование характера состояний равновесия. Для дифференциальных уравнений мультивибратора проведено исследование “медленных” и “быстрых” движений. Приведен пример численного решения.

Список литературы

  1. Китаев А.Е. // Радиотехника. 2020. № 1. С. 74.

  2. Китаев А.Е. // Радиотехника. 2017. № 10. С. 189.

  3. Дубинов А.Е., Дубинова И.Д., Сайков С.К. W-функция Ламберта и ее применение в математических задачах физики. Саров: РФЯЦ, 2006.

  4. Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний. М.: Физматлит, 1959.

  5. Китаев А.Е. // Нелинейный мир. 2018. № 5. С. 16.

  6. Бонч-Бруевич В.Л., Калашников С.Г. Физика полупроводников. М.: Наука, 1990.

  7. Степаненко И.П. Основы теории транзисторов и транзисторных схем. М.: Энергия, 1977.

Дополнительные материалы отсутствуют.