9400

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

комплексное число α + iβ. Рассмотрим для начала случай действительной характеристики, т.е. когда правая часть имеет вид

Будем искать решение в таком же виде, что и правая часть уравнения, т.е.

где Qn (x) – многочлен, коэффициенты которого подлежат определению

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

соотношения между характеристикой правой части α и корнями характеристического уравнения

λ 2+ a1λ + a2 = 0.

∙ Пусть α не совпадает с корнями характеристического уравнения. Тогда в равенстве (46.4) слева многочлен тоже степени n , как и справа. Приравнивая коэффициенты при одинаковых степенях, получим систему

41

уравнения, то в левой части равенства (46.4) будет многочлен степени n −1. В этом случае умножим искомый многочлен на x . Степень многочлена увеличится, а число неизвестных коэффициентов не изменится. Задача свелась к предыдущей.

∙ И, наконец, если характеристика α совпадёт с двукратным корнем характеристического уравнения, то искомый многочлен умножаем на x2 .

Итак, если правая часть неоднородного уравнения имеет вид (46.2), то его частное решение отыскиваем в виде

y ( x) = xmQn ( x)eαx ,

где m = 0,1, 2 число совпадений корней характеристического уравнения с характеристикой правой части α , а Qn (x) – многочлен степени n с неопределёнными коэффициентами, которые определяются подстановкой

решение ищем в виде

y (x) = x( A1x + A2 ) .

Найдем

y′ = 2 A1x + A2 , y′′ = 2 A1

и подставив в уравнение, получим тождество

−6 A1x + (2 A1 − 3A2 ) ≡ x .

Приравнивая теперь коэффициенты многочленов справа и слева, получим A1 = −1/ 6 , A2 = 1/ 9 и в результате

y (x) = −(1/ 6)x2 + (1/ 9) x .

В случае комплексной характеристики частное решение неоднородного уравнения ищем в виде

y ( x) = xmeα x (Qn ( x) cos βx + Rn (x) sin βx) ,

степени n , коэффициенты которых подлежат определению. Заметим, что в решение входит «полный» тригонометрический многочлен, содержащий как синус, так и косинус, и что перед каждой тригонометрической функцией в качестве множителя находится «свой» многочлен.

В заключение этого раздела подчеркнем, что если правая часть уравнения (46.1) представляет собой сумму двух функций

y (x) = y1 (x) + y2 (x) .

Это легко проверяется подстановкой y (x) в уравнение

( y1 + y2 )′′ + a1 ( y1 + y2 )′ + a2 ( y1 + y2 ) =

=y1′′ + a1 y1′ + a2 y1 + y2′′ + a1 y2′ + a2 y2 = f1 + f2 .

46.2.Метод вариаций произвольных постоянных – это метод

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

y′′ + a1 y′ + a2 y = f (x)

с помощью известного общего решения однородного уравнения

y′′ + a1 y′ + a2 y = 0 .

Если метод неопределённых коэффициентов применим для правых частей, имеющих специальный вид, то рассматриваемый метод является

и запишем частное решение неоднородного уравнения

y(x) = ∫j1(x)dx × y1(x) + ∫j2 (x)dx × y2 (x) .

λ1,2 = ±i . Поэтому общее решение однородного уравнения имеет вид

yодн = C1 cos x + C2 sin x .

Применяя метод вариаций, будем искать частное решение неоднородного уравнения в виде

y(x) = u1(x)cos x + u2 (x)sin x .

Для производных искомых функций составляем систему

Умножим первое из уравнений системы на cos x , второе – на ( − sin x ) и сложим. Тогда получим

и после интегрирования

u1 (x) = −x , u2 ( x) = ln sin x .

Итак, общее решение уравнения

y = C1 cos x + C2 sin x - x cos x + sin x × ln sin x .

Лекция 47. Биения и резонанс

45

Рассмотрим теперь случай, когда на колебательную систему воздействует периодическая внешняя сила. Этот случай встречается наиболее часто: давление газа на поршень в двигателе внутреннего сгорания, удары рельса на стыке по колесу железнодорожного вагона, действие струи газа или пара на лопатку турбины, удары волн о корпус судна и т.д. Для простоты анализа будем предполагать, что сопротивление среды отсутствует. В этом случае уравнение линейного осциллятора y′′ + a1 y′ + a2 y = f (t) , примет вид

где ω1 – частота колебаний вынуждающей силы.

Решение однородного уравнения приведено (см.45.4). Это

Сначала пусть частота вынуждающей силы не совпадает с частотой

и общее решение уравнения (47.1) представляет собой сумму двух гармонических колебаний с разными частотами

Посмотрим, что произойдет, когда частота вынуждающей силы и собственная частота близки, т.е. ω ≈ ω1 . Вернемся к записи общего решения уравнения (47.1) в виде

46

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

ωsin ω1t − ω1sin ω t ≈ ω (sinω1t − sinω t) =

=w 2cos ω1+ ω t ×sin ω1− ω t » 2w sin ω1− ω t × cos w t .

второй – cos ωt быстро изменяющаяся функция. В приближённом решении уравнения

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

Пример «сложения» таких колебаний приведён на рисунке 47.1. Теперь пусть частота внешней вынуждающей силы и собственная

частота совпадают, т.е. ω = ω1 . Найдем решение уравнения

y′′ + ω2 y = Asin ωt .

Характеристика правой части уравнения совпадает с корнем характеристического уравнения, поэтому частное решение неоднородного уравнения ищем в виде

y(t) = t(M cosωt + N sinωt) .

47

Рис. 47.2

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

48

y(t) = e−0,1t ( 451sin t + 1000 cost) + 50 sin t − 1000 cost 401 401 401 401

или

y(t) ≈ e−0,1t (9 /8sin t + 5/ 2cost) +1/8sin t − 5/ 2cost .

49

Как видим, хотя амплитуда внешнего воздействия на колебательную систему всего 0,5, тем не менее эта внешняя сила «раскачала» систему до колебаний с амплитудой примерно 2,5. В дальнейшем система будет колебаться с частотой вынуждающей силы, а амплитуда колебаний расти не будет (см. рис. 47.4).

50