русс | укр

Языки программирования

ПаскальСиАссемблерJavaMatlabPhpHtmlJavaScriptCSSC#DelphiТурбо Пролог

Компьютерные сетиСистемное программное обеспечениеИнформационные технологииПрограммирование

Все о программировании


Linux Unix Алгоритмические языки Аналоговые и гибридные вычислительные устройства Архитектура микроконтроллеров Введение в разработку распределенных информационных систем Введение в численные методы Дискретная математика Информационное обслуживание пользователей Информация и моделирование в управлении производством Компьютерная графика Математическое и компьютерное моделирование Моделирование Нейрокомпьютеры Проектирование программ диагностики компьютерных систем и сетей Проектирование системных программ Системы счисления Теория статистики Теория оптимизации Уроки AutoCAD 3D Уроки базы данных Access Уроки Orcad Цифровые автоматы Шпаргалки по компьютеру Шпаргалки по программированию Экспертные системы Элементы теории информации

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


Дата добавления: 2015-01-16; просмотров: 2171; Нарушение авторских прав


К этому классу относятся ДУ вида

y// +py/+qy = f(x), (13)

где p,q – вещественные числа. Это же уравнение с нулевой правой частью, т.е. когда f(x)=0, называется линейным однородным ДУ. Вначале решают линейное однородное ДУ, а затем находят общее решение линейного неоднородного ДУ.

 

а) Линейные однородные уравнения дифференциальные уравнения с

постоянными коэффициентами.

y// +py/+qy = 0 . (14)

Они попадают в класс ДУ из раздела 3.3. Однако для них имеется более простой метод решения.

Пусть y1(x) и y2(x) - частные решения этого уравнения.

Определение. Два решения называются линейно зависимыми, если

они пропорциональны друг другу: y2(x) = α · y1(x)или y2(x) = b · y1(x) . Понятия линейной зависимости и линейной независимости применимы к любой паре функций. Например, и при линейно независимы, так как .Две функции ekx и z(x)ekx также линейно независимы, если z(x) – многочлен степени ≥ 1.

Для двух линейно независимых частных решений y1и y2ДУ(14) и произвольных постоянных C1и C2 линейная комбинация

y=C1 y1+ C2 y2 (15)

является решением ДУ (14). Эта формула в действительности даёт все решения линейного однородного ДУ (14),так как число (два) произвольных постоянных совпадает с порядком ДУ, и разным значениям произвольных постоянных соответствуют разные решения. Если решение дифференциального уравнения есть комплексная функция y(x)=u(x)+iv(x), то отдельно действительная часть u(x) и мнимая часть v(x) являются решениями.

Будем искать решение уравнения (14) в виде y(x)= ekx, где k - пока неизвестное число. Тогда y/=kek x, y//= k2ekx ,и уравнение примет вид:

k 2 ekx+ pkekx+qekx = 0.

Поделив обе части уравнения на ekx, получим алгебраическое, а именно, квадратное уравнение относительно неизвестного k:

k2+pk+q = 0

(это уравнение называется характеристическим), решение которого есть



.

Числа k1 и k2 называются корнями характеристического уравнения.

Рассмотрим три возможности ( >0, =0, <0) для подкоренного выражения. (Напомним, что p2 - 4 q называется дискриминантом характеристического уравнения.)

1. >0, тогда , и - два линейно независимые решения уравнения (14). Общее решение имеет вид

y(x)=C1 +C2 . (16)

2. =0, тогда . Формула y(x)= ekx дает только одно из двух линейно независимых решений: y1(x)= ekx. Второе решение y2(x) ищем в виде y2(x)= z(xekx. Подставляя это в уравнение (14), получаем : z//=0, одним из решений которого является z = x. В результате получаем пару линейно независимых решений y1(x)= ekx, y2(x) =xekx иобщее решение ДУ (14):

y(x)= ekx + ekx. (17)

3. <0, k1=α +βi, k2=α – βi, где . Общее решение

y(x)= e( α + βi)x+ e( α - βi)x = eα x ( e βix+ e- βix).

Согласно формулам Эйлера e βix=cos βx+isin βx, e- βix= cos βx-isin βx.

Как указывалось выше, в случае комплексного решения отдельно eα x cos βx и eα x sin βx являются решениями. Следовательно, общее решение (вещественное) ДУ (14) имеет вид

y(x)= eα x1cos βx+С2 sin βx). (18)

