Уравнение гармонического осциллятора с затуханием и его решение

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

Взяв за основу ту же модель, добавим в нее силу вязкого трения. Сила вязкого трения направлена против скорости движения груза относительно среды и пропорциональна этой скорости. Тогда полная сила, действующая на груз, записывается так:

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

Здесь введено обозначение: . Коэффициент γ носит название постоянной затухания. Он тоже имеет размерность частоты.

Решение же распадается на три случая.

При малом трении (γ − γt

• При сильном же трении γ > ω₀ решение выглядит следующим образом:

, где

Критическое затухание примечательно тем, что именно при критическом затухании осциллятор быстрее всего стремится в положение равновесия. Если трение меньше критического, он дойдет до положения равновесия быстрее, однако «проскочит» его по инерции, и будет совершать колебания. Если трение больше критического, то осциллятор будет экспоненциально стремиться к положению равновесия, но тем медленнее, чем больше трение.

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

Затухание осциллятора также часто характеризуют безразмерным параметром, называемым добротностью. Добротность обычно обозначают буквой Q. По определению, добротность равна:

Чем больше добротность, тем медленнее затухают колебания осциллятора.

У осциллятора с критическим затуханием добротность равна 0,5. Соответственно, добротность указывает характер поведения осциллятора. Если добротность больше 0,5, то свободное движение осциллятора представляет собой колебания; со временем он пересечёт положение равновесия неограниченное количество раз. Добротность меньше или равная 0,5 соответствует неколебательному движению осциллятора; в свободном движении он пересечёт положение равновесия не более одного раза.

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

Также добротность примерно равна количеству колебательных циклов, за которое амплитуда колебаний уменьшается в e раз, умноженному на π.

В случае колебательного движения затухание еще характеризуют такими параметрами, как:

• Время жизни колебаний, оно же время затухания, оно же время релаксации. τ — время, за которое амплитуда колебаний уменьшится в e раз.

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

• Логарифмический декремент затухания. Определяется как логарифм отношения двух последовательных максимальных отклонений в одну сторону. . Величина, обратная d, есть количество колебаний, которое пройдёт за время затухания τ.

Коэффициент затухания

величина, характеризующая скорость затухания колебаний

Закон сохранения энергии — основной закон природы, заключающийся в том, что энергия изолированной (замкнутой) системы сохраняется во времени. Другими словами, энергия не может возникнуть из ничего и не может в никуда исчезнуть, она может только переходить из одной формы в другую. Закон сохранения энергии встречается в различных разделах физики и проявляется в сохранении различных видов энергии. Например, в классической механике закон проявляется в ]] закон сохранения энергии называется первым началом термодинамики и говорит

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

Добро́тность — характеристика колебательной системы, определяющая остроту резонанса и показывающая, во сколько раз запасы энергии в реактивных элементах контура больше, чем потери энергии на активных за один период колебаний.

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

Общая формула для добротности любой колебательной системы:

,

• f — частота колебаний
• W — энергия, запасённая в колебательной системе
• P_d — рассеиваемая мощность.

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

Обратимый процесс можно в любой момент заставить протекать в обратном направлении, изменив какую-либо независимую переменную на бесконечно малую величину.

Обратимые процессы дают наибольшую работу. Большую работу от системы вообще получить невозможно. Это придает обратимым процессам теоретическую важность. На практике обратимый процесс реализовать невозможно. Он протекает бесконечно медленно, и можно только приблизиться к нему.

Следует отметить, что термодинамическая обратимость процесса отличается от химической обратимости. Химическая обратимость характеризует направление процесса, а термодинамическая — способ его проведения.

понятие равновесного состояния и обратимого процесса играют большую роль в термодинамике. Все количественные выводы термодинамики применимы только к равновесным состояниям и обратимым процессам.

Вынужденные колебания — колебания, происходящие под действием внешней силы, меняющейся во времени.

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

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

