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

126 yx2 = yx2 + yy(i) *xx(i) ^2; yx3 = yx3 + yy(i) *xx(i) ^3; end A = [x6 x5 x4 x3; x5 x4 x3 x2; x4 x3 x2 x1; x3 x2 x1 n] B = [yx3; yx2; yx1; y1] c = A\B for i = 1: n f(i) = c(1) *xx(i) ^3 + c(2) *xx(i) ^2 + c(3) *xx(i) + c(4) end plot (xx, yy, 'ko'); hold on plot (xx, f, 'r-'); for i = 1: n d(i) = abs(yy(i) -f(i)) end Приложение 3 Нахождение коэффициентов по МНМ D = [21952 784 28 1 19683 729 27 1 17576 676 26 1 13824 576 24 1]; H = [9.85; 9.87; 9.84; 9.92]; c = D\H; xx = [1 2 3 4 5 7 8 9 10 11 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28]; yy = [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]