数值计算_第6章 曲线拟合的最小二乘法

最小二乘拟合的概念-概述说明以及解释1.引言1.1 概述最小二乘拟合是一种常用的数据分析方法，通过最小化观测值与拟合值之间的残差平方和来求取最优拟合曲线或平面，从而描述数据的模式和趋势。

该方法被广泛应用于统计建模、机器学习、信号处理、金融分析等领域。

最小二乘法的核心思想是寻找一条曲线或平面，使得该曲线或平面与数据点的残差之和最小。

通过最小二乘法，我们可以得到最佳拟合曲线或平面，从而对数据进行更准确的描述和预测。

因此，最小二乘拟合在数据分析中具有重要的意义。

本文将详细介绍最小二乘拟合的定义、原理和应用，从而帮助读者更好地理解和运用这一重要的数据分析方法。

1.2 文章结构文章结构部分的内容如下：文章结构部分将介绍整篇文章的组织结构和主要内容安排，以便读者对文章的整体框架有一个清晰的认识。

在本文中，主要分为引言、正文和结论三个部分。

- 引言部分包括对最小二乘拟合的概念进行简要介绍，阐述本文撰写的目的和重要性。

- 正文部分将详细讨论最小二乘拟合的定义、原理和应用，以便读者全面了解这一重要的数据分析方法。

- 结论部分将对最小二乘拟合的重要性进行总结，探讨最小二乘法在数据分析中的价值，并展望最小二乘拟合在未来的发展趋势。

通过这样的结构安排，读者可以清晰地了解本文的主要内容和章节布局，有助于他们更好地理解和掌握最小二乘拟合的相关知识。

1.3 目的本文的主要目的是介绍最小二乘拟合这一重要的数学方法。

通过对最小二乘拟合的定义、原理和应用进行详细讨论，希望读者能够深入了解这一方法在数据分析和模型拟合中的重要性。

此外，本文还将探讨最小二乘法在实际问题中的应用，以及展望未来最小二乘拟合在数据分析领域的发展趋势。

通过阐述这些内容，旨在让读者更加深入地理解和应用最小二乘拟合方法，为其在数据分析和模型拟合中提供有效的工具和思路。

2.正文2.1 最小二乘拟合的定义最小二乘拟合是一种常用的数学方法，用于通过调整参数来拟合一个数学模型以最小化观测数据和模型之间的残差平方和。

曲线拟合最小二乘法
最小二乘法是统计学中最常用的数据拟合方法，也被称为**最小平方法**。

该方法在数学和统计学中已经有很长的历史，广泛应用于各种学科的科学研究和实际应用。

最小二乘法的主要思想是最小化所给数据点与目标曲线之间的误差平方和，以此来确定目标曲线的参数。

具体而言，最小二乘法是根据**基函数**与参数之间的函数关系，采用多元函数去拟合所给数据点，旨在最小化拟合数据点和多元函数之间误差平方和的拟合方法。

最小二乘法可以用来拟合任何形式的曲线，在各种应用中都大量应用。

比如在政治学、经济学和心理学中，研究者通过最小二乘法来拟合某种结果与输入变量之间的联系，以更好地理解呈现结果的背景机制；在数值计算中，最小二乘法可用来拟合数值计算数据，从而精确地求解各种方程；而在工程学中，最小二乘法常用来拟合统计数据，估计影响工作效率的各种自变量。

总之，最小二乘法是一种统计学中经久不衰的拟合方法，可以用来拟合任何形式的曲线，在广泛的应用领域有着重要地位。

标准曲线的最小二乘法拟合和相关系数（合肥工业大学控释药物研究室尹情胜）1 目的用最小二乘法拟合一组变量（，，i＝1-n）之间的线性方程（y＝ax+b），表示两变量间的函数关系；（开创者：德国数学家高斯）一组数据（，，i＝1-n）中，两变量之间的相关性用相关系数（R）来表示。

