我们发现近似表达式
x
̃
(
k
,
n
,
一个
)和
y
̃
(
k
,
n
,
一个
)的实部和虚部
kth零
z
k
=
x
k
+
我
y
k贝塞尔的多项式
y
n
(
x
;
一个
)。获得这些封闭公式我们使用定义良好的曲线在复平面上的点是零的极限点的归一化贝塞尔多项式。因此,这些零通过静电的实现首先计算数值解释公式,然后,一个合适的函数的实部和虚部
k,
n和
一个是获得。结果表明,由此产生的复数
x
̃
(
k
,
n
,
一个
)
+
我
y
̃
(
k
,
n
,
一个
)是
O
(
1
/
n
2
)收敛到
z
k固定
k。
1。介绍
的多项式解微分方程
(1.1)
z
2
y
′′
(
z
)
+
(
一个
z
+
2
)
y
′
(
z
)
- - - - - -
n
(
n
+
一个
- - - - - -
1
)
y
(
z
)
=
0
,
一个
>
0
,
z
∈
ℂ
,系统地研究了在
1通过第一次)。他们命名(广义)贝塞尔多项式,给出明确
(1.2)
y
n
(
z
;
一个
)
=
∑
k
=
0
n
n
!
(
n
+
一个
- - - - - -
1
)
k
(
n
- - - - - -
k
)
!
k
!
(
z
2
)
k
,因为它可以显示在
2]。在这里,
(
x
)
kPochhammer符号和吗
n
=
0 1
,
…。许多属性以及相关的应用程序这个方程;旅行波在径向方向上的解决方案,在球坐标中波动方程可以写成多项式的解决方案(
1.1)。同时,这个方程在网络应用程序和滤波器的设计,各向同性湍流场等等(见专著[
2]或[
3- - - - - -
14)和引用其中一些其他的结果)。其中,一些结果的重要的问题关于取得了零的位置(
8- - - - - -
11)和(
12),显式表达式和规则和同质产品和对称函数的零多项式。另一方面,这些零平衡配置的静电解释在复杂平面对数电势和偶极子在原点已在
13],在[
14]表明,这种均衡配置是不稳定的。因此,这些情况表明,需要获得新的分析了解贝塞尔的零多项式的位置。
在本文中,我们给出近似显式公式的实部和虚部
kth零
z
k
=
x
k
+
我
y
k的
y
n
(
z
;
一个
)和显示这些新公式的近似阶贝塞尔的确切零多项式
O
(
1
/
n
2
)固定
k。
让
z
k
=
x
k
+
我
y
k,
k
=
1、2
,
…
,
n贝塞尔的零多项式
y
n
(
z
;
一个
)下令根据虚部。然后,从(
1.1),(程序获取这种非线性方程组的零多项式满足第二并给出高阶微分方程(
19]。)
(2.1)
∑
k
=
1
n
1
z
j
- - - - - -
z
k
+
(
一个
z
j
+
2
)
2
z
j
2
=
0
,在哪里
j
=
1、2
,
…
,
n,即零的实部和虚部应该满足静电方程
(2.2)
∑
k
=
1
n
x
j
- - - - - -
x
k
(
x
j
- - - - - -
x
k
)
2
+
(
y
j
- - - - - -
y
k
)
2
+
一个
x
j
3
+
2
x
j
2
+
一个
x
j
y
j
2
- - - - - -
2
y
j
2
2
(
x
j
2
+
y
j
2
)
2
=
0
,
∑
k
=
1
n
y
j
- - - - - -
y
k
(
x
j
- - - - - -
x
k
)
2
+
(
y
j
- - - - - -
y
k
)
2
+
y
j
(
一个
x
j
2
+
4
x
j
+
一个
y
j
2
)
2
(
x
j
2
+
y
j
2
)
2
=
0
。这一组非线性方程可以通过标准方法来解决。我们使用一个牛顿法来解决这些问题
n
=
500年和
一个
=
One hundred.。
让
ω
k
=
μ
k
+
我
ν
k是
k规范化的贝塞尔多项式零点
y
n
(
2
z
/
(
2
n
+
一个
- - - - - -
2
)
;
一个
),也就是说,
μ
k
=
(
2
n
+
一个
- - - - - -
2
)
x
k
/
2和
ν
k
=
(
2
n
+
一个
- - - - - -
2
)
x
k
/
2。因为它是图所示
1分段线性插值的实部和虚部
ω
k可以通过多项式拟合指数的第二个和第三个学位吗
k。
实部和虚部的零规范化贝塞尔多项式
y
n
(
2
z
/
(
2
n
+
一个
- - - - - -
2
)
;
一个
)为
一个
=
2,40100。他们绘制的功能
k为
n
=
100200300400500年在灰度强度,从低到高,根据的价值
n。
因此,我们提出以下表达式
(2.3)
μ
~
(
k
,
n
,
一个
)
=
一个
2
(
n
,
一个
)
k
2
+
一个
1
(
n
,
一个
)
k
+
一个
0
(
n
,
一个
)
,
ν
~
(
k
,
n
,
一个
)
=
b
3
(
n
,
一个
)
k
3
+
b
2
(
n
,
一个
)
k
2
+
b
1
(
n
,
一个
)
k
+
b
0
(
n
)
,的近似为零
ω
~
k
=
μ
~
(
k
,
n
,
一个
)
+
我
ν
~
(
k
,
n
,
一个
)符合我们的数据。找到这些多项式的系数的依赖
n和
一个我们考虑到数值行为数据的中间和结束点。
我们首先发现的系数二阶多项式拟合值的实部
μ
~
(
k
,
n
,
一个
)在
k
=
1和
k
=
n
/
2。在图
2我们的依赖性
μ
~
(
1
,
n
,
一个
)和
μ
~
(
n
/
2
,
n
,
一个
)在
n一些值的
一个。
的依赖关系
μ
~
(
1
,
n
,
一个
)和
μ
~
(
n
/
2
,
n
,
一个
)在
n为
一个
=
10年,20年,30,40100年。的情节所示灰度强度,从低到高,根据的价值
一个。
这些数据的模型
- - - - - -
一个
/
(
n
+
B
)和
- - - - - -
一个
/
(
n
+
B
)
- - - - - -
3
/
2收益率
(2.4)
μ
~
(
1
,
n
,
一个
)
=
- - - - - -
54
一个
+
860年
One hundred.
n
+
50
一个
+
715年
,
μ
~
(
n
2
,
n
,
一个
)
=
- - - - - -
75年
n
- - - - - -
2
一个
+
400年
50
n
+
11
一个
+
220年
。这些条件的对称
μ
~
(
k
,
n
,
一个
)对中间点导致以下系数:
(2.5)
一个
2
(
n
,
一个
)
=
p
2
(
n
,
一个
)
r
(
n
,
一个
)
,
一个
1
(
n
,
一个
)
=
p
1
(
n
,
一个
)
r
(
n
,
一个
)
,
一个
0
(
n
,
一个
)
=
p
0
(
n
,
一个
)
r
(
n
,
一个
)
,在哪里
(2.6)
p
2
(
n
,
一个
)
=
4
(
7500年
n
2
+
50625年
n
+
850年
一个
n
- - - - - -
694年
一个
2
- - - - - -
2770年
一个
+
96800年
)
,
p
1
(
n
,
一个
)
=
- - - - - -
4
(
n
+
1
)
(
7500年
n
2
+
50625年
n
+
850年
一个
n
- - - - - -
694年
一个
2
- - - - - -
2770年
一个
+
96800年
)
,
p
0
(
n
,
一个
)
=
- - - - - -
One hundred.
(
130年
n
2
- - - - - -
993年
n
- - - - - -
7656年
)
- - - - - -
20.
(
135年
n
2
+
627年
n
- - - - - -
1580年
)
一个
- - - - - -
2
(
297年
n
+
794年
)
一个
2
,
r
(
n
,
一个
)
=
5
n
(
n
- - - - - -
2
)
(
50
n
+
11
一个
+
220年
)
(
20.
n
+
10
一个
+
143年
)
。现在,找到三阶多项式的系数虚部我们遵循类似的过程。的实部的依赖
ν
~
(
1
,
n
,
一个
)和
∂
ν
~
(
k
,
n
,
一个
)
/
∂
k
|
k
=
n
/
2在
n一些值的
一个如图
3。
的依赖关系
ν
~
(
1
,
n
,
一个
)和
∂
ν
~
(
k
,
n
,
一个
)
/
∂
k
|
k
=
n
/
2在
n为
一个
=
10年,20年,30,40100年。的情节所示灰度强度,从低到高,根据的价值
一个。
再一次,这些数据的模型
一个
/
(
n
+
B
)
- - - - - -
1和
一个
/
n收益率
(2.7)
ν
~
(
1
,
n
,
一个
)
=
- - - - - -
25
n
- - - - - -
一个
- - - - - -
50
25
(
n
- - - - - -
1
)
,
∂
ν
~
(
k
,
n
,
一个
)
∂
k
|
k
=
n
/
2
=
96年
n
(
一个
+
25
)
。此外,我们有
(2.8)
ν
~
(
n
2
,
n
,
2
)
=
0
,
ν
~
(
n
,
n
,
2
)
=
- - - - - -
ν
~
(
1
,
n
,
2
)
,因此,给出了系数
(2.9)
b
3
(
n
,
一个
)
=
问
3
(
n
,
一个
)
年代
(
n
,
一个
)
,
b
2
(
n
,
一个
)
=
问
2
(
n
,
一个
)
年代
(
n
,
一个
)
,
b
1
(
n
,
一个
)
=
问
1
(
n
,
一个
)
年代
(
n
,
一个
)
,
b
0
(
n
,
一个
)
=
问
0
(
n
,
一个
)
年代
(
n
,
一个
)
,在哪里
(2.10)
问
3
(
n
,
一个
)
=
- - - - - -
200年
(
23
n
3
- - - - - -
92年
n
2
+
90年
n
+
4
)
+
200年
(
n
3
- - - - - -
5
n
2
+
8
n
- - - - - -
6
)
一个
- - - - - -
8
(
n
2
- - - - - -
2
n
+
2
)
一个
2
,
问
2
(
n
,
一个
)
=
One hundred.
(
69年
n
4
- - - - - -
255年
n
3
+
186年
n
2
+
92年
n
+
8
)
- - - - - -
One hundred.
(
3
n
4
- - - - - -
12
n
3
+
9
n
2
+
4
n
- - - - - -
12
)
一个
+
4
(
3
n
3
- - - - - -
3
n
2
+
4
)
一个
2
,
问
1
(
n
,
一个
)
=
50
(
21
n
4
+
54
n
3
- - - - - -
522年
n
2
+
748年
n
- - - - - -
176年
)
- - - - - -
50
(
3
n
4
- - - - - -
9
n
3
- - - - - -
6
n
2
+
26
n
- - - - - -
24
)
一个
+
2
(
3
n
3
- - - - - -
6
n
+
8
)
一个
2
,
问
0
(
n
,
一个
)
=
- - - - - -
25
(
25
n
4
- - - - - -
217年
n
3
+
618年
n
2
- - - - - -
660年
n
+
184年
)
- - - - - -
25
(
n
4
- - - - - -
2
n
3
- - - - - -
7
n
2
+
16
n
- - - - - -
12
)
一个
+
(
n
3
+
n
2
- - - - - -
4
n
+
4
)
一个
2
年代
(
n
,
一个
)
=
25
(
n
- - - - - -
1
)
2
(
n
- - - - - -
2
)
2
(
一个
+
25
)
。因此,替换(
2.10),(
2。9),(
2。6)和(
2。5),分别在(
2。3)收益率近似封闭表达式
(2.11)
z
~
k
=
x
~
(
k
,
n
,
一个
)
+
我
y
~
(
k
,
n
,
一个
)
=
2
μ
~
(
k
,
n
,
一个
)
2
n
+
一个
- - - - - -
2
+
我
2
ν
~
(
k
,
n
,
一个
)
2
n
+
一个
- - - - - -
2
,在哪里
k
=
1、2
,
…
,
n非规范贝塞尔的零多项式
y
n
(
z
;
一个
)。中给出的表达式(
2.11收敛到零
z
k这些多项式,当我们将显示在下面。
3所示。收敛
后(
9),我们定义
(3.1)
W
(
z
)
=
e
1
+
1
/
z
2
z
(
1
+
1
+
1
/
z
2
)
,和表示
Γ定义的曲线
(3.2)
Γ
=
{
z
∈
ℂ
:
|
W
(
z
)
|
=
1
,
|
参数
z
|
≥
π
2
}
,它包含的极限点
ω
^
k0的规范化贝塞尔多项式
y
n
(
2
z
/
(
2
n
+
一个
- - - - - -
2
)
;
一个
)。然后,事实证明在
8),零
ω
k的
y
n
(
2
z
/
(
2
n
+
一个
- - - - - -
2
)
;
一个
)订单的方法
O
(
1
/
n
)极限的值
ω
^
k,也就是说,
(3.3)
|
ω
k
- - - - - -
ω
^
k
|
=
O
(
1
n
)
,作为
n
→
∞。
因此,如果我们显示
|
ω
~
k
- - - - - -
ω
^
k
|
=
O
(
1
/
n
),我们将证明
(3.4)
|
ω
k
- - - - - -
ω
~
k
|
=
O
(
1
n
)
,因此,考虑到这一点
ω
k
=
(
2
n
+
一个
- - - - - -
2
)
z
k
/
2显式表达式(
2.11)方法
O
(
1
/
n
2
)零
z
k贝塞尔的多项式
y
n
(
z
;
一个
)。
为了这个目的,我们只是用表达式代替
ω
~
k
=
μ
~
(
k
,
n
,
一个
)
+
我
ν
~
(
k
,
n
,
一个
)由(
2。3)(
3所示。1)获得,经过长时间的计算,这的扩张
W
(
ω
~
k
)而言,
1
/
n是
(3.5)
W
(
ω
~
k
)
=
1
+
(
2
一个
2
- - - - - -
One hundred.
一个
(
3
k
- - - - - -
2
)
- - - - - -
50
(
42
k
- - - - - -
67年
)
25
(
一个
+
25
)
+
h
(
k
,
一个
)
25
(
一个
+
25
)
(
因为
2
(
t
(
k
,
一个
)
)
- - - - - -
罪
2
(
t
(
k
,
一个
)
)
)
]
1
n
+
O
(
1
n
3
/
2
)
,在哪里
(3.6)
h
(
k
,
一个
)
=
(
25
+
一个
)
2
(
130年
+
27
一个
+
300年
k
)
2
+
(
2
一个
2
- - - - - -
One hundred.
一个
(
3
k
- - - - - -
2
)
- - - - - -
50
(
42
k
- - - - - -
67年
)
)
2
,
t
(
k
,
一个
)
=
(
25
+
一个
)
(
130年
+
27
一个
+
300年
k
)
4
(
1675年
- - - - - -
One hundred.
一个
+
一个
2
- - - - - -
1050年
k
+
150年
一个
k
)
。这意味着
(3.7)
|
W
(
ω
~
k
)
|
=
1
+
O
(
1
n
)
,固定
k和
一个。因此,
ω
~
k订单的方法
O
(
1
/
n
)的
Γ曲线和(
3所示。4)。从这里我们有
(3.8)
|
z
k
- - - - - -
z
~
k
|
=
O
(
1
n
2
)
,作为
n
→
∞。数值计算证实和扩展这个结果。图
4显示了最大值的行为
|
z
k
- - - - - -
z
~
k
|在
k因为他们依靠
n的具体情况
一个
=
2。数字计算(
2。2)作为确切的0
z
k。这些数据的
1
/
n
一个与
一个
=
1.7。
块的值
马克斯
k
=
1
n
|
z
k
- - - - - -
z
~
k
|对
n的情况下
一个
=
2。
4所示。一些测试
为了让示例应用的近似表达式(
2.11),我们考虑下面的情况。
4.1。真正的零
一个真正独特的封闭公式为零
α
n
(
一个
)贝塞尔的多项式
y
n
(
z
;
一个
)可以获得的替代吗
k
=
(
n
+
1
)
/
2在现实的一部分(
2.11),
x
~
(
k
,
n
)。这给了
(4.1)
α
~
n
(
一个
)
=
- - - - - -
4
n
3
+
1
75年
(
50
- - - - - -
71年
一个
)
+
(
- - - - - -
707年
一个
2
+
5900年
一个
+
35000年
)
7500年
n
+
O
(
1
n
2
)
,作为我们的新结果。在[
2,
11非常精确的表达式
α
n
(
一个
)给出了。特别是,以下公式:
(4.2)
2
α
n
(
一个
)
=
- - - - - -
1.325486838
n
- - - - - -
1.00628995
一个
+
1.349836480
+
O
(
1
2
n
+
一个
- - - - - -
2
)
,给出了在
11]。扩大
2
/
α
~
n
(
一个
)的权力
n我们发现
(4.3)
2
α
~
n
(
一个
)
≃
- - - - - -
1.33333
n
- - - - - -
0.946667
一个
+
0.666667
,表明良好的相对两个结果之间的协议。
4.2。权力和
这里我们进行相应的乘法和使用某些情况下Faulhaber的公式。然后我们比较我们的结果和精确的。
(1)
的和零。简单的总和
z
~
k(cf。
2.11)给出了一个复杂的实部的表达式。然而,扩大这个总和的实部和虚部
(4.4)
年代
~
1
(
n
)
=
∑
k
=
1
n
z
~
k
=
- - - - - -
1
+
(
53
一个
One hundred.
- - - - - -
71年
30.
+
我
32
一个
+
25
)
1
n
+
O
(
1
n
2
)
。确切的结果
年代
1
(
n
)
=
- - - - - -
1我们可以看到从(
1.2)。
(2)
零的平方之和。在这种情况下我们的具体情况
一个
=
1。为我们获得这个值
(4.5)
年代
~
2
(
n
)
=
∑
k
=
1
n
z
~
k
2
=
3469年
5915年
n
+
O
(
1
n
2
)
。确切的数目
(4.6)
年代
2
(
n
)
=
1
2
n
- - - - - -
1
=
1
2
n
+
O
(
1
n
2
)可以找到其他地方(
12]。