Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

РГР / prakticheskaya_rabota_metod_naimenshih_kvadratov

.doc
Скачиваний:
116
Добавлен:
22.02.2014
Размер:
81.41 Кб
Скачать

Практична робота №1

Тема: Використання методу найменших квадратів для створення моделі і її дослідження.

Мета: Закріпити навички використання методу найменших квадратів.

Теоретичні відомості.

Для досягнення поставленої мети необхідно визначити функціональну залежність за допомогою точкової апроксимації функції f(x).

Дана вибірка значень функції f(x) де xi (i=1,n).

Необхідно підібрати функцію (x)=c00+ c11+…+ cmm , що зветься загальним поліномом, таким чином щоб :

Відомо, що якщо система функцій — лінійно незалежна, то існує многочлен Фm і він є єдиним. Так як узагальнений поліном m(х) залежить від параметрів С1, С2,…Сm то рівняння (1) можна подати у вигляді :

Функціонал S є позитивним, тому min S існує і причому єдиний. Використовуючи необхідні умови екстремуму функцій кількох змінних отримаємо нормальну систему рівнянь для визначення коефіцієнту С:

Хід роботи.

  1. Визначимо місце для застосування даного методу аналізу:

Розглянемо задачу прогнозування, яка передбачає визначення урожайності сільськогосподарської культури на 1997 рік, маючи дані про попередні урожайності з 1990 по 1996 роки.

  1. Побудуємо таблицю, яка описує задачу:

Рік

Урожайність картоплі

1990

1991

1992

1993

1994

1995

1996

Урожайність

117

95

119

137

105

96

119

Дані наведені в таблиці приведено у додатку А.

Треба визначити скільки буде зібрано урожаю у 1997 році (на основі даних, поданих в таблиці).

Задаємо функція залежності:

1. Записуємо функцію S.

  1. Обрахуємо вираз.

3.Знайдемо частинні похідні:

Знаходимо частинні похідні по .

  1. Складаємо систему і знаходимо коефіцієнти:

>> a=[7 28 140;28 140 784;140 784 4676]

a =

7 28 140

28 140 784

140 784 4676

>> b=[788; 3146; 15672]

>> x=a\b

x =

107.5749

3.6882

-0.4878

С0=107.7143

С1= 3.5952

С2= -0.4762

Таким чином модель набуває вигляду:

Q = 107. 7143+ 3. 5952* t -0. 4762* t2 ;

  1. Розрахуємо значення Q в точці

Q1=107.7143+ 3.5952* 1 -0.4762*(1)2 =110.8333;

Q2=107.7143+ 3.5952* 2 -0.4762*(2)2 =113.0000;

Q3=107.7143+ 3.5952* 3 -0.4762*(3)2 =114.2143;

Q4=107.7143+ 3.5952* 4 -0.4762*(4)2 =114.4762;

Q5=107.7143+ 3.5952* 5 -0.4762*(5)2 =113.7857;

Q6=107.7143+ 3.5952* 6 -0.4762*(6)2 =112.1429;

Q7=107.7143+ 3.5952* 7 -0.4762*(7)2 =109.5476;

5. Знайдемо Q за 1997 рік:

Обрахуємо похибку(рівень апрокимації) за формулою:

, де

Q = 107. 7143+ 3. 5952* 8 -0. 4762* (8)2 = 106;

Отже врожайність на 1997 рік буде дорівнювати 106. Якщо звірити з даними поданими у додатку А, то можна побачити, що модель обрана правильно і обрахунок отримано точно.

Висновок:

У процесі виконання роботи були використані реальні статистичні дані і отримана на основі них таблиця деяких залежностей. Було поставлено задачу на основі табличних даних побудувати одно кроковий прогноз за допомогою методу найменших квадратів.

Необхідно було знайти таку формулу, яка виражала б залежність аналітично. Одним із підходів до вирішення задачі є побудова інтерполяційного многочлена, значення якого в точках х1, х2, .. , хn будуть співпадати з відповідними значеннями функції. Проте беззаперечне спів падання значень є неоправданим, якщо значення функції відомі з деякою похибкою. Тому з самого початку треба враховувати характер початкової функції, а це значить, що функція має бути заданого виду: y=F(x), яка в точках х1, х2, .. , хn приймає значення як можна найближчі до табличних y1, y2, .. , yn. Від того яка була вибрана функція і буде залежати точність прогнозу.

Даний метод використовується для обробки статистичних даних.