数值分析函数逼近和曲线拟合
（5）设{ pn( x)}0 为以( x)为权函数的[a,b]上的正交多项式
序列. 则pn( x)(n 1)的n个根都是在(a,b)内的单重实根;并有
pn1( x) (n x n ) pn( x) n pn1( x), n 1,2，, (2.4)
其中
p1( x)
0，n
an1 an
,n
an1 an
( 三角不等式)
则称 || || 为线性空间S上的范数，S与 || || 一起称为赋范 线性空间，记为X .
例如，对Rn上的向量x ( x1,, xn )T，有 三种常用范数：
|| x || max | xi |， 称为 范数或最大范数，
1 i n n
|| x ||1 | xi ，| 称为1 范数,
正交化手续立得正交正交多项式序列 :
p0( x) 0,
pn( x)
xn
n1( xn,
j0 ( pj,
p p
j j
) )
p
j
,
n 1,2,.
(2.3)
性质：
(1) pn( x)的首项系数为1. (2)Qn( x) Hn均可表为p0( x), p1( x),, pn( x)的线性组合. (3)当k j时，( p j , pk ) 0,且pk ( x)与任一次数小于k的多
若令x cos，则Tn( x) cos(n ),0 .
T0( x) cos(0) 1, T1( x) cos(arccos x) x, T2( x) cos(2arccos x) 2x2 1, T3( x) 4x3 3x,
切比雪夫多项式的性质:
(1) 递推关系
TTn0(1x( x)
(2.2)
例如，三角函数族 1,cos x,sin x,cos 2x,sin 2x,,
为[ , ]上的正交函数族, (1,1) 2 ,(cos kx,cos kx) (sin kx,sin kx) ,其他内积 0.
定义6 设pn( x)是[a,b]上首项系数an 0的n次多项式, ( x)
n
i
i 1
xi2
1/2
.
若x, y Cn，则定义加权内积
n
( x, y) i xi yi .
i 1
可以有限或
定义4 设( x)是区间[a,b]上的非负无限函区数间, 如果满足条件
(1) ab xk ( x)dx存在, k 0,1,2,; (2) 对于[a,b]上的非负连续函数g( x),若abg( x)( x)dx 0,
( xpn, pn ) , ( pn , pn )
n
an1an1
a
2 n
(
( pn , pn 1 ,
pn ) . pn1 )
二、勒让德多项式
区间[1,1]上带权(x) 1的正交多项式
Pn
( x)
1 2n n!
dn dx n
[(x2
1)n ],
(n 0,1,2,)
称为n次 Legendre多项式 .
例1 考察Rn与Cn的内积和范数.
设x ( x1,, xn )T , y ( y1,, yn )T Rn，则定义
内积
( x,
y)
n
xi
i 1
yi；范数
||
x
ቤተ መጻሕፍቲ ባይዱ
||2
n
xi2
i 1
1/2
.
若给定i 0(i 1,, n)为权系数,则定义
内积
( x,
y)
n
i xi
i 1
yi；范数
||
x
||2
称为Gram矩阵，则G非奇异的充要条件是u1, u2,, un线性
无关.
证明：1) G非奇异 以G为系数矩阵的齐次线性方程组
n
n
( ju j , uk ) (u j , uk ) j 0,
j1
j1
只有零解。
k 1,,n.
n
n
n
2) juj 0 ( juj , juj ) 0
三、内积与内积空间
R n中向量x及y定义内积 : (x, y) x1 y1 xn yn .
定义3 设X是数域K(R或C)上的线性空间，对u,v X， 有K中一个数与之对应，记为( u, v )，并满足条件：
(1) (u,v) (v,u), u,v X;
(2) (u,v) (u,v), R;
项式正交.
（4）有递推关系
pn1( x) ( x n ) pn( x) n pn1( x), n 0,1,, (2.4)
其中 p0( x) 1，p1( x) 0，
n ( xpn, pn ) /( pn, pn ), n ( pn , pn ) /( pn1, pn1)，n 1,2,,
j1
j1
j1
n
( ju j ,uk ) 0, k 1,,n.
j1
G非奇异 u1, u2,, un线性无关(反证法)；反之亦然.
在内积空间X上可以由内积导出一种范数,即对u X ,记
|| u || (u,u)，
(1.10)
易证它满足范数定义的正定性和齐次性, 而三角不等式由
Cauchy Schwarz不等式得出.
1
1
x 2 Tm
(
x
)Tn
(
x
)dx
/
2, ,
m n 0, m n 0.
则在[a,b]上g( x) 0；
就称( x)为[a,b]上的权函数.
例2 设f ( x), g( x) C[a,b], ( x)为[a,b]上的权函数,则可
定义内积
( f , g) ab( x) f ( x)g( x)dx.
1,( f , g) ab f (x)g(x)dx.
容易验证内积定义中的四个性质，并导出范数
||
f ( x) ||2
ab( x)
f
2( x)dx
1/ 2
.
||
f
(x) ||2
ab
f
2
(
x)dx
1/
2
.
设0,,n C[a,b],则Gram矩阵为
G G(0 ,,n )
(0 ,0 ) (0 ,1 )
(1
,0
)
(1,1 )
(n ,0 ) (n ,1 )
(0 ,n )
(1
,
n
)
(n ,n )
根据定理3，0,,n线性无关 det(G) 0.
§2 正交多项式
一、正交函数族与正交多项式
定义5 若f ( x), g( x) C[a,b], ( x)为[a,b]上的权函数, 且
( f , g) ab( x) f ( x)g( x)dx 0，
(2.1)
则称f ( x)与g( x)在[a,b]上带权ρ(x)正交 .
, 1
m n, m n.
(2.7)
(2) 奇偶性 Pn( x) (1)n Pn( x).
(2.8)
(3) Pn( x)在(1,1)内部有n个互异的实零点.
(4) 递推关系
Pn1(
P0( x) 1,
x)
2n 1 n1
xPn
(
P1( x) x,
x)
n
n
1
Pn1
(
x),
(n
1,2,)
(2.9)
1, )2
T1( x) xTn( x)
x, Tn1(
x).
(2.11)
Tn( x)的最高次幂xn的系数为2n1,(n 1).
事实上，只需由
cos(n 1) 2cos cos n cos(n 1) , n 1.
代入 x cos , 即得递推关系式.
(2) 正交性
0, m n,
11
i 1
1
||
x
||2
n
xi2
i 1
2
，
称为2 范数.
类似地，对C[a,b]上的f ( x)，可定义三种常用范数：
|| f || max | f ( x) |， 称为 范数，
a xb
|| f ||1 ab| f ( x) | dx， 称为1 范数,
1
|| f ||2 ab f 2( x)dx 2， 称为2 范数.
第3章 函数逼近与曲线拟合
§1 函数逼近的基本概念
一、函数逼近与函数空间
在数值计算中经常要计算函数值，当函数只 在有限点集上给定函数值，要在包含该点集的区 间上用公式给出函数的简单表达式，这些都涉及 在区间[a,b]上用简单函数逼近已知复杂函数的问 题，这就是函数逼近问题。本章讨论的函数逼近， 是指“对函数类A中给定的函数f(x),记作f(x)∈A， 要求在另一类简单的便于计算的函数类B中求函 数p(x)∈B,使p(x)与f(x)的误差在某种度量意义下 最小。
如果存在不全为零的数1,,n P,使得
1x1 n xn 0，
（1 .1）
则称x1,, xn线性相关. 否则,若(1.1)只对1 n 0成
立，则称x1,, xn线性无关.
定理1(Weierstrass) 如果f (x) C[a,b], 那么 0,
多项式p(x),使得
max | f (x) p(x) | , 在a, b上一致成立。
为[a,b]上的权函数, 若多项式序列{ pn( x)}0 ,满足正交性
(2.2)，则称{ pn( x)}0 为以( x)为权函数的[a,b]上的正交 多项式序列. 称pn( x)为以( x)为权函数的[a,b]上的n次正
交多项式.
只要给定[a,b]上的权函数( x), 由{1, x, xn,}利用逐个
(3) (u v, w) (u,w) (v,w), u,v,w X; (4) (u,u) 0,当且仅当u 0时，(u,u) 0. 则称(u,v)为X上的u与v的内积. 定义了内积的线性空间称 为内积空间. (v,u)为(u,v)的共轭，当K R时 (v,u) (u,v).
常用函数的逼近和曲线拟合
常用函数的逼近和曲线拟合在数学中，函数逼近和曲线拟合都是常见的问题。
函数逼近是指找到一个已知函数，尽可能地接近另一个函数。
而曲线拟合则是给定一组数据点，找到一条曲线来描述这些数据点的分布。
本文将讨论常用的函数逼近和曲线拟合方法。
一、函数逼近1. 插值法插值法是最简单的函数逼近方法之一。
它的基本思想是：给定一组已知点，通过构造一个多项式，使得该多项式在这些点处的函数值与已知函数值相等。
插值法的优点是精度高，缺点是易产生龙格现象。
常用的插值多项式有拉格朗日插值多项式和牛顿插值多项式。
拉格朗日插值多项式的形式为：$f(x)=\sum_{i=0}^{n}y_{i}\prod_{j=i,j\neq i}^{n}\frac{x-x_{j}}{x_{i}-x_{j}}$其中，$x_{i}$是已知点的横坐标，$y_{i}$是已知点的纵坐标，$n$是已知点的数量。
牛顿插值多项式的形式为：$f(x)=\sum_{i=0}^{n}f[x_{0},x_{1},...,x_{i}]\prod_{j=0}^{i-1}(x-x_{j})$其中，$f[x_{0},x_{1},...,x_{i}]$是已知点$(x_{0},y_{0}),(x_{1},y_{1}),...,(x_{i},y_{i})$的差商。
2. 最小二乘法最小二乘法是一种常用的函数逼近方法。
它的基本思想是：给定一组数据点，找到一个函数，在这些数据点上的误差平方和最小。
通常采用线性模型，例如多项式模型、指数模型等。
最小二乘法的优点是适用性广泛，缺点是对于非线性模型要求比较高。
最小二乘法的一般形式为：$F(x)=\sum_{i=0}^{n}a_{i}\varphi_{i}(x)$其中，$a_{i}$是待求的系数，$\varphi_{i}(x)$是一组已知的基函数，$n$是基函数的数量。
最小二乘法的目标是使得$\sum_{i=1}^{m}[f(x_{i})-F(x_{i})]^{2}$最小，其中$m$是数据点的数量。
数值分析Ch3函数逼近与曲线拟合
正交，这就需要引进范数与赋范线性空间，内积
3．1 函数逼近的基本概念
• 定义 设集合 S 是数域 P 上的线性空间，元 素 x1 , x2 , , xn S ，若存在不全为零的数 1 , 2 , , n P ，使得 1 x1 2 x2 n xn 0 则称 x1 , x2 , , xn 线性相关，否则，若仅对
数 值 分 析
Computational Method
Chapter 3 函数逼近
第三章 函数逼近与曲线拟合 设函数 y f x 的离散数据（有误差）为
x y
，
x0 y0
x1 y1
x2 y2
xn yn
希望找到简单函数 Px 整体上有 是某度量， 0 是指定精度。
f x Px
1 x1
2 x2 x 2 , 1 1 1 , 1 x , x , 3 2 2 3 x3 3 1 1 2 , 2 1 , 1
xn , 1 xn , 2 xn , n1 1 2 n1 n xn 1 , 1 2 , 2 n1 , n1 k 1 xk , i i ( k 1,2,, n) 简写为: k x xk i 1 i , i
，
x
2
。
（连续） f x Ca, b
b
常见范数：
f x 1 f x dx • 1范数： a ，
• 2-范数：
f x 2
2 f x dx a b
1 2
f x max f x • 范数： ， a ,b
曲线拟合实验报告[优秀范文5篇]
![曲线拟合实验报告[优秀范文5篇]](https://360docs-1393072327.cos.ap-beijing.myqcloud.com/ed56c4cd0342a8956bec0975f46527d3240ca6e1.webp)
曲线拟合实验报告[优秀范文5篇]第一篇：曲线拟合实验报告数值分析课程设计报告学生姓名学生学号所在班级指导教师一、课程设计名称函数逼近与曲线拟合二、课程设计目的及要求实验目的: ⑴学会用最小二乘法求拟合数据的多项式,并应用算法于实际问题。
⑵学会基本的矩阵运算,注意点乘与叉乘的区别。
实验要求: ⑴编写程序用最小二乘法求拟合数据的多项式,并求平方误差,做出离散函数与拟合函数的图形;⑵用MATLAB 的内部函数polyfit 求解上面最小二乘法曲线拟合多项式的系数及平方误差,并用MATLAB的内部函数plot作出其图形,并与(1)结果进行比较。
三、课程设计中的算法描述用最小二乘法多项式曲线拟合,根据给定的数据点,并不要求这条曲线精确的经过这些点,而就是拟合曲线无限逼近离散点所形成的数据曲线。
思路分析 : 从整体上考虑近似函数)(x p 同所给数据点）（i iy x , 误差i i iy x p r -=)(的大小,常用的方法有三种:一就是误差i i iy x p r -=)(绝对值的最大值im ir≤≤ 0max ,即误差向量的无穷范数;二就是误差绝对值的与∑=miir0,即误差向量的 1成绩评定范数;三就是误差平方与∑=miir02的算术平方根,即类似于误差向量的 2 范数。
前两种方法简单、自然,但不便于微分运算,后一种方法相当于考虑 2 范数的平方,此次采用第三种误差分析方案。
算法的具体推导过程: 1、设拟合多项式为:2、给点到这条曲线的距离之与,即偏差平方与:3、为了求得到符合条件的 a 的值,对等式右边求偏导数,因而我们得到了:4、将等式左边进行一次简化,然后应该可以得到下面的等式5、把这些等式表示成矩阵的形式,就可以得到下面的矩阵:⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∑∑∑∑∑∑∑∑∑∑∑=====+==+====niininiiknikinikinikinikiniiniinikiniiyyyaax x xx x xx x11i11012111111211 1an MMΛM O M MΛΛ 6.将这个范德蒙得矩阵化简后得到⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡n kkn nkkyyyaaax xx xx x M MΛM O M MΛΛ21102 21 1111 7、因为 Y A X = * ,那么 X Y A / = ,计算得到系数矩阵,同时就得到了拟合曲线。
实验二函数逼近与曲线拟合
《数值分析》课程设计实验报告实验二函数逼近与曲线拟合一、问题提出从随机的数据中找出其规律性，给出其近似表达式的问题，在生产实践和科学实验中大量存在，通常利用数据的最小二乘法求得拟合曲线。
在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量与时间t 的拟合曲线。
二、实验步骤先写出线性最小二乘法的M文件function c=lspoly(x,y,m)% x是数据点的横坐标组成的向量，y是纵坐标组成的向量% m是要构成的多项式的次数，c是多项式由高到低次的系数所组成的向量n=length(x);b=zeros(1:m+1);f=zeros(n,m+1);for k=1:m+1f(:,k)=x.^(k-1);enda=f'*f;b=f'*y';c=a\b;c=flipud(c);方法一：近似解析表达式为：y（t）=a1t+a2t2+a3t3第二步在命令窗口输入：lspoly([0,5,10,15,20,25,30,35,40,45,50,55],[0,1.27,2.16,2.86,3.44 ,3.87,4.15,4.37,4.51,4.58,4.02,4.64],2)回车得到：ans =0.0000-0.00520.26340.0178即所求的拟合曲线为y=-0.0052t2+0.2634t+0.0178在编辑窗口输入如下命令：>>x=[0,5,10,15,20,25,30,35,40,45,50,55];y=[0,1.27,2.16,2.86,3.44, 3.87,4.15,4.37,4.51,4.58,4.02,4.64];>> t=0:0.1:55;>> z=-0.0052*t.^2+0.2634*t+0.0178;>> plot(x,y,'ro',t,z);grid命令执行得到如下图（图2-1）0102030405060拟合多项式与数据点的关系方法二：假设近似表达式为：y（t）=c0+c1t+c2t2第一步在命令窗口输入：>>lspoly([0,5,10,15,20,25,30,35,40,45,50,55],[0,1.27,2.16,2.86,3. 44,3.87,4.15,4.37,4.51,4.58,4.02,4.64],2)回车得到：ans =-0.00240.20370.2305即所求的拟合曲线为y=-0.0024t2+0.2037t+0.2305在编辑窗口输入如下命令：>>x=[0,5,10,15,20,25,30,35,40,45,50,55];y=[0,1.27,2.16,2.86,3.44,3.87,4.15,4.37,4.51,4.58,4.02,4.64];>> t=0:0.1:55;>> z=-0.0024*t.^2+0.2037*t+0.2305;>> plot(x,y,'ro',t,z);grid命令执行得到如下图（图2-2）拟合多项式与数据点的关系三、实验结论在利用数据的最小二乘法求拟合曲线时，选取合适的近似表达式很重要，应通过不断的试验找出较为合适的近似表达式，这样才能尽可能的提高拟合精度。
数值分析 第3章 函数逼近与曲线拟合)
在[a, b]上一致成立 。
定理：设X为一个内积空间,u1,u2,…,un∈X,矩阵
(u1, u1) (u2 , u1)
G
(u1, u2
(u1, un
) )
(u2 , u2 )
(u2 , un )
(un , u1)
(un , u2 )
(un
, un
)
称为格拉姆矩阵,则G非奇异的充分必要条件是 u1,u2,…,un线性无关 。
n1(x) (x an )n (x) n n1(x)
(n 0,1,...)
其中 0 (x) 1, －1(x) 0, n (xn (x),n (x)) /(n (x),n (x)), n (n (x),n (x)) /(n1(x),(n1(x))
(n 1,2,.....)
并且(
中找一个元素 * (x) 使 f (x) *(x) 在某种意义下
最小.
3、 范数的定义
设S为线性空间,x∈S,若存在唯一实数 || || 满足条件：
(1)‖x‖≥0;当且仅当x＝0时,‖x‖＝0; (正定性)
(2)‖αx‖＝|α|‖x‖,α∈R; (齐次性)
(3)‖x＋y‖≤‖x‖＋‖y‖,x,y∈S. (三角不等式)
类较简单的便于计算的函数类B中，求函数 P(x) B , 使P(x)与f(x)
之差在某种度量意义下最小” . 函数类A通常是区间[a，b]上的连续 函数，记作C[a，b]；函数类B通常是代数多项式，分式有理函数或 三角多项式．
2、函数空间 数学上常把在各种集合中引入某些不同的确定关系称为赋予
集合以某种空间结构，并将这样的集合称为空间.
1 2n n!
dn dxn
{(
数值分析实验报告中北大学
实验类别：数值分析专业：信息与计算科学班级：学号：姓名：中北大学理学院实验二 函数逼近与曲线拟合【实验内容】从随机的数据中找出其规律性，给出其近似表达式的问题，在生产实践和科学实验中大量存在，通常利用数据的最小二乘法求得拟合曲线。
在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量y 与时间t 的拟合曲线。
【实验方法或步骤】1、用最小二乘法进行曲线拟合；2、近似解析表达式为；33221)(t a t a t a t ++=ϕ3、打印出拟合函数)(t ϕ，并打印出)(j t ϕ与)(j t y 的误差，12,,2,1 =j ；4、另外选取一个近似表达式，尝试拟合效果的比较；5、* 绘制出曲线拟合图。
#include "stdio.h" #include "conio.h" #include "stdlib.h" #include "math.h"#define N 12//N 个节点 #define M 2//M 次拟合 #define K 2*Mvoid zhuyuan (int k,int n,float a[M+1][M+2]) {int t,i,j;float x,y;x=fabs(a[k][k]);t=k;for (i=k+1;i<=n;i++)if (fabs(a[i][k])>x){x=fabs(a[i][k]);t=i;}for (j=k;j<=n+1;j++){y=a[k][j];a[k][j]=a[t][j];a[t][j]=y;}}void xiaoyuan(int n,float a[M+1][M+2]){int k,i,j;for(i=0;i<n;i++){zhuyuan(i,n,a);for (j=i+1;j<=n;j++)for (k=i+1;k<=n+1;k++)a[j][k]=a[j][k]-a[j][i]*a[i][k]/a[i][i];}}void huidai(int n,float a[M+1][M+2],float x[M+1]){int i,j;x[n]=a[n][n+1]/a[n][n];for (i=n-1;i>=0;i--){ x[i]=a[i][n+1];for (j=i+1;j<=n;j++)x[i]=x[i]-a[i][j]*x[j];x[i]=x[i]/a[i][i];}}void main(){float x_y[N][2],A[N][K+1],B[N][M+1],AA[K+1],BB[M+1],a[M+1][M+2],m[M+1]; int i,j,n;printf("请输入%d个已知点：\n",N);for(i=0;i<N;i++){printf("(x%d y%d):",i,i);scanf("%f %f",&x_y[i][0],&x_y[i][1]);}for(i=0;i<N;i++){A[i][0]=1;for(j=1;j<=K;j++)A[i][j]=A[i][j-1]*x_y[i][0];for(j=0;j<=M;j++)B[i][j]=A[i][j]*x_y[i][1];}for(j=0;j<=K;j++)for(AA[j]=0,i=0;i<N;i++)AA[j]+=A[i][j];for(j=0;j<=M;j++)for(BB[j]=0,i=0;i<N;i++)BB[j]+=B[i][j];for(i=0;i<M+1;i++){a[i][M+1]=BB[i];for(j=0;j<=M;j++)a[i][j]=AA[i+j];}n=M;printf("正规系数矩阵为:\n");for(i=0;i<=n;i++){for(j=0;j<=n+1;j++)printf("%f ",a[i][j]);printf("\n");}xiaoyuan(n,a);huidai(n,a,m);printf("拟合曲线方程为:\ny(x)=%g",m[0]); for(i=1;i<=n;i++){printf(" + %g",m[i]);for(j=0;j<i;j++){printf("*X");}}}p3=polyfit(x,y,3);y3=polyval(p3,x);e3=norm(y-y3);t=0:5:60;pt3=polyval(p3,t);plot(t,pt3);plot(t,pt3);title('3次拟合函数')>>3次拟合函数实验三 数值积分与数值微分【实验内容】选用复合梯形公式，复合Simpson 公式，Romberg 算法高斯算法计算（1） )5343916.1(sin 44102≈-=⎰I dx x I（2） )9460831.0,1)0((sin 10≈==⎰I f dx xxI （3） dx x e I x⎰+=1024 ;（4） dx x x I ⎰++=1021)1ln( 【实验前的预备知识】1、 深刻认识数值积分法的意义；2、 明确数值积分精度与步长的关系；3、 根据定积分的计算方法，可以考虑二重积分的计算问题。
实验三 函数逼近与曲线拟合
实验三 函数逼近与曲线拟合实验3.1（曲线逼近方法的比较） 问题提出：曲线的拟合和插值，是逼近函数的基本方法，每种方法具有各自的特点和特定的适用范围，实际工作中合理选择方法是重要的。
实验内容：考虑实验2.1中的著名问题。
下面的MATLAB 程序给出了该函数的二次和三次拟合多项式。
x=-1:0.2:1;y=1/(1+25*x.*x); xx=-1:0.02:1; p2=polyfit(x,y,2); yy=polyval(p2,xx); plot(x,y,’o’,xx,yy); xlabel(‘x’); ylabel(‘y’); hold on;p3=polyfit(x,y,3); yy=polyval(p3,xx); plot(x,y,’o’,xx,yy); hold off;适当修改上述MATLAB 程序，也可以拟合其他你有兴趣的函数。
实验要求：（1）将拟合的结果与拉格朗日插值及样条插值的结果比较。
（2）归纳总结数值实验结果，试定性地说明函数逼近各种方法的适用范围，及实际应用中选择方法应注意的问题。
实验3.2：（最小二乘拟合的经验公式和模型）1.（已知经验公式）：某类疾病发病率为y ‰和年龄段x (每五年为一段，例如0~5岁为第一段，6~10岁为第二段……)之间有形如bxae y =的经验关系，观测得（1）用最小二乘法确定模型bxae y =中的参数a 和b (提示函数：lsqcurvefit ，lsqnonlin )。
（2）利用MATLAB 画出离散数据及拟合函数bxae y =图形。
（3）利用MATLAB 画出离散点处的误差图，并计算相应的均方误差。
2．（最小二乘拟合模型未知） 某年美国轿车价格的调查资料如表，其中i x 表示轿车的使用年数，i y 表示相应的平均价格，实验要求：试分析用什么形式的曲线来拟合表中的数据，并预测实验3.3（研究最佳平方逼近多项式的收敛性质）实验内容和要求：取函数xe xf =)(，在[-1，1]上以勒让德多项式为基函数，对于10,,1,0 =n 构造最佳平方逼近多项式)(x p n ，令)()()(x p x f x n n -=ε，将x x n ~)(ε的曲线画在一个图上。
数值分析实验报告--实验3--函数逼近与曲线拟合
数值分析实验三：函数逼近与曲线拟合1曲线逼近方法的比较1.1问题描述曲线的拟合和插值，是逼近函数的基本方法，每种方法具有各自的特点和特定的适用范围，实际工作中合理选择方法是重要的。
考虑实验2.1中的著名问题。
下面的MATLAB程序给出了该函数的二次和三次拟合多项式。
x=-1:0.2:1;y=1./(1+25*x.*x);xx=-1:0.02:1;p2=polyfit(x,y,2);yy=polyval(p2,xx);plot(x,y,’o’,xx,yy);xlabel(‘x’);ylabel(‘y’);hold on;p3=polyfit(x,y,3);yy=polyval(p3,xx);plot(x,y,’o’,xx,yy);hold off;实验要求：(1) 将拟合的结果与拉格朗日插值及样条插值的结果比较。
(2) 归纳总结数值实验结果，试定性地说明函数逼近各种方法的适用范围，及实际应用中选择方法应注意的问题。
1.2算法设计对于曲线拟合，这里主要使用了多项式拟合，使用Matlab的polyfit函数，可以根据需要选用不同的拟合次数。
然后将拟合的结果和插值法进行比较即可。
本实验的算法比较简单，此处不再详述，可以参见给出的Matlab脚本文件。
1.3实验结果1.3.1多项式拟合1.3.1.1多项式拟合函数polyfit和拟合次数N的关系1 / 13首先使用polyfit函数对f(x)进行拟合。
为了便于和实验2.1相比较，这里采取相同的参数，即将拟合区间[-1,1]等分为10段，使用每一段区间端点作为拟合的数据点。
分别画出拟合多项式的次数N=2、3、4、6、8、10时，f(x)和多项式函数的图像，如图1所示。
Matlab 脚本文件为Experiment3_1_1.m。
Figure 1 多项式拟合与拟合次数N的关系可以看出，拟合次数N=2和3时，拟合效果很差。
增大拟合次数，N=4、6、8时，拟合效果有明显提高，但是N太大时，在区间两端附近会出现和高次拉格朗日插值函数类似的龙格现象。
数值分析—第3章函数逼近与数据拟合法
称为广义多项式。
数值分析
三、函数的最佳平方逼近 对于给定的函数 f ( x) C[a, b] 如果存在 使
* ( x) Span 0 , 1 , , n } {
b
a
( x) f ( x) ( x) dx min
* 2
( x ) a
mn mn0 mn0
(2) 递推关系
相邻的三个切比雪夫多项式具有三项递推关系式： T0 ( x ) 1, T1 ( x ) x (n 1, 2, ) Tn1 ( x ) 2 x Tn ( x ) Tn1 ( x ) Tn (x) 的最高次项系数为 2n-1 (n = 1, 2, …)。
连续函数在[a, b]上线性无关的充分必要条件是它们 的Gramer行列式Gn 0，其中
( 0 , 0 ) ( 0 , 1 ) ( 0 , n ) G n G n ( 0 , 1 , , n ) (1 , 0 ) (1 , 1 ) (1 , n ) ( n , 0 ) ( n , 1 ) ( n , n )
(n 1, 2, )
(3) 奇偶性： 当n为偶数时，Pn (x)为偶函数； 当n为奇数时，Pn (x)为奇函数。 (4) Pn (x)的n个零点都是实的、相异的，且全
部在区间[-1, 1]内部。
数值分析
2．切比雪夫(Tchebyshev)多项式 称多项式
Tn ( x) cos(narc cos x)
Span{ 0 , 1 , , n }
并称 0 ( x), 1 ( x), , n ( x) 是生成集合的一个基底。 设函数系{
0 ( x), 1 ( x), , n ( x) ，…}线性无关，
数值分析第3讲
n
n
n
( j u j , uk ) 0, k 1,, n.
j 1
G非奇异 u1 , u2 ,, un线性无关(反证法)；反之亦然 .
在内积空间X上可以由内积导出一种 范数, 即对u X , 记 || u || ( u, u)， (1.10) 易证它满足范数定义的 正定性和齐次性, 而三角不等式由 Cauchy Schwarz不等式立得.
设0 ,, n C[a, b], 则Gram矩阵为
G G ( 0 ,, n ) ( 0 , 0 ) ( , ) 1 0 ( , ) n 0 ( 0 ,1 ) (1 ,1 ) ( n ,1 ) ( 0 , n ) (1 , n ) ( n , n )
2 内积 ( x , y ) i xi yi；范数 || x ||2 i xi i 1 i 1
n n
1/2
.
若x , y Cn，则定义加权内积 ( x , y ) i xi y i .
i 1
n
定义4 设 ( x )是区间[a , b]上的非负函数, 如果满足条件 (1) (2)
则称 || || 为线性空间S上的范数，S与 || || 一起称为赋范 线性空间，记为X .
例如，对R n上的向量x ( x1 ,, xn )T ，有 三种常用范数： || x || max | xi | ， 称为 范数或最大范数，
|| x ||1 | xi ， |
例如，三角函数族 1, cos x , sin x , cos 2 x , sin 2 x ,, 为[ , ]上的正交函数族,