非线性拟合

其中
定理：当RTR可逆时，超定方程组（3）存在最小二乘解， 且即为方程组
RTRa=RTy
的解：a=(RTR)-1RTy
24
线性最小二乘拟合 f(x)=a1r1(x)+ …+amrm(x)中 函数{r1(x), …rm(x)}的选取 1. 通过机理分析建立数学模型来确定 f(x)； 2. 将数据 (xi,yi) i=1, …n 作图，通过直观判断确定 f(x)： f=a1+a2x + + + + + f=a1+a2x+a3x2 + + + + + f=a1+a2x+a3x2 + + + + +
线性最小二乘法的求解：预备知识 超定方程组：方程个数大于未知量个数的方程组
r11a1 r12 a2 r1m am y1 ( n m) r a r a r a y nm m n n1 1 n 2 2
r11 R 其中 rn1 r 12 rn 2
14
例
已知飞机下轮廓线上数据如下，求x每改变0.1时的y值。
X Y
0 0
3 1.2
5 1.7
7 2.0
9 2.1
y
11 2.0
12 1.8
13 1.2
14 1.0
15 1.6
机翼下 轮廓线







x
返回
15
二、拟 合
1. 拟合问题引例 2.拟合的基本原理
16
拟 合 问 题 引 例 1
18
10
0
曲 线 拟 合 问 题 的 提 法
已知一组（二维）数据，即平面上 n个点（xi,yi) i=1,…n, 寻求一个函数（曲线）y=f(x), 使 f(x) 在某种准则下与所 有数据点最为接近，即曲线拟合得最好。 y + +
+
+
+ i (x+ i,yi)
+ +
+
y=f(x)
x
i 为点（xi,yi) 与曲线 y=f(x) 的距离
‘nearest’ ：最邻近插值 ‘linear’ ： 线性插值； ‘spline’ ： 三次样条插 值； ‘cubic’ ： 立方插值。 缺省时： 分段线性插值。 注意：所有的插值方法都要求 x是单调的，并且 xi不
13
能够超过x的范围。
例：在1-12的11小时内，每隔1小时测量一次温 度，测得的温度依次为：5，8，9，15，25，29，31， 30，22，25，27，24。试估计每隔1/10小时的温度 值。
f=a1+a2/x + + +
f=aebx +
+
-bx f=ae + +
Pn ( x ) L i ( x ) y i
i 0
n
其中Li(x) 为n次多项式：
( x x 0 )( x x 1 ) ( x x i 1 )( x x i 1 ) ( x x n ) L i (x) ( x i x 0 )( x i x 1 ) ( x i x i 1 )( x i x i 1 ) ( x i x n )
si ( xi ) si 1 ( xi ), si( xi ) si1 ( xi ), si( xi ) si1 ( xi ) (i 1, , n 1)
4) S ( x0 ) S ( xn ) 0 ( 自然边界条件） 2) 3) 4) ai , bi , ci , di S ( x)
求任一插值点
x ( x j ) 处的插值 y * .
*
节点可视为由
y1 y0
y
*





y g ( x)产生， g表达式复杂,
或无封闭形式, 或未知.。
3
x0 x1 x*
xn
构造一个(相对简单的)函数 y f ( x), 通过全部节点, 即
f ( x j ) y j ( j 0,1,n)
hours=1:12; temps=[5 8 9 15 25 29 31 30 22 25 27 24]; h=1:0.1:12; t=interp1(hours,temps,h,'spline'); (直接输出数据将是很多的) plot(hours,temps,'+',h,t,hours,temps,'r:') %作图 xlabel('Hour'),ylabel('Degrees Celsius’)
使n个点（xi,yi) 与曲线 y=f(x) 的距离i 的平方和最小 。
记
J (a1 , a2 , am ) i2 [ f ( xi ) yi ]2
i 1 n i 1
n
n
[ ak rk ( xi ) yi ]2
i 1 k 1
m
ቤተ መጻሕፍቲ ባይዱ
(2)
22
问题归结为，求 a1,a2, …am 使 J(a1,a2, …am) 最小。
n
即 Ra=y

r1m a1 y1 , y , a rnm am yn
超定方程一般是不存在解的矛盾方程组。
2 ( r a r a r a y ) 如果有向量a使得 i1 1 i 2 2 达到最小， im m i i 1
计算量与n无关;
xn
例
用分段线性插值法求插值,并观察插值误差.
1 g ( x) , 6 x 6 2 1 x
1.在[-6,6]中平均选取5个点作插值(xch11)
2.在[-6,6]中平均选取11个点作插值(xch12)
3.在[-6,6]中平均选取21个点作插值(xch13)
4.在[-6,6]中平均选取41个点作插值(xch14)
再用
f ( x) 计算插值，即 y f ( x ).
* *
y1 y0
y
*





x0 x1 x*
xn
返回
4
拉格朗日(Lagrange)插值
已知函数f(x)在n+1个点x0,x1,…,xn处的函数值为 y0,y1,…,yn 。求一n次多项式函数Pn(x)，使其满足： Pn(xi)=yi,i=0,1,…,n. 解决此问题的拉格朗日插值多项式公式如下
c (g/ml) 19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01
求血药浓度随时间的变化规律c(t). 作半对数坐标系(semilogy)下的图形
10
2
MATLAB(aa1)
10
1
c(t ) c0 e
kt
c, k为待定系数
0 2 4 6 8
实例：下面数据是某次实验所得，希望得到X和 f之间的关系？
x f 1 1.5 2 3.9 4 6.6 7 11.7 9 15.6 12 13 18.8 19.6 15 20.6 17 21.1
MATLAB(cn)
20
最临近插值、线性插值、样条插值与曲线拟合结果：
25
0 0 2 4 6 8 10 12 14 16 18
0C) 20.5 32.7 51.0 73.0 95.7 温度 t( 已知热敏电阻数据：
电阻R() 765 求600C时的电阻R。
1100 1000 900 800 700 20
826
873
942 1032
设 R=at+b a,b为待定系数
40
60
80
100
17
拟 合 问 题 引 例 2 已知一室模型快速静脉注射下的血药浓度数据(t=0注射300mg) t (h) 0.25 0.5 1 1.5 2 3 4 6 8
称为拉格朗日插值基函数。
5
拉格朗日(Lagrange)插值
特别地: 两点一次(线性)插值多项式:
x x0 x x1 L1 x y0 y1 x0 x1 x1 x0
三点二次(抛物)插值多项式:
x x0 x x2 x x0 x x1 x x1 x x2 L2 x y0 y1 y2 x0 x1 x0 x2 x1 x0 x1 x2 x2 x0 x2 x1
10
三次样条插值
S ( x) {si ( x), x [ xi 1, xi ],i 1,n}
1) si ( x ) ai x 3 bi x 2 ci x d i (i 1, n) 2) S ( xi ) yi (i 0,1, n) 3) S ( x ) C 2 [ x0 , xn ]
lim S ( x ) g ( x )
n
g(x)为被插值函数。
11
例
1 g ( x) , 6 x 6 2 1 x
用三次样条插值选取11个基点计算插值(ych)
返回
12
用MATLAB作插值计算
一维插值函数：
yi=interp1(x，y，xi，'method')
xi处的插 值结果 插值节点 被插值点 插值方法
ÒÑÖªÊý¾Ýµã 20
5
15
10 ÒÑÖªÊý¾Ýµã
linest Èý´Î¶àÏîÊ ½²åÖµ
10
15 nearest Èý´Î¶àÏîÊ ½²åÖµ