曲线拟合的数值计算方法实验

曲线拟合的数值计算方法实验

【摘要】实际工作中，变量间未必都有线性关系，如服药后血药浓度与时间的关系；疾病疗效与疗程长短的关系；毒物剂量与致死率的关系等常呈曲线关系。曲线拟合（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 =的加权平方和达到最小，此时所求曲线称作在加权最小二乘意义下对数据的拟合曲线。有许多求解拟合曲线的成功方法，对于线性模型一般通过建立和求解方程组来确定参数，从而求得拟合曲线。至于非线性模型，则要借助求解非线性方程组或用最优化方法求得所需参数才能得到拟合曲线，有时称之为非线性最小二乘拟合。

曲线拟合：贝塞尔曲线与路径转化时的误差。值越大，误差越大；值越小，越精确。

2.最小二乘法拟合：

最小二乘法（又称最小平方法）是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。其他一些优化问题也可通过最小化能量或最大化熵用最小二乘法来表达。

函数曲线为： y=Ax+B

其中系数满足下列的正规方程：

∑∑∑====⎪⎭⎫

⎝⎛+⎪⎭⎫ ⎝⎛N

k k

k N k k N k k y x B x A x 1112

∑∑===+⎪⎭⎫

⎝⎛N

k k

N k k y NB A x 11

3.幂函数拟合：

函数曲线为：

设(){}Ν

k k ，k y x 1=有N 个点，其中横坐标是确定的。最小二乘幂函数拟合曲线

的系数A 为：

)

/()(1

21∑∑===N

k M

k N

k k M

k x y x A 、

4.对数函数拟合：

对数函数（lograrithmic function ）的标准式形式为 )0(ln a Y >+=X X b

b>0时，Y随X增大而增大，先快后慢；b<0时，Y随X增大而减少，先快后慢，见图(c)、(d)。当以Y和lnX绘制的散点图呈直线趋势时，可考虑采用对数函数描述Y与X之间的非线性关系，式中的b和a分别为斜率和截距。

更一般的对数函数

Y=a+bln(X+k)

式中k为一常量，往往未知。

(a)lnY=lna+bX

(b)lnY=lna-bX

(c)Y=a+blnX

(d)Y=a-blnX

5.线性插值：

在代数插值中,为了提高插值多项式对函数的逼近程度一般是增加节点的个数，即提高多项式的次数，但这样做往往不能达到预想的效果。

如下图所示：

=5*(k/3-1) f(x) = 1 / (1 + x2)如果在区间[-5，5]上取7个等距节点:x

k

(k=0,1,2,...,6),由lagrange插值公式可得到f(x)的次L7(x)。如图所示：L7(x)仅在区间的中部能较好的逼近函数f(x), 在其它部位差异较大，而且越接近端点，逼近效果越差。可以证明，当节点无限加密时，Ln(x)也只能在很小的范围内收敛，这一现象称为Runge现象。它表明通过增加节点来提高逼近程度是不适宜的，因而不采用高次多项式插值。

如果我们把以上的点用直线连接起来，显然比L7(x)要更逼近f(x)。这就是分段线性插值。

而在实际应用中通常采用分段低次插值来提高近似程度，比如可用分段线性插值或分段三次埃尔米特插值来逼近已知函数，但它们的总体光滑性较差。为了克服这一缺点，一种全局化的分段插值方法——三次样条插值成为比较理想的工具。

6.三次样条插值：

设(){}N k k k ，y

x 0=有N+1个点，其中b x ...x x x a N =<<<<=210。如果存在N 个三次

多项式x ）S k （，系数为3

2

1

0k ，k ，k ，k ，

，S

，S

，S S 满足如下性质：

3

3,22,1,0,)()()()()(k k k k k k k k x x s x x s x x s s x S x S -+-+-+==

1,....,1,0],,[1-=∈+N k x x x k k

)('')('')(')(')()()(111111111+++++++++====k k k k k k k k k k

k x S x S x S x S x S x S y x S 2,....1,02,....1,02,....1,0,....1,0-=-=-==N k N k N k N k

则成函数S(x)为三次样条函数。

7.端点约束：

紧压样条：存在唯一的三次样条曲线，其一阶导数的边界条件是：

N d b S d a S ==)(',)('0

natural 样条：存在唯一的三次样条曲线，它的自由边界条件是：

0)('',0)(''==b S a S

外推样条：存在唯一的三次样条曲线，其中通过对点x1和x2进行外推得到

()a S ''，同时通过对点X(n-1)和X(N-2)进行外推得到()b S ''。

端点曲率调整：存在唯一的三次样条曲线，其中二阶导数的边界条件()a S ''和

()b S ''是确定的。

抛物线终结样条：存在唯一的三次样条曲线，其中二阶在区间[X 0，X 1]内0x S ≡'''）（，而在[Xn-1,Xn]内0x S ≡'''）（。

三、实验内容