Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Презентации по математике / Презентации / 18. Диф.ур.!-го порядка.ppt
Скачиваний:
48
Добавлен:
11.02.2015
Размер:
378.37 Кб
Скачать

План лекции.

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

 

 

 

 

 

 

 

Нетрудно видеть, что ошибки выбора

 

 

 

формулы достаточно малы (могут быть

 

 

 

порядка ошибок измерения).