Алгоритм решения линейного однородного уравнения 2-го порядка с постоянными коэффициентами



 

- По заданному уравнению y''(x) + p*y'(x) + q*y(x) =0 составляется алгебраическое квадратное уравнение k2+p*k+q=0 и находятся его корни k1 и k2.       

 

- Полученные значения k1 , k2 используются при построении формул для решения ДУ (1):

В зависимости от знака дискриминанта D = p2/4 - q и корней k1 , k2 имеется три случая:

1. D>0 (т.е. k1≠k2 - действительные неравные числа)

Имеются два линейно независимых решения: y1(x) = ek1x,   y2(x) = ek2x

2. D=0 (т.е. k1=k2 - действительные равные числа)

Имеются два линейно независимых решения: y1(x)=ek1x, y2(x)=x*ek1x

3. D<0 (т.е. k1 ≠ k2 - комплексные неравные числа)

k1= a + i*b, k2= a - i*b, где a = -p/2, b = √|D|, i=√(-1)

Имеются два линейно независимых решения:      y1(x)=eаxcos(b*x),

                                                                        y2(x)= eаxsin(b*x)

Общее решение: у(х) = С1*y1(x)+С2* y2(x), где С1, С2 - произвольные константы.

Примеры.

1) D>0     

2) Решим уравнение

      

1a)           k2+4k+13=0

 k1= -4/2 + Ö(42/4-13)= -2 + Ö-9 = -2+√((-1)9)=-2+√(-1)√9)= -2+ i*3; k2= -2- i*3

 a=-2, b=3

Ответ:               y1(x)=exp(-2x)*cos(3*x), y2(x)=exp(-2x)*sin(3*x)

у = С1*exp(-2x)*cos(3*x) + С2*exp(-2x)*sin(3*x)

 

2) D < 0

 

 

                  

 

4)

      

 

Определение общего решения (общего интеграла):

3) D < 0 

Решение:  

 

4) D=0 

 

Упражнения (задачи с ответами)

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

 y'' + py' + qy = f(x)          (1)

Например: y(x)'' + 3*y(x)' + 2*y(x) =xe2x  

Общим решением дифференциального. уравнения 2-го порядка называется функция которая зависит от двух произвольных постоянных С1, С2 и удовлетворяет дифференциальному уравнению при любом конкретных значениях постоянных.  

 

 

Отыскание частного решения дифференциального уравнения, удовлетворяющего начальным условиям вида: у=у0 при х=х0 и у'=у1 при х=х0, называется задачей Коши.

 

Теорема. Если в правой части уравнения (1) f(x)=f1(x)+f2(x),то частное решение (1) равно сумме частного решения уравнения с f1(x) и решения уравнения с f2(x). 

                                                       * * *

Вид частного решения в зависимости от вида f(x) и корней характеристического уравнения.

0. Пусть f(x) = Pn(x), где Pn-многочлен n-й степени.

Pn(x)= a0*x0 +a1*x1+a2*x2+...+an*xn.

Например, P2(x)= a0*x0 +a1*x1+a2*x2.

Тогда:

1) Если ЧИСЛО 0 не равно ни одному из корней характеристического уравнения

k2+pk+q=0 (т.е. 0 ¹ k1 , 0 ¹ k2),

то решение имеет вид y*=Qn(x)-многочлен той же степени n.  

Алгоритм решения. Неизвестные множители многочлена Qn(x) определяются следующим образом. Выписываем выражение для у* (например, у*= (ax+b), где a и b - неизвестные множители). Находим формулы для y*' и y*''. (В нашем примере у*'=a, y*''=0 ).   Подставляем выражения для у*, y*' и y*'' в уравнение (1), приводим все соответствующие члены в правой и левой частях уравнения (например, -x+3x заменяем на 2х). Затем приравниваем множители при одинаковой степени x (x0,x1,x2,...,xn) в левой и правой части полученного тождества (например, если получено тождество (a+b)x+2b=3x-4, то оно равносильно системе двух уравнений: a+b=3 и 2b= - 4 ). Таким образом, мы получаем n+1 уравнений для определения n+1 неизвестных коэффициентов в многочлене Qn(x).

