定义于椭球面上的多元切触插值问题研究

Research on the Multivariate Contact Interpolation Problem Defined on an Ellipsoid

摘要: 以定义于平面代数曲线上的切触插值研究结果为基础，对定义于椭球面上的切触插值问题进行了研究。给出了定义于椭球面上的切触插值和插值正则性条件组问题提法，对插值条件组的拓扑结构进行了较为深入的研究，得到了定义于椭球面上的切触插值正则条件组的判定定理以及迭加构造方法，最后给出了实验算例验证了算法的有效性。

Abstract: Based on the research results of tangent interpolation defined on planar algebraic curves, the problem of tangent interpolation defined on an ellipsoid was studied. The formulation of the tan-gent interpolation problem defined on an ellipsoid and the regularization condition set for inter-polation were given. The topological structure of the interpolation condition set was thoroughly studied, and the decision theorem and superposition method for the tangent interpolation regu-larization condition set defined on an ellipsoid were obtained. Finally, experimental examples were provided to verify the effectiveness of the algorithm.

1. 引言

近年来，将多元函数插值 [1] 应用到实际的生产生活中一直是计算数学专业研究领域的一个极其重要的研究内容 [2] ，而多元多项式切触插值 [3] 则是多元多项式插值研究的重点方向。在进行多元切触插值问题时，首先要解决的就是插值适定性问题，1965年，梁学章研究了二元多元多项式插值问题，提出了多元插值适定结点组的基本概念，并且给出两种全新得构造方法：一种是添加相关直线法，另一种是用圆锥曲线表达的相关方法 [1] 。1977年，Chung和Yao也提出了类似构造方法 [4] 。1982年M. Gasca对直线迭加插值法进行了扩展与补充 [5] 。目前，有关在整个空间进行插值以及关于定义于空间中一般代数流形切触插值的研究结果相对系统与完整，关于定义于一些具体流形上的Lagrange插值也有一部分详细研究 [6] [7] [8] ，而关于有着重要生产生活实用价值的具体流形上的插值结果却相对较少。

椭球面是一类重要的二次代数曲面，其在地理测绘和导航、工程测量以及通信和卫星技术等方面有着重要实际应用。比如卫星通信和导航系统通常需要考虑地球的椭球面形状。通过使用椭球面模型，可以更精确地计算卫星轨道、通信传播路径等参数，提高通信系统的性能。除此之外，在石油和矿产资源勘探中，对地球表面的准确建模对于找到资源的位置和规模至关重要。椭球面模型可以提供更真实的地球形状，有助于资源勘探的精确性。因此，椭球面上的多元切触插值的研究有着十分重要的实际应用价值。

2. 基本定义和基本定理

本文主要研究沿三维欧式空间 $R^3$ 中的椭球面上 $F=\left\{\left(x,y,z\right)\in R^3|\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1=0\right\}$ 进行切触插值的正则性问题，为此本文首先引入若干基本定义。

设k为正整数， $r\in N$ 。令 $p_n^{\left(3\right)}$ 表示全部全次数小于等于n的三元代数多项式所构成的集合， $P_{n,r}^{\left(3\right)}\left[q\left(k\right)\right]$ 表示定义在k次代数曲面 $q\left(x,y,z\right)=0$ 上的全次数小于等于nx且有n阶方向导数的三元代数多项式空间。定义 $d_{n\cdot u}\left(k\right)$ 如下：

