- •Аналіз та обробка експериментальних даних засобами MathCAD в проектуванні СОС
- •1 Робота з файлами даних
- •1.1 Запис даних y файл
- •1.2 Зчитування даних з файла
- •1.3 Імпортування даних з файла
- •1.4 Експортування даних y файл
- •1.6 Обмін даних з іншими програмами
- •1.6.1 Обмін даними з Excel.
- •1.7 Підготовка даних до аналізу
- •1.8 Створення Масивів, використовуючи індекси
- •1.8.1 Створення елементів масиву математично
- •1.8.2 Створення масиву прямим введенням елементів
- •1.8.3 Створення масиву, використовуючи таблицю
- •1.8.4 Інші шляхи створення масиву даних
- •1.9 Сортування даних
- •1.9.1 Звертання до підмножин даних
- •1.9.2 Розділення даних на групи, використовуючи програму
- •2 Регресія
- •2.1 Використання line і medfit в лінійній регресії
- •2.1.1 Функції лінійної регресії
- •2.1.2 Функція line
- •2.1.3 Експоненційні дані
- •2.1.4 Medfit
- •2.2 Багатовимірна регресія
- •2.2.1 Багатовимірна поліноміальна регресія
- •2.2.2 Результати функцій regress і loess
- •2.2.3 Прогляд результатів interp
- •2.2.4 Зауваження
- •2.3 Спеціалізовані підгоночні функції
- •2.3.1 Використання спеціалізованих підгоночних функцій
- •2.4 Лінійна регресія з використанням linfit
- •2.4.1 Функція linfit
- •2.4.2 Представлення результатів функції
- •2.5 Узагальнена нелінійна регресія (genfit)
- •2.5.1 Функція Genfit
- •3 Інтерполяція
- •3.1.1 cspline інтерполяція багатьох змінних
- •3.1.2 Функція cspline
- •3.1.3 Функція interp
- •3.1.4 Прогляд результатів interp
- •4 Згладжування даних
- •4.1 Вбудовані функції згладжування
- •4.1.1 Функція medsmooth
- •4.1.2 Функція ksmooth
- •4.1.3 Функція supsmooth
- •5 КОНТРОЛЬНІ ЗАПИТАННЯ
- •6 ЛАБОРАТОРНЕ ЗАВДАННЯ
- •8 СПИСОК ЛІТЕРАТУРИ
2Регресія
2.1Використання line і medfit в лінійній регресії
2.1.1Функції лінійної регресії
Функція line повертає параметри лінії найкращого
наближення згідно методу найменших квадратів. Для вхідних значень (x, y), що знаходяться у векторах vx і vy відповідно, згідно методу найменша квадратів, лінія найкращого наближення задається формулою:
y
line(vx,vy)0 + line(vx,vy)1 x
В подальшому для демонстрації лінійних методів аналізу будемо використовувати наступну матрицю з (x, y) даними (x – в 1-му стовпці, y – в 2-му) :
data :=
0 0.3
-1 -0.26
1 0.38
2.1.2Функція line
line повертає параметри (кут нахилу і зміщення) прямої, яка дає найкраще наближення для набору даних в vx і vy.
line отримує два параметри: line(vx, vy)
1. vx є вектор n-елементів дійсних х-значень даних:
vx := DATA 0 |
n := length (vx) |
n = 18 |
|
|
2. vy є вектор n-елементів дійсних у-значень даних:
vy := DATA 1 |
fit := line(vx,vy) |
|||
Функція line повертає: |
|
|
|
|
Зміщення: |
Нахил: |
|||
fit0 = 0.033 |
|
fit1 = 0.235 |
|
|
|
|
|||
|
|
|||
Тоді, найкраще наближення задається : f(x) := fit0 + fit1 x
15
Для |
оцінки |
та |
1.209 |
2 |
|
|
|
||
|
|
|
|
||||||
|
|
|
|
|
|||||
графічної |
візуалізації |
|
1 |
|
|
|
|||
наближення |
розпо- |
vy |
|
|
|
|
|||
ділу |
точок |
з |
2-х |
f(x) |
0 |
|
|
|
|
масивів vx і vy та |
|
|
|
|
|||||
|
|
|
|
|
|||||
отриманої |
|
прямої, |
|
1 |
|
|
|
||
можна |
скористатись |
−1.143 |
|
|
|
|
|||
графіком: |
|
|
|
2 |
5 |
0 |
5 |
||
|
|
|
|
||||||
|
|
|
|
|
|
|
min( vx) |
vx,x |
max( vx) |
|
Існує також вбудована функція співвідношення, яка |
||||||||
повертає коефіцієнт кореляції даних : corr(vx,vy) = 0.986 |
|
||||||||
|
Чим ближче це значення до 1, тим кращим є наближення. |
||||||||
2.1.3Експоненційні дані
Кут нахилу та зміщення є також корисними для експонентних даних. Особливо, якщо x і y зв'язані наступною
функціональною залежністю – y
A ekx .
Функція line може бути застосована і до логарифмічних
значень за допомогою виразу log(y) |
|
log(A) + k x, де коефіцієнти |
||||||
|
||||||||
|
||||||||
( 0) |
і |
1 . |
|
|
||||
A |
|
exp fit |
|
k |
|
fit |
||
2.1.4Medfit
Якщо розраховані параметри, отримані функцією line, не відповідають експериментальним даним, можна спробувати використати функцію medfit, що використовує серединну регресію. Така ситуація може виникнути, якщо дані мають екстремальні викиди. В цьому випадку серединна регресія може повернути лінію, яка відповідає таким даним більш точно.
Функція medfit повертає параметри (кут нахилу і зміщення) згідно методу найменших квадратів, які містяться в vx і vy – medfit(vx, vy):
1.vx є вектор n-елемента дійсних значень даних. Значення в vx відповідають x-значенням.
vx := DATA 0 |
n := length (vx) |
n = 18 |
|
|
16