（开创者：英国统计学家卡尔·皮尔逊）2 最小二乘法原理用最小二乘法拟合线性方程时，其目标是使拟合值（）与实测值（）差值的平方和（Q）最小。

式（1）3 拟合方程的计算公式与推导当Q最小时，；得到式（2）、式（3）：式（2）式（3）由式（3）和式（4），得出式（4）和式（5）：式（4）式（5）式（4）乘以n，式（5）乘以，两式相减并整理得斜率a：斜率（k＝xy／xx，n*积和-和积）式（6）截距b的计算公式为公式（5），也即：截距b＝（y-x）／n，差平均差）式（7）4 相关系数的意义与计算公式相关系数（相关系数的平方称为判定系数）是用以反映变量之间相关关系密切程度的统计指标。

相关系数（也称积差相关系数）是按积差方法计算，同样以两变量与各自平均值的离差为基础，通过两个离差相乘来反映两变量之间相关程度；着重研究线性的单相关系数。

相关系数r xy取值在-1到1之间。

r xy = 0时，称x,y不相关；| r xy | = 1时，称x,y完全相关，此时，x,y之间具有线性函数关系；| r xy | < 1时，X的变动引起Y的部分变动，r xy的绝对值越大，x的变动引起y的变动就越大，|r xy | > 0.8时称为高度相关，当0.5< | r xy|<0.8时称为显著相关，当0.3<| r xy |<0.5时，成为低度相关，当| r xy | < 0.3时，称为无相关。

(式（7）5 临界相关系数的意义5.1 临界相关系数中显著性水平（α）与置信度（P）的关系显著性水平取0.05，表示置信度为95%；取0.01，置信度就是99%。

曲线拟合最小二乘法
曲线拟合是指通过已知数据点来推导出一条函数曲线，使得该曲线尽
可能地贴近这些数据点。

而最小二乘法（Least Squares Method）是求解
这种拟合问题的一种常用方法。

最小二乘法的核心思想是尽量减小误差平方和。

假设已知的数据点为$(x_i, y_i)$，曲线函数为 $y=f(x)$，我们希望找到一组参数 $\theta$，使得 $f(x_i;\theta)$ 与 $y_i$ 的差距最小，即：
$$\min_{\theta}\sum_{i=1}^n [y_i - f(x_i;\theta)]^2$$。

这个式子被称为目标函数，也叫做残差平方和（RSS）。

通过对目标
函数进行求导，可以得到最优参数 $\theta^*$ 的解析解：
$$\theta^* = (\mathbf{X}^T \mathbf{X})^{-1} \mathbf{X}^T
\mathbf{y}$$。

其中，$\mathbf{X}$ 是一个 $n \times p$ 的矩阵，每一行代表一
个数据点的特征向量，$p$ 是曲线函数的参数个数。

$\mathbf{y}$ 是一
个 $n \times 1$ 的列向量，代表数据点的真实输出值。

最小二乘法在实际应用中有很广泛的应用。

例如，可以用它来构建多
项式回归模型、高斯过程回归模型等。

此外，在机器学习领域，最小二乘
法也被用于求解线性回归模型、岭回归模型等。

曲线拟合的数值计算方法实验Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】曲线拟合的数值计算方法实验【摘要】实际工作中，变量间未必都有线性关系，如服药后血药浓度与时间的关系；疾病疗效与疗程长短的关系；毒物剂量与致死率的关系等常呈曲线关系。

曲线拟合（curve fitting）是指选择适当的曲线类型来拟合观测数据，并用拟合的分析两变量间的关系。

曲线直线化是曲线拟合的重要手段之一。

对于某些非线性的资料可以通过简单的变量变换使之直线化，这样就可以按原理求出变换后变量的，在实际工作中常利用此直线方程绘制资料的标准工作曲线，同时根据需要可将此直线方程还原为，实现对资料的曲线拟合。

常用的曲线拟合有最小二乘法拟合、幂函数拟合、对数函数拟合、线性插值、三次样条插值、端点约束。

关键词曲线拟合、最小二乘法拟合、幂函数拟合、对数函数拟合、线性插值、三次样条插值、端点约束一、实验目的1.掌握曲线拟合方式及其常用函数指数函数、幂函数、对数函数的拟合。

2.掌握最小二乘法、线性插值、三次样条插值、端点约束等。

3.掌握实现曲线拟合的编程技巧。

二、实验原理1.曲线拟合曲线拟合是平面上离散点组所表示的坐标之间的函数关系的一种数据处理方法。

用解析表达式逼近的一种方法。

在或社会活动中，通过实验或观测得到量x 与y 的一组数据对（X i ，Y i )（i=1，2，...m ），其中各X i 是彼此不同的 。

