Дифференциальные уравнения RLC-цепей
В электрической цепи, содержащей
сопротивление \(R\),
индуктивность \(L\) и
емкость \(C\), могут возбуждаться электрические колебания. С точки зрения топологии чаще всего рассматриваются два вида электрических цепей:
последовательная \(RLC\)-цепь (рисунок \(1\)) и
параллельная \(RLC\)-цепь (рисунок \(2\)).
Выведем дифференциальное уравнение, описывающее закон изменения тока в
последовательной \(RLC\)-цепи.
Напряжения \({V_R},{V_C},{V_L},\) соответственно, на резисторе \(R,\) конденсаторе \(C\) и катушке индуктивности \(L\) выражаются формулами \[ {{V_R}\left( t \right) = RI\left( t \right),}\;\; {{V_C}\left( t \right) = \frac{1}{C}\int\limits_0^t {I\left( \tau \right)d\tau } ,}\;\; {{V_L}\left( t \right) = L\frac{{dI}}{{dt}}.} \] Из второго закона Кирхгофа следует, что \[{V_R}\left( t \right) + {V_C}\left( t \right) + {V_L}\left( t \right) = E\left( t \right),\] где \(E\left( t \right)\) − электродвижущая сила (э.д.с.) источника питания.
В случае постоянной э.д.с. \(E\) после подстановки выражений для \({V_R},{V_C}\) и \({V_L},\) и последующего дифференцирования получаем следующее дифференциальное уравнение: \[\frac{{{d^2}I\left( t \right)}}{{d{t^2}}} + \frac{R}{L}\frac{{dI\left( t \right)}}{{dt}} + \frac{1}{{LC}}I\left( t \right) = 0.\] Если ввести обозначения \(2\beta = {\large\frac{R}{L}\normalsize},\;\omega _0^2 = {\large\frac{1}{{LC}}\normalsize},\) то уравнение записывается в виде \[\frac{{{d^2}I}}{{d{t^2}}} + 2\beta \frac{{dI}}{{dt}} + \omega _0^2I = 0.\] Данное дифференциальное уравнение совпадает с уравнением, описывающим
затухающие колебания грузика на пружинке. Следовательно, в последовательной \(RLC\)-цепи при определенных значениях параметров также могут возникать затухающие колебания.
Теперь рассмотрим
параллельную \(RLC\)-цепь и выведем для нее аналогичное дифференциальное уравнение.
По первому закону Кирхгофа полный ток будет равен сумме токов через сопротивление \(R,\) катушку индуктивности \(L\) и конденсатор \(C\) (рисунок \(2\)): \[{I_R}\left( t \right) + {I_L}\left( t \right) + {I_C}\left( t \right) = I\left( t \right).\] Учитывая, что \[ {{I_R} = \frac{V}{R},}\;\;\; {{I_L} = \frac{1}{L}\int\limits_0^t {Vd\tau } ,}\;\;\; {{I_C} = C\frac{{dV}}{{dt}},} \] для случая постоянного полного тока \(I\left( t \right) = {I_0}\) получаем следующее дифференциальное уравнение \(2\)-го порядка относительно переменной \(V:\) \[ {\frac{V}{R} + \frac{1}{L}\int\limits_0^t {Vd\tau } + C\frac{{dV}}{{dt}} = {I_0},}\;\; {\Rightarrow C\frac{{{d^2}V}}{{d{t^2}}} + \frac{1}{R}\frac{{dV}}{{dt}} + \frac{1}{L}V = 0.} \] Как видно, мы снова приходим к уравнению, описывающему затухающие колебания. Таким образом, колебательный режим может возникать и в
параллельных \(RLC\)-цепях.
Простейший колебательный контур. Формула Томсона
В простейшем случае, когда омическое сопротивление равно нулю (\(R = 0\)) и источник э.д.с. отсутствует (\(E = 0\)), колебательный контур состоит лишь из конденсатора \(C\) и катушки индуктивности \(L\) и описывается дифференциальным уравнением \[\frac{{{d^2}I}}{{d{t^2}}} + \omega _0^2I = 0,\;\; \text{где}\;\;\omega _0^2 = \frac{1}{{LC}}.\] В таком контуре будут происходить
незатухающие электрические колебания с периодом \[{T_0} = \frac{{2\pi }}{{{\omega _0}}} = 2\pi \sqrt {LC} .\] Данная формула называется
формулой Томсона в честь английского физика Уильяма Томсона (\(1824-1907\)), который теоретически вывел ее в \(1853\) году.
Затухающие колебания в последовательной \(RLC\)-цепи
Выше мы получили дифференциальное уравнение второго порядка, описывающее
затухающие колебания в последовательном \(RLC\)-контуре, которое записывается как \[\frac{{{d^2}I}}{{d{t^2}}} + \frac{R}{L}\frac{{dI}}{{dt}} + \frac{1}{{LC}}I = 0.\] Соответствующее характеристическое уравнение имеет вид \[{\lambda ^2} + \frac{R}{L}\lambda + \frac{1}{{LC}} = 0.\] Его корни вычисляются по формулам: \[ {{\lambda _{1,2}} = \frac{{ - \frac{R}{L} \pm \sqrt {\frac{{{R^2}}}{{{L^2}}} - \frac{4}{{LC}}} }}{2} } = { - \frac{R}{{2L}} \pm \sqrt {{{\left( {\frac{R}{{2L}}} \right)}^2} - \frac{1}{{LC}}} } = { - \beta \pm \sqrt {{\beta ^2} - \omega _0^2} ,} \] где величина \(\beta = \large\frac{R}{{2L}}\normalsize\) называется
коэффициентом затухания, а \({\omega_0}\) −
резонансной частотой колебательного контура.
В зависимости от значений параметров \(R, L, C\) могут возникнуть три режима.
Случай 1. \({R^2} > \large\frac{{4L}}{C}\normalsize\)
В этом случае оба корня характеристического уравнения \({\lambda_1}\) и \({\lambda_2}\) действительны, различны и отрицательны. Общее решение дифференциального уравнения определяется формулой \[I\left( t \right) = {C_1}{e^{{\lambda _1}t}} + {C_2}{e^{{\lambda _2}t}}.\] В этом режиме ток монотонно уменьшается, приближаясь к нулю (рисунок \(3\)).
Случай 2. \({R^2} = \large\frac{{4L}}{C}\normalsize\)
Данный режим можно назвать
граничным или
критическим. Здесь оба корня характеристического уравнения совпадают, но при этом являются действительными и отрицательными. Общее решение уравнения выражается функцией \[ {I\left( t \right) = \left( {{C_1}t + {C_2}} \right){e^{ - \beta t}} } = {\left( {{C_1}t + {C_2}} \right){e^{ - {\large\frac{R}{{2L}}\normalsize} t}}.} \] В начале процесса ток может даже возрастать, но в дальнейшем он быстро уменьшается по экспоненциальному закону.
Случай 3. \({R^2} < \large\frac{{4L}}{C}\normalsize\)
В этом случае корни характеристического уравнения являются комплексно-сопряженными. В электрической цепи возникают затухающие колебания. Закон изменения тока имеет вид \[I\left( t \right) = {e^{ - \beta t}}\left( {A\cos \omega t + B\sin \omega t} \right),\] где величина \(\beta = \large\frac{R}{{2L}}\normalsize\) − как и выше,
коэффициент затухания, \(\omega = \sqrt {{\large\frac{1}{{LC}}\normalsize} - {{\left( {\large\frac{R}{{2L}}\normalsize} \right)}^2}} \) − частота колебаний, \(A, B\) − постоянные интегрирования, зависящие от начальных условий. Заметим, что частота затухающих колебаний \(\omega\) меньше резонансной частоты \({\omega_0}\) колебательного контура. Типичный вид кривой \(I\left( t \right)\) в этом режиме также представлен на рисунке \(3.\)
Вынужденные колебания и резонанс
Если колебательный контур содержит генератор с периодически изменяющейся э.д.с., то в нем устанавливаются
вынужденные колебания. Если э.д.с. \(E\) источника тока изменяется по закону \[E\left( t \right) = {E_0}\cos \omega t,\] то дифференциальное уравнение вынужденных колебаний в последовательной \(RLC\)-цепи записывается в виде \[ {\frac{{{d^2}q\left( t \right)}}{{d{t^2}}} + \frac{R}{L}\frac{{dq\left( t \right)}}{{dt}} + \frac{1}{{LC}}q\left( t \right) = \frac{1}{L}{E_0}\cos \omega t}\;\; {\text{или}\;\;\frac{{{d^2}q}}{{d{t^2}}} + 2\beta \frac{{dq}}{{dt}} + \omega _0^2q = \frac{{{E_0}}}{L}\cos \omega t,} \] где \(q\) − заряд конденсатора, \(2\beta = \frac{R}{L},\;\omega _0^2 = \frac{1}{{LC}}.\)
Данное уравнение аналогично уравнению вынужденных колебаний пружинного маятника, рассмотренного на странице
Механические колебания. Его общее решение представляет собой сумму двух слагаемых − общего решения однородного уравнения и частного решения неоднородного уравнения. При этом общее решение однородного уравнения описывает затухающий переходный процесс, по истечении которого в системе устанавливаются
вынужденные колебания. Эти вынужденные колебания будут происходить по закону \[ {q\left( t \right) } = {\frac{{{E_0}}}{{L\sqrt {{{\left( {\omega _0^2 - {\omega ^2}} \right)}^2} + 4{\beta ^2}{\omega ^2}} }}\cos \left( {\omega t + \varphi } \right) } = {\frac{{{E_0}}}{{\omega \sqrt {{R^2} + {{\left( {\omega L - \frac{1}{{\omega C}}} \right)}^2}} }}\cos \left( {\omega t + \varphi } \right),} \] где фаза \(\varphi\) определяется формулой \[ {\varphi = \arctan \left( { - \frac{{2\beta \omega }}{{\omega _0^2 - {\omega ^2}}}} \right) } = {\arctan \frac{R}{{\omega L - \frac{1}{{\omega C}}}}.} \] Зная закон изменения заряда \(q\left( t \right),\) легко найти закон изменения тока \(I\left( t \right):\) \[ {I\left( t \right) = \frac{{dq\left( t \right)}}{{dt}} } = { - \frac{{{E_0}}}{{\sqrt {{R^2} + {{\left( {\omega L - \frac{1}{{\omega C}}} \right)}^2}} }}\sin\left( {\omega t + \varphi } \right) } = {\frac{{{E_0}}}{{\sqrt {{R^2} + {{\left( {\omega L - \frac{1}{{\omega C}}} \right)}^2}} }}\cos\left( {\omega t - \theta } \right),} \] где введен угол \(\theta,\) равный \(\theta = - \left( {\varphi + \frac{\pi }{2}} \right).\) Угол \(\theta\) показывает отставание колебаний тока \(I\left( t \right)\) по отношению к колебаниям напряжения источника питания \(E\left( t \right) = {E_0}\cos \omega t.\)
Амплитуда тока \({I_0}\) и сдвиг фаз \(\theta\) определяются формулами \[ {{I_0} = \frac{{{E_0}}}{{\sqrt {{R^2} + {{\left( {\omega L - \frac{1}{{\omega C}}} \right)}^2}} }} = \frac{{{E_0}}}{Z},}\;\;\; {\theta = \arctan \frac{{\omega L - \frac{1}{{\omega C}}}}{R}.} \] Величина \(Z = \sqrt {{R^2} + {{\left( {\omega L - \large\frac{1}{{\omega C}}\normalsize} \right)}^2}} \) называется
полным сопротивлением или
импедансом контура. Она состоит из омического сопротивления \(R\) и
реактивного сопротивления \({\omega L - \large\frac{1}{{\omega C}}}\normalsize\) Импеданс колебательного контура в комплексной форме записывается как \[Z = R + i\left( {\omega L - \frac{1}{{\omega C}}} \right).\] Из полученных формул видно, что амплитуда установившихся колебаний тока будет максимальной когда \[\omega L = \frac{1}{{\omega C}}\;\;\text{или}\;\;\omega = {\omega _0} = \frac{1}{{\sqrt {LC} }}.\] При этом условии в колебательном контуре наступает
резонанс. Резонансная частота \({\omega_0}\) равна частоте свободных колебаний в контуре и не зависит от сопротивления \(R.\)
Формулу для амплитуды тока вынужденных колебаний можно преобразовать, выделив в явном виде зависимость от отношения частот \(\large\frac{\omega }{{{\omega _0}}}\normalsize,\) где \({\omega_0}\) − резонансная частота. В результате получаем \[\require{cancel} {{I_0} = \frac{{{E_0}}}{{\sqrt {{R^2} + {{\left( {\omega L - \frac{1}{{\omega C}}} \right)}^2}} }} } = {\frac{{\frac{{{E_0}}}{{{\omega _0}}}}}{{\frac{{\sqrt {{R^2} + {{\left( {\omega L - \frac{1}{{\omega C}}} \right)}^2}} }}{{{\omega _0}}}}} } = {\frac{{\frac{{{E_0}}}{{{\omega _0}}}}}{{\sqrt {\frac{{{R^2}}}{{\omega _0^2}} + {{\left( {\frac{\omega }{{{\omega _0}}}L - \frac{1}{{\omega {\omega _0}C}}} \right)}^2}} }} } = {\frac{{{E_0}\sqrt {LC} }}{{\sqrt {{R^2}LC + {{\left( {\frac{\omega }{{{\omega _0}}}L - \frac{1}{{\frac{\omega }{{{\omega _0}}}\frac{\cancel{C}}{{L\cancel{C}}}}}} \right)}^2}} }} } = {\frac{{{E_0}\sqrt {LC} }}{{\sqrt {{R^2}LC + {{\left( {\frac{\omega }{{{\omega _0}}}L - \frac{L}{{\frac{\omega }{{{\omega _0}}}}}} \right)}^2}} }} } = {\frac{{{E_0}\sqrt C }}{{\sqrt {{R^2}C + {{\left( {\frac{\omega }{{{\omega _0}}} - \frac{1}{{\frac{\omega }{{{\omega _0}}}}}} \right)}^2}} }}.} \] Типичные зависимости амплитуды тока от отношения частот \(\large\frac{\omega }{{{\omega _0}}}\normalsize\) при различных значениях \(R\) и \(C\) показаны на рисунках \(5\) и \(6.\) Данные графики построены при \(E = 100\;\text{В},\) \(L = 1\;\text{мГн},\) \(C = 10\;\text{мкФ}\) (на рисунке \(5\)), \(R = 10\;\text{Ом}\) (на рисунке \(6\)).
Резонансные свойства колебательного контура характеризуются
добротностью \(Q,\) которая численно равна отношению резонансной частоты \({\omega_0}\) к ширине резонансной кривой \(\Delta\omega\) на уровне убывания амплитуды в \(\sqrt 2\) раз (см. выше рисунок \(4\)).
В
последовательном колебательном контуре добротность вычисляется по формуле \[Q = \frac{1}{R}\sqrt {\frac{L}{C}} .\] Для
параллельной \(RLC\)-цепи добротность определяется обратным выражением: \[Q = R\sqrt {\frac{C}{L}} .\]