Резона́нс (фр. resonance, от лат. resono — откликаюсь) — явление резкого возрастания амплитуды вынужденных колебаний, которое наступает при приближении частоты внешнего воздействия к некоторым значениям (резонансным частотам), определяемым свойствами системы. При помощи явления резонанса можно выделить и/или усилить даже весьма слабые периодические колебания. Резонанс — явление, заключающееся в том, что при некоторой частоте вынуждающей силы колебательная система оказывается особенно отзывчивой на действие этой силы.

Но это далеко не полное определение явления резонанса. Для более детального восприятия этой категории необходимы некоторые факты из теории дифференциальных уравнений и математического анализа. В теории обыкновенных дифференциальных уравнений известна проблема собственных векторов и собственных значений. Резонанс в динамической системе, описываемой дифференциальными уравнениями (и не только ими), формально наступает, когда проблема собственных значений приводит к кратным собственным числам. При этом в математическом аспекте не очень существенно, являются ли собственные числа комплексными или действительными. В физическом аспекте явление резонанса обычно связывают только с колебательными динамическими системами. Наиболее ярко понятие явления резонанса развито в современной теории динамических систем. Примером является известная теория Колмогорова-Арнольда-Мозера. Центральная проблема этой теории — вопрос сохранения квазипериодического или условно-периодического движения на торе (теорема КАМ). Эта теорема дала мощный толчок к развитию современной теории нелинейных колебаний и волн. В частности, стало ясно, что резонанс может и не наступить, хоть собственные числа совпадают или близки. Напротив, резонанс может проявиться в системе, где никакие собственные числа не совпадают, а удовлетворяют лишь определенным резонансным соотношениям или условиям синхронизма.

теплова́я маши́на

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

Цикл Карно́ — идеальный термодинамический цикл. Тепловая машина Карно, работающая по этому циклу, обладает максимальным КПД из всех машин, у которых максимальная и минимальная температуры осуществляемого цикла совпадает соответственно с максимальной и минимальной температурами цикла Карно.

КПД тепловой машины Карно

Количество теплоты, полученное рабочим телом от нагревателя при изотермическом расширении, равно

.

Аналогично, при изотермическом сжатии рабочее тело отдало холодильнику

.

Отсюда коэффициент полезного действия тепловой машины Карно равен

.

Из последнего выражения видно, что КПД тепловой машины Карно зависит только от температур нагревателя и холодильника. Кроме того, из него следует, что КПД может составлять 100 % только в том случае, если температура холодильника равна абсолютному нулю. Это невозможно, но не из-за недостижимости абсолютного нуля (этот вопрос решается только третьим началом термодинамики, учитывать которое здесь нет необходимости), а из-за того, что такой цикл или нельзя замкнуть, или он вырождается в совокупность двух совпадающих адиабат и изотерм.

Можно показать, что КПД любой тепловой машины, работающей по циклу, отличному от цикла Карно, будет меньше КПД тепловой машины Карно, работающей при тех же температурах нагревателя и холодильника.

Моме́нт и́мпульса (кинетический момент, угловой момент, орбитальный момент, момент количества движения) характеризует количество вращательного движения. Величина, зависящая от того, сколько массы вращается, как она распределена относительно оси вращения и с какой скоростью происходит вращение.

Замечание: момент импульса относительно точки — это псевдовектор, а момент импульса относительно оси — скалярная величина.

Следует учесть, что вращение здесь понимается в широком смысле, не только как регулярное вращение вокруг оси. Например, даже при прямолинейном движении тела мимо произвольной воображаемой точки, оно также обладает моментом импульса. Наибольшую роль момент импульса играет при описании собственно вращательного движения.

Момент силы (синонимы: крутящий момент; вращательный момент; вращающий момент) — физическая величина, характеризующая вращательное действие силы на твёрдое тело.