$\begin{array}{c}{d}_{n\cdot u}\left(k\right)=\left(\begin{array}{c}n+3\\ 3\end{array}\right)-\left(\begin{array}{c}n-\left(\mu +1\right)k+3\\ 3\end{array}\right)\\ =\{\begin{array}{l}\frac{1}{6}\left(n+1\right)\left(n+2\right)\left(n+3\right),n<\left(\mu +1\right)k\\ \frac{1}{6}\left(\mu +1\right)k\left(3n\left(n-\left(\mu +1\right)k\right)+12n+{\left(\mu +1\right)}^2{k}^2-6\left(\mu +1\right)k+11\right),n\ge \left(\mu +1\right)k\end{array}\end{array}$

$l_{n,r}\left(k\right)=\frac{1}{2}\left(n-rk\right)k\left(n-\left(r+1\right)k+4\right)+\left(\begin{array}{c}k-1\\ 3\end{array}\right)+1$

而且有 $\dim P_{n,r}^3\left[q\left(k\right)\right]=d_{n,u}\left(k\right)$ 。

假设 $q\left(x,y,z\right)=0$ 是一个k次没有重复分量的代数曲面， $B=\left\{Q_i^{\left(r\right)}|r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(k\right)\right\}$ 是这个曲面上的一个切触插值条件组。对于任意给定的一组实数组 $\left\{f_i^{\left(r\right)}|r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(k\right)\right\}$ ，要寻找一个在 $p_n^{\left(3\right)}$ 中的多项式 $p\left(x,y,z\right)$ ，并使其满足下述切触插值条件：

$\frac{{\partial }^r}{\partial n^r}p\left(Q_i^{\left(r\right)}\right)=f_i^{\left(r\right)},r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(k\right)$ (1)

式中 $\frac{{\partial }^r}{\partial n^r}p\left(Q_i^{\left(r\right)}\right)$ 代表多项式 $p\left(x,y,z\right)$ 在曲面 $q\left(x,y,z\right)=0$ 上条件组 $B$ 中点 $Q_i\left(i=1,\cdots ,l_{n-r}\left(k\right)\right)$ 处沿这个曲面的r阶法向导数。(即沿曲面位于该点处的切平面的法方向的r阶方向导数)

定义1：( $R^3$ 中的全次数型多元切触插值)

设 $z_q\in Z$ 且 $z_q\ge 0$ ， $q=1,\cdots ,m;m\ge 0$ 及一个由m个结点组成的集合 $\mathcal{A}={\left\{Q_q\right\}}_{q=1}^m$ 。 $R^3$ 中的全次数型多元切触插值问题就是如果对于给定的数组 $C_{q,\alpha }$ ， $q=1,\cdots ,m$ 及 $\left|\alpha \right|\le z_q$ ，在 $p_n^{\left(3\right)}$ 中寻找一个多项式 $p\left(x,y,z\right)$ 满足：

$\frac{{\partial }^{{\alpha }_1+\cdots +{\alpha }_d}}{\partial x_1^{{\alpha }_2}\cdots \partial x_d^{{\alpha }_d}}p\left(Q_q\right)=C_{q,\alpha },q=1,\cdots ,m;\left|\alpha \right|\le z_q$

假定 $z_q$ 及n满足 $\left(\begin{array}{c}n+d\\ d\end{array}\right)={\displaystyle \underset{q=1}{\overset{m}{\sum }}\left(\begin{array}{c}{z}_q+d\\ d\end{array}\right)}$ ，则称 $\mathcal{A}$ 及求导条件构成的插值条件组为 $p_n^{\left(3\right)}$ 的一个正则切触插值正则条件组。如果 $z_q$ ， $q=1,\cdots ,m$ 都相同，则该切触插值问题被称为一致切触插值问题。

定义2：(椭球面 $F\left(x,y,z\right)$ 上的切触插值正则条件组)

设 $F\left(x,y,z\right)$ 为如上所定义的椭球面， $B=\left\{Q_i^{\left(r\right)}|r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(2\right)\right\}$ 为沿二次代数曲面 $\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1=0$ 的一个n次r阶切触插值正则条件组，并简记为 $B\in I_{n,r}^3\left(q\right)$ (这里 $I_{n,r}^3\left(q\right)$ 表示全部沿该曲面的一个n次r阶切触插值正则条件组的集)，如果对于任意给定的一组数组

$\left\{f_i^{\left(r\right)},r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(2\right)\right\}$ ，在 $p_n^{\left(3\right)}$ 中恒存在多项式 $p\left(x,y,z\right)$ ，满足 $\frac{{\partial }^r}{\partial n^r}p\left(Q_i^{\left(r\right)}\right)=f_i^{\left(r\right)},r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(2\right)$ 。

注记：设 $F\left(x,y,z\right)$ 为如上所定义的锥面，是一个二次没有重复分量的代数曲面，

$B=\left\{Q_i^{\left(r\right)}|r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(2\right)\right\}$ 为这个曲面上的一个切触插值条件组。如果对于每个任意给定的数组 $\left\{f_i^{\left(r\right)},r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(2\right)\right\}$

方程组(1)总是存在一组解等价于若 $p_n^{\left(3\right)}$ 中存在 $p\left(x,y,z\right)$ 满足齐次切触插值条件： $\frac{{\partial }^r}{\partial n^r}p\left(Q_i^{\left(r\right)}\right)=f_i^{\left(r\right)},r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(2\right)$ ，可推出在曲面 $F\left(x,y,z\right)$ 上总是有 $p\left(x,y,z\right)$ 恒等于零。

我们同时能够得到结论：定义于锥面上的切触插值正则条件组所包含的条件总数与插值空间的维数是要相等的。

定义3：(理想)

一个子集 $I\subset K\left[x_1,\cdots ,x_n\right]$ ，如果其满足下述三个条件：

1) $0\in I$ ；

