ЭВМ_Семестр4_МетодПособие
.pdf
Пример 8.2. Используя программу на VBA, найти приближенное решение задачи Коши по методу Рунге-Кутта и Эйлера модифицированного с шагом h=0.1
Sub main()
Dim t As Single, ye As Single, yr As Single, yt As Single t = 0
yr = 4: ye = 4
Call Runge_1(3, 30, t, yr) Call eiler_1(3, 30, t, ye)
t = 3
yt = 2 * t - 2 + 6 * Exp(-t)
Debug.Print t, ye, yr, yt
Debug.PrintAbs(yr - yt), Abs(ye - yt) Debug.PrintAbs(yr - yt) / yt, Abs(ye - yt) / yt
End Sub
Function Rigth(t As Single, y As Single) Rigth = 2 * t - y
End Function
Sub eiler _1(h As Single, m, t As Single, y As Single)
Dim x As Single x = t
dx = h / m
For k = 1 To m k1 = Rigth(x, y)
y = y + dx * Rigth(x + dx / 2, y + dx / 2 * k1) x = x + dx
Next
End Sub
71
Sub Runge _1(h As Single, m, t As Single, y As Single)
Dim x As Single x = t
dx = h / m
For k = 1 To m k1 = Rigth(x, y)
k2 = Rigth(x + dx / 2, y + dx * k1 / 2) k3 = Rigth(x + dx / 2, y + dx * k2 / 2) k4 = Rigth(x + dx, y + dx * k3)
x = x + dx
y = y + dx * (k1 + 2 * k2 + 2 * k3 + k4) / 6
Next EndSub
Результаты расчетов в конечной точке
3 |
4,300336 |
4,298723 |
4,298722 |
t |
Эйлера |
Рунге |
аналитическое |
Задание к расчетно-графической работе №8
Для обыкновенного дифференциального уравнения (ОДУ) 1 порядка
y '(t) f (t, y(t)), t [t0,T ], y(t0) y0
исходные данные: функцию f правой части, начальное значение y0 взять из таблицы 8.1 выполнить следующие задания:
–используя функцию dsolve пакета MATLAB, найти аналитическое решение дифференциального уравнения.
–используя функцию ode23 пакета MATLAB, найти приближенное решение задачи Коши по методу Рунге-Кутта
–используя программу на VBA, найти приближенное решение задачи Коши по методу Рунге-Кутта с шагом h=0.1
72
Таблица 8.1 Варианты заданий
№ |
Функция f(t,y) |
t0 |
T |
y0 |
|
вар. |
|||||
|
|
|
|
||
|
|
|
|
|
|
1. |
y /t t2 |
1 |
2 |
0 |
2. |
|
|
yctgt 2t sint |
|
|
+1 |
0 |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
3. |
y cost |
sin(2t) |
0 |
|
1 |
0 |
|
||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
||
4. |
ytgt cos2 t |
|
|
+1 |
0.5 |
||||||||||||||||||||
4 |
4 |
||||||||||||||||||||||||
5. |
|
|
|
y |
t |
2 |
2t |
-1 |
|
0 |
1.5 |
||||||||||||||
|
t 2 |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
6. |
|
|
|
|
|
|
e |
(t 1) |
0 |
|
1 |
1 |
|
||||||||||||
|
|
t 1 |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
7. |
|
|
y / t t sin t |
|
|
+1 |
1 |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
8. |
|
|
y / t sint |
|
+1 |
1 |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
9. |
|
|
|
|
|
y |
|
t 2 |
|
|
1 |
|
2 |
1 |
|
||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
2t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
2t |
|
|
|
|
|
|
|
2t2 |
|
|
|
|
|
2 |
|
|||||
10. |
|
|
1 t2 y 1 t2 |
0 |
|
1 |
|
|
|||||||||||||||||
|
|
3 |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
11. |
|
|
|
2t 5 |
|
y 5 |
|
|
2 |
|
3 |
4 |
|
||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
t 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
12. |
|
|
y / t |
|
t 1 |
e |
t |
1 |
|
2 |
|
e |
|||||||||||||
|
|
|
|
t |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
№ |
|
|
Функция f(t,y) |
t0 |
|
T |
y0 |
||||||||||||||||||
вар. |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
13. |
|
|
y / t 2ln t / t |
1 |
|
2 |
1 |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
14. |
|
|
|
y / t 12 / t3 |
|
|
1 |
|
2 |
4 |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
73
15. |
|
2 y / t t3 |
1 |
2 |
- |
5 |
|
||||||||||||||||
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
16. |
|
|
y /t 3t |
|
|
|
1 |
2 |
1 |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
17. |
|
|
|
2ty |
1 t2 |
|
|
1 |
2 |
3 |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
1 t2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
18. |
|
|
|
2t 1 |
y 1 |
|
|
|
1 |
2 |
1 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
t2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
19. |
|
|
|
3y |
|
2 |
|
|
|
|
1 |
2 |
1 |
|
|||||||||
|
|
|
|
t3 |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
||||||
20. |
|
2ty 2t3 |
|
1 |
2 |
e 1 |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
21. |
|
y / t 2 / t |
2 |
1 |
1 |
1 |
|
||||||||||||||||
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
22. |
|
|
|
ty t3 |
|
|
|
0 |
1 |
3 |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
23. |
|
2 |
y et (t 1)2 |
0 |
1 |
1 |
|
||||||||||||||||
|
|
|
|||||||||||||||||||||
|
|
t 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
24. |
2ty te |
t2 |
sint |
0 |
1 |
1 |
|
||||||||||||||||
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 y |
|
|
|
|
|
|
|
3 |
|
|
|
|
|
||||||||
25. |
|
|
(t |
1) |
|
0 |
1 |
0.5 |
|||||||||||||||
|
t 1 |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
26. |
|
y cost sin 2t |
0 |
1 |
3 |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
27. |
4ty |
4t |
3 |
0 |
1 |
-0.5 |
|||||||||||||||||
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
28. |
|
y / t ln t / t |
1 |
2 |
1 |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
29. |
3t 2 y t 2 (1 t3 ) / 3 |
0 |
1 |
0 |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
30. |
y cost sin 2t |
0 |
1 |
-1 |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
74
РАСЧЕТНО-ГРАФИЧЕСКАЯ РАБОТА № 9
Численные методы безусловной оптимизации Краткие теоретические сведения
Необходимые и достаточные условия экстремума в задачах без-
условной оптимизации |
|
|
|
|
|
|
|
||||
|
Пусть |
задано |
множество |
|
и |
функция |
|||||
( ) |
|
( |
|
) определенная на множестве X. Точка |
|
|
|||||
называется точкой локального минимума функции |
( |
) на множе- |
|||||||||
стве |
, |
если |
существует |
окрестность |
или шар |
( |
) |
* |
‖ |
||
‖ |
+ такой, что для любого |
( |
) выполняется неравенство |
||||||||
|
|
|
|
( |
) |
( ) |
|
|
|
(9.1) |
|
|
Если неравенство (10.1) выполняется как строгое (при |
|
), |
||||||||
то говорят, что |
- точка строгого локального минимума. Точка |
||||||||||
|
называется точкой глобального минимума функции |
( |
) на |
||||||||
множестве , если неравенство (10.1) выполняется для любого |
из |
||||||||||
множества . |
|
|
|
|
|
|
|
|
|
||
|
Аналогично определяются точки локального и глобального |
||||||||||
максимума функции |
( ) на множестве . |
|
|
|
|
||||||
|
Точки локального минимума и максимума функции ( |
) назы- |
|||||||||
вают точками экстремума этой функции. |
|
|
|
|
|||||||
|
Задача отыскания всех локальных минимумов (максиму-мов) |
||||||||||
функции |
( ), |
если множество |
совпадает со всем -мерным про- |
||||||||
странством, т.е. |
, называется задачей безусловной оптими- |
||||||||||
зации, а функция ( |
) – целевой функцией. |
|
|
|
|
||||||
|
Задачу отыскания точек локального минимума целевой функ- |
||||||||||
ции |
( ) символически записывают так: |
|
|
|
|
|
|||||
|
|
|
( |
) |
|
|
|
|
|
(9.2) |
|
Аналогично, задачу отыскания точек локального максиму-ма функции ( ) символически записывают следующим образом:
75
( ) |
(9.3) |
Задача (3) эквивалентна задаче
( )
в том смысле, что множества локальных и глобальных решений этих задач соответственно совпадают.
Теорема (необходимые условия локального минимума).
Пусть |
– точка локального минимума функции ( ), которая имеет |
||
в этой |
точке непрерывные частные производные |
|
( |
|
|||
) тогда частные производные функции ( ) в этой точке равны нулю, т.е.
( ) |
( |
) |
Иначе говоря, в этой точке градиент функции
( |
)=. |
|
|
|
|
|
|
/ |
|
|
|
|
|
|
|
|
|
|
|
|
|||
равен нулевому |
вектору, т.е. |
( ) |
. |
|
|
||||||
Точкa |
, удовлетворяющая условию |
|
|
||||||||
|
|
|
( |
) |
|
|
|
||||
называется стационарной точкой функции ( |
). |
|
|||||||||
Не всякая стационарная точка функции |
( ) является решени- |
||||||||||
ем задачи (2). |
|
|
|
|
|
|
|
|
|
||
Пусть даны квадратная матрица |
( |
) порядка |
и вектор |
||||||||
. Умножение матрицы |
на вектор |
дает вектор |
, коор- |
||||||||
динаты которого определяются соотношениями |
|
||||||||||
= ∑ |
( |
|
|
) |
|
|
|
||||
76
|
Квадратная |
матрица |
называется |
симметричной, если |
||||
|
|
. |
|
|
|
|
|
|
|
Симметричная матрица |
называется неотрицательно опреде- |
||||||
ленной, если для любого вектора |
скалярное произведение век- |
|||||||
торов |
|
и |
неотрицательно, т.е. |
( |
) |
; положительно |
||
определенной, если ( |
) |
, |
; неположительно опреде- |
|||||
ленной, если ( |
) |
|
; отрицательно определенной, если |
|||||
( |
) |
, |
. |
|
|
|
|
|
Теорема (критерий Сильвестра). Симметричная матрица неотрицательно (положительно) определена тогда и только тогда, когда все главные (угловые) миноры неотрицательны (положительны):
= |
( |
) |
| |
| |
( )
| |
| |
( |
) |
( |
) |
Симметричная матрица является неположительно (отрицательво) определенной тогда и только тогда, когда знаки последовательных главных миноров чередуются, причем
= |
( ) |
( ) |
( ) |
|
|
Матрица вида |
|
|
( )
( |
|
|
|
|
|
) |
77
называется матрицей Гессе функции ( ). |
||
Теорема. Если точка |
– локальное решение задачи (2) и в |
|
этой точке функция |
( ) имеет непрерывные частные производные |
|
до второго порядка включительно, то матрица Гессе функции ( ) в |
||
точке является неотрицательно определенной, т.е. |
||
( ( |
) |
) |
|
Теорема (о достаточных условиях локального экстрему- |
|||
ма). |
Если точка |
является стационарной точкой функции ( ), |
||
т.е |
( ) |
., и матрица Гессе функции ( ) в тоже |
положи- |
|
тельно определена, то |
- строгое локальное решение задачи (2); ес- |
|
ли _в точка |
являемся стационарной для функции ( ) и матрица |
|
Гессе в ней отрицательно определена, то ─ строгое локальное решение задачи (3).
Пример 1. Дана функция
( )
Решением этой системы уравнений являются две точки
. / |
. / |
при |
этом |
( ) ( |
* |
( |
) |
. |
/ |
( ) |
. |
/ |
|
|
|
|
|
|
|
|
В силу критерия Сильвестра матрица |
|
( ) не является неот- |
|||||
рицательно определенной, а матрица |
( ) положительно определе- |
|||||||
на. Следовательно, на основании теоремы 3 |
точка |
не является ре- |
||||||
шением задачи; согласно теореме 4, в точке |
|
данная целевая функ- |
||||||
ция достигает строгого локального минимума. |
|
|
|
|||||
78
Пример 2. Исследовать на экстремум функцию
( ) |
|
( |
) . |
|
|
|
/ |
|
|
|
Найдем стационарные точки:
( )
( )
или
|
|
{ |
|
|
|
|
|
|
|
|
||||
Решением этой системы является точка |
. /. |
|||||||||||||
При этом |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( ) ( |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|||||||||||
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Следовательно, в точке |
имеет место строгий локальный мак- |
|
симум. |
|
|
Выпуклые множества и выпуклые функции |
||
Множество |
называется выпуклым, если вместе с лю- |
|
быми двумя точками |
и |
ему целиком принадлежит отрезок, со- |
единяющий эти точки.
Аналитически условие выпуклости множества записывается
так:
79
( |
) |
|
|
|
|
|
|
, |
- |
(9.4) |
|
Функция |
( ) |
определенная |
на |
выпуклом |
множестве |
||||||
|
, называется выпуклой, если |
|
|
, |
- |
справед- |
|||||
ливо неравенство |
|
|
|
|
|
|
|
|
|
||
|
(( |
) |
) ( |
|
) ( |
) |
( ) |
|
|
(9.5) |
|
Если неравенство (5) выполняется как строгое, то функция |
( ) |
||||||||||
называется строго выпуклой на множестве . |
|
|
|
|
|||||||
Теорема. Если функция |
( ) |
выпукла на множестве |
и |
яв- |
|||||||
ляется стационарной точкой функции |
( |
), т.е. |
( ) |
|
, то |
|
|||||
строгое локальное решение задачи |
|
|
|
|
|
|
|
||||
|
|
|
( ) |
|
|
|
|
|
|
(9.6) |
|
Теорема. |
Если функция |
( ) |
и множество |
выпуклы, то лю- |
|||||||
бое локальное решение задачи (6) является также глобальным решением на множестве . Пусть функция ( ) выпукла на множестве ; тогда из двух предыдущих теорем следует, что отыскание стационарной точки задачи (6) означает нахождение глобального решения этой задачи.
Достаточные условия выпуклости функции содержит следующая теорема.
Теорема. Если функция ( ) имеет непрерывные производные
до второго порядка включительно и матрица Гессе функции ( ) |
по- |
ложительно определена в любой точке х выпуклого множества |
, то |
функция ( ) является выпуклой на множестве . |
|
Пример 1. Дана функция |
|
( )
Показать, что стационарная точка функции ( ) является глобальным решением задачи ( ) .
80