В физике момент силы можно понимать как «вращающая сила». В системе СИ единицами измерения для момента силы является ньютон-метр, хотя сантиньютон-метр (cN•m), футо-фунт (ft•lbf), дюйм-фунт (lbf•in) и дюйм-унция (ozf•in) также часто используются для выражения момента силы. Символ момента силы τ (тау). Момент силы иногда называют моментом пары сил, это понятие возникло в трудах Архимеда над рычагами. Вращающиеся аналоги силы, массы и ускорения есть момент силы, момент инерции и угловое ускорение соответственно. Сила, приложенная к рычагу, умноженная на расстояние до оси рычага, есть момент силы. Например, сила в 3 ньютона, приложенная к рычагу, расстояние до оси которого 2 метра, это то же самое, что 1 ньютон, приложенный к рычагу, расстояние до оси которого 6 метров. Более точно, момент силы частицы определяется как векторное произведение:

где — сила, действующая на частицу, а — радиус-вектор частицы!

Теплоемкость идеального газа — это отношение тепла, сообщенного газу, к изменению температуры δТ, которое при этом произошло.

Идеальный гармонический осциллятор. Уравнение идеального осциллятора и его решение. Амплитуда, частота и фаза колебаний

КОЛЕБАНИЯ. ВОЛНЫ. ОПТИКА

КОЛЕБАНИЯ

Лекция 1

ГАРМОНИЧЕСКИЕ КОЛЕБАНИЯ

Идеальный гармонический осциллятор. Уравнение идеального осциллятора и его решение. Амплитуда, частота и фаза колебаний

Колебание – один из самых распространённых процессов в природе и технике. Колебания – это процессы, повторяющиеся во времени. Колеблются высотные здания и высоковольтные провода под действием ветра, маятник заведённых часов и автомобиль на рессорах во время движения, уровень реки в течение года и температура человеческого тела при болезни. Звук – это колебания давления воздуха, радиоволны – периодические изменения напряжённости электрического и магнитного поля, свет – это тоже электромагнитные колебания. Землетрясения – колебания почвы, приливы и отливы – изменение уровней морей и океанов, вызываемые притяжением луны и т.д.

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

Первыми учёными, изучавшими колебания, были Галилео Галилей и Христиан Гюйгенс. Галилей установил независимость периода колебаний от амплитуды. Гюйгенс изобрёл часы с маятником.

Любая система, которая, будучи слегка выведена из положения равновесия, совершает устойчивые колебания, называется гармоническим осциллятором. В классической физике такими системами являются математический маятник в пределах малых углов отклонения, груз в пределах малых амплитуд колебаний, электрический контур, состоящий из линейных элементов ёмкости и индуктивности.

Гармонический осциллятор можно считать линейным, если смещение от положения равновесия прямо пропорционально возмущающей силе. Частота колебаний гармонического осциллятора не зависит от амплитуды. Для осциллятора выполняется принцип суперпозиции — если действуют несколько возмущающих сил, то эффект их суммарного действия может быть получен как результат сложения эффектов от действующих сил в отдельности.

Гармонические колебания описываются уравнением (рис.1.1.1)

(1.1.1)

где х -смещение колеблющейся величины от положения равновесия, А – амплитуда колебаний, равная величине максимального смещения, — фаза колебаний, определяющая смещение в момент времени , — начальная фаза, определяющая величину смещения в начальный момент времени, — циклическая частота колебаний.

Время одного полного колебания называется периодом, , где — число колебаний, совершенных за время .

Частота колебаний определяет число колебаний, совершаемых в единицу времени, она связана с циклической частотой соотношением , тогда период .

Скорость колеблющейся материальной точки

,

. (1.1.2)

Таким образом, скорость и ускорение гармонического осциллятора также изменяются по гармоническому закону с амплитудами и соответственно. При этом скорость опережает по фазе смещение на , а ускорение – на (рис.1.1.2).

Из сопоставления уравнений движения гармонического осциллятора (1.1.1) и (1.1.2) следует, что , или

. (1.1.3)

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

.

Отсюда .

