matlab n次多项式拟合代码实现

最新推荐文章于 2023-04-26 14:59:15 发布

卡洛斯曹

最新推荐文章于 2023-04-26 14:59:15 发布

阅读量1.3k

点赞数

分类专栏： MATLAB 文章标签： matlab 算法

本文链接：https://blog.csdn.net/qq_35395845/article/details/117701460

版权

本文介绍如何在MATLAB中实现n次多项式拟合，使用最小二乘法求解拟合曲线的系数。通过化简矩阵和解线性方程组，得到拟合方程的详细步骤，最终计算出拟合方程的系数。

摘要由CSDN通过智能技术生成

MATLAB下n次多项式拟合，且Robust为OFF时的代码实现

% f(x) = a(n)*x^n + a(n-1)*x^(n-1) + ...+ a1*x +a0

    % 最小二乘法求拟合曲线系数
    % 拟合n次多项式
    n = 3;
    for k=1:1:(n+1)
        for g=1:1:(n+2)
            sumdata1 = 0;
            for j=1:1:length(x)
                if g <= (n+1)
                    sumdata1 = sumdata1 + (x^(k-1))*(x^(g-1));
                else
                    sumdata1 = sumdata1 + y(j)*(x^(k-1));
                end
            end
            A(k,g) = sumdata1;
        end
    end

    % 多次多项式拟合
    % f(x) = a(n)*x^n + a(n-1)*x^(n-1) + ...+ a1*x +a0
    % A(1,1)*a0 + A(1,2)*a1 + A(1,3)*a2 + ...+ A(1,n+1)*an= A(1,n+2)
    % A(2,1)*a0 + A(2,2)*a1 + A(2,3)*a2 + ...+ A(2,n+1)*an= A(2,n+2)