2) 如果 $f,g\in I$ ，则有 $f+g\in I$ ；

3) 如果 $f\in I$ 且 $h\in K\left[x_1,\cdots ,x_n\right]$ ，则有 $hf\in I$ 。那么I被称为是一个理想。

定义4：(生成理想)

令 $f_1,\cdots ,f_s$ 为n元多项式环中的s个多项式，则定义：

$\langle f_1,\cdots ,f_s\rangle =\left\{{\displaystyle \underset{i=1}{\overset{s}{\sum }}{h}_i{f}_i},h_1,\cdots ,h_s\in K\left[x_1,\cdots ,x_n\right]\right\}$

显然， $\langle f_1,\cdots ,f_s\rangle$ 是一个理想，称 $\langle f_1,\cdots ,f_s\rangle$ 是 $f_1,\cdots ,f_s$ 的生成理想。

定义5：(根理想)

设 $I\subset K\left[x_1,\cdots ,x_n\right]$ 是一个理想，用 $\sqrt{I}$ 表示I的根理想，它是集合： $\left\{f|f^m\in I,m\in \mathbb{Z},m\ge 1\right\}$

定义6：(关于r阶理想I的一个强H-基)

设 $f_1,\cdots ,f_s$ 为 $K\left[x_1,\cdots ,x_s\right]$ 中的s个多项式， $\mathrm{deg}\left(f_i\right)=l_i$ ，且理想 $I=\langle f_1,\cdots ,f_s\rangle$ 。若对于每个给定的多项式 $p\in I\cap P_m^{\left(n\right)}$ ，在 $K\left[x_1,\cdots ,x_s\right]$ 中恒存在多项式 $p_1,\cdots ,p_s$ ，使得

$p={\displaystyle \underset{i=1}{\overset{s}{\sum }}{p}_i{f}_i^{r+1},\mathrm{deg}\left(p_i\right)\le n-\left(r+1\right)l_i,i=1,\cdots ,s}$ 称多项式集 $\left\{f_1,\cdots ,f_s\right\}$ 是关于r阶理想I的一个强H-基。

命题1：

设I为一个理想且令 $V_1,V_2$ 为两个仿射簇，则有 $V_1\subset V_2$ 蕴含 $I\left(V_1\right)\supset I\left(V_2\right)$ 。

命题2：

若f是一个 $K\left[x_1,\cdots ,x_s\right]$ 中的s元多项式， $I=\langle f\rangle$ 是由f生成的素理想， $f=f_1^{{\alpha }_1}\cdots f_n^{{\alpha }_n}$ 是f的不可约多项式的分解，则有I的根理想为 $\sqrt{\langle f\rangle }=\langle f_1,\cdots ,f_n\rangle$ 。

特别，假若f是一个没有重复分量的代数多项式，则有I的根理想与I相等。(注：素理想是由一个元素生成的理想)

本文所获得的研究结果如下：

定理1：(椭球面 $F\left(x,y,z\right)$ 上的n次r阶切触插值正则条件组判定定理)

由如上定义的椭球面 $F\left(x,y,z\right)$ 上的切触插值条件组 $B=\left\{Q_i^{\left(r\right)}|r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(2\right)\right\}$ 能够做成定义在该曲面的一个n次r阶切触插值正则条件组的充分必要条件是：若 $p_n^{\left(3\right)}$ 中存在多项式 $p\left(x,y,z\right)$ ，满足下述齐次切触插值条件：

$\frac{{\partial }^r}{\partial n^r}p\left(Q_i^{\left(r\right)}\right)=f_i^{\left(r\right)},r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(2\right)$ 则一定在 $P_{n-2\left(\mu +1\right)}$ 中存在多项式 $p_1\left(x,y,z\right)$ ，使得

$p\left(x,y,z\right)={\left(\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1\right)}^{\mu +1}\cdot p_1\left(x,y,z\right)$

若 $n<2\left(\mu +1\right)$ ， $p_1\left(x,y,z\right)$ 恒等于零。

定理2：(构造 $p_n^{\left(3\right)}$ 上的切触插值正则条件组的添加椭球面法)

