VB计算程序课程设计报告

测试与光电工程学院

课程设计任务书

测控技术与仪器系100813班学号10081329 姓名吴辉

课程名称：用最小二乘法求拟合曲线

课题要求：

利用VB语言编程实现对给定离散点的拟合（不小于10个）的拟合

用最小二乘法求数据的拟合曲线。要求有良好的输入、输出界面，输出应包含直线方程并图形显示拟合效果。完成软件的整体设计。

课题进程：

1）熟悉VB编程语言、最小二乘法算法分析3天

2）编写程序实现以上功能3天

3）软件调试、测试2天

4）撰写课程设计报告2天

指导老师：杨琳瑜

目录

摘要 ------------------------------------------------------------------ 2第一章最小二乘法---------------------------------------------------- 2

1) 理论依据------------------------------------- 错误!未定义书签。

2) 线性拟合分析 -------------------------------- 错误!未定义书签。

3) 非线性拟合分析------------------------------------------------ 5第二章系统设计------------------------------------------------------ 5

1) 采用的软件及开发平台 ----------------------------------------- 5

2) 项目的总体方案------------------------------------------------ 5

3) 项目的详细设计------------------------------------------------ 6第三章设计实现----------------------------------------------------- 10

1) 主要功能模块的具体实现-------------------------------------- 10

2) 主要技术问题或难题的解决方法------------------------------- 10

3) 亮点或创新点的实现 ------------------------------------------ 11第四章结束语 ------------------------------------------------------- 11参考文献------------------------------------------------------------- 12附录 ----------------------------------------------------------------- 13

摘要

最小二乘法最早是由高斯提出的，这是数据处理的一种很有效的统计方法。高斯用这种方法解决了天文学方面的问题，特别是确定了某些行星和彗星的天体轨迹。这类天体的椭圆轨迹由5个参数确定，原则上，只要对它的位置做5次测量就足以确定它的整个轨迹。但由于存在测量误差，由5次测量所确定的运行轨迹极不可靠，相反，要进行多次测量，用最小二乘法消除测量误差，得到有关轨迹参数的更精确的值。最小二乘法近似将几十次甚至上百次的观察所产生的高维空间问题降到了椭圆轨迹模型的五维参数空间。

最小二乘法普遍适用于各个科学领域，它在解决实际问题中发挥了重要的作用。它在生产实践、科学实验及经济活动中均有广泛应用。比如说，我们引入等效时间的概念,根据

Arrhenius 函数和指数函数研究水化热化学反应速率随温度的变化,最后采用最小二乘法回归分析试验数据,确定绝热温升和等效时间的关系式。

第一章 最小二乘法

1) 理论依据

最小二乘的思想就是要使得观测点和估计点的距离的平方和达到最小.这里的“二乘”指的是用平方来度量观测点与估计点的远近（在古汉语中“平方”称为“二乘”），“最小”指的是参数的估计值要保证各个观测点与估计点的距离的平方和达到最小.

对于回归模型y=S(x),若（i i y x ,）（i=1,2,3···m ）为收集到的观测数据，则应该用来估计，这里是（)(,i i x S x ）（i=1,2,3···m ）的估计值。它们之间距离的平方和就是21])([i m

i y x S -∑。

进而最小二乘估计量就是使

2

2

δ

=∑=m i i

2

δ=2

])(*[i m

i i y x S -∑==2

1

)(]

)([min 1

i

m

i

x S y x S i -∑=∈ϕ (1)

达到最小值的参数。

2) 线性拟合分析

给定一组测量数据{（i i y x ,），i=0,1,2,…,m}，基于最小二乘原理，求得变量x 和y 之间的函数关系f(x,A),使它最佳地逼近已知数据。其中A=（n a a a ,...,,10）是一些待定参数。

通常把最小二乘法中的2

2δ都考虑为加权平方和，即

2

2δ=20

))()((i i m

i i y x f x -∑=ω

其中，0)(>=i x ω是[a,b]上的权函数，它表示反应数据（i i y x ,）在实验中所占数据的

比重。

选择参数A 使得加权平方和最小，即求满足

0)(,)

)()((m in

)

)(*)((2

2

>=-=-∑∑==i i

i

m

i i

i

i

m

i i

x y x f x y x f x ωωω (2)

的f*(x)。

要使（2）最小，它转换为求多元函数

∑∑==-=m i n

j i i j j i n x f x a x a a a I 0

20

10])()()[(),,(ϕωΛ，

的极小点),(*

*

1*

0n a a a Λ问题。

由求多远函数极值的必要条件，有

∑∑====-=∂∂m i i k n

j i i j j i k n k x x f x a x a I

00

).,,1,0(0)(])()()[(2Λϕϕω 若记∑==m

i i k i j i k j x x x 0

)()()(),(ϕϕωϕϕ，

则

∑==≡=m

i k i k i i k n k d x x f x f 0),,,1,0()()()(),(Λϕωϕ

可改写为

∑===n

j k j k j

n k d a 0

).,,1,0(),(Λϕϕ

（3）

此方程成为法方程。它也可以写成矩阵形式

⎥

⎥

⎥⎥⎦

⎤⎢

⎢⎢⎢⎣⎡====),(),(),(),(),(),(),(),(),(,),,,(,),,(,

1011

10

1010001010n n n n n n T

n T n G d d d d a a a a d Ga ϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕΛ

M M M Λ

ΛΛΛ，其中

由于0ϕ，1ϕn ϕΛ线性无关，故0≠G ,方程组（3）存在唯一解