Вы можете написать быстрый код, чтобы сделать это.Решение задачи методом наименьших квадратов является b = (G ^ (T) G) ^ (- 1) G ^ (T) m, где b - вектор коэффициентов линейного приближения, m - вектор всех yзначения и G является ядром матрицы размерности (длина (м), 2), и каждая строка G является (1, x), где x является значением x, соответствующим значению yy.G ^ T - это транспонирование G.
Транспонирование и инверсия могут быть выполнены с использованием библиотеки LAPACK (или любой другой библиотеки линейной алгебры на ваш выбор) или путем извлечения соответствующего кода из книги «Числовые рецепты».