Пример 11. y// -2y/-8y = 0.

Решение: k2-2k-8=0, k1=4, k2=-2. Согласно формуле (16) пишем

Ответ: y=C1e4x+C2e-2x.

Пример 12. y// - 6y/ + 9y = 0.

Решение: k0-6k+9=0, k1= k2=k =3. Согласно формуле (16)

Ответ: y = C1e3x+C2xe3x.

Пример 13. y// - 6y /+ 13y = 0

Решение: k2-6k+13=0, k1=3-2i, k2=3+2i. Согласно формуле (18) имеем

Ответ: y = e3x(C1 cos2x+C2 sin2x).

б) Линейные неоднородные дифференциальные уравнения с постоянными коэффициентами.

Общее решение неоднородного линейного ДУ (13) равно сумме

y=C1y1+C2y2+y*

общего решения C1y1+C2y2 соответствующего однородного уравнения (14) и частного решения y* данного неоднородного уравнения (13).Частное решение y* определяется видом правой части f(x) в (13),и его можно явно определить для функций f(x) специального вида – произведений функции еlх и многочленов Pn(x) некоторой степени ≤ n . (При комплексном l берется вещественная часть такой такого произведения; см. ниже).

При вещественном l ,когда

f(x) = Pn(xelx ,

частное решение y* можно взять в виде

y* = xm·Qn(xelx,

где показатель m означает, со сколькими корнями k1 и k2 характеристического уравнения совпадает l ( показатель m может быть 0, или 1, или 2).

Если l=λ+μi- комплексное число, то правая часть f(x) cпециального вида есть вещественная часть произведения еlх на многочлен, т.е.

f(x) = eλx · (Pn(x)·cos μx+ n(x)·sin μx).

Тогда частное решение ищется в виде

y* = xm ·eλx ·( Qn(x)·cos μx + n(x)·sin μx),

где m=0. если λ+μi не совпадает ни с одним из корней k1 и k2 ,

и m=1, если λ+μi совпадает с одним из корней k1 и k2 .

Перечислим типичные случаи.

(а) f(x)=Pn(x) (Pn(x) - многочлен n -й степени), тогда y* ищем в виде:

y*=Qn(x), если k1≠ k 2, k10, k 2 0 ;

y*=xQn(x), если k1=0, k 2 0;

y*=x2Qn(x), если k1= k 2=0.

(b) f(x)=aelx с вещественным l, тогда y* ищем в виде:

y*=Aelx, если l k1, l≠k2;

y*=Axelx, если k1≠ k 2, l=k1либо l=k2;

y*=Ax2elx, если l=k1=k2.

(c) f(x)= eλx(a ·sinµ+b·cosµx) с комплексным l= λ+iμ, тогда y* ищем в виде:

y*= eλx· (A ·sin µx+B ·cos µx),если λ+iμ≠k 1 и λ+iμ≠k2;

y*=x·eλx (A ·sin µx+B · cosµx), если λ+iμ=k 1 или λ+iμ=k2.

Для нахождения частного решения y* применяется метод неопределённых коэффициентов : подставляем y* в уравнение (13) и находим неизвестные коэффициенты многочлена Qn(x) или неизвестные числа A,B.

Пример 14. y// -2y/-8y =2ex.

Решение: общее решение однородного уравнения есть С1e4x2e-2x, Частное ищем в виде y*=Аex. Подставляя y* в исходное уравнение, получим

A ex -2A ex -8A e x =2 ex . откуда находим А=-2/9 и y*=- 2/9 ex.

Ответ: y(x) = С1e4x2e-2x - 2/9 ex.

Пример 15. y// -2y/-8y = e4x .

Решение: общее решение однородного уравнения: С1e4x2e -2x; частное решение ищем в виде y*=Aхe4x . Подставляя в исходное уравнение y* и приравнивая коэффициенты при e4x , получим A = 1/6

и y*=1/6 e4x .

Ответ: y(x) = С1e4x2e-2x +1/6 xe4x.

Пример 16. y// -6y/+9y = e3x.

Решение: общее решение однородного уравнения есть С1e3x2 xe3x, тогда частное решение ищем в виде y*=Aх2e3x . Подставляя y* в исходное уравнение и приравнивая коэффициенты при е3х ,получим: A =1/2,

y*=1/ x2e3x

Ответ: y(x) = С1e3x2 xe3x + 1/ x2e3x.

Пример 17. y// - 6y /+ 25y =50 x.

Решение: общее решение однородного уравнения есть e3x(C1 ·cos 4x+ C2·sin 4x); частное решение ищем в виде y*=Ax+B. Подставляя это в исходное уравнение и приравнивая коэффициенты при одинаковых степенях х,получим: A=2, B=12/25 , y*=2x+12/25.

Ответ: y = e3x(C1 cos 4x+ C2sin 4x)+2x+12/25.

Пример 18. y// -6y /+ 13y =sin x.

Решение: общее решение однородного уравнения естьe3x(C1 cos 2x+

C2sin 2x).Частное решение ищем в виде y*=A sin x +B cos x. Подставляя y* в исходное и приравнивая коэффициенты при sin x и cos x, получим А=1/15, В=1/30 , y*=1/15·sin x +1/30 ·cos x.

Ответ: y = e3x(C1cos2x+ C2sin2x)+ 1/15·sin x +1/30 ·cos x.

4. Задачи для самостоятельного решения.

 

1. (1+y)dx - (1-x)dy=0;

Ответ: у=С/(х-1) - 1.

2. ;

Ответ: у=С·е sin х.

3. ;

Ответ: у= .

4.(1+x2)y/= (1+y2);

Ответ: .

5. ;

Ответ: .

6. x2 y/- 2xy=3;

Ответ: у=1/х+ Сх2 .

7. xy/+y=x+1;

Ответ: у = х/2 + 1 +С/х .

8. (1+x2)y/-2xy=(1+x2)2;

Ответ: у=(1+х2)·(х+С) .

9. y//= x+2;

Ответ: у = х3/3 +х2 +С1 х +С2 .

10. xy//=2x-y/;

Ответ: у= х2/2 +С1·ln х +С2 .

11. 8y/ y//=9;

Ответ: у= (х+С1)3/2 +C2 .

12. ; начальные условия у(0)=0, у/(0)=1.

Ответ: у= tg x.

13. y//+4y=e -2x;

Ответ: у = С1 sin 2x+C2 cos 2x - 1/4 ·e -2x.

14. y// + 5y/ + 6y = 52 cos2x;

Ответ: y= C1e -2x+C2e -3x +5 sin 2x + cos 2x .

15. y// - 6y/ + 9y =9/2· x2+21 x +1;

Ответ: y=C1e 3x + C2 x·e3x +x2/2 +3x +2.

 



<== предыдущая лекция | следующая лекция ==>
Предназначено для студентов всех специальностей. | Применение дифференциала в приближенных вычислениях


Карта сайта Карта сайта укр


Уроки php mysql Программирование

Онлайн система счисления Калькулятор онлайн обычный Инженерный калькулятор онлайн Замена русских букв на английские для вебмастеров Замена русских букв на английские

Аппаратное и программное обеспечение Графика и компьютерная сфера Интегрированная геоинформационная система Интернет Компьютер Комплектующие компьютера Лекции Методы и средства измерений неэлектрических величин Обслуживание компьютерных и периферийных устройств Операционные системы Параллельное программирование Проектирование электронных средств Периферийные устройства Полезные ресурсы для программистов Программы для программистов Статьи для программистов Cтруктура и организация данных


 


Не нашли то, что искали? Google вам в помощь!

 
 

© life-prog.ru При использовании материалов прямая ссылка на сайт обязательна.

Генерация страницы за: 0.009 сек.