Если периодически повторяющийся процесс описывается уравнениями, не совпадающими с (1.1.1), он н6азывается ангармоническим. Система, совершающая ангармонические колебания, называется ангармоническим осциллятором.

1.1.2. Свободные колебания систем с одной степенью свободы. Комплексная форма представления гармонических колебаний

В природе очень распространены малые колебания, которые система совершает вблизи своего положения равновесия. Если система, выведенная из положения равновесия, предоставлена себе, то есть на неё не действуют внешние силы, то такая система будет совершать свободные незатухающие колебания. Рассмотрим систему с одной степенью свободы.

Устойчивому равновесию соответствует такое положение системы, в котором её потенциальная энергия имеет минимум (q – обобщённая координата системы). Отклонение системы от положения равновесия приводит к возникновению силы , которая стремится вернуть систему обратно. Значение обобщённой координаты, соответствующей положению равновесия, обозначим , тогда отклонение от положения равновесия

Будем отсчитывать потенциальную энергию от минимального значения . Примем Полученную функцию разложим в ряд Маклорена и оставим первый член разложения, имеем: о

,

где . Тогда с учётом введённых обозначений:

, (1.1.4)

С учётом выражения (1.1.4) для силы, действующей на систему, получаем:

Согласно второму закону Ньютона, уравнение движения системы имеет вид: ,

, (1.1.5)

Выражений (1.1.5) совпадает с уравнением (1.1.3) свободных гармонических колебаний при условии, что

, (1.1.6)

и имеет два независимых решения: и , так что общее решение:

, или

,

где

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

Зависимость координаты колеблющейся системы от времени можно определить в виде вещественной части комплексного выражения , где A=Xe-iα– комплексная амплитуда, её модуль совпадает с обычной амплитудой, а аргумент – с начальной фазой.

1.1.3. Примеры колебательных движений различной физической природы

Колебания груза на пружине

Рассмотрим колебания груза на пружине, при условии, что пружина не деформирована за пределы упругости. Покажем, что такой груз будет совершать гармонические колебания относительно положения равновесия (рис.1.1.3). Действительно, согласно закону Гука, сжатая или растянутая пружина создаёт гармоническую силу:

где – коэффициент жёсткости пружины, – координата положения равновесия, х – координата груза (материальной точки) в момент времени , — смещение от положения равновесия.

Поместим начало отсчета координаты в положение равновесия системы. В этом случае .

Если пружину растянуть на величину х, после чего отпустить в момент времени t=0, то уравнение движения груза согласно второму закону Ньютона примет вид -kx =ma, или , и

(1.1.6)

Это уравнение совпадает по виду с уравнением движения (1.1.3) системы, совершающей гармонические колебания, его решение будем искать в виде:

. (1.1.7)

Подставим (1.17) в (1.1.6), имеем: то есть выражение (1.1.7) является решением уравнения (1.1.6) при условии, что

Если в начальный момент времени положение груза было произвольным, то уравнение движения примет вид:

.

Рассмотрим, как меняется энергия груза, совершающего гармонические колебания в отсутствие внешних сил (рис.1.14). Если в момент времени t=0 грузу сообщить смещение х=А, то его полная энергия станет равной потенциальной энергии деформированной пружины , кинетическая энергия равна нулю (точка 1).

На груз действует сила F= -kx, стремящаяся вернуть его в положение равновесия, поэтому груз движется с ускорением и увеличивает свою скорость, а, следовательно, и кинетическую энергию. Эта сила сокращает смещение груза х, потенциальная энергия груза убывает, переходя в кинетическую. Система «груз — пружина» замкнутая, поэтому её полная энергия сохраняется, то есть:

. (1.1.8)

В момент времени груз находится в положении равновесия (точка 2), его потенциальная энергия равна нулю, а кинетическая максимальна . Максимальную скорость груза найдём из закона сохранения энергии (1.1.8):

