数值计算方法三次样条插值.pptVIP

4．5曲线拟合的最小二乘法插值法是用多项式近似的表示函数,并要求在他们的某些点处的值相拟合.同样也可以用级数的部分和作为函数的近似表达式.无论用那种近似表达式,在实际应用中都要考虑精度,所以我们给出最佳逼近的讨论.第30页，讲稿共69页，2023年5月2日，星期三4.5.1最佳平方逼近定义4.5.1设称为函数在区间[a,b]上的内积.其中为区间[a,b]上的权函数,且满足下面两个条件:第31页，讲稿共69页，2023年5月2日，星期三容易验证,上述定义的函数内积满足一般内积概念中四条基本性质.第32页，讲稿共69页，2023年5月2日，星期三内积的性质第33页，讲稿共69页，2023年5月2日，星期三函数的欧几里得范数定义4.5.2设称为函数f(x)的欧几里得范数,或2范数.第34页，讲稿共69页，2023年5月2日，星期三函数的欧几里得范数性质第35页，讲稿共69页，2023年5月2日，星期三线性相关的函数系定义4.5.3设函数,如果存在一组不全为零的数使成立,则称函数系是线性相关的,否则称是线性无关的.第36页，讲稿共69页，2023年5月2日，星期三线性相关的函数系的判定定理4.5.1函数在区间[a,b]上线性相关的充分必要条件是Gramer行列式第37页，讲稿共69页，2023年5月2日，星期三不难证明在R上线性无关.定理4.5.1的等价说法是:函数系线性无关的充分必要条件是Gramer行列式.第38页，讲稿共69页，2023年5月2日，星期三最佳平方逼近定义4.5.4设函数及函数系且线性无关.记为连续函数空C[a,b]的子空间,如果存在元素满足第39页，讲稿共69页，2023年5月2日，星期三则称为f(x)在上的最佳平方逼近函数.且其中是法方程唯一的一组解.第40页，讲稿共69页，2023年5月2日，星期三令则误差为第41页，讲稿共69页，2023年5月2日，星期三特例取则法方程为其中第42页，讲稿共69页，2023年5月2日，星期三例题例4.5.1设求f(x)在区间[0,1]上的一次最佳平方逼近多项式.解设由于第43页，讲稿共69页，2023年5月2日，星期三故法方程为解得第44页，讲稿共69页，2023年5月2日，星期三平方误差为第45页，讲稿共69页，2023年5月2日，星期三4.5.2对离散数据的曲线拟合最小二乘法曲线拟合问题对于f(x)插值问题,要想提高精度,就要增加节点,因此多项式的次数也就太高,计算量过大,而节点少,多项式的次数低,但误差精度不能保证,为了消除误差干扰,取多一些节点利用最小二乘法确定低次多项式近似表示f(x),这就是曲线拟合问题.第46页，讲稿共69页，2023年5月2日，星期三在科学实验中,得到函数y=f(x)的一组实验数据:,求曲线与实验数据误差在某种度量意义下最小.第47页，讲稿共69页，2023年5月2日，星期三设是[a,b]上一组线性无关的连续函数系,令记误差.为寻求我们常以误差加权平方和最小为度量标准,即第48页，讲稿共69页，2023年5月2日，星期三达到极小值,这里是[a,b]上的权函数.类似前述最佳平方逼近方法,有多元函数极