Математика. В 2. Ч. 1
.pdf
Пример 4.1. Найти решение уравнения ex |
|
, удовле- |
||
творяющее начальному условию y(0) 1 |
(задача Коши). Изобразите график |
|||
решения (интегральную кривую, проходящую через точку (0, 1)). |
||||
Решение. |
|
|
|
|
1. Разделите в уравнении переменные: |
ydy |
ex |
|
dx . |
ex |
|
|||
|
|
1 |
||
2.Установите режим автоматических вычислений.
3.Установите режим отображения результатов символьных вычислений по горизонтали, установив метку Horizontally в окне диалога строки Evaluation Style меню Symbolics.
4. |
Введите начальные условия y(x0 ) y0 . |
|||
5. |
Определите подынтегральные функции Y ( y) и X (x) : |
|||
x0 : 0; |
y0 : 1; |
|
||
|
|
ex |
|
|
X x : |
|
; |
Y y : y . |
|
1 ex |
||||
6. |
Решите уравнение, задающее неявно y(x) как функцию переменной x : |
|||
Given
x |
y |
|
= |
x0 |
y0 |
Find y |
. |
Здесь Given – ключевое слово.
Логический знак жирного равенства вводится с палитры Booleans. Завершает блок встроенная функция Find y , после которой следует набрать сим-
вольный знак равенства из меню Symbolic (подменю Evaluation).
Для решения этого уравнения используется символьный процессор MathCAD (аналогично применению вычислительного блока для численного решения нелинейных уравнений).
7. Решение, удовлетворяющее условию y(0) 1:
yx :
.
8.Постройте график найденного решения.
51
10 |
10 |
|
|
|
|
|
|
|
5 |
|
|
y(x) |
|
|
|
5 |
0 |
5 |
10 |
5 |
5 |
|
|
|
|
|
|
1.54 |
|
x |
10 |
Численный метод Рунге–Кутта решения задачи Коши для ОДУ первого порядка. Нахождение точного решения задачи Коши для многих типов ОДУ затруднительно, иногда невозможно, поэтому было создано множество приближенных методов. Один из самых популярных – метод Рунге–Кутта четвертого порядка, дающий приближенное решение в виде таблицы значений этого решения в отдельных точках заданного интервала для независимой переменной x . В пакете MathCAD есть встроенная функция, реализующая этот метод:
rkfixed y, x1, x2, npoints, D ,
где |
y |
– вектор начальных условий (для ОДУ первого порядка – |
это точка |
|
y0 |
y |
x1 ; |
|
|
|
x1, x2 – границы интервала, на котором ищется решение ДУ; |
|
||
|
npoints |
– число точек (не считая начальной), в которых ищется прибли- |
||
женное решение; |
|
|||
|
D x, y |
– для системы ОДУ первого порядка вектор, состоящий из первых |
||
производных неизвестных функций (для ДУ первого порядка y |
– |
|||
это |
f |
x, y ). |
|
|
|
Пример 4.2. Решите на отрезке [0, 3] задачу Коши y |
, y(0) 1 |
||
методом Рунге–Кутта с постоянным шагом в 20 40, 100
равностоящих точках отрезка 0, 3 .
Решение.
1.Установите режим автоматических вычислений.
2.Присвойте переменной ORIGIN значение, равное 1: ORIGIN : 1.
3.Введите начальное условие:
y : 1.
52
4. Определите правую часть уравнения f (x, y) : |
|
|
|
|
|||||||||
f x, y : |
|
|
в 20 |
40, 100 |
равностоящих точках отрезка 0, 3 . |
||||||||
5. Вычислите решение |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
Y1: |
rkfixed |
y,0,3,19, f |
|
|
|
|
||||
|
|
|
Y 2: |
rkfixed |
y,0,3,39, f |
|
|
|
|
||||
|
|
|
Y 3: |
rkfixed |
y,0,3,99, f |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
|
|
|
1 |
2 |
|
|
1 |
2 |
|
|
1 |
0 |
1 |
|
|
1 |
0 |
1 |
|
1 |
0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
0.15 |
1.011 |
|
2 |
0.075 |
1.003 |
|
2 |
0.03 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
0.3 |
1.046 |
|
3 |
0.15 |
1.011 |
|
3 |
0.06 |
1.002 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
0.45 |
1.104 |
|
4 |
0.225 |
1.026 |
|
4 |
0.09 |
1.004 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
0.6 |
1.189 |
|
5 |
0.3 |
1.046 |
|
5 |
0.12 |
1.007 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
0.75 |
1.301 |
|
6 |
0.375 |
1.072 |
|
6 |
0.15 |
1.011 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
7 |
0.9 |
1.436 |
|
7 |
0.45 |
1.104 |
|
7 |
0.18 |
1.016 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y1 |
8 |
1.05 |
1.584 |
Y2 |
8 |
0.525 |
1.143 |
Y3 |
8 |
0.21 |
1.022 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
1.2 |
1.727 |
|
9 |
0.6 |
1.189 |
|
9 |
0.24 |
1.029 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
1.35 |
1.84 |
|
10 |
0.675 |
1.242 |
|
10 |
0.27 |
1.037 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
11 |
1.5 |
1.907 |
|
11 |
0.75 |
1.301 |
|
11 |
0.3 |
1.046 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12 |
1.65 |
1.924 |
|
12 |
0.825 |
1.366 |
|
12 |
0.33 |
1.055 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
13 |
1.8 |
1.9 |
|
13 |
0.9 |
1.436 |
|
13 |
0.36 |
1.066 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
14 |
1.95 |
1.846 |
|
14 |
0.975 |
1.509 |
|
14 |
0.39 |
1.078 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
15 |
2.1 |
1.771 |
|
15 |
1.05 |
1.584 |
|
15 |
0.42 |
1.091 |
|
|
16 |
2.25 |
1.684 |
16 |
1.125 |
1.658 |
|
16 |
0.45 |
1.104 |
|||
6. Постройте на одном графике найденные решения. |
|
|
||||||||||
|
|
1.924 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y2 |
2 |
1.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Y3 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
2 |
|
3 |
|
|
|
|
|
0 |
|
Y1 |
1 |
Y2 |
1 |
Y3 |
1 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|||
Дифференциальные уравнения второго порядка. ОДУ второго порядка |
||||||||||||
имеет вид |
|
|
|
|
|
|
|
|
|
|
|
|
F x, y, y |
или y |
. |
53
Найти |
такое |
решение, которое |
удовлетворяет начальным условиям |
y x0 y0 , |
y |
. |
|
Решение ДУ второго порядка методом Рунге–Кутта. С помощью заме- |
|||
ны переменных y |
x |
вместо одного ДУ второго порядка |
|
y
f x, y, y
получим систему двух ДУ первого порядка:
|
1,Y2 |
|
с вектором начальных условий |
|
|
Y 0 |
. |
|
Дальше применяем метод Рунге–Кутта с помощью встроенной функции |
||
rkfixed Y 0, x1, x2,npoints,D , где |
|
|
D x,Y |
. |
|
Пример 4.3. Решите задачу Коши методом Рунге–Кутта и постройте гра- |
||
фик приближенного решения ДУ y |
, y 0 1, y |
по 60 точ- |
кам отрезка 0, 6 . |
|
|
С помощью замены y x Y1 x , y |
x вместо заданного ДУ полу- |
|
чим систему ДУ |
|
|
с вектором начальных условий Y 0
.
Решение.
1. ORIGIN : 1.
2. Y 0 : |
. |
54
3. |
D x,Y : |
. |
4. |
Z: rkfixed |
Y 0, 0, 6, 59, D , |
где Z
– решение в форме матрицы, I-й столбец которой состоит из значений x , II-й – значений y , III-й – y
.
5. Построить график решения.
|
|
|
|
|
1.109 |
1.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Z |
2 |
0.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-0.12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-0.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
2 |
4 |
6 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
Z 1 |
|
|
|
|
|
|
Линейные однородные дифференциальные уравнения второго порядка |
|||||||||||||||
с |
постоянными |
коэффициентами. |
Структура |
общего решения |
ДУ |
|||||||||||
y |
a1 y |
a2 y |
|
0 |
зависит от характера корней соответствующего характери- |
|||||||||||
стического уравнения |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
Если характеристическое уравнение имеет два различных действительных |
|||||||||||||||
корня |
и |
, |
то фундаментальная система решений имеет вид |
y1(x) |
e |
|||||||||||
и y2(x) |
e . Если характеристическое уравнение имеет два равных действи- |
|||||||||||||||
тельных |
корня |
|
|
, |
то |
фундаментальная система |
решений |
имеет |
вид |
|||||||
y1(x) e |
и |
|
y2(x) xe |
. |
Если |
характеристическое |
уравнение |
имеет |
два |
|||||||
комплексных корня |
и |
, то фундаментальная система ре- |
|
шений имеет вид y1(x) |
и y2(x) |
. |
|
Общее решение ДУ |
|
|
|
y |
x |
. |
|
Пример 4.4. Найдите общее решение уравнения y |
. Решите |
||
задачу Коши с начальными условиями y(0) 1, y |
, решенную прибли- |
||
55
женно методом Рунге–Кутта в примере 4.3. Проверьте правильность решения. Изобразите его график.
Решение.
1.Установите режим автоматических вычислений.
2.ORIGIN: 1.
3.Найдите корни характеристического уравнения.
|
Given |
|
|
|
|
|
|
л2 |
2 л |
3 |
|
|
|
Find |
л |
|
|
4. |
y1 x : |
, |
y2 |
x : |
– функции фундамен- |
тальной системы решений. |
|
|
|
||
5. |
Запишите общее решение уравнения (как функцию переменных x , c1 |
||||
и c2 ). |
|
|
|
|
|
6. |
Для определения значений c1 и c2 , при которых выполняются началь- |
||||
ные условия, найдите y 0 |
и y |
0 : |
|
|
|
y 0,c1,c2
dy1 x,c1,c2 : ddx y x,c1,c2
dy1 0,c1,c2
Упростим это выражение
dy1 0,c1,c2 simplify |
. |
7. Используйте вычислительный блок для нахождения c1 и c2 , учитывая,
что 
Given
c1= 1
c1 c2 
2 = 1
Find c1,c2
1 
2 
56
Таким образом, искомое решение:
yc x :
yc x
8. Проверьте решение подстановкой в уравнение.
dyc1 x : |
d |
d2 |
dx |
dx |
|
dyc2 x |
|
|
9. Проверьте выполнение начальных условий
yc 0
1
dyc1 0 1
10. Постройте график решения:
y x,c1,c2 :
y x,c1,c2
50
50
yc(x)
0
25
10 |
0 |
10 |
10 |
x |
10 |
11. Сравните графики точного и приближенного решений.
Линейные неоднородные ДУ второго порядка с постоянными коэффициентами и правой частью специального вида. Общее решение линейного неоднородного уравнения записывается как сумма общего решения однородного уравнения и любого частного решения неоднородного уравнения.
Вид частного решения устанавливается по виду правой части уравнения
y
,
57
при этом |
|
|
|
|
|
|
1) для |
f |
x |
|
|
, |
|
если |
|
|
|
, то частное решение ДУ имеет вид |
||
|
|
|
|
ycn x,b2,b1,b0 |
; |
|
если |
|
|
или |
, то частное решение ДУ |
||
|
|
|
|
ycn x,b2,b1,b0 |
; |
|
если |
|
|
|
то частное решение ДУ |
||
|
|
|
|
ycn x,b2,b1,b0 |
; |
|
2) для |
f |
x |
Ae , |
|
|
|
если |
|
|
|
, то частное решение ДУ имеет вид |
||
|
|
|
|
ycn |
x, B |
; |
если |
|
|
или |
, то частное решение ДУ |
||
|
|
|
|
ycn |
x, B |
; |
если |
|
|
|
, то частное решение ДУ |
||
|
|
|
|
ycn x, B |
; |
|
3) для |
f |
x |
|
|
|
(в частности, при A1 0 или B1 0 ), |
если |
|
|
|
, то частное решение ДУ имеет вид |
||
|
|
|
ycn x, A2, B2 |
|
; |
|
если |
|
|
|
или |
, то частное решение ДУ имеет вид |
|
|
|
|
ycn x, A2, B2 |
|
; |
|
Пример 4.5. |
Найдите общее |
решение неоднородного уравнения |
||||
y |
|
|
. Проверьте правильность решения. |
|||
58
Решение.
1. ORIGIN: 1.
2–3. Найдите общее решение соответствующего однородного уравнения вашего задания (см. пункты 2–5 решения задачи 4.4).
4.Запишите выражение для частного решения как функцию переменной х
инеизвестных коэффициентов по виду правой части неоднородного уравнения:
ycn x, A2, A1, A0 : |
. |
5. Подставьте выражение частного решения в левую часть уравнения:
dycn1 x, A2, A1, A0 |
: |
|
d |
ycn x, A2, A1, A0 |
|
|
|
||||
|
|
dx |
|||
|
|
|
d2 |
||
dycn2 x, A2, A1, A0 |
: |
|
|
ycn x, A2, A1, A0 |
|
|
dx2 |
||||
dycn2 x, A2, A1, A0 |
|
|
|
|
|
.
6. В полученном выражении приведите подобные относительно степеней х, для чего выделите переменную х и щелкните по строке Collect в меню
Symbolics:
2 
3 
7. Приравняв коэффициенты при степенях x полученного выражения левой части уравнения и выражения правой части, запишите и решите систему относительно параметров a2 , a1, a0 :
3
Given
3 |
a2 =1 |
3 |
=0 |
2 |
=1 |
|
3 |
|
Find a2, a1, a0 |
4 |
|
9 |
||
|
||
|
11 |
|
|
|
59
8. |
Запишите частное решение с найденными коэффициентами a2 , a1, a0 : |
||||
|
ycn x : 1 |
4 |
11 . |
|
|
|
3 |
9 |
27 |
|
|
9. |
Запишите общее решение неоднородного уравнения: |
|
|
||
ycno x,c1,c2 : |
|
1 |
4 |
11 . |
|
|
|
|
3 |
9 |
27 |
10. Проверьте решение подстановкой:
dycno1 x,c1,c2 : ddx ycno x,c1,c2
d2
dycno2 x,c1,c2 : dx2 ycno x,c1,c2
dycno2 x,c1,c2 |
|
. |
|
Замечание. |
В примере 4.5 правая |
часть имеет |
вид многочлена |
f x |
. Если в задании f (x) |
или f x |
, |
соответственно измените в решении пункты 4–9.
Порядок выполнения работы
Студенту рекомендуется внимательно изучить теоретический материал, проделать все примеры, в нем встречающиеся, и после этого приступать к выполнению своего варианта задания.
Содержание отчета
1.Краткий обзор по теоретической части.
2.Файл MathCAD с выполненными заданиями своего варианта.
Варианты заданий |
|
Вариант 1 |
|
1. Найти решение уравнения с разделяющимися переменными ex |
3ydy xdx , |
удовлетворяющее начальному условию y(0) 1 (задача Коши). |
Изобразить |
график решения (интегральную кривую, проходящую через точку (0,1)).
60