За счёт запаса кинетической энергии груз совершает работу против упругой силы – и пролетает положение равновесия. Кинетическая энергия постепенно переходит в потенциальную. При груз имеет максимальное отрицательное смещение –А, кинетическая энергия Wk=0, груз останавливается и начинает движение к положению равновесия под действием упругой силы F= -kx. Далее движение происходит аналогично.

Маятники

Под маятником понимают твёрдое тело, которое совершает под действием силы тяжести колебания вокруг неподвижной точки или оси. Различают физический и математический маятники.

Математический маятник – это идеализированная система, состоящая из невесомой нерастяжимой нити, на которой подвешена масса, сосредоточенная в одной материальной точке.

Математическим маятником, например, является шарик на длинной тонкой нити.

Отклонение маятника от положения равновесия характеризуется углом φ, который образует нить с вертикалью (рис.1.15). При отклонении маятника от положения равновесия возникает момент внешних сил (силы тяжести) : , где m – масса, – длина маятника

Этот момент стремится вернуть маятник в положение равновесия (аналогично квазиупругой силе) и направлен противоположно смещению φ, поэтому в формуле стоит знак «минус».

Уравнение динамики вращательного движения для маятника имеет вид: Iε= ,

.

Будем рассматривать случай малых колебаний, поэтому sin φ ≈φ, обозначим ,

имеем: , или , и окончательно

.

Это уравнение гармонических колебаний, его решение:

.

Частота колебаний математического маятника определяется только его длиной и ускорением силы тяжести, и не зависит от массы маятника. Период равен:

.

Если колеблющееся тело нельзя представить, как материальную точку, то маятник называют физическим (рис.1.1.6). Уравнение его движения запишем в виде:

.

В случае малых колебаний , или =0 , где . Это уравнение движения тела, совершающего гармонические колебания. Частота колебаний физического маятника зависит от его массы, длины и момента инерции относительно оси, проходящей через точку подвеса.

Обозначим . Величина называется приведённой длинной физического маятника. Это длина математического маятника, период колебаний которого совпадает с периодом данного физического маятника. Точка на прямой, соединяющей точку подвеса с центром масс, лежащая на расстоянии приведённой длины от оси вращения, называется центром качания физического маятника (О’). Если маятник подвесить в центре качания, то приведённая длина и период колебаний будут теми же, что и в точке О. Таким образом, точка подвеса и центр качания обладают свойствами взаимности: при переносе точки подвеса в центр качения прежняя точка подвеса становится новым центром качения.

Математический маятник, который качается с таким же периодом, как и рассматриваемый физический, называется изохронным данному физическому маятнику.

1.1.4. Сложение колебаний (биения, фигуры Лиссажу). Векторное описание сложения колебаний

Сложение одинаково направленных колебаний можно производить методом векторных диаграмм. Любое гармоническое колебание можно представить в виде вектора следующим образом. Выберем ось х с началом отсчета в точке О (рис.1.1.7)

Из точки О построим вектор , который составляет угол с осью х. Пусть этот вектор поворачивается с угловой скоростью . Проекция вектора на ось Х равна:

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

Рассмотрим два гармонических колебания одинакового направления и одинаковой циклической малой , заданные векторами и . Смещения по оси Х равны:

результирующий вектор имеет проекцию и представляет собой результирующее колебание (рис.1.1.8), по теореме косинусов Таким образом, сложение гармонических колебаний производится сложением векторов.

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

.

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

Из уравнения движения следует: ,

. (1.1.9)

Из уравнения (1.1.9) можно получить уравнение эллипса (рис.1.1.9):

Рассмотрим частные случаи этого уравнения:

1. Разность фаз колебаний α= 0. При этом т.е. или Это уравнение прямой, и результирующее колебание происходит вдоль этой прямой с амплитудой (рис.1.1.10).

2. Если разность фаз то уравнение (1.1.9) переходит в уравнение эллипса, приведенного к координатным осям, При материальная точка движется по окружности, уравнение которой (рис.1.1.11).

3. Если частоты колебаний неодинаковы, то материальная точка описывает фигуры Лиссажу (рис.1112).

Рассмотрим сложение колебаний одного направления, частоты которых мало отличаются друг от друга. В этом случае результирующее движение можно рассматривать как гармоническое колебание с пульсирующей амплитудой. Такие колебания называются биениями.

Пусть частота одного колебания , второго . Амплитуды обоих колебаний одинаковы и равны а. Начальные фазы равны нулю. В таком случае уравнения колебаний имеют вид:

Сложим эти выражения:

(1.1.10)

График функции х(t) представлен на рис. 1.1.13. Множитель меняется гораздо медленнее, чем , поэтому (1.1.10) можно рассматривать как гармоническое колебание частоты , амплитуда которого меняется по некоторому периодическому закону

Частота изменения амплитуды – частота биений – равна разности частот складываемых колебаний .

Энергия колебаний

Смещение колеблющейся точки от положения равновесия, описывается уравнением:

ее ускорение равно второй производной от смещения по времени тогда сила, действующая на колеблющуюся точку, по второму закону Ньютона равна

— то есть сила пропорциональна смещению х и направлена против смещения к положению равновесия. Эта сила называется возвращающей силой. В случае груза на пружине возвращающей силой является сила упругости, в случае математического маятника – составляющая силы тяжести.

Возвращающая сила по характеру подчиняется закону Гука F= -kx, где

– коэффициент возвращающей силы. Тогда потенциальная энергия колеблющейся точки равна:

(постоянную интегрирования выбирают равной нулю, чтобы при х=0 энергия Wn=0).

где , тогда

Полная механическая энергия равна сумме кинетической и потенциальной энергий, и в случае свободных колебаний без трения сохраняется (рис.1.1.15). Когда материальная точка совершает колебания, кинетическая энергия переходит в потенциальную, и наоборот. В крайних точках (х = ±А) скорость , кинетическая энергия равна нулю, и полная энергия равна потенциальной:

Таким образом, полная механическая энергия гармонического осциллятора пропорциональна квадрату амплитуды колебаний.

В положении равновесия (х=0) потенциальная энергия переходит в кинетическую:

В промежуточных точках полная энергия равна

а скорость

На рисунке 1.1.16 приведена кривая потенциальной энергии , горизонтальная линия соответствует полной энергии. Расстояние от этой линии до кривой равно кинетической энергии. Движение ограничено значениями х, заключёнными в пределах от –А до +А.

Средние за период значения кинетической и потенциальной энергии одинаковы и равны , так что средняя полная энергия системы равна полной энергии системы ( средние значения ).

АНГАРМОНИЧЕСКИЙ ОСЦИЛЛЯТОР

Обыкновенные дифференциальные уравнения

Дифференциальные уравнения – это язык физики

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

ДУ второго порядка часто встречаются в физике: закон Ньютона, уравнение Шредингера, диффузия.

• restart:
• E2:=diff(y(t),t$2)=-g;
• S2:=dsolve(,y(t));

Это очевидный результат.

Заметьте несколько деталей.

Для ввода начальных условий для производной от y применена D(y) , что означает dy/dt ( D – это общий символ производной для Maple), затем добавили (0) , чтобы сказать, что производная вычислена при t=0 . Это означает, что diff и D – это производные, но между ними есть тонкое различие, которое надо понимать. Оно обсуждено в разделе о производных (см. D и diff ). Поэтому или запоминайте синтаксис команды dsolve (выше), или запомните, где его легко найти. Кроме того, учтите, что применяемая в пакете Physics команда diff может давать другой результат – см. описание команд пакета Physics.

Это второе знаменитое ДУ второго порядка.

Maple решает его в общем виде с произвольными константами:

