Методы наименьших квадратов и наименьших модулей в научно-технических расчетах

124 ПРИЛОЖЕНИЕ Приложение 1 Интерполяция по методу наименьших квадратов x = [1 2 3 4 5 7 8 9 10 11 14 15 16 17 18 19 20 21 22 23 24 26 26 27 28]; y = [12 11 11 10.8 10.8 10.8 10.6 10.8 10.3 10.3 10.3 10.4 10.5 10.2 ... 10 9.95 10.14 10.06 10.25 9.99 9.92 9.96 9.84 9.87 9.85]; xx = 29; yy = 9.69; % вывод графика табличной функции маркерами plot(x,y,'ko',xx,yy,'go'); % вычисление коэффициентов полиномов разных степеней p2 = polyfit(x,y,2); p3 = polyfit(x,y,3); p4 = polyfit(x,y,4); % построение графиков полиномов t = [1:0.01:29]; P2 = polyval(p2,t); P3 = polyval(p3,t); P4 = polyval(p4,t); hold on grid on plot(t,P2,'k-',t,P3,'r-',t,P4,'b-')