|
|
| Тема |
 Още в началото щях да ти дам Excel file [re: Sargonlll] |
|
| Автор |
Fire_Fury (Огнен) |
|
| Публикувано | 26.06.12 12:16 |
|
|
|
в който с помощта на матрици можеш да сметнеш коефициентите на регресия за уравненията от вида
a*x1^2 + b*x1 + c = y1
a*x2^2 + b*x2 + c = y2
a*x3^2 + b*x3 + c = y3
. . .
В твоя случай тези уравнения са от вида A*X = B, където
А е 4 на 3 матрица
x1^2, x1, 1
x2^2, x2, 1
x3^2, x3, 1
x4^2, x4, 1
X е 3 на 1 матрица от неизвестните коефициенти a, b, c, а
матрицата B е 4 на 1 от измерванията на Y: y1, y2, y3, y4
Тогава имаме X = (AT*A)-1*AT*B.
Виж и Excel файла parabola.xls на следния адрес
Изчислените коефициенти са в региона B24:B26.
Като умножиш A*X получаваш Bcalc, което в твоя случай е с много лошо RMSD.
Въобще на Excel можеш да решиш по този начин многопроменлива регресия и всяка една друга еднопроменлива, стига коефициентите да са линейни във формулата (т.е. полином). Трябва само да знаеш как се прилагат в Excel умножение на матрици, обръщане и транспониране, с функциите =mmult(), =minverse() и =transpose(). За тези функции и матриците виж материала, който съм написал за студентите по Компютърна химия от ПУ в дисциплината Хемометрика
AT и A-1 означават транспонирана матрица и обратна матрица.
| |
| |
| 
|
|
