| Похоже, что один из основных авторов этой статьи тесная связь со своим предметом. Может потребоваться очистка для соответствия политике содержания Википедии, в частности нейтральная точка зрения. Пожалуйста, обсудите подробнее страница обсуждения. (Январь 2017 г.) (Узнайте, как и когда удалить этот шаблон сообщения) |
Синусоидальные колебания F = 0.01
В Уравнение Бирюкова (или же Осциллятор Бирюкова), названная в честь Вадима Бирюкова (1946), является нелинейным дифференциальное уравнение используется для модели с демпфированием генераторы.[1]
Уравнение имеет вид
![{ displaystyle { frac {d ^ {2} y} {dt ^ {2}}} + f (y) { frac {dy} {dt}} + y = 0, qquad qquad (1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aa20b16d7154b1e7306dc4abc6abc2a937ddd764)
куда ƒ(у) - кусочно-постоянная функция, положительная, за исключением малых у так как
![{ displaystyle f (y) = { begin {cases} -F, & | y | leq Y_ {0}; F, & | y |> Y_ {0}. end {cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/da32e7215432296790a5d2b635bf1b5ee2676a38)
![{ displaystyle F = { text {constant}}> 0, quad Y_ {0} = { text {constant}}> 0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/501778127c1d4809680f3c9cec456e3b007224f7)
Уравнение (1) является частным случаем Уравнение Лиенара; он описывает автоколебания.
Решение (1) через отдельные промежутки времени, когда f (y) постоянна, имеет вид[2]
![{ Displaystyle y_ {k} (t) = A_ {1, k} exp (s_ {1, k} t) + A_ {2, k} exp (s_ {2, k} t) qquad qquad (2)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9470b37a9347feb2fda1339bc87d1ea6d461a25a)
Здесь
, в
и
иначе. Выражение (2) можно использовать для действительных и комплексных значений
.
Решение первого полупериода при
является
Релаксационные колебания F = 4
![{ displaystyle y (t) = { begin {cases} y_ {1} (t), & 0 leq t <T_ {0}; y_ {2} (t), & T_ {0} leq t < Т / 2. End {case}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/53942e47a3c6e629e431d6b77f4f7bae2ea1e4f9)
![{ Displaystyle y_ {1} (t) = A_ {1, k} cdot exp (s_ {1, k} t) + A_ {2, k} cdot exp (s_ {2, k} t) ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8060695a8be6edf7b75881fbd3e1a4af44914578)
![{ displaystyle y_ {2} (t) = A_ {3, k} cdot exp (s_ {3, k} t) + A_ {4, k} cdot exp (s_ {4, k} t) .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f4714be4be9af48f1eb77922cf39ea9cefa8d4d3)
Решение второго полупериода:
![{ displaystyle y (t) = { begin {cases} -y_ {1} (tT / 2), & T / 2 leq t <T / 2 + T_ {0}; - y_ {2} (tT / 2), & T / 2 + T_ {0} leq t <T. end {cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0d5b8f21337375d7bb949c64efddd1c4cbbcdeb1)
Решение содержит четыре константы интегрирования
,
,
,
, Период
и граница
между
и
нужно найти. Граничное условие выводится из непрерывности
) и
.[3]
Таким образом, решение (1) в стационарном режиме получается путем решения системы алгебраических уравнений в виде
;
;
;
;
;
.
Константы интегрирования получаются Алгоритм Левенберга-Марквардта. С участием
,
, Уравнение (1) названный Генератор Ван дер Поля. Ее решение нельзя выразить элементарными функциями в замкнутой форме.
Рекомендации
- ^ Х. П. Гэвин, Метод Левенберга-Марквардта для нелинейных задач аппроксимации кривой наименьших квадратов (включая реализацию MATLAB)
- ^ Эроусмит Д. К., Плейс С. М. Динамические системы. Дифференциальные уравнения, карты и хаотическое поведение. Чепмен и Холл (1992)
- ^ Пилипенко А. М., Бирюков В. Н. «Исследование современных методов численного анализа эффективности автоколебательных схем», Радиоэлектроника, № 9, (2013). http://jre.cplire.ru/jre/aug13/9/text-engl.html