1. 引言
熟知,多元函数插值长期以来一直是计算数学研究领域的一个主要研究内容(详见文献 [1] ),有关多元函数插值基本理论和方法研究中一个基本问题是多元插值函数的唯一存在性问题,也就是插值的唯一可解性问题。由文献 [1] 可知,国内外学者对这一问题的研究主要有两个判别,一种是给定插值空间,去构造相应唯一正解结点组;另一种是给定结点组,去构造相应正则插值空间,而且要求空间的次数尽可能低 [2] [3] [4] ,对于某一类问题,目前,有关在整个空间进行插值以及关于定义于空间中一般代数流形插值的研究结果相对完备 [5] ,而关于有着重要实用价值的具体流形上的插值结果并不多见。梁学章等人 [6] [7] 讨论了单位球面上的纬线组选取插值唯一正解结点组的方法。Castell等人 [8] 利用球面上偶数条纬线上等距点组构造了球面上的唯一正解结点组。
抛物柱面是除球面外的另一类主要的二次代数曲面,其在军事上,天文方面上以及生活上均有广泛应用。例如,在电子对抗微波毫米波系统中,往往需要把大功率的定向电磁辐射波束辐射到作战空间中,传递己方的有用信息或者干扰同频的敌方信息,抛物柱面天线就是通常采用多个天线组成的线来进行空间功率的合成,用抛物柱面作为反射面来增加天线的口径,合成高功率的尖锐波束。星载抛物柱面天线是各类型星载天线的一种,由于其具有方向性强、增益高、易于光束自动扫描等特点,星载抛物柱面天线已经成为星载天线新的发展方向。抛物柱面聚光器,这种新型的太阳能聚光器,它生产工艺简单,生产成本低。因此,对抛物柱面上的插值问题研究意义重大。
2. 基本定义和基本定理
本文主要研究沿三维欧式空间
中的抛物柱面
上进行Lagrange插值的唯一可解性问题,为此本文首先引入若干基本概念。
设n为非负整数,令
表示所有全次数为n的三元代数多项式构成的集,即
,而且有
。
定义1:(
的插值唯一可解结点组)设
,令
为
中m个互异点构成的点集,如果对于任意给定的数组
,恒存在唯一多项式
,满足:
,则称
为
的一个唯一可解结点组。
定义2:(
上的插值唯一可解结点组)设
为如上所定义的旋转抛物面,
为
在
上的限制,
,称
为定义于
上的一个n次插值唯一可解结点组,如果对于任意给定的数组
,恒存在多项式
,满足
。
本文所获得研究结果如下:
定理1:(构造
插值唯一可解结点组的添加旋转抛物面法)
设
为如上所定义,
,结点组
关于
的一个唯一可解结点组,而
是定义于
的一个
次唯一可解结点组,则
必定构成
的唯一可解结点组。
由定理1和文献 [2] 、 [3] 、 [4] 中的结果,我们得到如下推论:
推论1:(构造
上插值正则结点组添加圆锥曲线法)
设
为
上的
次插值唯一可解结点组,平面
与
横截相交于圆锥曲线
,
是定义于圆锥曲线
上的一个
次唯一可解结点组,则
必定构成定义于
上的一个
次唯一可解结点组。
定理2:(判定定理)
由抛物柱面
上的
个互异点构成的点组
能够做成该抛物柱面
上的
次插值唯一可解结点组的充分必要条件是,若
,满足
,蕴含着如此的
沿抛物柱面
恒为零。
3. 定理的证明
定理1的证明:
证明:设
。因为
为定义于
上的
次唯一可解结点组,由定义2,对任意给定数组
恒存在多项式
使得
。
又因为
关于
的一个唯一可解结点组,由定义1知,对任意的数组
,恒存在多项式
,使得
成立,其中
为
的三维坐标。构造一个多项式
显然有
,且满足:
则由定义1知,
为
的唯一可解结点组。
定理2的证明:
证明:首先证明充分性。只需证明命题“若有多项式
,满足
,
,则必有
。”成立即可。设
满足
,
,由条件可知,
,
。从而,对于定义于抛物面
的一个
次唯一可解结点组
,亦有,
,即,
。又因为
为
的唯一可解结点组,故此有
。
必要性:令
,取
为关于
的唯一可解结点组,可以断言:
构成了
的唯一可解结点组。事实上,对于任意给定的实数
,由于
为抛物柱面
上的
次唯一可解结点组,故存在多项式
满足
又因为
且为
的唯一可解结点组,则存在多项式
满足条件
其中
为
的三维坐标。则多项式
(2.1)
满足
,即
是
的唯一可解结点组。此外,在上述证明过程中特别取
,则(2.1)式中的
,此时
为满足定理中的插值条件的多项式,故由
空间中满足相同插值条件的多项式的唯一存在性,有
即
在抛物柱面
上恒取值为零。
4. 实验算例
设被插值函数为
,抛物柱面方程:
,取该抛物柱面外的一点
并且在此抛物柱面上取互异的9个点
,
,
,
,
,
,
,
,
由上述定理2知,点组
构成定义于抛物柱面
上的2次插值唯一可解结点组,如图1所示。设被插值函数在这些点上的二次插值多项式为
将插值条件带入
,得到方程组为
,其中
,
,
Figure 1. Parabolic cylinder point taking
图1. 旋转抛物面取点图
解方程组得到:
代入得到插值多项式为:
被插值多项式在
,
的插值结果分别为4.7996,9.5429,而插值函数在这两点处的精确值分别为
,
,误差分别为
,
。