2) Если число 0 равно одному из неравных (при D>0)корней характеристического уравнения ( 0 = k1 , 0 ¹ k2), то решение имеет вид y*=x*Qn(x)  (вводится дополнительный множитель х). Неизвестные множители многочлена Qn(x) определяются аналогично.

3) Если 0 является двукратным корнем характеристического уравнения (0 =k1=k2), то решение имеет вид   y*=x2*Qn(x)  (вводится дополнительный множитель х2).

Неизвестные множители многочлена Qn(x) определяются аналогично.

 

Следует подчеркнуть, что рассмотренный случай (f(x)=Pn(x)) является частным по отношению к случаю  I., рассмотренному ниже.

I. Пусть f(x) = Pn(x)*exp(a*x), где Pn-многочлен n-й степени,exp(ax)= eax.

Pn(x)= a0*x0 +a1*x1+a2*x2+...+an*xn   a - действительное число.

Тогда:

1) Если a не равно ни одному из корней характеристического уравнения

k2+pk+q=0 ( a ¹ k1 , a ¹ k2),

то решение имеет вид y*=Qn(x)exp(ax).  

Неизвестные множители многочлена Qn(x) определяются следующим образом.

Выписываем выражение для у* (например, у*= (ax+b)e2x). Находим формулы для y*' и y*''.Подставляем выражения для у*, y*' и y*'' в уравнение (1), делим обе части получившегося уравнения на eαx , приводим все подобные члены в правой и левой частях уравнения. Затем приравниваем множители при одинаковой степени x (x0,x1,x2,...,xn) в левой и правой части полученного тождества. Таким образом, мы получаем n+1 уравнений для определения n+1 неизвестных коэффициентов в многочлене Qn(x).

2) Если a равно одному из неравных (при D>0) корней характеристического уравнения ( a = k1 , a ¹ k2), то решение имеет вид y*=x*Qn(x)exp(ax). Неизвестные множители многочлена Qn(x) определяются аналогично - (вводится дополнительный множитель х).

3) Если α является двукратным корнем характеристического уравнения (α =k1=k2), то решение имеет вид y*=x2*Qn(x)exp(α x) - (вводится дополнительный множитель х2).

Неизвестные множители многочлена Qn(x) определяются аналогично.

4) Рассмотренный выше случай f(x) = Pn(x) можно рассматривать как частный случай f(x)=Pn(x)*exp(a*x) при α=0  (т.к. exp(0)=1).

Мы можем полагать f(x)= Pn(x)*1= Pn(x)*exp(0*x), т.е. α=0 и , таким образом, сравнивать k1,k2  с  нулем .

Пример.

В приведенном ниже примере 1 правая часть f(x) = x. Представляем ее в виде: x=x*1=x*e0*x. Таким образом, aв выражении Pn(x)*exp(ax) равно 0.Это значение aсравнивается с корнями k1,k2.

Пример. 

Пример.    

Решение: в правой части уравнения нет множителя exp(ax)(формально, он есть : exp(0*x)=1 с a = 0 ). Итак a =0 и (a ¹ k1 , a ¹ k2).

 

 

II.  Правая часть уравнения содержит sin и cos

Пусть f(x) = P(x)*exp(ax)cos(bx) + Q(x)*exp(ax)sin(bx), где P, Q -

многочлены. Тогда:

 

1) Если комплексное число α + i*β не является корнем характеристического

уравнения( α + i*β ≠ k1 , α + i*β ≠ k2) , то частное решение следует искать в виде

y*= U(x)*exp(α x)cos(βx) + V(x)*exp(α x)sin(βx)

где U и V - многочлены, степень которых равна наивысшей степени многочленов P(x) и Q(x).

 

2) Если комплексное число α + iβ   является корнем характеристического

