Главная | Обратная связь | Поможем написать вашу работу!

Эмпирические зависимости. Метод наименьших квадратов.

ЛАБОРАТОРНАЯ РАБОТА № 1

1. Теоретические сведения

При исследовании многочисленных физических, химических, технологических и других процессов очень часто сталкиваются со следующей задачей: в итоге опыта получен ряд значений переменных и , требуется по полученным данным найти аналитическое выражение зависимости между и . Такая зависимость называется эмпирической.

Пусть заданы результаты наблюдений

…

…

(1)

и нужно получить такую эмпирическую зависимость

, (2)

где ─ параметры, чтобы значения в точках , мало отличались от опытных данных

Задача определения эмпирической зависимости состоит из двух этапов:

1) определение вида функциональной зависимости (выбор класса функций, которому должна принадлежать искомая функция ;

2) определение параметров эмпирической зависимости. Определение вида зависимости может быть произведено на основе теоретических представлений о характере изучаемой зависимости или из геометрических соображений. На плоскости строят точки , и по характеру их расположения выбирают вид функциональной зависимости. Например, расположение экспериментальных точек может навести на мысль о линейной, или квадратичной, или экспоненциальной или другой зависимости. Однако общего метода для нахождения наилучшего типа эмпирической зависимости, соответствующей опытным данным, указать нельзя.

После того как определен класс, которому должна принадлежать искомая эмпирическая зависимость, встает вопрос о нахождения конкретных значений параметров .

Одним из основных методов нахождения параметров эмпирической зависимости является метод наименьших квадратов. Этот метод не решает вопроса о выборе общего вида аналитической функции, а дает возможность при заданном типе аналитической функции подобрать наиболее вероятные значения для ее параметров.

Сущность метода наименьших квадратов состоит в том, что параметры выбираются из условия минимума суммы квадратов уклонений

(3)

Если рассматривать как непрерывную функцию её параметров , то в точке минимума её частные производные должны равняться нулю, т.е. искомые значения параметров должны удовлетворять системе уравнений

, , …, (4)

Известно, что если в качестве функции берется многочлен , то функция имеет одну точку экстремума, в которой достигается минимум. Следовательно, в этом случае искомые значения представляют собой решение системы (4). В этом случае

а систему (4) можно преобразовать к виду

(5)

Полученная система решается любым известным методом (Гаусса, Крамера, матричным) или с помощью стандартных программ на компьютере.

Отметим два частных случая наиболее часто встречающихся на практике — случаи линейной и квадратичной зависимости от . Если зависит от линейно , то система (5) примет вид

(6)

Чтобы получить систему (6), нужно в системе (5) положить .

В случае квадратичной зависимости система (5) преобразуется к виду

(7)

Для этого в системе (5) нужно положить .

Метод наименьших квадратов обладает тем преимуществом, что если сумма квадратов уклонений мала, то сами эти уклонения также малы по абсолютной величине. Кроме того, этот метод допускает веское обоснование с вероятностной точки зрения. Недостатком метода наименьших квадратов является громоздкость вычислений. Поэтому к нему прибегают обычно при обработке наблюдений высокой точности, когда нужно получить также весьма точные значения параметров. Заметим, что в этом случае промежуточные вычисления нужно проводить с достаточным количеством десятичных знаков, так как в противном случае искомые коэффициенты могут иметь мало верных знаков.

2. Примеры решения задач.

Задание 1. Найти эмпирическую зависимость времени валки дерева от его диаметра, по результатам испытаний, приведенных в таблице

d, cм

t, c

Решение.

Обозначим и . На плоскости построим точки ,

Рис. 1. Точечная диаграмма

Из точечной диаграммы (рис.1) видно, что точки расположены в близи некоторой прямой, поэтому можно считать, что зависимость от будет линейной, т.е. имеет вид . Для вычисления коэффициентов системы (6) составим таблицу:

Таблица 1

x_i y_i x_i² x_iy_i

∑

Система (6) для нашего примера имеет вид:

