数值分析课程设计题目非线性曲线拟合最小二乘法学生杨国连.DOCVIP

《数值分析》课程设计 题 目 非线性曲线拟合最小二乘法 学 生 杨国连 杨明洁 郑露 指导教师 郭阁阳 天津工程师范学院 课 程 设 计 任 务 书 数理 系 数学0701 班 学生 杨国连 杨明洁 郑露 课程设计课题： 非线性曲线拟合最小二乘法 一、课程设计工作日自 2009 年 6 月 22 日至 2009 年 6 月 28 日 二、同组学生： 杨国连 杨明洁 郑露 三、课程设计任务要求（包括课题来源、类型、目的和意义、基本要求、完成时间、主要参考资料等）： 【来源与意义】 本课题来源于教材第三章非线性曲线拟合最小二乘法，目的是通过实例，考察用迭代法求解非线性最小二乘拟合曲线。另一方面重在考查学生对用计算机与Matlab解决实际问题的能力。 【基本要求】 要求自编程序；掌握编程思想，学会一门编程语言；报告要有较强的理论分析；有较强说服力的数据表或图像；对结果进行分析；给出相应结论；鼓励创新； 【参考资料】 数值分析，李庆扬，王能超，易大义，2001，清华大学出版社（第四版）。 数值方法，关治，陆金甫，2006，清华大学出版社。 数值分析与实验学习指导，蔡大用，2001，清华大学出版社。 数值分析与实验，薛毅，2005，北京工业大学出版社. 指导教师签字： 教研室主任签字： 天津工程师范学院 课 程 设 计 评 审 表 数理 系 数学0701 班 学生 杨国连 杨明洁 郑露 设计任务完成情况及指导教师评语 答辩情况 评定成绩 成绩： 指导教师签字： 日期： 教研室主任： 主任签字： 日期： 日期： 一、问题提出 设数据（）,(i=0,1,2,3,4).由表3-1给出，表中第四行为,可以看出数学模型为,用最小二乘法确定a及b。 i 0 1 2 3 4 1.00 1.25 1.50 1.75 2.00 5.10 5.79 6.53 7.45 8.46 1.629 1.756 1.876 2.008 2.135 二、理论基础 根据最小二乘拟合的定义:在函数的最佳平方逼近中，如果f(x)只在一组离散点集{,i=0,1,…,m}，上给定，这就是科学实验中经常见到的实验数据{（）, i=0,1,…,m}的曲线拟合，这里，i=0,1,…,m,要求一个函数与所给数据{（）,i=0,1,…,m}拟合，若记误差,i=0,1,…,m,,设是上线性无关函数族，在中找一函数，使误差平方和 , 这里 (nm) 这就是一般的最小二乘逼近，用几何语言来说，就称为曲线拟合的最小二乘法。 在建模的过程中应用到了求和命令（sum）、求偏导命令(diff)、化简函数命令（simple）、用迭代方法解二元非线性方程组的命令（fsolve），画图命令（plot）等。 三、实验内容 用最小二乘法求拟合曲线时，首先要确定S(x)的形式。这不单纯是数学问题，还与所研究问题的运动规律及所得观测数据（）有关；通常要从问题的运动规律及给定数据描图，确定S(x)的形式，并通过实际计算选出较好的结果。S(x)的一般表达式为线性形式，若是k次多项式，S(x)就是n次多项式，为了使问题的提法更有一般性，通常在最小二乘法中都考虑为加权平方和 . 这里是[a,b]上的权函数，它表示不同点处的数据比重不同。 用最小二乘法求拟合曲线的问题，就是求形如S(x)的一个函数，使取得最小。它转化为求多元函数 的极小点问题。再由求多元函数极值的必要条件，有 (k=0,1,…,n) 此题中假设，由已知所给数据点（）画出图形，根据离散点的位置观察出它们所拟合的曲线图形应类似于指数函数的曲线图形，故设拟合曲线的函数为。 本题编程过程中，令f=y,z1=a,z2=b, 令拟合曲线中对应的函数值与的差的平方和为J,即J=sum(fy.^2);分别求J关于z1,z2的偏导，简化后并令其分别为0得一关于z1,z2的二元非线性方程组，最后利用fsolve命令求得z1,z2的值分别为 z1=3.0751 z2=0.5052 故得到拟