Затухающие колебания. Характеристики гармонических колебаний

Вернемся к основной теме — переходным решениям. Пе­реходными решениями называются решения дифференциаль­ного уравнения, соответствующие ситуации, когда внешняя сила не действует, но система тем не менее не находится в покое. (Конечно, лучше всего решать задачу, когда сила не действует, а система покоится, покоится — ну и пусть покоится!) Соответ­ствующие переходным решениям колебания можно вызвать так: заставить силу поработать, а потом выключить ее. Что тогда случится с осциллятором? Сначала подумаем, как будет вести себя система с очень большой Q. Если сила действовала долго, то запасенная энергия была постоянной и работа тратилась лишь для того, чтобы поддержать ее. Предположим теперь, что мы выключили силу, тогда трению, которое раньше поглощало энергию поставщика, питаться больше нечем — кормильца-то нет. И трение начинает пожирать запасенную осциллятором энергию. Пусть добротность системы Q/2p=1000. Это значит, что работа, произведенная за цикл, равна 1/1000 запасенной энергии. Пожалуй, разумно предположить, что при не поддерживае­мых внешней силой колебаниях за каждый цикл будет теряться одна тысячная часть имеющейся к началу цикла энергии. Будем считать, что при больших Q изменение энергии описывается угаданным нами приближенным уравнением (мы еще вернемся к этому уравнению и сделаем его совсем верным!)

Уравнение это приближенное, потому что оно справедливо только для больших Q. За каждый радиан система теряет 1 /Q часть запасенной энергии Е. Значит, за промежуток времени dt энергия уменьшится в (wdt/Q раз (частота появляется при переводе радианов в настоящие секунды). А какая это частота? Предположим, что система устроена очень жестко, поэтому даже при действии силы она сколько-нибудь заметно колеблется толь­ко со своей собственной частотой. Поэтому будем считать, что w — это резонансная частота w₀. Таким образом, из уравнения (24.8) следует, что запасенная энергия меняется
следующим образом:

Теперь нам известно значение энергии в любой момент. Какой будет приближенная формула, определяющая амплитуду коле­баний как функцию времени? Той же самой? Нет! Потенциаль­ная энергия пружины изменяется как квадрат смещения, кинетическая энергия — как квадрат скорости; это приводит к тому, что полная энергия пропорциональна квадрату сме­щения. Таким образом, смещение (амплитуда колебаний) будет уменьшаться с половинной скоростью. Иначе говоря, мы ожидаем, что решение в случае затухающего переходного дви­жения будет выглядеть как колебание с частотой, близкой к ре­зонансной частоте w₀; амплитуда этого колебания будет умень­шаться как ехр(-g t/2)

Эта формула и рис. 24.1 дают представление о том, чего следует ожидать, а теперь приступим к точному анализу движе­ния, т. е. к решению дифференциального уравнения движения.

Рис. 24.1. Затухающие колебания.

Как же решить уравнение (24.1), если выкинуть из него внешнюю силу? Будучи физиками, мы интересуемся не столько методом, сколько самим решением. Поскольку мы люди уже опытные, попытаемся представить решение в виде экспоненци­альной кривой, х=Аexp (iat). (Почему мы так поступили? Оттого, что экспоненту легче всего дифференцировать!) Подставим это выражение в (24.1), помня о том, что каждое дифференцирование х по времени сводится к умножению на ia [напомним, что F(t)=0]. Сделать это очень легко, и наше уравнение примет вид

(-a²+i g a+ w ²₀)Ае^iat=0. (24.11)

Левая часть равенства должна быть равна нулю все время, но это возможно только в двух случаях: а) А=0, однако это даже и не решение: ведь тогда все покоится, или б)

Если мы сможем решить это уравнение и найти a, то мы найдем и решение, амплитуда которого А не обязательно равна нулю!

Чтобы не думать о том, как извлечь квадратный корень, предположим, что g меньше w₀, и поэтому w²₀-g²/4 — положи­тельная величина. Беспокоит другое: почему мы получили два решения! Им соответствуют

и

Займемся пока первым решением, предположив, что мы ничего не знаем о том, что квадратный корень принимает два значе­ния. В этом случае смещение х равно x₁=Aexp(i a ₁t), где А — произвольная постоянная. Чтобы сократить запись, введем специальное обозначение для входящего в a _t квадратного корня:

Так, и , или, если воспользоваться замечательным свойством экспоненты,

Итак, система осциллирует с частотой w_g, которая в точности не равна частоте w₀, но практически близка к ней, если система достаточно добротна. Кроме того, амплитуда колебаний экспо­ненциально затухает! Если взять действительную часть (24.16), то мы получим

Это решение очень напоминает угаданное нами решение (24.10), вот только частота немного другая, w_g. Но это лишь небольшая поправка, значит, первоначальная идея была правильной.

И все-таки не все благополучно! А не благополучно то, что су­ществует второе решение.

Этому решению соответствует a₂, и оно отличается от пер­вого лишь знаком w_g

Что все это значит? Скоро мы докажем, что если x₁ и х₂— воз­можные решения (24.1) при F(t)=0, то х₁+х₂—тоже решение этого уравнения! Таким образом, общее решение имеет вид