设关于 $p_n^{\left(3\right)}$ 的一个切触插值正则条件组为 $\Phi =\left\{D_{\mathfrak{A}}\left(Q_i|\mathfrak{A}\in {\mathcal{A}}_i,i=0,1,\cdots ,t\right)\right\}$ ，且 $\Phi$ 中任何点都不在椭球面 $F\left(x,y,z\right)$ 上，则对 $F\left(x,y,z\right)$ 上的一个 $n+2\left(\mu +1\right)$ 次r阶切触插值正则条件组 $B=\left\{Q_i^{\left(r\right)}|r=0,\cdots ,\mu ;i=1,\cdots ,l_{n+2r\left(\mu +1\right)}\left(2\right)\right\}$ ，则 $B\cup \Phi$ 一定构成一个关于 $P_{n+2\left(\mu +1\right)}^3$ 的切触插值正则条件组。

定理1和定理2中的结果是在前人研究的基础上，构造了椭球面切触插值的正则条件组，并将二元切触插值的通过“添加直线法”来构造插值适定泛函组，推广到三元情况的添加椭球面法构造切触插值正则条件组。

由定理1和文献3中的结果，我们得到如下推论：

推论1：(构造椭球面 $F\left(x,y,z\right)$ 上切触插值正则条件组的添加圆锥曲线法)

设二次代数曲线 $s\left(x,y\right)=p^2+c^2\left(\frac{x^2}{a^2}+\frac{y^2}{b}-1\right)=0$ ，p为正的常数，与另一条 $l\left(l\ge 2\right)$ 次代数曲线 $\omega \left(x,y\right)=0$ 正好交于2l个互不相同的点，由此确定一个插值条件组，将其记为 $\mathcal{A}=\left\{Q_i^{\left(r\right)}|r=0,\cdots ,\mu ;i=1,\cdots ,2l\right\}$ (求导方向既有沿曲线 $p^2+c^2\left(\frac{x^2}{a^2}+\frac{y^2}{b}-1\right)=0$ ，p为正的常数的法方向又有沿曲线 $\omega \left(x,y\right)=0$ 的法方向)，而 $\left\{s,w\right\}$ 正好做成关于理想 $\langle s,w\rangle$ 的强H-基。若 $B=\left\{Q_i^{\left(r\right)}|r=0,\cdots ,\mu ;i=1,\cdots ,e_{n-r}\left(2\right)\right\}\in I_{n-r}^{\left(2\right)}\left(w\right)$ ，其中 $e_{n-r}\left(2\right)=2\left(n-2r\right)-\left(\begin{array}{c}2-1\\ 2\end{array}\right)+1$ ，且满足 $B\cap \mathcal{A}=\varnothing$ ，则有 $B\cap \mathcal{A}\in I_{n+\left(\mu +1\right)l\cdot r}^2\left(w\right)$ 。

这样我们就可以通过垂直于z轴的平面与锥面 $F\left(x,y,z\right)$ 一次次相交，即把添加圆锥曲线法构造锥面 $F\left(x,y,z\right)$ 上切触插值正则条件组的问题转化为添加曲线交点法构造平面代数曲线上切触插值正则条件组的问题。

3. 定理的证明

定理1的证明：

由注记，充分性得证，以下只需证明定理的必要性。

设 $B=\left\{Q_i^{\left(r\right)}|r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(2\right)\right\}$

是定义在椭球面 $F\left(x,y,z\right)$ 上的一个n次r阶切触插值正则条件组，且在 $p_n^{\left(3\right)}$ 中存在多项式 $p\left(x,y,z\right)$ 满足

$\frac{{\partial }^r}{\partial n^r}p\left(Q_i^{\left(r\right)}\right)=f_i^{\left(r\right)},r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}(\; 2\; )$

则若 $r=0$ ，由注记可得，沿椭球面 $F\left(x,y,z\right)=\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1=0$ 恒有 $p\left(x,y,z\right)=0$ ，将 $\langle F\rangle$ 记为 $I_1$ ，将 $\langle p\rangle$ 记为 $I_2$ ，则有 $V\left(I_1\right)\subset V\left(I_2\right)$ 蕴含 $I\left(V\left(I_1\right)\right)\supset I\left(V\left(I_2\right)\right)$ ，而椭球面 $F\left(x,y,z\right)=\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1=0$

是一个二次没有重复分量的代数曲面，所以有

$I\left(V\left(I_1\right)\right)=\sqrt{I_1}=I_1$ 又由 $I\left(V\left(I_2\right)\right)=\sqrt{I_2}=I_2$ 故有 $I_2\subset I_1$ ，则根据理想的定义可知，存在 $P_{n-2}^{\left(3\right)}$ 中的 $p_1\left(x,y,z\right)$ ，使得

$p\left(x,y,z\right)=\left(\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1\right)\cdot p_1\left(x,y,z\right)$

假定结论对 $r=s,r\in Z$ 且 $r>o$ 成立，即有 $B=\left\{Q_i^{\left(r\right)}|r=0,\cdots ,\mu ;i=1,\cdots ,l_{n-r}\left(2\right)\right\}$ 为定义在曲面 $F\left(x,y,z\right)=\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1=0$ 的n次s阶切触插值正则条件组，同时满足条件

$\frac{{\partial }^r}{\partial n^r}p\left(Q_i^{\left(r\right)}\right)=0,\forall Q_i^{\left(s\right)}\in B$ (2)

有

$p\left(x,y,z\right)=\left(\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1\right)\cdot p_1\left(x,y,z\right)$ (3)

则若 $r=s+1$

我们对公式(2)两端求法向导数直到 $s+1$ 阶，并使用Leibniz公式，有

$r\left(Q_i^{\left(r\right)}\right)=0,\forall Q_i^{\left(r\right)}\in B$ 因为 $p_1\left(x,y,z\right)\in P_{n-2}^{(\; 3\; )}$

同时经过唯一且确定的 $F\left(x,y,z\right)=\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1=0$ 的所有条件点，所以沿 $F\left(x,y,z\right)$ 上总有 $p_1\left(x,y,z\right)\equiv 0$ 。由 $r=0$

情况同理可得

$p_1\left(x,y,z\right)=\left(\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1\right)\cdot {p^{\prime }}_1\left(x,y,z\right)$ (4)

将公式(4)代入到公式(3)中得：

$p\left(x,y,z\right)={\left(\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1\right)}^{s+2}\cdot {p^{\prime }}_1\left(x,y,z\right)$

由归纳法知，结论成立。

定理2的证明：

只需证明仅存在零多项式满足所给出的齐次切触插值条件即得出结论成立。可知所有的条件数是：

$\begin{array}{l}\left(\begin{array}{c}n+3\\ 3\end{array}\right)+{\displaystyle {\sum }_{r=0}^{\mu }{l}_{n+2\left(\mu +1\right),r}\left(2\right)}\\ =\left(\begin{array}{c}n+3\\ 3\end{array}\right)+{\displaystyle {\sum }_{r=0}^{\mu }\left(2\cdot \frac{1}{2}\left(n+2\left(\mu +1\right)-2r\right)\left(n+2\left(\mu +1\right)-2\left(r+1\right)+4\right)+\left(\begin{array}{c}2-1\\ 3\end{array}\right)+1\right)}\\ =\left(\begin{array}{c}n+3\\ 3\end{array}\right)+2\cdot \frac{1}{6}\left(\mu +1\right)\left(3\cdot 2n\left(n+\left(\mu +1\right)\right)+12n+4{\left(\mu +1\right)}^2-6\cdot 2\left(\mu +1\right)+11\right)\\ =\left(\begin{array}{c}n+3\\ 3\end{array}\right)+\left(\begin{array}{c}n+2\left(\mu +1\right)+3\\ 3\end{array}\right)-\left(\begin{array}{c}n+3\\ 3\end{array}\right)\\ =\left(\begin{array}{c}n+2\left(\mu +1\right)+3\\ 3\end{array}\right)\end{array}$

显然这与多项式空间 $P_{n+2\left(\mu +1\right)}^3$ 的维数相等。

假设 $P_{n+2\left(\mu +1\right)}^3$ 中存在满足齐次切触插值条件的多项式 $p\left(x,y,z\right)$ ，即有：

$D_{\mathfrak{A}}p\left(Q_i\right)=0,\forall Q_i\in \mathcal{N}$

$\frac{{\partial }^r}{\partial n^r}p\left(Q_i^{\left(r\right)}\right)=0,\forall Q_i^{\left(s\right)}\in B$

因为 $\frac{{\partial }^r}{\partial n^r}p\left(Q_i^{\left(r\right)}\right)=0,\forall Q_i^{\left(s\right)}\in B$ ，同时有 $B\in I_{n+2\left(\mu +1\right),r}^3\left(F\right)$ ， $p\left(x,y,z\right)\in P_{n+2\left(\mu +1\right)}^3$ 。

则据定理1可得，在 $p_n^{\left(3\right)}$ 中存在多项式 $p_1\left(x,y,z\right)$ ，使得

$p\left(x,y,z\right)=\left(\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1\right)\cdot p_1\left(x,y,z\right)$ (5)

对公式(5)两端求直到 $\mu +1$ 阶导数，得到

$0=D_{\mathfrak{A}}p\left(Q_2\right)=D_{\mathfrak{A}}\left\{{\left(\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}-1\right)}^{s+2}\cdot {p^{\prime }}_1\left(x,y,z\right)\right\},\forall Q_i\in \mathcal{N}$

因为对 $\forall Q_i\in \mathcal{N}$ 在 $F\left(x,y,z\right)$ 上取值都不为零，则根据Leibniz公式得对任意 $\forall Q_i\in \mathcal{N}$ ， $D_{\mathfrak{A}}r\left(Q_i\right)=0$ 。但 $\Phi$ 是关于 $p_n^{\left(3\right)}$ 的正则切触插值条件组，则有 $p_1\left(x,y,z\right)\equiv 0$ ，即有 $p\left(x,y,z\right)\equiv 0$ 。

4. 实验算例

设被插值函数为 $f\left(x,y,z\right)=\sqrt{x^2+y^2+z^2}$ ，椭球面方程： $\frac{x^2}{8}+\frac{y^2}{4}+\frac{z^2}{2}=1$ ，取该椭球面外的一点 $Q_0^0\left(0,0,0\right)$ ，则该点是 $P_0^{\left(3\right)}$ 的一个正则结点组，并且在此椭球面上取点 $Q_1^0\left(2,0,1\right)$ ，此时求出在这点处一阶法向导数 $Q_1^1\left(\frac{2\sqrt{5}}{5},0,\frac{\sqrt{5}}{5}\right)$ 由上述定理1知，条件组 $\left\{Q_0^0,Q_1^0,Q_1^1\right\}$ 构成沿椭球面 $\frac{x^2}{8}+\frac{y^2}{4}+\frac{z^2}{2}=1$ 上的1次插值正则条件组，如图1所示。设被插值函数在这些条件下的一次插值多项式为

$p\left(x,y,z\right)=a_{1}x+a_{2}y+a_{3}z+a_4$

Figure 1. The picture of ellipsoidal point taking (Including nodes)

图1. 椭球面取点图(含结点)

插值条件： $\{\begin{array}{l}p\left(Q_0^0\right)=f_0^{\left(0\right)}\\ p\left(Q_1^0\right)=f_1^{\left(0\right)}\\ \frac{{\partial }^1}{\partial n^1}p\left(Q_1^1\right)=f_1^{\left(1\right)},n=1\end{array}$ 。

代入条件得： $\{\begin{array}{l}{a}_4=0\\ 2a_1+a_3+a_4=\sqrt{5}\\ a_1=\frac{2\sqrt{5}}{5}\\ a_2=0\\ a_3=\frac{\sqrt{5}}{5}\end{array}$ 。

解方程组得到： $\{\begin{array}{l}{a}_1=\frac{2\sqrt{5}}{5}\\ a_2=0\\ a_3=\frac{\sqrt{5}}{5}\\ a_4=0\end{array}$ 。

代入得到插值多项式为： $p\left(x,y,z\right)=\frac{2\sqrt{5}}{5}x+\frac{\sqrt{5}}{5}z$ 。

经计算求得被插值函数和插值多项式在点 $\left(\frac{1}{2},0,0\right)$

处的值分别是： $f\left(\frac{1}{2},0,0\right)=\frac{1}{2}$ ， $p\left(\frac{1}{2},0,0\right)=\frac{\sqrt{5}}{5}$ 。

误差为 $m=\left|\frac{1}{2}-\frac{\sqrt{5}}{5}\right|\approx 0.0527$ 。

5. 结论

本文首先介绍了多元切触插值的相关定义与基本定理，同时重点研究了定义于椭球面上的多元切触插值正则条件组，提出构造椭球面上多元切触插值正则条件组的判定定理，最后给出实验算例说明并验证有关结论。本文创新点为给出定义于一类具体流形椭球面上的切触插值结果，其对生产生活有着重要的实用价值。在此研究基础上，多元切触插值的研究领域还需我们进行不断的探索与完善，从而使其应用范围更加广泛。

NOTES

^*通讯作者。