最小二乘曲线拟合及MATLAB实现_测绘专业论文(51页).doc
《最小二乘曲线拟合及MATLAB实现_测绘专业论文(51页).doc》由会员分享,可在线阅读,更多相关《最小二乘曲线拟合及MATLAB实现_测绘专业论文(51页).doc(49页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、-最小二乘曲线拟合及MATLAB实现_测绘专业论文-第 42 页内蒙古科技大学本科生毕业设计说明书(毕业论文)题 目:最小二乘曲线拟合及MATLAB实现学生姓名:李亚学 号:0972143230专 业:测绘工程班 级:2009测绘2班指导教师:孙同贺 讲师最小二乘曲线拟合及MATLAB实现摘 要介绍曲线拟合的基本理论,对最小二乘原理进行了全方位的理论阐述,同时也阐述了曲线拟合的基本原理及多项式曲线拟合模型的建立。详细的解答了曲线拟合中的最小二乘法,并介绍了部分的正交最小二乘法理论。重点讲解多项式拟合的具体步骤,同时也介绍了非线性方程的最小二乘拟合,在建立理论的基础上对最小二乘曲线拟合法的MAT
2、LAB实现方法进行研究,利用MATLAB2012b的平台对测量数据进行最小二乘曲线拟合,介绍MATLAB的具体构造和曲线拟合工具。利用MATLAB中的ployfit函数对实测数据进行多项式曲线拟合,并给出曲线拟合MATLAB实现的源程序,给出拟合曲线,并评定拟合的精度证明该方法是行之有效的。关键词:最小二乘法,曲线拟合,MATLAB,测量数据Curve Fitting in Least-Square Methodand Its Realization with MatlabAbstract To introduce the basic theory of curve fitting and d
3、iscuss the least squares principle in this paper, whats more, we also discuss the basic principle of curve fitting and the establishment of polynomial curve fitting model. Meanwhile, we also introduce the least-square method of curve fitting in detail and part of the theory of orthogonal least squar
4、e method. We mainly discuss the specific steps of polynomial fitting, and also introduces the nonlinear equation of the least squares fitting at the same time, which established on the theory of least squares curve fitting in MATLAB in order to realize the method to do research. Using MATLAB2012b pl
5、atform to achieve the goal of measuring data and introducing the special structure of MATLAB and curve fitting tool. We can use ployfit function in MATLAB to polynomial curve fitting of experimental data, and get the MATLAB source program about curve fitting and the fitting curve. Finally, we need t
6、o prove the method of assessing the precision of the fitting is effective.Key words: least square method; curve fitting; MATLAB, metrical data最小二乘曲线拟合及MATLAB实现I摘 要ICURVE FITTING IN LEAST-SQUARE METHOD AND ITS REALIZATION WITH MATLABIIABSTRACTII第一章 引 言11.1研究背景11.1.1 历史理论原理11.1.2 现代研究11.2 问题定义21.2.1 曲
7、线拟合的思想21.2.2 多项式拟合31.2.3 利用Matlab的polyfit函数进行多项式拟合31.3 论文结构3第二章 数据曲线拟合42.1测量数据42.2拟合模型42.3最小二乘原理52.3.1最小二乘法52.3.2最小二乘估计与极大似然估计72.4数据拟合92.4.1曲线拟合理论92.4.2最小二乘法线性拟合原理102.4.3最小二乘非线性拟合122.4.4正交多项式132.4.5正交最小二乘曲线拟合152.5曲线拟合精度评定17第三章MATLAB193.1MATLAB概述193.1.1MATLAB简介193.1.2MATLAB的主要组成部分213.2MATLAB2012b的运行简
8、介233.2.1启动和退出MATLAB2012b233.2.2MATLAB2012b桌面系统243.2.3MATLAB函数调用系统263.2.4MATLAB2012b的帮助系统273.2.5附件管理系统283.2.6数据交换系统283.2.7MATLAB 中的其他系统293.3最小二乘曲线拟合法的MATLAB实现30第四章 最小二乘法曲线拟合的MATLAB实现324.1 使用polyfit函数实现多项式拟合324.2 二次多项式的曲线拟合334.3三次多项式的曲线拟合344.4 四次多项式曲线拟合354.5数据处理和精度评定36第五章 总结40参考文献41附录1:43MATLAB语言编程源代码
9、43附录2:45各次拟合的拟合曲线方程45致谢46外文翻译47外文部分47翻译部分54第一章 引 言1.1研究背景1.1.1 历史理论原理Weierstrass第一逼近定理1 对任意函数和任意给定的,都存在n次代数多项式,满足 (1-1-1)Bernstein多项式(bernstein polynomial)1前苏联数学家Bernstein曾经给出这样的多项式序列: (1-1-2)在整体上一致逼近,但它的收敛缓慢,要达到一定的精度,则n要取很大,计算量大,所以研究如何在给定的精度下,对进行整体逼近,成为逼近论中的一个重要问题。1.1.2 现代研究曲线拟合问题是诸多试验和工程实际中广泛应用的数据
10、处理方法。试验数据的正确处理,关系到是否能达到试验目的,得出明确结论。传统的数据处理方法,很难得到一条很好适应所有点的曲线,同时也无法估计所得曲线的精度,由此所确定的特征值就可能有较大的误差,且没有建立起由这些点构成曲线的数学模型,直接影响利用数学方法进行解析分析。在进行试验数据的分析时,通常可采用曲线拟合法寻找一条光滑曲线,曲线在某种准则下最佳的拟合数据。测量工作中,通常根据测定的一系列坐标点,选取一定的数学模型拟合直线、二次曲线或者其他高次曲线。拟合的目的是根据测量点寻求曲线的特征,求解曲线的相关参数,为工程建设管理提供必要的基础信息。如在既有铁路工程、又有公路工程测量中,通常根据一系列的
11、测量点和线路工程的特点求解线路工程的线性特征,为线路维护养护、二线工程建设、行车安全分析等提供必要的基础信息。在GIS数据获取中,通常根据一系列的实际测量点或者是地图数字化点拟合道路、水系、等高线、等曲线。这类问题的做法通常是根据线形的特点选取一定的数学模型,以待求的线形参数作为未知参数,以测点的纵坐标或者横坐标为观测值,采用最小二乘法处理。在测量中获取的数据均为随机数据,它们是由一些离散的数据组成,单就获得的原始数据本身来说根本反映不出事物的本质。如何从这些离散的数据中找出观测数据的变化规律?在实际中传统的数据处理方法,很难得到一条很好地适应所有点的曲线,同时也无法估计所得曲线的精度,且没有
12、建立起由这些点构成曲线的数学模型,直接影响到利用数学方法进行解析分析。用Matlab进行数据拟合可以形象直观地发现所有数据体现出来的规律性。在进行分析时,通常可采用曲线拟合法,曲线拟合法的目的是寻找一条光滑曲线,即对观测的几个变量进行多次观测从而求出反映变量之间的相对函数关系,它在某种准则下最佳的拟合数据。1.2 问题定义本文介绍最小二乘曲线拟合法的基本原理,就其MATLAB的实现方法进行研究,给出曲线拟合MATLAB 的实现方法进行研究,给出曲线拟合MATLAB实现的源程序,并进行仿真测试,对测试误差进行分析。1.2.1 曲线拟合的思想如果不要求所构造的函数精确的通过所有由离散数据所确定的离
13、散点,而只要求是相对与同一函数类H中的其他函数而言达到最优的。即我们希望找到一条曲线,既能反映给定数据的一般趋势又不至于出现局部较大波动。在这种逼近方式下,只要构造的近似函数与被逼近函数在区间a,b上的偏差满足某种要求即可。1.2.2 多项式拟合有时所给的数据点的分布并不一定近似的呈一条直线,这时若仍用直线拟合显然是不合适的,对于这种情况可以考虑用多项式拟合。多项式方程的一般形式是: (1-2-1)解出多项式系数,可得到函数模型。1.2.3 利用Matlab的polyfit函数进行多项式拟合在Matlab中曲线拟合的形式非常简单,他的形式是:,该拟合函数的结果将保证在数据点上的拟合值与数据值之
14、差的平方和最小,满足最小二乘法则标准的最小二乘曲线拟合。 1.3 论文结构 本文主要分为五章,第一章介绍本文的主旨和需要解决的问题的介绍,第二章介绍最小二乘法曲线拟合的基本理论和具体步骤,第三章通过MATLAB2012b的平台介绍MATLAB实现最小二乘曲线拟合的具体方法和步骤,第四章利用MATLAB的ployfit函数对一组矿山沉陷数据进行多项式曲线拟合,并对多项式拟合的精度进行分析,最后第五章对全文进行一个总结。第二章 数据曲线拟合2.1测量数据测量数据或观测数据是指用一定的仪器、工具、传感器或其他手段获取的反映地球与其他实体的空间分布有关信息的数据。观测数据可以是直接测量的结果,也可以是
15、经过某种变换后的结果。任何观测数据总是包含信息和干扰两部分,采集数据就是为了获取有用的信息,干扰也称为误差,是除了信息以外的部分,要设法予以排除或减弱其影响4。2.2拟合模型拟合模型是测量平差中常遇到的一种特殊的函数模型。拟合模型是一种函数逼近型或是统计回归模型。用一个函数去逼近所给定的一组数据,或者利用变量与变量之间统计相关性质给定的回归模型都属于拟合模型4。拟合模型误差方程的组成举例:1.在地图数字化中,已知圆上m个点的数字化观测值(i=1,2,m),设为等权独立观测试求该圆的曲线方程。由于数字化观测值有误差,m个点并不在同一圆线上,需要在这些观测点上拟合一条最佳圆曲线,这就是拟合模型问题
16、。圆曲线的参数方程以平差值表示为 (2-2-1)公式中为圆心坐标平差值,和分别为半径和矢径方位角的平差值,它们为平差的未知参数,故此例n=2m,t=3+m。令 =+,=+.将上式线性化,最后得误差方程为 . (2-2-2)式中,2.在摄影测量学中,数字高程模型、GPS水准的高程异常拟合模型等,常采用多项式拟合模型。已知m个点的数据是(i=1,2,m),其中是点的高程或高程异常(GPS水准拟合模型),为点的坐标,视为无误差,并认为Z是坐标的函数,即可取拟合函数为, (2-2-3)式中=,未知参数为.为常数,则其误差方程为 . (2-2-4)2.3最小二乘原理2.3.1最小二乘法在生产实践中,经常
17、会遇到利用一组观测数据来估计某些未知参数的问题。例如,一个做匀速运动的质点在时刻的位置是,可以用如下的线性函数来表达描述: (2-3-1)式中,是质点在时刻的初始位置,是平均速度,它们是待估计的未知数参数,可见这类问题为线性参数的估计问题。对于这一问题,如果观测没有误差,则只要在两个不同时刻观测出质点的相应位置,由上式分别建立两个方程,就可以解出的值。但是在实际的观测时,考虑到观测值带有偶然误差,所以总是作多余观测。在这种情况下,为了求得,就需要在不同时刻来测定其位置,得出一组观测值,这时,由上式可以得到 (2-3-2)若令则 (2-3-3)这就是间接平差的模型4。如果我们将对应的用图解表示,
18、从图2.1看出由于存在观测误差的缘故,由观测数据绘出的点观测点,描绘不成直线,而有些摆动。图2.1根据观测点确定直线这里就产生了一个问题:用什么准则来对参数进行估计,从而使估计直线“最佳”的拟合于各观测值点。通常的做法有以下几种: (1) (2-3-4)(2) (2-3-5)(3) (2-3-6)其中第一种较复杂,第二种不可导,求解困难,所以目前采用较多的方法是第三种方法,这种方法就叫做最小二乘法。所谓的最小二乘原理,就是要在满足 (2-3-7)的条件下解出参数的估值,也可以表达为 (2-3-8)式中,表示未知参数的估计向量,在上述例子中,。满足上式的估计值称为的最小二乘估计,这种求估计量的方
19、法就叫做最小二乘法1。2.3.2最小二乘估计与极大似然估计 测量中的观测值是服从正态分布的随机变量,最小二乘原理可用数理统计中最大似然估计来解释,两种估计准则的估值相同。设观测向量为L,L为随机正态向量,其数学期望和方差分别为由极大似然估计准则知道,其似然函数(即L的正态密度函数)为 (2-3-9)按最大似然估计的要求,应选取能使取得极大值的作为的估计量,考虑到,为的估计量也就是以改正数V作为真误差的估计量。由于上式中右边第一项为常量第二项前为负号,所以只有当第二项取得极小值时,似然函数才能取得极大值。因此由极大似然估计求得的V值必须满足条件考虑到为常量,则上式等价于 (2-3-10)此方程即
20、为最小二乘原理。由此可见,当观测值为正态随机变量时,最小二乘估计可由最大似然估计导出,由以上两个准则出发,平差结果完全一致。最小二乘原理中的P阵,称为权阵,定义是。设为独立观测值,其权为,则有式中,为的权倒数或协因数,权阵及协因数阵为如果为相关的观测值,则有协因数Q与协方差D统计含义相同,数值的表达式形式上仅差一个常量,如果=1,则D=Q。因为权阵由于为的权倒数,但是,所以权阵P中的主对角线不具有权的意义,P仅表示,但在运算中起着权的作用。特别时,当观测为同精度观测时,P=I,则最小二乘原理是 (2-3-11)2.4数据拟合2.4.1曲线拟合理论在测量学上,常常使用一组测定的数据,i=1,2,
21、n,求得一个近似的函数关系y=f。由于y值来自观测或者实验,数据不可避免地带有一定程度的误差,因此不能像插值那样要求曲线严格通过数据点,只能是y=最优地靠近这些数据点,这样,在某种意义下的偏差为最小,部分抵消数据误差,进而反应数据的一般趋势。假定有n对实验数据,其中(=1,2,n)。设由这些点得到的数据关系为。在一般情况下,有线性模型= (2-4-1)假设令,则有非线性模型为 (2-4-2)利用线性模型与非线性模型,基于最小二乘法,可以计算出经验公式和参数。2.4.2最小二乘法线性拟合原理前面我们已经介绍了最小二乘法,现在我们就最小二乘法拟合直线和曲线(多项式)做一个详细的原理分析,着重介绍多
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 最小 曲线拟合 MATLAB 实现 测绘 专业 论文 51
限制150内