matlab householder变换,数据拟合线性最小二乘法及其应用(householder变换).docx

数据拟合线性最小二乘法及其应用(householder变换)

摘要本文介绍了数据拟合中最小二乘法的基本原理，并根据最小二乘数据拟合方法，针对线性最小二乘拟合及其快速算法(运用HOUSEHOLDER变换及QR分解)进行了分析和研究，并通过数学工具Matlab实现了对相关问题的举例应用，使相关数据拟合理论更加生动易懂. 关键词：数据拟合，最小二乘法， householder变换，应用ABSTRACTThis thesis introduces the basic principle of the least square data fitting. According to the least square data fitting ,the paper studies and analyzes the linear least square fitting using the Householder matrix transformations and QR decompositions of the matrix. Using Matlab makes the theory of data fitting easier to be understood. The method of operation and compose of Matlab is quite easy and the result is accurate. It is very convenient for solving practical problems.Key Words:data fitting least square method householder transformation application第一章引言实验数据的正确处理关系到能否达到实验目的、得出明确结论. 在进行实验数据的处理分析时，数据拟合是经常采用的方法之一. 数据拟合的目标是寻找一条光滑数据，使之在某种准则下最佳地拟合数据. 本文论述了最小二乘数据拟合法的若干原理，介绍了如何对数据应用Matlab软件实现数据拟合，并给出了若干实例帮助理解相关理论及Matlab实现最小二乘问题在数据拟合、参数估计和函数逼近等方面中有着广泛的应用. 第二章最小二乘法2.1最小二乘拟合给出m组观测数据(ti,yi),i=1,2,…,m.可以认为它是某一模型得到的，其中是“真正”的光滑函数，是误差函数，y(x)是经干扰所得到的，因此借用物理学名词，又把成为白噪音过程，,,,而是观测的随机误差. 设是带有n个参数的数据拟合函数，记决定中的参数，使达到极小，这就是最小二乘拟合问题，也常称为最小二乘法或最小二乘逼近，所求得的解称最小二乘解. 2.2 线性最小二乘拟合k次代数多项式是的线性组合由于由k+1个系数唯一确定，所以又称为k+1阶多项式. 一般地，设是x的函数，并且它们是线性无关组，则其线性组合称为n阶广义多项式，取广义多项式为拟合函数的最小二乘拟合问题称为线性最小二乘问题. 对于m组观测数据，计算广义多项式P(x)在处所得的值与之差的平方和为目的是要求出(即确定参数)使达到最小值. 由于是的函数. 因此问题化为求多元函数的极小值. 由多元函数极值的必要条件可得n元线代数方程组其中n元线代数方程组称为线性最小二乘问题的正规方程组. 以为元素的n元线代数方程组的系数行列式不为零时，可以唯一地解出广义多项式的参数特别地，当时，相应的n元线性代数方程组为其中2.3 小结最小二乘问题的解归结为求正则方程组. 由于正则方程系数矩阵的条件数是原矛盾方程组系数条件数的平方，正则方程的病态程度大大增加. 因此，有必要采用有较高数值稳定性的计算方法. 一般使用H-变换，平方根法和Gram-Schmidt正交化法. 这三种方法用于一次性处理时，H-变换所需的运算量较少，而得到最广泛的使用. 平方根法便于形成递推算法而首先被用于最小二乘估计的递推算法中，但与直接求解正则方程的递推算法相比较，运算量有所增加. 第三章快速算法及基本原理3.1H-变换Gauss消去法是用左乘一系列初等下三角阵来约化一个矩阵为上三角阵. 而Householder变换(注释一)是应用非常广泛的另外一种三角化方法. 定义 Householder矩阵是指形式为的矩阵，其中，通常记为H. 根据定义可以看出，Householder矩阵H具有良好的性质：对称性()；正交性()；对称性()；反射性：对任意的应用中经常是，确定一个Householder矩阵H的u并不总是单位向量，把它规范化需要计算，因而要求平方根，这是比较费时的. 为此，用以下定理的方法：定理1：设，令则是一个Householder矩阵. 证明：设，那么，且Householder矩阵的一个关键是能把一个给定向量的若干个指定的分量变为零，