11 мая 2011

gnuplot : аппроксимация полиномом

Исходные данные: файл test.dat c экспериментальными точками (x, y, dy).

Задача: построить график по исходным точкам и вписать в них полином. 

Решение:
итоговый код gnuplot

#Определяем функцию аппроксимации
f(x)=a*x*x+b*x+c

#Находим параметры аппроскимаци
fit f(x) 'test.dat' using 1:2:3 via a,b,c

#Задаём заметки, которые будут выведены на график
set label 1 "Y=AX^2+BX+C" at 11,1225
set label 2 "A = %g", a at 12,1200
set label 3 "B = %g", b at 12,1180
set label 4 "C = %g", c at 12,1160

#Выводим график аппроксимации и исходных точек одной командой plot
#с (with) отображением ошибок по Y (yerr)
plot f(x),'test.dat' using 1:2:3 with yerr

В результате получаем нечто подобное:


Код и рисунок взят с сайта "not so Frequently Asked Questions".

Комментариев нет:

Отправить комментарий