- •Содержание
- •Вводная лабораторная работа
- •Раздел 1 Лабораторная работа №1
- •Лабораторная работа №3
- •Лабораторная работа №4
- •Раздел 2 Лабораторная работа №1
- •Задания для лабораторной работы №1
- •Лабораторная работа №2
- •Задания для лабораторной работы №2
- •Лабораторная работа №3
- •Задания для лабораторной работы №3
- •Лабораторная работа №4
- •Задания для лабораторной работы №4
- •Лабораторная работа №5
- •Задания для лабораторной работы №5
- •Лабораторная работа №6
- •Задания для лабораторной работы №6
- •Литература
Лабораторная работа №3
Дифференциальные уравнения
Задача 1. Найти решение уравнения с разделенными переменными ydy=(exp(x)/1+exp(x))dx, удовлетворяющее начальному условию y(0)=1 (задача Коши). Изобразите график решения (интегральную кривую, проходящую через точку (0,1)).
Решение:
1) Установите режим автоматических вычислений.
2) Установите режим отображения результатов символьных вычислений по горизонтали, установив метку Horizontaly в окне диалога строки Evaluation Style меню Symbolics.
3) Введите начальные условия y(x0)=y0:
4) Если уравнение имеет вид Y(y)dy=X(x)dx, определите подынтегральные функции Y(y) и X(x):
5) Вычислите символьно интегралы с переменными верхними пределами и нижними пределами, равными начальным условиям x0 и y0:
6) Запишите уравнение , задающее неявно y(x) как функцию x, и решите его символьно относительно переменной y:
7) Выбираете решение, удовлетворяющее условию y(0)=1, и определите как функцию переменной x:
8) Постройте график найденного решения:
Задача 2. Решите на отрезке [0,3] задачу Коши y'=sin(xy), y(0)=1 методом Рунге-Кутты с постоянным шагом на сетке из 20 (40, 100) равнооотстоящих узлов.
Решение:
1) Установите режим автоматических вычислений.
2) Присвойте переменной ORIGIN значение, равное 1:
3) Присвойте начальное значение решения переменной y1:
4) Определите правую часть уравнения f(x,y):
5) Вычислите решение, используя функцию rkfixed(y,x1,x2,npoints,f), где у - вектор начальных условий, х1 и х2 - концы отрезка интегрирования, npoints- число узлов на отрезке интегрирования, f - правая часть уравнения. В результате получите матрицу размерности (npoints, 2), в первом столбце которой содержатся значения х, во втором - значения у.
6) Постройте на одном графике найденные решения:
Задача 3. Решите на отрезке [0,3] задачу Коши у"=ехр(-ху), у(0)=1, у'(0)=1 методом Рунге-Кутты с постоянным шагом на сетке из 20 (40 , 100) равноотстоящих узлов.
Решение:
Сведите решение задачи для уравнения к задаче для системы. Обозначьте у1(х)=у(х) и у2(х)=у'(х). Поскольку у"(х)=(у'(x))'=y2'(x), то получим
у1'=y2 y1(0)=1
y2'=exp(-xy1) y2(0)=1
1) Установите режим автоматических вычислений.
2) Присвойте переменной ORIGIN значение, равное 1:
3) Присвойте начальное значение решения вектору-столбцу с именем у:
4) Определите правую часть уравнения, присвойте соответствующие выражения элементам вектора-столбца с именем f(x,y):
5) Вычислите решение, используя функцию rkfixed:
6) Постройте на одном графике найденные решения:
Задача 4. Найдите общее решение однородного уравнения y''+2y'+3y=0. Решите задачу Коши с начальными условиями у(0)=1, y'(0)=1. Проверьте правильность решения. Изобразите его график.
Решение:
1) Установите режим автоматических вычислений.
2) Присвойте переменной ORIGIN значение, равное 1:
3) Запишите характеристический многочлен уравнения и найдите его корни:
4) Если характеристическое уравнение имеет 2 различных действительных корня 1 и 2, то фундаментальная система решений имеет вид и . Если характеристическое уравнение имеет 2 равных действительных корня 1= 2= , то фундаментальная система решений имеет вид и . Если характеристическое уравнение имеет 2 комплексных корня 1= +i и 2= -i , то фундаментальная система решений имеет вид и .
Запишите функции фундаментальной системы решений:
5) Запишите общее решение уравнения (как функцию переменных х, с1 и с2):
6) Найдите значения констант с1 и с2, при которых выполняются заданные начальные условия у(0)=1 и у'(0)=1:
7) Запишите решение задачи Коши:
8) Проверьте решение подстановкой в уравнение:
9) Проверьте выполнение начальных условий:
10) Постройте график решения:
11) Решите задачу Коши методом Рунге-Кутты и постройте график приближенного решения:
12) Сравните графики.
Задача 5. Найдите общее решение неоднородного уравнения y''+2y'+3y=x*x+1. Проверьте правильность решения.
Решение:
Общее решение линейного неоднородного уравнения записывается как сумма общего решения однородного уравнения и любого частного решения неоднородного уравнения.
1) Установите режим автоматических вычислений.
2) Присвойте переменной ORIGIN значение, равное 1:
3) Найдите общее решение однородного уравнения у"+2y'+3y=0:
4) Запишите выражение для частного решения как функцию переменной х и неизвестных коэффициентов. Вид частного решения устанавливается по виду правой части уравнения.
5) Подставьте выражение частного решения в левую часть уравнения:
6) В полученном выражении приведите подобные отностительно степеней х, для чего выделите переменную х и щелкните по строке Collect в меню Symbolics:
7) Приравняв коэффициенты при степенях х полученного выражения левой части уравнения и выражения правой части, запишите и решите систему относительно параметров а1, а2, а0:
8) Запишите частное решение с найденными коэффициентами а2,а1,а0:
9) Запишите общее решение:
10) Проверьте решение подстановкой: