线性插值函数的基函数构造

栀彧

已于 2022-07-12 11:38:27 修改

阅读量1.4k

点赞数 2

分类专栏：笔记文章标签：算法 matlab

于 2022-07-12 11:31:44 首次发布

本文链接：https://blog.csdn.net/qq_42158869/article/details/125739650

版权

笔记专栏收录该内容

7 篇文章 4 订阅

订阅专栏

线性插值函数的基函数构造

这里介绍的基函数为线性插值对应的基函数，也就是通过不在同一直线上的三个点构造的平面。这里我们直接通过三个点构造一个平面，再将其写法简单的变形就可以找到对应的基函数。

设 $A:(x_1,y_1,z_1),B:(x_2,y_2,z_2),C:(x_3,y_3,z_2)$ 是空间中不在同一直线上的三个点，且 $A, B, C$ 三点在平面 $\pi:ax+by+cz+d=0$ 上。

这里，利用点法式确定平面 $\pi$ 。首先，根据向量 $\vec{AB}=[x_2-x_1,y_2-y_1,z_2-z_1]$ ， $\vec{AC}=[x_3-x_1,y_3-y_1,z_3-z_1]$ 构造垂直于 $\vec{AB}，\vec{AC}$ 的向量 $\vec{n}$ ， $\vec{n}$ 即为所求平面的法向量，

$\vec{n} = \vec{AB}\times\vec{AC}\\\ \ \ = [(y_2-y_1)(z_3-z_1)-(y_3-y_1)(z_2-z_1), \\\ \ \ \ \ \ \ \ \ (z_2-z_1)(x_3-x_1)-(z_3-z_1)(x_2-x_1), \\\ \ \ \ \ \ \ \ \ (x_2-x_1)(y_3-y_1)-(x_3-x_1)(y_2-y_1)]$

进而，设点 $P (x, y, z)$ 在平面 $\pi$ 上，则有 $\vec{PA}\perp\vec{n}$ ，即
$0 = (x-x_1)((y_2-y_1)(z_3-z_1)-(y_3-y_1)(z_2-z_1))\\\ \ \ +(y-y_1)((z_2-z_1)(x_3-x_1)-(z_3-z_1)(x_2-x_1))\\\ \ \ +(z-z_1)((x_2-x_1)(y_3-y_1)-(x_3-x_1)(y_2-y_1))$

整理可得：

$x_2-x_1)(y_3-y_1)-(x_3-x_1)(y_2-y_1))z \\= -(z_1((x-x_1)(y_3-y_1-(y_2-y_1))\\\ \ \ +z_2((y-y_1)(x_3-x_1)-(x-x_1)(y_3-y_1))\\\ \ \ +z_3((x-x_1)(y_2-y_1)-(y-y_1)(x_1-x_1))\\\ \ \ +z_1(y-y_1)(x_2-x_1-(x_3-x_1)))\\\ \ \ +z_1((x_2-x_1)(y_3-y_1)-(x_3-x_1)(y_2-y_1)) \\=z_1((x-x_2)(y_2-y_3)+(y-y_2)(x_3-x_2))\\+z_2((x-x_3)(y_3-y_1)+(y-y_3)(x_1-x_3))\\+z_3((x-x_1)(y_1-y_2)+(y-y_1)(x_2-x_1))$

由三角形面积 $S_{\triangle ABC}=1/2((x_2-x_1)(y_3-y_1)-(x_3-x_1)(y_2-y_1))$ 因此有

$[z_1((x-x_2)(y_2-y_3)+(y-y_2)(x_3-x_2))\\\ \ +z_2((x-x_3)(y_3-y_1)+(y-y_3)(x_1-x_3))\\\ \ +z_3((x-x_1)(y_1-y_2)+(y-y_1)(x_2-x_1))]/(2*S_{\triangle ABC})$

因此线性插值基函数可以写为

$N_A = ((x-x_2)(y_2-y_3)+(y-y_2)(x_3-x_2))/(2*S_{\triangle ABC})$

$N_B=((x-x_3)(y_3-y_1)+(y-y_3)(x_1-x_3))/(2*S_{\triangle ABC})$

$N_C=((x-x_1)(y_1-y_2)+(y-y_1)(x_2-x_1))/(2*S_{\triangle ABC})$

上述基函数也可以写为

$N_A = (x(y_2-y_3)+y(x_3-x_2)+x_2y_3-y_2x_3)/(2*S_{\triangle ABC})$

$N_B=(x(y_3-y_1)+y(x_1-x_3)+x_3y_1-y_3x_1)/(2*S_{\triangle ABC})$

$N_C=(x(y_1-y_2)+y(x_2-x_1)+x_1y_2-x_2y_1)/(2*S_{\triangle ABC})$

注意：这里利用空间中不在同一条直线上的三个点确定唯一平面，通过这种方式构造平面方程，进一步改写方程样式，得到基函数。

实验结果

我们绘制过点 $A : (1, 2, 1), B : (2, 3, 6), C : (1, 6, 3)$ 的平面 $\pi$ ，结果如下
在这里插入图片描述

实例代码如下：

clc,clear
x = [1,2,4];
y = [2,3,1];
z = [1,6,3];
p = (x(1)-x(2))*(y(1)-y(3)) - (y(1)-y(2))*(x(1)-x(3));
% K1 = @(t,s)(s-y(2))*(x(3)-x(2)) - (t-x(2))*(y(3)-y(2));
% K2 = @(t,s)(s-y(3))*(x(1)-x(3)) - (t-x(3))*(y(1)-y(3));
% K3 = @(t,s)(s-y(1))*(x(2)-x(1)) - (t-x(1))*(y(2)-y(1));
K1 = @(t,s)((y(2)-y(3))*t - (x(2)-x(3))*s + x(2)*y(3)-x(3)*y(2));
K2 = @(t,s)((y(3)-y(1))*t - (x(3)-x(1))*s + x(3)*y(1)-x(1)*y(3));
K3 = @(t,s)((y(1)-y(2))*t - (x(1)-x(2))*s + x(1)*y(2)-x(2)*y(1));
disp('下面这个矩阵为基函数在对应点处的值，即克罗内克尔符号矩阵')
disp([K1(x(1),y(1)),K1(x(2),y(2)),K1(x(3),y(3))
 K2(x(1),y(1)),K2(x(2),y(2)),K2(x(3),y(3))
 K3(x(1),y(1)),K3(x(2),y(2)),K3(x(3),y(3))]/p)
F = @(t,s) (z(1) * K1(t,s) + z(2) * K2(t,s) + z(3) * K3(t,s))/p;
% F = triangle_2D(x,y,z);
xx = 0:0.1:5;
for i = 1:length(xx)
    for j = 1:length(xx)
        Z(j,i) = F(xx(i),xx(j));
    end
end
% 注意到surf绘图中，Z的列数对应的是X，行数对应的是Y，在制作Z矩阵时，应该注意这一点。
[X,Y] = meshgrid(xx);
surf(X,Y,Z)
hold on
plot3(x,y,z,'r*')