曲面拟合原理与实例.docxVIP

0 0 问题: 给定一组坐标(x 给定一组坐标(xg’yg’Zg) , g 1,2,…,n , 表示有 n个点。要求用以下二元 多项式函数对所给的坐标进行拟合: f(x,y) i 1 j 3j x y ij 1,1 1 ai i 1 j 1 i 1 j 1 iX y f (x, y) a11 a12y 2 a13y L q 1 a21x a22xy 2 a23xy L q 1 a2qxy M 1 i 1 i 1 2 L i 1 q 1 ai1x y 33x y aiqX yq M p 1 p 1 p 1 2 L p 1 q ap1X ap2X y ap3X y apqX y p,q p q 即 1 x 2 x x M xp ,y y 2 y M yq ,A a12 L a1q a22 L a2q M O M ap2 L apq aii a21 M ap1 则函数又可表示为 f (x, y) xTAy ,拟合的目标就是求出系数矩阵 A。 最小二乘法: 构造关于系数aj的多元函数: p q/ p q / i 1 j 1 g( ajX y i 1 j 1 zg zg)2 s( ai1，L，apq) g[f(xg，yg) Zg ] g 1 点（311，…，apq）是多元函数s（a11,L ,apq）的极小点，其中 g为权函数，默 认为1,所以点（811，…，apq ）必须满足方程组 s 3ij 在g 1的情况下，有 aij aij aij 2 [f (Xg,yg) Zg] g i 2[f(Xg,yg) 叩石If)] 2[f (Xg,yg)  i 1 j 1 Zg]Xg yg g i xg xg1yg 1f(xg,yg) xg+gz 2 g i 因此可得 u (i, j)n(Xgg 11 u (i, j) n (Xg g 1 1 yg 1 i 1 Xg yg1), v(i,j) n i 1 j 1 Xg yg Zg g 1 n n i 1 j 1 i 1 j 1 Xg yg f(Xg,yg) Xg yg Zg g 1 g 1 n p q n i Xg 1yg 1 1 a Xg yg 1 i 1 j Xg Yg 1 Zg g 1 1 1 g 1 n p,q n i Xg 1yg 1 1 1 a Xg yg i 1 j 1 Xg yg zg g 1 1,1 g 1 p,q n n a i 1 j (xg yg Xg yg 1、 i 1 ) Xg yg 1zg 1,1 g 1 g 1 p,q a u (i, j) v(i, j) (i, j) (1,1),…,(p,q) 1,1 上式实际共有p q个等式，可将这 比1(1,1) L Upq(1,1) an M O M M Un(p,q) L Upq(p,q) apq 也就是U*a=V的形式，其中 Un(1,1) L Upq (1,1) U M O M ， Un(P,q) L Upq(p,q)  p q个等式写成矩阵的形式有: v(1,1) M v( p, q) an v(1,1) a M ，V M apq v(p,q) U为pq pq阶矩阵，实现函数为function A=leftmatrix（x,p,y,q） ； V为长pq的列 向量，实现函数为function B=rightmatrix（x,p,y,q,z）。这样就可以算出列矩阵 a, 然后转化成A。 例子： 某地区有一煤矿，为估计其储量以便于开采，先在该地区进行勘探。假设 该地区是一长方形区域，长为 4公里，宽为5公里。经勘探得到如下数据: 煤矿勘探数据表 编号 1 2 3 4 5 6 7 8 9 10 横坐标（公里） 1 1 1 1 1 2 2 2 2 2 纵坐标（公里） 1 2 3 4 5 1 2 3 4 5 煤层厚度（米） 13.72 25.80 8.47 25.27 22.32 15.47 21.33 14.49 24.83 26.19 编号 11 12 13 14 15 16 17 18 19 20 横坐标（公里） 3 3 3 3 3 4 4 4 4 4 纵坐标（公里） 1 2 3 4 5 1 2 3 4 5 煤层厚度（米） 23.28 26.48 29.14 12.04 14.58 19.95 23.73 15.35 18.01 16.29 请你估计出此地区内（2x4, 1 y 5 ）煤的储量，单位用立方米表示, 并用电脑画出该煤矿的三维图象。 如果直接画出三维曲面图形: clear;x=1:4;y=1:5;[X,Y]=meshgrid(x,y) Z=[13.72 25.80 8.47 25.27 22.32; 15.47 21.3314.49 24.83 26.19;23.28 26.48 29.14 12.0414