НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ "ЛЬВІВСЬКА ПОЛІТЕХНІКА"




АПРОКСИМАЦІЯ ЕКСПЕРИМЕНТАЛЬНИХ ЗАЛЕЖНОСТЕЙ МЕТОДОМ НАЙМЕНШИХ КВАДРАТІВ

І н с т р у к ц і я
до лабораторної роботи
з курсу "Числові методи і моделювання на ЕОМ, ч.1"





Львів 2009
АПРОКСИМАЦІЯ ЕКСПЕРИМЕНТАЛЬНИХ ЗАЛЕЖНОСТЕЙ
МЕТОДОМ НАЙМЕНШИХ КВАДРАТІВ
Мета роботи: вивчення і програмна реалізація алгоритмів побудови апроксимуючих залежностей методом найменших квадратів. Одержання навиків програмування і роботи в середовищах Сі.
ОСНОВНІ ВІДОМОСТІ
Нехай в результаті вимірювань в процесі експериментальних досліджень одержана таблиця залежності функції f від аргументу х:
Таблиця №1
Треба знайти формулу, яка виражає цю залежність аналітично.
Поставимо задачу так: знайти функцію заданого вигляду:
y= f(x), (1)
яка в точках х1, х2, ..., хn приймає значення якнайближче до табличних значень y1, y2,..., yn, але яка одночасно враховує характер експериментально знайденої функції.
Практично вигляд наближеної функції f(x) можна визначити наступним способом. За таблицею 1 будують точковий графік функції f(x), а відтак проводять плавну криву, яка відображає характер розташування експериментальних точок (див. рис. 1).
Рис.1. Вигляд апроксимуючої експериментальні точки кривої.
За одержаною таким чином кривою встановлюють (вибирають) вигляд наближеної функції (часто з числа простих за виглядом аналітичних функцій).
Процес наближення для експериментально встановленої функціональної залежності y= f(x) за методом НК складається з двох етапів: спочатку вибирають вид формули і вже після цього визначають числові значення параметрів, для яких наближення буде найкращим.
Вид формули наближення вибирають виходячи з теоретичних міркувань або з числа елементарних функцій за виглядом точкового графіку.
Розглянемо метод знаходження параметрів наближеної функції в загальному вигляді на прикладі наближеної функції з трьома параметрами:
EMBED Equation.3 (1)
Сума квадратів різниць відповідних значень функцій f(x) і EMBED Equation.3 має вигляд EMBED Equation.3.
Ця сума є функцією трьох змінних (параметрів EMBED Equation.3). Задача зводиться до знаходження її мінімуму. Застосуємо необхідну умову
екстремуму: EMBED Equation.3; EMBED Equation.3; EMBED Equation.3.
З врахуванням того, що EMBED Equation.3; EMBED Equation.3; EMBED Equation.3,
одержимо систему EMBED Equation.3
Розв'язавши цю систему трьох рівнянь з трьома невідомими відносно параметрів EMBED Equation.3одержимо конкретний вигляд шуканої функції EMBED Equation.3.
Як видно з розглянутого прикладу, зміна кількості параметрів не міняє суті даного методу, а виразиться лише у зміні кількості рівнянь у системі (5).
Природно очікувати, що значення знайденої функції EMBED Equation.3 у точках х1, х2, ..., хn відрізнятимуться від табличних значень y1, y2,..., yn. Значення різниць
EMBED Equation.3; EMBED Equation.3, (6)
називають відхиленнями виміряних значень EMBED Equation.3 від обчислених за формулою (4).
Для знайденої емпіричної формули (4) по відношенню до таблиці 1 можна знайти суму квадратів відхилень
EMBED Equation.3
яка згідно з принципом найменших квадратів для заданого виду наближеної функції має бути найменшою.
З двох різних наближень однієї і тієї ж табличної функції ліпшим є те для якого сума (7) має менше значення. Отже величину EMBED Equation.3 можна використовувати для вибору найліпшої з розглядуваних функцій EMBED Equation.3, які наближують задану функцію EMBED Equation.3. При цьому слід пам'ятати, що параметр EMBED Equation.3 при наближенні експериментальних даних є, по-перше, розмірною величиною і, по-друге, не дозволяє оцінити якість кожного окремого наближення.
Для оцінки якості наближення можна використовувати також відносне середнє квадратичне відхилення:
EMBED Equation.3
При виборі вигляду функції y= f(x) за характером точкового графіку слід розглядати такі випадки:
1) вигляд точкового графіку вказує на те, що функція є монотонно-зростаючою (спадною) і не має точок перегину;
2) функція має екстремум;
3) функція має точки перегину, або точки перегину і екстремуми.
У першому випадку достатню точність забезпечить наближення елементарними функціями (лінійною, степеневою, показниковою і т.д.). Зокрема, якщо значення абсциси значно більші від значень ординати, то застосовують логарифмічну, або дробово-логарифмічну залежності.
Для випадку 2) доцільно застосовувати наближення квадратним тричленом.
У випадку 3) потрібно проводити наближення складнішими функціями: поліномом вищого порядку, тригонометричними функціями і іншими.
Якщо характер точкового графіку відповідає вимогам першого випадку, то часто застосовують елементарні двопараметричні функції наведені в таблиці № 2.
Параметри EMBED Equation.3 і EMBED Equation.3 кожного з наведених рівнянь можна одержати складанням і розв'язуванням системи рівнянь типу (5). Саме таким шляхом одержані формули для параметрів EMBED Equation.3 і EMBED Equation.3 лінійної регресії. Формули для розрахунку параметрів всіх інших вищенаведених рівнянь одержані нижче через формули лінійної регресії шляхом перетворення вихідного рівняння регресії до лінійного виду. Формули для визначення параметрів EMBED Equation.3 і EMBED Equation.3 вищевказаних наближуючих функцій зведені в таблицю №2.
Таблиця №2
Формули розрахунку параметрів елементарних рівнянь регресії.
EMBED Equation.3; EMBED Equation.3, EMBED Equation.3; (8)
EMBED Equation.3; (9)
EMBED Equation.3; (10)
EMBED Equation.3; EMBED Equation.3; EMBED Equation.3
ПРИКЛАД ПОБУДОВИ РОЗРАХУНКОВИХ ФОРМУЛ
Нехай як апроксимуючу залежність вибрано показникову функцію EMBED Equation.3. Для розрахунку коефіцієнтів a і b в таблиці 2 дано формули: a=T3/G1; b=exp(H3/G1).
Вираз для знаходження G1 будується наступним чином:
згідно з (8) EMBED Equation.3, де EMBED Equation.3 - перший елемент масиву EMBED Equation.3; EMBED Equation.3. Тоді EMBED Equation.3, а так як сумування потрібно здійснити для EMBED Equation.3 пар значень, EMBED Equation.3.
Щоб побудувати розрахункові формули для T3, H3, потрібно вияснити зміст елементів двомірного масиву Z - EMBED Equation.3, EMBED Equation.3. Для T3, G3 j=3, отже Z =Z =x, Z =Z =v. Тоді формули для T3, G3 приймуть такий вигляд :
T3= n S x v - S x S v ; H3= S v S x - S x S x v .
Формули розрахунку параметрів квадратного тричлена
EMBED Equation.3
EMBED Equation.3;
EMBED Equation.3; EMBED Equation.3;
EMBED Equation.3; EMBED Equation.3
EMBED Equation.3; EMBED Equation.3;
EMBED Equation.3; EMBED Equation.3.
ПОРЯДОК ВИКОНАННЯ РОБОТИ
1. Побудувати точковий графік експериментальної залежності згідно з отриманим завданням. За виглядом графіку вибрати формули функцій для апроксимації.
2. Скласти блок-схему алгоритму розрахунку коефіцієнтів апроксимаційних формул та похибок апроксимації.
3. Розробити програми мовами TurboС++. Після перевірки їх викладачем відредагувати програми у відповідному середовищі і отримати результати обчислень.
4. Побудувати в одній системі координат графіки апроксимуючих функцій.
5. Оформити звіт згідно з поданою схемою.
РЕКОМЕНДОВАНА ЛІТЕРАТУРА
1. Заварыкин В.М. и др. Численные методы. - М.: Просвещение, 1990.
2. Брановицкая С.В. Вычислительная математика в химии и химической технологии. К: Высшая школа, 1986.
4. Гаврилюк М.А. и др. Прикладные программы и лабораторный практикум для персонального компьютера. - Киев: УМК ВО, 1988.
ЗАВДАННЯ НА САМОСТІЙНУ РОБОТУ
Додаток
СХЕМА ОФОРМЛЕННЯ ЗВІТУ
1. Завдання.
2. Точковий графік залежності. Обґрунтування вибору функцій для наближення.
3. Формули для розрахунку коефіцієнтів вибраних функцій.
4. Блок-схема алгоритму.
5. Програмна реалізація алгоритму.
6. Результати обчислень. Графіки апроксимаційних функцій.
7. Порівняння точності апроксимації різними функціями.
Висновки.