Из нее находим коэффициенты и : , .

Искомая эмпирическая функция будет . На рис.1 построим график полученной прямой. Точечная диаграмма и график эмпирической зависимости подтверждает их соответствие.

Решение при помощи Excel.

Нахождение эмпирических формул методом наименьших квадратов можно осуществлять в среде Excel, используя:

· функции СТЕПЕНЬ, ПРОИЗВЕД, СУММ, СУММКВ, СУММПРОИЗВ,

· формулы массива МОБР, МУМНОЖ,

· графический анализ с помощью Линий тренда.

Рассмотрим на конкретных примерах, как работают эти функции.

1. Внесите данные задачи в ячейки B4:B13 и C4:C13.

2. Вычислите , введя в ячейку D4 формулу

=B4*B4

и затем скопировав ее в ячейки D5:D13 (направив указатель мыши на черный квадрат в правом нижнем углу ячейки, чтобы он принял вид черного крестика, и протащив в ячейку D5 и т.д.).

Аналогично найдите произведения и результат запишите в ячейках E4:E13.

Эти же результаты можно получить, используя функции СТЕПЕНЬ и ПРОИЗВЕД в категории Математические Мастера функций.

3. Найдите сумму значений в ячейках B4:B13 и результат запишите в ячейку B14. Для этого выделим соответствующие ячейки и нажмем пиктограмму суммирование на Панели инструментов. Для суммирования элементов массива можно воспользоваться функциями СУММ, СУММКВ, СУММПРОИЗВ в категории Математические Мастера функций.

3. Система (6) для нашего примера имеет вид:

Обозначим − матрица системы, − столбец свободных членов, − столбец неизвестных. Систему решите матричным методом по формуле: , выполнив действия, записанные ниже.

· В ячейках L3:M4 запишите матрицу А.

· Найдите обратную матрицу . Для этого выделите область соответствующего размера диапазон L6:M7 для записи обратной матрицы. Вызовите формулу массива

→ Математические → МОБР → ОК.

Задайте адреса матрицы А L3:M4 и нажмите Ctrl+Shift+Enter для выполнения действия.

· В ячейках M9:M10 запишите столбец свободных членов B.

· Умножьте матрицу на В и результат запишите в ячейках M12:M13. Для этого выделите диапазон ячеек M12:M13. Вызовите формулу массива

→ Математические → МУМНОЖ → ОК.

Задайте адреса перемножаемых матриц L6:M7 и M9:M10 и нажмите Ctrl+Shift+Enter для выполнения действия.

Таким образом, искомая эмпирическая функция имеет вид .

Рис.2. Рабочий лист 1

4. Постройте точечную диаграмму зависимости t от d. используя Мастер диаграмм, который вызывается с помощью команды Вставка → Диаграмма или при нажатии соответствующей кнопки на панели инструментов.

· Выделите данные, которые будут включены в диаграмму, т.е. ячейки B4:C13.

· Выберите команду Вставка → Диаграмма (можно щелкнуть кнопку Мастер диаграмм на панели инструментов).

· На первом шаге работы Мастера диаграмм выберите Тип диаграммы вариант Точечная и щелкните Далее. Можно пропустить остальные шаги и сразу нажать кнопку Готово. На экране появится встроенный график. При необходимости его можно перетащить в более удобное место, зацепив мышью за край области диаграммы.

· Добавьте к исходной точечной диаграмме линию тренда. Для этого используйте команду меню Диаграмма → Добавить линию тренда, которая появляется, если диаграмма выделена.

· В окне Линия тренда во вкладке Тип и выберите шаблон кривой. В нашем случае Линейная . Во вкладке Параметры поставьте флажок: Показывать уравнение на диаграмме и нажимите ОК.

На диаграмме появится линия тренда и ее уравнение (рис 2).

Задание 2. Установить вид эмпирической зависимости y от x и методом наименьших квадратов найти ее параметры для следующих значений и :

0,78 1,56 2,34 3,12 3,81

