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

PAGE

PAGE4

《数值方法》实验报告

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

【摘要】实际工作中，变量间未必都有线性关系，如服药后血药浓度与时间的关系；疾病疗效与疗程长短的关系；毒物剂量与致死率的关系等常呈曲线关系。曲线拟合（curvefitting）是指选择适当的曲线类型来拟合观测数据，并用拟合的\t/_blank曲线方程分析两变量间的关系。曲线直线化是曲线拟合的重要手段之一。对于某些非线性的资料可以通过简单的变量变换使之直线化，这样就可以按\t/_blank最小二乘法原理求出变换后变量的\t/_blank直线方程，在实际工作中常利用此直线方程绘制资料的标准工作曲线，同时根据需要可将此直线方程还原为\t/_blank曲线方程，实现对资料的曲线拟合。常用的曲线拟合有最小二乘法拟合、幂函数拟合、对数函数拟合、线性插值、三次样条插值、端点约束。

关键词曲线拟合、最小二乘法拟合、幂函数拟合、对数函数拟合、线性插值、三次样条插值、端点约束

一、实验目的

掌握曲线拟合方式及其常用函数指数函数、幂函数、对数函数的拟合。

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

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

二、实验原理

曲线拟合

曲线拟合是平面上离散点组所表示的坐标之间的函数关系的一种数据处理方法。用解析表达式逼近\t/_blank离散数据的一种方法。在\t/_blank科学实验或社会活动中，通过实验或观测得到量x与y的一组数据对（Xi，Yi)（i=1，2，...m），其中各Xi是彼此不同的。人们希望用一类与数据的背景材料规律相适应的解析表达式，y=f(x，c）来反映量x与y之间的依赖关系，即在一定意义下“最佳”地逼近或拟合已知数据。f(x，c)常称作拟合模型，式中c=(c1，c2，…cn)是一些待定参数。当c在f中\t/_blank线性出现时，称为线性模型，否则称为\t/_blank非线性模型。有许多衡量拟合优度的标准，最常用的一种做法是选择参数c使得拟合模型与实际\t/_blank观测值在各点的\t/_blank残差(或\t/_blank离差)的\t/_blank加权平方和达到最小，此时所求曲线称作在加权最小二乘意义下对数据的拟合曲线。有许多求解拟合曲线的成功方法，对于线性模型一般通过建立和求解\t/_blank方程组来确定参数，从而求得拟合曲线。至于\t/_blank非线性模型，则要借助求解非线性方程组或用最优化方法求得所需参数才能得到拟合曲线，有时称之为非线性\t/_blank最小二乘拟合。

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

最小二乘法拟合：

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

函数曲线为： y=Ax+B

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

幂函数拟合：

函数曲线为：

设有N个点，其中横坐标是确定的。最小二乘幂函数拟合曲线的系数A为：

、

对数函数拟合：

对数函数（lograrithmicfunction）的标准式形式为

P2021

胡克定律指出F=kx，其中F是拉伸弹簧的拉力（单位为盎司），x为拉伸的长度（单位为英寸）。根据下列试验数据，求解拉伸常量k的近似值。

xk

Fk

0.2

3.6

0.4

7.3

0.6

10.9

0.8

14.5

1.0

18.2

Xk

Fk

0.2

5.3

0.4

10.6

0.6

15.9

0.8

21.2

1.0

26.4

（a） （b）

P2151

洛杉矶（美国城市）郊区11月8日的温度记录入下表所示，其中共有24个数据点。

（a）根据例5.5中的处理过程（使用fmins命令），对给定的数据求解最小二乘曲线

。

（b）求

（c）在同一坐标系下画出这些点集和(a)得出的最小二乘曲线。

时间，p.m.

温度

时间，a.m.

温度

1

66

1

58

2

6