Предупреждение: эта форма решения опасна, поскольку нет гарантии, что при каждом запуске рабочего листа с командой dsolve в таком виде получим на выходе _C1 с синусом и _C2 с косинусом. В Maple есть собственная логика такого выбора, которая даже для более простых ДУ кажется случайной. Поэтому не давайте Maple выбирать, делайте это сами. Для этого скопируйте решение мышкой в новую строку и выберите неизвестные коэффициенты, например избавьтесь от неудобного y(t)= , заменив его естественным присвоением:

Теперь каждый раз расчет величин A и B будет срабатывать одинаково, потому что «насильственно» установлено, что A – с синусом, а B – с косинусом.

Если константы определяются начальными условиями, есть форма команды dsolve , которая будет определять их автоматически. Например, нужно решить ДУ E3 с начальными условиями y(0) = 1 и dy/dt(0) = 2. Тогда пишем:

Применим assign , чтобы Maple дал выражение для y(t) :

Теперь построим график решения:

Не работает. Рамка есть, но нет функции. И нет догадок, что не так. Для отладки надо бы подставить несколько значений t в y(t) . Но вспомним, что y(t) выглядит как выражение, но не совсем им является, и что подстановка значений для t запустит его. Рассмотрим y(t) с этой точки зрения и прикинем, можно ли увидеть, что не так:

Команда plot даст численные значения t , но ведь просто включено в выражение и не имеет численного значения. Это и есть ошибка: мы не задали численные значения !

OК, теперь сделаем так:

Решение этой задачи было достаточно легким. Но так же легко Maple решает и более сложные задачи.

Усложним условия и рассмотрим

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

и мы получаем общее решение:

Где же тогда появятся синусы и косинусы? Проблема в том, что Maple не знает, насколько велики и τ. Подумаем с точки зрения физического смысла. Пусть естественная частота порядка 2π, тогда период N = 2π/ = 1 с. Теперь пусть осциллятор (подразумеваем – маятник) находится в моторном масле при 50 градусах ниже нуля. Тогда характеристическое время затухания 0.05 с (трение столь велико, что движение остановится за 0.05 с). Интуиция должна подсказать, что если толкнуть маятник назад и освободить его, то ожидается некоторое покачивание, т. е. маятник медленно проплывет к своему вертикальному положению и останется в нем. Посмотрим, что скажет Maple:

OK, Maple согласен с интуицией.

Теперь изменим затухание и представим, что произойдет, если нагреть масло или использовать вместо него WD-40 (жидкая смазка) или, может быть, воздух. Маятник начнет колебаться с малым затуханием. Например, если τ = 20, то Maple дает:

Теперь он колеблется дольше и медленно затухает. Но какое чудо в решении превратило экспоненту в синусы и косинусы? Вспомним формулу Эйлера:
т. е. в экспоненте константы стали комплексными, поэтому решение перешло от распада к колебаниям.

Это поднимает вопрос о величине τ, при которой происходит переход от чистого затухания к затухающим колебаниям. Посмотрите на общий вид решения и увидите, что в экспонентах есть квадратный корень Если аргумент под корнем положителен, то оба фундаментальных решения в сумме есть затухающие экспоненты, и получаем затухание. Если аргумент под корнем отрицателен, то результат комплексный и экспонента комплексная, т. е. это синусы и косинусы, тогда получаем колебания. Переход между двумя режимами происходит, если квадратный корень = 0, т. е. когда τ = 1/(2). Этот особый случай называется критическим затуханием.

Пусть Maple решит задачу о гармоническом осцилляторе с затуханием для случая критического затухания при = 1 и нарисует график функции y(t) для t = 0..20 при начальных условиях y(0) = 1 и dy/dt(0) = 0. График будет выглядеть в точности как для осциллятора с затуханием при полном отсутствии намеков на то, что это точно граница между затуханием и колебаниями. Для проверки этого увеличьте τ на 2 % и снова постройте график. Растяните окно вниз, чтобы можно было видеть, что y(t) становится немного отрицательным, показывая начало колебаний (нужны очень серьезные размеры окна).

Примените Maple, чтобы получить общие решения ДУ второго порядка:

(a) ,

(b) ,

(c) .