2,50 1,20 1,12 2,25 4,28

Решение. На плоскости построим точки , .

Рис. 3. Точечная диаграмма

Из точечной диаграммы (рис.3) видно, что точки расположены в близи некоторой параболы, поэтому можно считать, что зависимость от будет квадратичной, т.е. имеет вид . Для вычисления коэффициентов системы (7) составим таблицу:

0,78 0,6084 0,4746 0,3702 2,50 1,9500 1,5210 1,56 2,4336 3,7964 5,9224 1,20 1,8720 2,9203 2,34 5,4756 12,8129 29,9822 1,12 2,6208 6,1327 3,12 9,7344 30,3713 94,7585 2,25 7,0200 21,9024 3,81 14,5161 55,3063 210,7172 4,28 16,3068 62,1289

11,61 32,7681 102,7615 341,7505 11,35 29,7696 94,6053

Система (7) для нашего примера имеет вид:

Решая её, например, по формулам Крамера, находим, что , , . Следовательно, зависимость от имеет вид:

.

Решение при помощи Excel. Все исходные, промежуточные данные и результаты решения изображены на рабочем листе2 (рис.4).

Рис.4. Рабочий лист 2

Система (6) для нашего примера принимает вид:

Следовательно, зависимость от имеет вид:

.

3. Индивидуальные задания

1. Методом наименьших квадратов найти параметры эмпирической зависимость, если известно, что зависимость от — линейная .

2,1 2,7 3,3 3,8 4,2 4,9 5,6 6,1 6,8

1,2 1,6 2,1 2,4 2,5 2,8 3,4 3,8 4,0

1.

0,6 1,3 1,8 2,7 3,1 3,9 4,2 5,1

10,2 8,2 6,0 5,1 1,5 0,8 -1,6 -2,8 -5,5

2.

-4,2 -3,7 -3,3 -2,6 -1,8 -1,1 -0,8 -0,4

0,1 0,5 1,1 1,8 2,9 3,6 4,0 4,7 4,9

3.

-10,0 -8,5 -6,5 -5,0 -2,0 2,0 4,6 7,0 9,5

0,8 0,5 0,3 -0,5 -1,7 -1,8 -2,5 -2,8

4.

-4,0 -3,5 -2,4 -2,0 -0,6 0,6 2,1 3,8 5,5

-1,6 -1,0 -0,8 -1,1 -0,7 -0,4 -0,2 0,3

5.

-6,0 -4,9 -3,8 -2,0 -0,6 0,6 2,1 3,8 5,5

-0,5 0,1 0,7 1,0 1,3 1,4 2,0 2,5

6.

-5,0 -4,0 -3,1 -2,2 -1,0 0,4 1,0 2,2 3,1

9,9 7,8 6,8 5,0 3,9 1,5 -1,4 -2,1

7.

3,0 -2,0 -1,5 -1,0 0,8 0,9 2,0 2,5 3,2

-2,1 -1,3 -1,0 0,5 0,7 1,5 1,6 2,3 2,7

8.

-5,3 -3,7 -2,3 0,8 2,1 3,6 5,9 7,8 9,2

2,6 2,4 1,7 1,3 1,0 0,7 0,1 -0,5 -0,6

9.

-5,0 -4,3 -3,7 -1,0 0,9 2,0 2,8 3,4 4,1

-2,2 1,5 2,7 8,3 11,4 14,5 15,3 17,1 18,4

10.

0,2 0,9 1,6 1,7 2,4 2,5 3,1 3,4 4,1

-2,2 1,5 2,7 8,3 11,4 14,5 15,3 17,1 18,4

11.

2,1 2,8 3,5 3,7 4,2 4,4 5,7 5,8 6,2

16,7 13,4 11,3 10,7 8,4 7,2 2,7 1,6 0,6

12.

0,4 0,9 1,2 1,6 2,3 2,5 3,2 3,8 4,2

1,7 3,3 3,8 5,0 6,3 7,5 8,8 10,4 11,7

13.