уравнения ( α + iβ = k1) , то частное решение следует искать в виде

y*= x*U(x)*exp(α x)cos(βx) + x*V(x)*exp(α x)sin(βx)

Следует отметить, что даже если в правой части (1) присутствует только один член (например, нет слагаемого с sin(x)), то частное решение y* следует искать в указанном выше общем виде.

 

3) Частный случай.

 

 

Примеры для частного случая 3).

y''+y= sin(x)   

k2+1=0 k1=i, k2= -i

В правой части

 Примеры для общего случая 1),2):

 

1) y'' - y = 3e2xcos(x)

Степень многочлена в правой части равна нулю. k1=1, k2=-1.

Общее решение однородного уравнения: С1*ex + C2*e-x       

Т.к. α + i*β =2+i*1 не является корнем характеристического уравнения, то

: y* = e2x(A*cos(x) + Bsin(x)).

 Подставив это выражение в уравнение получим:    (2A+4B)e2xcos(x) + (-4A+2B)e2xsin(x) = 3e2xcos(x)

Приравнивая коэффициенты при cos(x) (и sin(x)) в правой и левой частях,

получим: 22A+4B=3; -4A+2B=0. Получаем A=3/4, B=3/5.

Следовательно: y* = e2x*[3/10 cos(x) + 3/5sin(x)]

Общее решение: С1*ex + C2*e-x + e2x*[3/10 cos(x) + 3/5sin(x)]

 

2)  

 

где 1-е слагаемое - решение однородного уравнения, а 2-е y* - решение неоднородного уравнения.

Упражнения.

 

 

Задачи с решениями          

1)                  

y(x) =

2)

Частное решение:

 

3)

Частное решение:

 

 

Приложение

Ликбез по мат. анализу

I       Функция

Определение: Переменная y называется функцией от переменной x (y=f(x)), если по некоторому правилу каждому значению х ставится в соответствие одно определенное значение y. (Читается: "y равен эф от икс"). Таким образом, правило f числу х сопоставляет число y. Число х называется значением аргумента функции, y –значением функции.

Переменная х называется независимой переменной, y - зависимой.

Способы задания функции: аналитический (формулой, например, y=x2), графиком, таблицей.

 

Предел функции

Это понятие лежит в основе математического анализа. Дадим пояснение.

Если значения переменной y приближаются к некоторому числу b (т.е., y - b стремится к 0) когда значения аргумента х приближаются к числу а, то число b называют пределом функции при стремлении х к а (говорят еще - в точке а).

 

 Принято писать: lim f(x) = b

                  при x--> a   

 

Пример: 1) Если y = x2, то: lim y = a2 при х → a        

Действительно: у -а2 = х22 = (х+а)*(х-а) --> 0, т.к. х-а --> 0 и х+а - ограничено.

 

       2) lim(x2-1)/(x-1)=lim(x-1)(x+1)/(x-1)=lim(x+1)= lim(x)+ lim(1)=1+1=2 при х -->1.

 

 

 

 

2) Степенная функция у = ха

 

 

                                                          

 

 

 

 

      

 

 

Логарифмическая функция y=loga(x) является обратной по отношению к показательной y=ax ( alog(x) =x; loga(ax)=x ).

Логарифм loge(x) по основанию a=e=2,71... называется натуральным логарифмом и обозначается через ln(x).

Свойства логарифма: loga(b) + loga(c)=loga(b*c); loga(b) - loga(c)= loga(b/c);

Важное свойство: eln(b)=b

Свойства корней и степеней:    √x = x1/2;   1/x=x-1;   1/x2=x-2;   1/√x= x-1/2;

                                                       √(x2)=| x|   (где | x| - модуль числа х, |-3|=3).

                                                       ab*ac = ab+c;    ab/ac = ab - c;

 

II Производная


Дата добавления: 2018-06-01; просмотров: 692; Мы поможем в написании вашей работы!

Поделиться с друзьями:






Мы поможем в написании ваших работ!