Теперь можно спросить: «А, собственно, зачем нам беспокоить себя еще одним решением, если нас вполне устраивало первое? К чему эти дополнительные решения, если мы все равно должны взять только действительную часть?» Мы знаем, что нужно взять действительную часть, но откуда математика знает, что мы хо­тим взять действительную часть? Когда у нас была внешняя сила F(t), то мы ее дополнили искусственной силой, и она каким-то образом управляла мнимой частью уравнения. Но когда мы по­ложили F(t)=0, то соглашение о том, что, каково бы ни было х, нужно взять только его действительную часть, стало нашим лич­ным делом, и математическое уравнение об этом ничего не знало. В мире физики есть только действительные решения, но реше­ние, которому мы так радовались, комплексно. Уравнению не из­вестно, что мы делаем совершенно неожиданный шаг и отбираем только действительную часть, и оно предлагает нам еще, так сказать, комплексно сопряженное решение, чтобы, сложив оба решения, мы получили настоящее действительное решение; вот для чего мы взяли еще и a₂. Чтобы х было действительным, Ввхр(-iw_gt) должно быть комплексно сопряженным к Aexp(iw_gt) числом, тогда мнимая часть исчезнет. Таким образом, В долж­но быть комплексно сопряжено с А, поэтому наше решение имеет вид

Значит, наши колебания — это колебания с фазовым сдвигом и, как полагается, с затуханием.

Посмотрим, как выглядят переходные колебания. Для этого соберем цепь, изображенную на рис. 24.2.

Рис. 24,2. Электрическая цепь для демонстраций переходных колебаний.

В этой цепи разность потенциалов между концами индуктивности L поступает в осцил­лоскоп. Неожиданное включение рубильника S включает допол­нительное напряжение и вызывает в осцилляторной цепи переходные колебания. Эти колебания аналогичны колебаниям механического осциллятора, вызванными неожиданным ударом. Сама цепь представляет собой электрический аналог механи­ческого осциллятора с затуханием, и мы можем наблюдать коле­бания при помощи осциллоскопа. Он покажет нам кривые, анализом которых мы и займемся. На рис. 24.3—24.6 представ­лены кривые затухающих колебаний, полученные на экране осциллоскопа. На рис. 24.3 показаны затухающие колебания в цепи с большой Q, т. е. с малым значением g.

Рис. 24.3. Затухающие коле­бания.

В такой цепи ко­лебания затухают не очень быстро; мы видим довольно длинную синусоиду с медленно убывающим размахом.

Теперь давайте посмотрим, что произойдет, если мы будем уменьшать Q, так что колебания должны затухать быстрее. Чтобы уменьшить Q, увеличим сопротивление цепи R. При повороте ручки сопротивления колебания действительно зату­хают скорее (рис. 24.4).

Рис. 24.4. Колебания затухают быстрее.

Если еще увеличить сопротивление, то колебания затухнут еще быстрее (рис. 24.5).

Фиг, 24.5. Колебания почти исчезли.

Но если сопротив­ление увеличить сверх некоторого предела, колебаний мы вооб­ще не увидим. А может быть, нам просто отказывают глаза? Увеличим еще сопротивление и получим тогда кривую, пред­ставленную на рис. 24.6; по ней можно лишь с натяжкой сказать, что в цепи произошли колебания, ну разве что одно.

Рис. 24.6. Колебаний нет.

Можем ли мы математически объяснить это явление?

Сопротивление механического осциллятора, конечно, про­порционально g. В нашем случае g — это R/L. Теперь, если уве­личивать g, то в столь приятных нам решениях (24.14) и (24.15) наступает беспорядок; когда g /2 становится больше w₀, реше­ния приходится записывать по-другому:

Это снова два решения, которые приводят нас к решениям exp(ia₁t) и ехр(ia₂t). Подставив теперь a₁, получим

Никаких колебаний. Чисто экспоненциальное убывание. То же самое дает и второе решение

Заметим, что квадратный корень не может превысить g/2; даже если w₀=0, оба члена равны. Если же w²₀отличается от g/²/4, то квадратный корень меньше g//2 и выражение в круглых скобках всегда положительно. Это очень хорошо! Почему? Да потому что если бы это выражение было отрицательным, то е пришлось бы возводить в положительную степень и мы по­лучили бы возрастающее со временем решение. Но при увели­чении в цепи сопротивления колебания не могут возрастать, зна­чит, мы избегли противоречия. Итак, мы получили два решения; оба решения экспоненциально затухают, но одно из них стре­мится «умереть» гораздо скорее. Общее решение, конечно, пред­ставляет собой комбинацию обоих решений, а значения коэф­фициентов А и В зависят от того, как начинаются колебания, каковы начальные условия. В нашей цепи случилось так, что А — отрицательное число, а В — положительное, поэтому на экране осциллоскопа мы увидели разность двух экспонент.

Давайте обсудим, как найти коэффициенты А и В (или А и A*), если известны начальные условия. Предположим, что в момент t=0 нам известны смещение х=х₀ и скорость dx/dt=v₀. Если в соотношения

подставить значения t=0, х=х₀, dx/dt=v₀ и воспользо­ваться тем, что е⁰=еⁱ⁰=1, то мы получим

x₀=A+A*=2A_R,

Значит,

Таким образом, зная начальные условия, мы полностью определили А и А*, а значит, и кривую переходного решения. Можно записать решение и по-другому. Вспомним, что

e^iq+e^-iq=2cosq и e^{i q} - e^{-i q} =2isin q, тогда

где wg=+Ö(w²₀-(g²/4). Мы получили формулу затухающих колебаний. Такая формула нам не понадобится, однако отметим ее особенности, справедливые и в более общих случаях.

Прежде всего поведение системы, на которую не действует внешняя сила, описывается суммой (суперпозицией) временных экспонент [мы записали их в виде exp(iat)]. Такое решение хорошо передает истинное положение вещей. В общем случае a — это комплексное число, и его мнимая часть соответствует затуханию колебаний. Наконец, тесная математическая связь синусоидальных и экспоненциальных функций, о которой го­ворилось в гл. 22, физически часто проявляется в переходе от колебаний к чисто экспоненциальному затуханию при крити­ческих значениях некоторых параметров системы (в нашем случае это было сопротивление g).