Решения дифференциальных уравнений

 

Аналитическое (символьное) решение дифференциальных уравнений осуществляется с помощью команды DSolve, которая имеет три формата DSolve[eqn,y,x] - решает обыкновенное дифференциальное уравнение eqn относительно неизвестной функции y[x] независимой переменной x;

DSolve[{eqn1,eqn2,…,},{y1,y2,…},x] - решает систему обыкновенных дифференциальных уравнений {eqn1,eqn2,…} относительно неизвестных функций {y1[x],y2[x],…} независимой переменной x;

DSolve[eqn,y,{x1,x2,…}] - решает дифференциальное уравнение в частных производных eqn относительно неизвестной функции y[x1,x2,…] независимых переменных x1,x2, … .

Решение представлено в виде подстановки, константа интегрирования обозначена как C[1]. Если желательно другое обозначение для этой константы, то следует использовать опцию DSolveConstants→K:

Для использования полученного решения можно, например, присвоить подстановке некоторое имя

и определить соответствующую функцию

Проверка показывает, что цель достигнута:

Если необходимо решить дифференциальное уравнение с определенными граничными (или начальными) условиями, то их надо определить в списке решаемых уравнений, например, следующая команда интегрирует уравнение y’[x]==ay[x] с условием y[0]=5:

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

то соответствующая команда имеет вид

Численное решение дифференциальных уравнений осуществляется с помощью команды NDSolve, которая имеет три формата:

NDSolve[eqns,y[x],{x,a,b}] - находит приближенное решение обыкновенного дифференциального уравнения eqns (с начальными условиями) относительно неизвестной функции y[x] независимой переменной x в интервале [a,b];

NDSolve[{eqn1,eqn2,…},{y1,y2,…},{x,a,b}] - численно решает систему обыкновенных дифференциальных уравнений {eqn1,eqn2,…) с начальными условиями относительно неизвестных функций {y1[x],y2[x],…} независимой переменной x на интервале [a,b];

NDSolve[eqn,y,{x1,x2,…}] - численно решает дифференциальное уравнение в частных производных eqn относительно неизвестной функции y[x1,x2,…] независимых переменных x1,x2, … .

Так, дифференциальное уравнение третьего порядка с начальными условиями можно численно проинтегрировать командой

Результат интегрирования Mathematica выводит в форме некоторой интерполяционной функции, для которой можно построить график:

Индивидуальные задания:

 

1. Решить уравнение символьно и численно уравнение:

Задание Задание

 

2. Решить систему уравнений:

Задание Задание

 

3. Решить символьно и численно дифференциальные уравнения:

Задание Задание
; ;
; ;
; ;
; ;
; ;
; ;
; ;

 

4. Решить дифференциальное уравнение:

Задание Задание