0,8 1,5 2,7 3,5 4,1 5,3 6,1 7,7 8,8

2,3 2,5 3,0 3,5 3,6 4,0 4,4 5,2 5,6

14.

2,0 1,5 1,1 0,5 -0,3 -1,1 -1,2 -2,0 -2,5

-1,1 1,0 2,1 2,3 2,9 4,0 5,1 5,8

15.

2. Установить вид эмпирической зависимости y от x и методом наименьших квадратов найти ее параметры для следующих значений и . Построить точечную диаграмму и график полученной эмпирической зависимости.

-1,9 -1,1 0,9 2,0 3,1 3,8 5,2 5,4

-19,5 -9,8 3,1 0,3 -1,0 -6,1 -14,4 -29,1 -34,5

1.

0,2 0,7 1,1 1,6 2,2 2,3 3,0 3,9 4,3

6,3 10,6 14,2 15,7 15,9 15,5 12,5 5,0 0,2

2.

-5,0 -4,2 -3,5 -2,8 -1,9 -1,2 -0,3 0,8 1,5

8,8 4,3 1,8 -0,2 -0,8 -0,5 1,8 7,6 12,2

3.

-6,0 -2,6 -2,1 -1,1 1,1 2,1 2,4 3,2

1,1 2,5 3,8 5,0 4,1 0,6 -4,6 -6,6 12,6

4.

-2,1 -1,5 -0,2 0,5 2,0 2,4 3,0 3,3 3,9

-10,2 -5,9 -3,2 2,1 6,3 5,6 4,4 3,4 1,1

5.

0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9

10,4 10,6 10,7 10,6 10,4 10,2 9,8 9,3 8,9

6.

-3,6 -3,1 -2,1 -1,1 0,5 1,1 2,1 3,1 3,5

-25,5 -16,4 -2,0 8,4 16,0 16,7 14,3 8,1 4,3

7.

-2,5 -2,1 -1,2 -0,3 1,1 2,0 3,1 3,2

38,1 30,8 15,2 13,6 6,8 2,0 5,5 15,2 17,5

8.

-3,2 -2,2 -1,6 -1,1 1,0 1,4 2,1 2,5

-16,4 -2,0 3,8 8,4 15,2 16,7 16,7 14,3 12,6

9.

-3,7 -2,5 -0,9 0,5 1,1 2,6 3,4 4,2 5,1

40,8 22,3 5,2 -0,5 -1,1 4,8 11,3 20,7 32,2

10.

-3,2 -2,6 -2,0 -0,6 0,5 0,9 1,7 2,0

0,9 1,6 2,7 3,7 4,2 4,0 3,7 3,5 2,3

11.

0,7 1,4 1,9 2,2 3,0 3,9 4,1 4,6

8,8 3,4 1,6 1,2 1,1 2,3 7,4 9,2 13,5

12.

-2,1 -1,6 -0,3 0,4 1,9 2,3 2,9 3,2 3,8

-9,9 -5,0 4,1 3,0 7,2 5,5 5,3 4,3 2,0

13.

-5,1 -4,3 -3,6 -2,9 -2,0 -1,3 -0,4 0,9 1,6

-7,7 -3,2 -0,7 1,3 1,9 1,6 -0,7 -6,5 -11,1

14.

0,4 0,9 1,3 1,8 2,3 2,6 3,3 4,1 4,5

4,4 8,7 12,2 13,7 13,9 13,5 10,5 3,0 -1,8

15.

3. Контрольные вопросы.

1. Какая зависимость называется эмпирической?

2. Из каких двух этапов состоит процесс получения эмпирической зависимости?

3. Из каких соображений выбирается вид эмпирической зависимости?

4. В чем состоит сущность метода наименьших квадратов?

5. Откуда находятся значения параметров эмпирической зависимости.

6. В чем состоит преимущество метода наименьших квадратов?

7. В чем недостаток метода наименьших квадратов?

8. Что дает анализ точечной диаграммы и графика эмпирической зависимости?