- •План лекции.
- •1. Метод наименьших
- •от соответствующих значений
- •Согласно необходимому условию экстремума должна выполнятся следующая система :
- •Пример.
- •Для решения этой системы составим следующую расширенную таблицу и заполним пустые
- •Найденные значения коэффициентов а и b подставляем в уравнение линейной функции:
- •2. Дифференциальные уравнения.
- •Определение.
- •*Решением дифференциального уравнения называется функцияy y( x) , удовлетворяющая этому уравнению.
- •найти решениеy y(x)
- •*Общее решение уравнения (1)- это решение в
- •Подставляем y’’ и y в уравнение:
- •Подставляем найденное С в (*):
- •Для решения уравнения разделим переменные x и y
План лекции.
1.Метод наименьших квадратов.
2.Дифференциальные уравнения.
1. Метод наименьших
квадратов.
В естествознании, в частности в физических и биологических науках, основным методом исследования являются наблюдения, опыты эксперименты. В связи с этим возникает необходимость в нахождении эмпирических формул, составленных на основании опыта и наблюдения. Одним из лучших методов получения таких формул является метод наименьших квадратов, который является эффективным приложением теории экстремумов функции нескольких переменных.
|
|
|
Итак, пусть дана таблица |
|
|
|
|
|||||||
|
|
|
|
измерений в некотором |
|
|
|
|
||||||
|
|
|
|
опыте, связывающая |
|
|
|
|
|
|
|
|||
|
|
|
|
переменные величины |
|
|
|
|
||||||
x |
|
|
x1X и Yx.2 |
… … |
xn |
|
|
|
|
|||||
yii |
|
y1 |
i |
y2 |
… |
yn |
|
|
|
|
||||
|
xi |
и |
|
|
… … |
|
|
|
|
|
|
|
||
Значения |
|
y будем считать также, как декартовые |
|
|||||||||||
|
|
|
||||||||||||
координаты точек на координатной плоскости XOY . |
|
|
||||||||||||
|
|
|
|
|
|
… |
|
|
, |
|
|
|||
|
|
|
|
|
|
|
y (x) |
|
|
|||||
Требуется найти аналитическую зависимость |
|
|
|
, |
||||||||||
|
|
|
|
|||||||||||
наилучшим образом отображающую опытную |
зависимость. |
|||||||||||||
Выберем “подходящую” функцию |
|
|
|
, где |
|
|||||||||
а,b… - параметры, так, чтобы |
y ( x, a, b,...) |
кривые |
|
|
||||||||||
соответствующие |
|
|
||||||||||||
|
для различных a, b, … проходили вблизи точек |
|
||||||||||||
|
из опыта |
(xi ; yi ) |
. |
|
|
|
|
|
|
|
|
|||
|
Найдём |
такой |
единственный набор значений |
|
|
|||||||||
|
параметров, чтобы соответствующая кривая |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
(x ; y ) |
|
||
|
распола- галась ближе всех других к точкамi |
изi |
|
|||||||||||
|
опыта |
|
, |
|
|
|
|
|
|
|
|
|
т.е. чтобы ошибки выбора формулы - отклонения
от соответствующих значений |
y |
( x ) |
из |
|||
i |
i |
|||||
Формулы были наименьшими |
|
|
||||
по абсолютной |
||||||
|
величине. |
|
|
|
|
|
|
Для этого составляется сумма |
|
||||
|
|
2 |
, |
|
|
|
|
n |
|
|
|
|
|
|
S(a, b,...) ( yi (xi , a, b,...)) |
|
|
|
|
|
|
i 1 |
|
|
|
|
|
где суммируются квадраты указанных ошибок выбора формулы.
Тогда ошибки выбора будет наименьшими (по абсолютной величине), если наименьшей будет сумма S .
Следовательно, нужно решить задачу на экстремум функции S (a, b, …): найти минимум функции нескольких переменных a, b, … .
Согласно необходимому условию экстремума должна выполнятся следующая система :
S 0a
S 0 (*)
b
...
...
Решение этой системы даст те значения параметров a, b, … , при которыхy ( xфункция, a, b,...)
будет наилучшей. (Можно доказать, что необходимое условие экстремума при решении таких задач будет и достаточным).
Пример.
Дана таблица измерений.
xi yi
1 |
3 |
5 |
7 |
Найти подходящую |
||
1 |
2 |
3 |
6 |
эмпирическую |
||
|
|
|
||||
формулу |
||||||
|
|
|
|
|
y ( x) |
|
Нанесем на координатную плоскость XOY точки ( xi ; yi ) из опыта:
y 6
3
2 1
1 3 5 7 |
x |
Все точки лежат вблизи некоторой прямой.
|
Найдем наилучшую из таких прямых, т.е. найдем |
||||||||||||||||
|
|
линейную зависимость |
|
|
|
|
|
, |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
( x) ax |
b |
|
|
|
||
|
|
наиболее точно описывающую опытную |
|||||||||||||||
|
|
зависимость. |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Для такой зависимости система (*) имеет вид: |
||||||||||||||||
|
|||||||||||||||||
|
|
n |
|
n |
|
2 |
|
|
n |
|
|
|
|
|
|
|
|
|
|
yi xi a xi |
|
b xi 0 |
(**) |
|
|
|
|||||||||
|
|
i 1 |
|
i 1 |
|
|
|
i 1 |
|
|
|
|
|
||||
|
|
n |
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
yi a xi |
bn 0 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|||||||||
|
|
i 1 |
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В нашем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
случае n=4 и система(**) |
|
|
|
|
|||||||||||||
|
|
|
|
4 |
y x a |
4 |
x2 |
b |
4 |
x 0 |
|
|
|||||
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
i |
|
i |
|
i |
|
|
|||||
|
|
перепишется :i |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
i 1 |
|
|
|
|
i 1 |
|
|
i 1 |
|
|
(***) |
||
|
|
|
|
4 |
|
|
|
4 |
|
|
|
|
|
|
|
||
|
|
|
|
yi |
a xi 4b 0 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
i 1 |
|
|
|
i 1 |
|
|
|
|
|
|
|
Для решения этой системы составим следующую расширенную таблицу и заполним пустые
клетки: |
4 |
4 |
|
|
xi 1 3 5 7 16 |
||
|
|||
|
i 1 |
i 1 |
xi |
1 |
3 |
5 7 16 |
yi |
1 2 3 6 12 |
|
|
||||||
yi |
1 |
2 |
3 |
6 |
12 |
4 |
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
i 1 |
|
|
|
|
|
|
|
2 |
1 |
9 |
25 49 84 |
xi2 1 9 25 49 84 |
|||||||||
xi |
i 1 |
|
|
|
|
|
|
|
|||||
xi yi |
1 |
6 |
15 42 64 |
xi yi 1 6 15 42 64 |
|||||||||
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
i 1 |
|
|
|
|
|
|
|
Найденные суммы подставляем в систему (***): |
|
|
|||||||||||
|
|
|
64 84a 16b 0 |: 4 |
|
16 21a 4b 0 |
|
|
||||||
|
|
|
|
16a 4b |
0 |: 4 |
|
|
3 4a b 0 |
|||||
|
|
|
12 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
16 21a 4(3 4a) 0 |
|
|
b |
5 |
|
|
||||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b 3 |
4a |
|
|
|
a |
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
Найденные значения коэффициентов а и b подставляем в уравнение линейной функции:
( x) 4 x |
1 |
По двум точкам строим |
|||
эту |
5 |
|
5 |
прямую на |
|
|
|
|
|
||
координатной |
|||||
|
|||||
x |
-1 |
4 |
|
плоскости , данной |
|
выше: |
3 |
|
|
||
y |
-1 |
|
|
y |
|
|
y 4 x |
1 |
|
|
|
|
|
||
|
|
|
|
5 |
5 |
|
6 |
|
|
|
-точки из опыта |
|
|
|
|
-точки для построения |
|
|
|
|
|
|
|
|
|
|
|
|
прямой |
|
3 |
|
|
|
-прямая (x) |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
1 |
|
|
|
|
-1 |
-1 |
1 |
4 5 |
7 |
x |
|
|
|
|
||
|
|
|
Нетрудно видеть, что ошибки выбора |
||
|
|
|
формулы достаточно малы (могут быть |
||
|
|
|
порядка ошибок измерения). |