人们希望用一类与数据的背景材料规律相适应的解析表达式，y=f(x ，c ）来反映量x 与y 之间的依赖关系，即在一定意义下“最佳”地逼近或拟合已知数据。

f(x ，c)常称作拟合模型 ，式中c=(c 1，c 2，…c n )是一些待定参数。

当c 在f 中出现时，称为线性模型，否则称为。

有许多衡量拟合优度的标准，最常用的一种做法是选择参数c 使得拟合模型与实际在各点的(或)，c）-f （f y e k k k 的平方和达到最小，此时所求曲线称作在加权最小二乘意义下对数据的拟合曲线。

第6章曲线拟合的最小二乘法6.1 拟合曲线通过观察或测量得到一组离散数据序列，当所得数据比较准确时，可构造插值函数逼近客观存在的函数，构造的原则是要求插值函数通过这些数据点，即。

此时，序列与是相等的。

如果数据序列，含有不可避免的误差（或称“噪音”），如图6.1所示；如果数据序列无法同时满足某特定函数，如图6.2所示，那么，只能要求所做逼近函数最优地靠近样点，即向量与的误差或距离最小。

按与之间误差最小原则作为“最优”标准构造的逼近函数，称为拟合函数。

图6.1 含有“噪声”的数据图6.2 一条直线公路与多个景点插值和拟合是构造逼近函数的两种方法。

插值的目标是要插值函数尽量靠近离散点；拟合的目标是要离散点尽量靠近拟合函数。

向量与之间的误差或距离有各种不同的定义方法。

例如：用各点误差绝对值的和表示：用各点误差按模的最大值表示：用各点误差的平方和表示：或（6.1）其中称为均方误差，由于计算均方误差的最小值的方法容易实现而被广泛采用。

按均方误差达到极小构造拟合曲线的方法称为最小二乘法。

本章主要讲述用最小二乘法构造拟合曲线的方法。

在运筹学、统计学、逼近论和控制论中，最小二乘法都是很重要的求解方法。

例如，它是统计学中估计回归参数的最基本方法。

关于最小二乘法的发明权，在数学史的研究中尚未定论。

有材料表明高斯和勒让德分别独立地提出这种方法。

勒让德是在1805年第一次公开发表关于最小二乘法的论文，这时高斯指出，他早在1795年之前就使用了这种方法。

但数学史研究者只找到了高斯约在1803年之前使用了这种方法的证据。

在实际问题中，怎样由测量的数据设计和确定“最贴近”的拟合曲线？关键在选择适当的拟合曲线类型，有时根据专业知识和工作经验即可确定拟合曲线类型；在对拟合曲线一无所知的情况下，不妨先绘制数据的粗略图形，或许从中观测出拟合曲线的类型；更一般地，对数据进行多种曲线类型的拟合，并计算均方误差，用数学实验的方法找出在最小二乘法意义下的误差最小的拟合函数。

数值分析中的最小二乘法与曲线拟合数值分析是现代理论与实践密切结合的一门交叉学科，其中最小二乘法和曲线拟合是其中两个非常重要的概念。

最小二乘法是一种数学运算方法，用于求解一组方程组的未知参数，使得每个方程的误差平方和最小。

在实际应用中，最小二乘法广泛应用于数据拟合、信号处理、回归分析等领域。

在数据拟合中，最小二乘法是一种常见的方法，它可以用于拟合曲线和函数。

它通过延伸曲线以获得局部数据之间的交点，并通过在它们上进行平均化的方法来尝试匹配数据。

最小二乘法的概念为我们提供了一个理论基础，以便在一定程度上预测新的数据中对象的行为或趋势。

但是，即使在相对简单的问题中，最小二乘法可能并不是最佳选择。

曲线拟合是对一系列数据进行插值的过程，以便获得与原始数据点更准确相匹配的曲线或函数。

曲线拟合可以通过在相邻数据点之间进行插值来完成。

在曲线拟合中，只有在数据有很好的统计关系或在相邻数据点
有很好的相关性时，才会产生准确的结果。

否则，结果可能并不
准确，因为这些结果取决于数据点的数量和分布。

需要注意的是，曲线拟合和最小二乘法并不是一个可以代替另
一个的工具。

它们的适用范围不同。

曲线拟合适用于对离散数据
点进行联合分析，而最小二乘法适用于求解连续数据的线性模型。

总之，数值分析中的最小二乘法和曲线拟合是非常实用的概念，可以应用于各种领域。

它们作为现代数据分析的主要工具之一，
不断吸引着越来越多的学者和工程师投入到其中，将继续发挥重
要作用。

§7 曲线拟合的最小二乘法7-1 一般的最小二乘逼近(曲线拟合的最小二乘法)最小二乘法的一般提法是:对给定的一组数据(,)(0,1,,)i i x y i m =L ,要求在函数类01{,,,}n ϕϕϕϕ=L 中找一个函数*()y S x =,使误差平方和22*222()001[()]min [()]m m m i i i i iS x i i i S x y S x y ϕδδ∈=====−=−∑∑∑其中 0011()()()()()n n S x a x a x a x n m ϕϕϕ=+++<L 带权的最小二乘法： 2220()[()()]mi i i i x S x f x δω==−∑其中()0x ω≥是［a,b ］上的权函数。

用最小二乘法求曲线拟合的问题，就是在()S x 中求一函数*()y S x =，使22δ取的最小。

它转化为求多元函数20100(,,,)()[()()]m n n i j j i i i j I a a a x a x f x ωϕ===−∑∑L 的极小点***01(,,,)n a a a L 问题。

由求多元函数极值的必要条件，有002()[()()]()0m n i j j i i k i i j k I x a x f x x a ωϕϕ==∂=−=∂∑∑ ),,1,0(n k L =若记 0(,)()()()m j k ij i k i i x x x ϕϕωϕϕ==∑ 0(,)()()()mk i i k i k i f x f x x d ϕωϕ==≡∑ ),,1,0(n k L = 则上式可改写为0(,)n kj j k j a d ϕϕ==∑ ),,1,0(n k L =这个方程称为法方程，矩阵形式.Ga d =其中 0101(,,,),(,,,)T T n n a a a a d d d d ==L L ，0001010111011(,)(,)(,)(,)(,)(,)(,)(,)(,)n n n n n n G ϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕ− = LL L L L L L 由于01,,,n ϕϕϕL 线性无关，故0G ≠，方程组存在唯一解*(0,1,,),k k a a k n ==L从而得到函数()f x 的最小二乘解为****0011()()()()n n S x a x a x a x ϕϕϕ=+++L 可证 *2200()[()()]()[()()]m mii i i i i i i x S x f x x S x f x ωω==−≤−∑∑ 故*()S x 使所求最小二乘解。

实验三函数逼近与曲线拟合、问题的提出：函数逼近是指“对函数类A中给定的函数f(x),记作f(x)・A,要求在另一类简的便于计算的函数类B中求函数p(x)・A,使p(x)与f (x)的误差在某中度量意义下最小”函数类A通常是区间[a,b]上的连续函数，记作C[a,b],称为连续函数空间，而函数类B通常为n次多项式，有理函数或分段低次多项式等，函数逼近是数值分析的基础。

主要内容有：(1)最佳一致逼近多项式(2)最佳平方逼近多项式(3 )曲线拟合的最小二乘法实验要求:1、构造正交多项式；2、构造最佳一致逼近；3、构造最佳平方逼近多项式；4、构造最小二乘法进行曲线拟合；5、求出近似解析表达式，打印出逼近曲线与拟合曲线，且打印出其在数据点上的偏差;6、探讨新的方法比较结果。

三、实验目的和意义:1、学习并掌握正交多项式的MATLAB编程；2、学习并掌握最佳一致逼近的MATLAB实验及精度比较;3、学习并掌握最佳平方逼近多项式的MATLAB实验及精度比较;4、掌握曲线拟合的最小二乘法；5、最小二乘法也可用于求解超定线形代数方程组；6、探索拟合函数的选择与拟合精度之间的关系；四、算法步骤：1、正交多项式序列的生成{ \ ( X)}o •:设\ ( X)是［a,b］上首项系数数，如果多项式序列{ \ ( X)}o：满足关系式则称多项式序列{ \(X)}o：为在［a,b］上带权的n次正交多项式。

1 )输入函数「(x)和数据a,b;2) 分别求(x n, j(x)),C j (x), j(x))的内积；. . n 2 (X n,®j(X)), ,3) 按公式①；：o(X)=1, -(X) =X n j j(X)计算；：n(X)，生成正交多项式;j鼻Wj(x),W j(x))流程图：开始a n=0的n次多项式，r(x)为［a,b］上权函；Q j秋A 0, jb(j, k)」(x) j(x) k(x)d(X> =a「(x)正交，称;:n (x)为［a,b］上带权「(x)cz>结束2、最佳一致逼近多项式f(x) C[a,b]，若存在 R*(x) H n 使得.：(f,P ；^E n ，则称 P ； (x)是 f (x)在[a,b]上的最佳一致逼近多项式或最小偏差逼近多项式，简称最佳逼近多项式。

Excel拟合曲线用的最小二乘法1. 介绍Excel作为一款常用的办公软件，被广泛应用于数据分析和处理，而拟合曲线是数据分析中常用的方法之一。

拟合曲线用的最小二乘法是一种常见的拟合方法，通过最小化数据点与拟合曲线之间的距离来找到最佳拟合曲线，从而对数据进行预测和分析。

在本文中，我将从深度和广度的角度来探讨Excel拟合曲线用的最小二乘法，带你深入探索这一主题。

2. 最小二乘法的原理在Excel中进行曲线拟合时，最小二乘法是一种常用的拟合方法。

其原理是通过最小化残差平方和来找到最佳拟合曲线。

残差是指每个数据点到拟合曲线的垂直距离，最小二乘法通过调整拟合曲线的参数，使得残差平方和最小化，从而得到最佳拟合曲线。

在Excel中，可以利用内置函数或插件来实现最小二乘法的曲线拟合，对于不同类型的曲线拟合，可以选择不同的拟合函数进行拟合。

3. Excel中的拟合曲线在Excel中进行拟合曲线时，首先需要将数据导入Excel，然后利用内置的数据分析工具或者插件来进行曲线拟合。

通过选择拟合函数、调整参数等操作，可以得到拟合曲线的相关信息，如拟合优度、参数估计值等。

可以根据拟合曲线的结果来对数据进行预测和分析，从而得到对应的结论和见解。

4. 个人观点与理解对于Excel拟合曲线用的最小二乘法，我认为这是一种简单而有效的数据分析方法。

它能够快速对数据进行拟合，并得到拟合曲线的相关信息，对于数据的预测和分析具有一定的帮助。

然而，也需要注意到拟合曲线并不一定能够准确描述数据的真实情况，需要结合实际背景和专业知识进行分析和判断。

在使用最小二乘法进行曲线拟合时，需要注意数据的可靠性和拟合结果的可信度，以避免出现不准确的结论和偏差的情况。

5. 总结通过本文的探讨，我们对Excel拟合曲线用的最小二乘法有了更深入的了解。

最小二乘法的原理、Excel中的实际操作以及个人观点与理解都得到了充分的展示和探讨。

在实际应用中，需要结合具体情况和专业知识来灵活运用最小二乘法进行曲线拟合，从而得到准确的分析和预测结果。

一、曲线拟合是什么？曲线拟合也就是求一条曲线,使数据点均在离此曲线的上方或下方不远处, 它既能反映数据的总体分布,又不至于出现局部较大的波动, 能反映被逼近函数的特性,使求得的逼近函数与已知函数从总体上来说其偏差按某种方法度量达到最小。

设函数y=f(x)在m个互异点的观测数据为求一个简单的近似函数φ(x)，使之“最好”地逼近f(x),而不必满足插值原则。

这时没必要取φ(xi) = yi, 而要使i=φ(xi)yi 总体上尽可能地小。

这种构造近似函数的方法称为曲线拟合,称函数y=φ(x)为经验公式或拟合曲线。

如下为一个曲线拟合示意图。

清楚什么是曲线拟合之后，我们还需要了解一个概念——残差。

曲线拟合不要求近似曲线严格过所有的数据点，但使求得的逼近函数与已知函数从总体上来说其偏差按某种方法度量达到总体上尽可能地小。

若令（1-1）则为残向量（残差）。

“使（1-1）尽可能地小”有不同的准则（1）残差最大值最小（2）残差绝对值和最小（绝对值的计算比较麻烦）（3）残差平方和最小（即最小二乘原则。

计算比较方便，对异常值非常敏感，并且得到的估计量具有优良特性。

）二、最小二乘法是什么？个人粗俗理解：按照最小二乘原则选取拟合曲线的方法,称为最小二乘法。

百度百科：最小二乘法（又称最小平方法）是一种数学优化技术。

它通过最小化误差的平方和寻找数据的最佳函数匹配。

利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。

最小二乘法还可用于曲线拟合。

其他一些优化问题也可通过最小化能量或最大化熵用最小二乘法来表达。

三、求解最小二乘法（包含数学推导过程）我们以最简单的线性模型来解释最小二乘法。

什么是线性模型呢？监督学习中，如果预测的变量是离散的，我们称其为分类（如决策树，支持向量机等），如果预测的变量是连续的，我们称其为回归。

回归分析中，n个自变量，且因变量和自变量之间是线性关系，则称为一/多元线性回归分析。

曲线拟合的最小二乘法姓名：学号：专业：材料工程学院：材料科学与工程学院科目：数值分析曲线拟合的最小二乘法一、目的和意义在物理实验中经常要观测两个有函数关系的物理量。

根据两个量的许多组观测数据来确定它们的函数曲线，这就是实验数据处理中的曲线拟合问题。

这类问题通常有两种情况：一种是两个观测量x 与y 之间的函数形式已知，但一些参数未知，需要确定未知参数的最佳估计值；另一种是x 与y 之间的函数形式还不知道，需要找出它们之间的经验公式。

后一种情况常假设x 与y 之间的关系是一个待定的多项式，多项式系数就是待定的未知参数，从而可采用类似于前一种情况的处理方法。

在两个观测量中，往往总有一个量精度比另一个高得多，为简单起见把精度较高的观测量看作没有误差，并把这个观测量选作 x，而把所有的误差只认为是y的误差。

设 x 和 y 的函数关系由理论公式y＝f（x；c1，c2，……cm）（0-0-1）给出，其中 c1，c2，……cm 是 m 个要通过实验确定的参数。

对于每组观测数据（xi，yi）i＝1，2，……，N。

都对应于xy 平面上一个点。

若不存在测量误差，则这些数据点都准确落在理论曲线上。

只要选取 m 组测量值代入式（0-0-1），便得到方程组yi ＝ f （x ；c1 ，c2 ，……cm）（0-0-2）式中 i＝1，2，……，m.求 m 个方程的联立解即得 m 个参数的数值。

显然N<m 时，参数不能确定。

y 2 y 在 N>m 的情况下，式（0-0-2）成为矛盾方程组，不能直接用解方程的方法求得 m 个参数值，只能用曲线拟合的方法来处理。

设测量中不存在着系统误差，或者说已经修正，则 y 的观测值 yi 围绕着期望值 <f （x ；c1，c2，……cm）> 摆 动，其分布为正态分布，则 yi 的概率密度为p y i1 exp,式中i是分布的标准误差。

为简便起见，下面用 C 代表（c1，c2，……cm ）。

曲线拟合--最小二乘法1：已知平面上四个点：（0，1）、（1，2.1）、（2，2.9）和（3，3.2），求出一条直线拟合这四个点，使得偏差平方和变为极小。

解：设直线方程为：0 1 0 01 2.1 1 2.12 2.9 4 5.83 3.2 9 9.6Sum=6 Sum=9.2 Sum=14 Sum=17.5 代入正规方程：，编程求解上方程组：>> eq1='14*A+6*B=17.5';>>eq2='6*A+4*B=9.2';>> [A,B]=solve(eq1,eq2,'A,B');>> disp(A)0.74>> disp(B)1.19所以直线方程为：2：已知数据如下表所示1 2 4 610 5 2 1试求（1）用抛物线拟合这些数据使得偏差平方和最小；（2）用型如的函数来拟合这些数据使得偏差平方和最小。

（3）比较这两种拟合结果。

解：（1）设抛物线方程为：1 10 1 1 1 10 102 5 4 8 16 10 20 4 2 16 64 256 8 326 1 36 216 1296 6 36 Sum=13 Sum=18 Sum=57 Sum=289 Sum=1569 Sum=34 Sum=98代入正规方程：得到系数A，B，C的方程组：编程求解上方程组：>>eq1='1569*A+289*B+57*C=98';>>eq2='289*A+57*B+13*C=34';>>eq3='57*A+13*B+4*C=18';>> [A,B,C]=solve(eq1,eq2,eq3,'A,B,C');>> disp(A); disp(B); disp(C)102/199-1048/1992848/199>> A=102/199; disp(A) 0.5126>> B=-1048/199; disp(B) -5.2663>> C=2848/199; disp(C) 14.3116所以得到抛物线的方程为：（2）设函数1 10 1 1 102 5 1/2 1/4 5/24 2 1/4 1/16 1/26 1 1/6 1/36 1/6Sum=13 Sum=18 Sum=23/12 Sum=193/144 Sum=79/6 得到系数A，B的方程组：编程求解上方程组：>> eq1='4*A+23*B/12=18';>>eq2='23*A/12+193*B/144=79/6';>> [A,B]=solve(eq1,eq2,'A,B');>> disp(A); disp(B)-160/243872/81>> A=-160/243; disp(A)-0.6584>> B=827/81; disp(B)10.2099所以得到的函数为：（3）比较（1）和（2）两种方法拟合的方程：编程画出抛物线的图像为：>> x=-2:0.1:12;>> y=0.5126*x.^2-5.2663*x+14.3116;plot(x,y);grid on（a）再编程画出的图像为：>> x=-2:0.1:12;>> y=-0.6584+10.2099*(x.^(-1));>> plot(x,y);grid on>> x=-1:0.01:1;>> y=-0.6584+10.2099*(x.^(-1));plot(x,y);grid on(b)比较两图像可知，图像（b）在点（0，0）处不连续。

曲线拟合的线性最小二乘法拟合是已知有限个数据点，求近似函数，不要求过已知数据点，只要求在某种意义下它在这些点上的总偏差最小。

线性最小二乘法曲线拟合问题的提法是，已知一组（二维）数据，即平面上的n 个点(,),i i x y 1,2,,i n =⋅⋅⋅，i x 互不相同，寻求一个函数（曲线）()y f x =，使()f x 在某种准则下与所有数据点最为接近，即曲线拟合的最好。

线性最小二乘法是解决曲线拟合最常用的方法，基本思路是，令1122()()()(),m m f x a r x a r x a r x =++⋅⋅⋅+其中：()k r x 是事先选定的一组线性无关的函数；k a 是待定系数(1,2,,;k m =⋅⋅⋅)m n <。

拟合准则是使(1,2,,)i y i n =⋅⋅⋅与()i f x 的距离i δ的平方和最小，称为最小二乘准则。

1.系数k a 的确定 记[]221211(,,,)()nnm i i i i i J a a a f x y δ====-∑∑为求12,,,m a a a ⋅⋅⋅使J 达到最小，只需利用极值的必要条件0jJa ∂=∂(1,2,,)j m =⋅⋅⋅，得到关于12,,,m a a a ⋅⋅⋅的线性方程组11()[()]0,1,2,,n mjik kiii k r x a r x y j m ==-==∑∑，即111[()()](),1,2,,.m n nkjikijiik i i a r x r x r x y j m =====∑∑∑ （1.1）记1111()()()()m n m n n mr x r x R r x r x ⨯⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ [][]TT1212,,,,,,,m n A a a a Y y y y =⋅⋅⋅=⋅⋅⋅方程（1.1）可表为T T .R RA R Y = （1.2） 当{}12(),(),,()m r x r x r x ⋅⋅⋅线性无关时，R 列满秩，T R R 可逆，于是方程组（1.2）有唯一解()1TT .A R R R Y -=2.函数()k r x 的选取面对一组数据(,),1,2,,i i x y i n =⋅⋅⋅，用线性最小二乘法作曲线拟合时，首要的也是关键的一步是恰当地选取12(),(),,()m r x r x r x ⋅⋅⋅。

第6章曲线拟合的最小二乘法6.1 拟合曲线通过观察或测量得到一组离散数据序列，当所得数据比较准确时，可构造插值函数逼近客观存在的函数，构造的原则是要求插值函数通过这些数据点，即。

此时，序列与是相等的。

如果数据序列，含有不可避免的误差（或称“噪音”），如图6.1所示；如果数据序列无法同时满足某特定函数，如图6.2所示，那么，只能要求所做逼近函数最优地靠近样点，即向量与的误差或距离最小。

按与之间误差最小原则作为“最优”标准构造的逼近函数，称为拟合函数。

图6.1 含有“噪声”的数据图6.2 一条直线公路与多个景点插值和拟合是构造逼近函数的两种方法。

插值的目标是要插值函数尽量靠近离散点；拟合的目标是要离散点尽量靠近拟合函数。

向量与之间的误差或距离有各种不同的定义方法。

例如：用各点误差绝对值的和表示：用各点误差按模的最大值表示：用各点误差的平方和表示：或（6.1）其中称为均方误差，由于计算均方误差的最小值的方法容易实现而被广泛采用。

按均方误差达到极小构造拟合曲线的方法称为最小二乘法。

本章主要讲述用最小二乘法构造拟合曲线的方法。

在运筹学、统计学、逼近论和控制论中，最小二乘法都是很重要的求解方法。

例如，它是统计学中估计回归参数的最基本方法。

关于最小二乘法的发明权，在数学史的研究中尚未定论。

有材料表明高斯和勒让德分别独立地提出这种方法。

勒让德是在1805年第一次公开发表关于最小二乘法的论文，这时高斯指出，他早在1795年之前就使用了这种方法。

但数学史研究者只找到了高斯约在1803年之前使用了这种方法的证据。

在实际问题中，怎样由测量的数据设计和确定“最贴近”的拟合曲线？关键在选择适当的拟合曲线类型，有时根据专业知识和工作经验即可确定拟合曲线类型；在对拟合曲线一无所知的情况下，不妨先绘制数据的粗略图形，或许从中观测出拟合曲线的类型；更一般地，对数据进行多种曲线类型的拟合，并计算均方误差，用数学实验的方法找出在最小二乘法意义下的误差最小的拟合函数。