2025-11-19T12:16:13.706625

Uniform Asymptotic approximation and numerical evaluation of the Reverse Generalized Bessel Polynomial zeros

Dunster, Gil, Ruiz-Antolin et al.
Uniform asymptotic expansions are derived for the zeros of the reverse generalized Bessel polynomials of large degree $n$ and real parameter $a$. It is assumed that $-Δ_{1} n+\frac{3}{2} \leq a \leq Δ_{2} n$ for fixed arbitrary $Δ_{1} \in (0,1)$ and bounded positive $Δ_{2}$. For this parameter range at most one of the zeros is real, with the rest being complex conjugates. The new expansions are uniformly valid for all the zeros, and are shown to be highly accurate for moderate or large values of $n$. They are consequently used as initial values in a very efficient numerical algorithm designed to obtain the remaining complex zeros using Taylor series.
academic

Равномерная асимптотическая аппроксимация и численное вычисление нулей обратного обобщённого полинома Бесселя

Основная информация

  • ID статьи: 2510.14295
  • Название: Uniform Asymptotic approximation and numerical evaluation of the Reverse Generalized Bessel Polynomial zeros
  • Авторы: T. M. Dunster, A. Gil, D. Ruiz-Antolin, J. Segura
  • Классификация: math.CA (Classical Analysis and ODEs)
  • Дата публикации: 16 октября 2025 г.
  • Ссылка на статью: https://arxiv.org/abs/2510.14295

Аннотация

В данной работе выведены равномерные асимптотические разложения нулей обратного обобщённого полинома Бесселя для больших степеней nn и вещественного параметра aa. Предполагается, что параметр находится в диапазоне −Δ1n+32≤a≤Δ2n-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n, где Δ1∈(0,1)\Delta_1 \in (0,1) — произвольная фиксированная константа, а Δ2\Delta_2 — ограниченная положительная константа. В этом диапазоне параметров может быть не более одного вещественного нуля, остальные представляют собой комплексно-сопряжённые пары. Новые разложения равномерно эффективны для всех нулей и демонстрируют высокую точность при средних и больших значениях nn. Эти разложения используются в качестве начальных приближений в эффективном численном алгоритме, который использует ряды Тейлора для вычисления оставшихся комплексных нулей.

Исследовательский контекст и мотивация

Описание проблемы

Обобщённый полином Бесселя определяется как: yn(z;a)=∑k=0n(nk)(n+a−1)k(12z)ky_n(z; a) = \sum_{k=0}^{n} \binom{n}{k} (n+a-1)_k \left(\frac{1}{2z}\right)^k

где (α)k=Γ(α+k)/Γ(α)(α)_k = Γ(α + k)/Γ(α) — символ Похгаммера. Нули этих полиномов обычно являются комплексными и имеют многочисленные приложения в прикладной математике и инженерии.

Мотивация исследования

  1. Практические потребности: нули обобщённых полиномов Бесселя имеют важные приложения в прикладной математике и инженерии
  2. Вычислительные вызовы: существующие методы вычисления этих нулей имеют ограничения в эффективности и точности
  3. Теоретическое совершенствование: необходимо установить равномерно эффективную асимптотическую теорию для описания поведения всех нулей
  4. Оптимизация численных алгоритмов: требуются эффективные численные методы для вычисления комплексных нулей

Ограничения существующих методов

  • Предыдущие работы 2, 13 имели ограничения в точности и эффективности вычислений
  • Отсутствовали асимптотические разложения, равномерно эффективные для всех нулей
  • Выбор начальных приближений для численных алгоритмов был неоптимальным

Основные вклады

  1. Выведены равномерные асимптотические разложения: установлены асимптотические разложения, равномерно эффективные для всех нулей обратного обобщённого полинома Бесселя
  2. Разработана теория точек поворота: на основе теории WKB и метода Лиувилля-Грина обработаны дифференциальные уравнения с точками поворота
  3. Построен эффективный численный алгоритм: разработан эффективный алгоритм вычисления комплексных нулей, объединяющий асимптотические разложения и ряды Тейлора
  4. Получены высокоточные приближения: разложения достигают относительной ошибки порядка 10−1510^{-15} при средних и больших значениях nn

Подробное описание методов

Постановка задачи

Поиск нулей tm(u,a)t_m(u, a) обратного обобщённого полинома Бесселя θn(z;a)=znyn(z−1;a)θ_n(z; a) = z^n y_n(z^{-1}; a), где u=n+12u = n + \frac{1}{2}, в диапазоне параметров −Δ1n+32≤a≤Δ2n-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n.

Основная теоретическая база

1. Преобразование дифференциального уравнения

Через масштабирующую функцию: wn(0)(z;a)=2−n−a+1z1−n−a/2e−zθn(z;a)w_n^{(0)}(z; a) = 2^{-n-a+1} z^{1-n-a/2} e^{-z} θ_n(z; a)

задача преобразуется в решение дифференциального уравнения: d2wdz2={1+a−2z+(2n+a)(2n+a−2)4z2}w\frac{d^2w}{dz^2} = \left\{1 + \frac{a-2}{z} + \frac{(2n+a)(2n+a-2)}{4z^2}\right\} w

2. Теория Лиувилля-Грина

Определяются параметры: u=n+12,α=a−2uu = n + \frac{1}{2}, \quad α = \frac{a-2}{u}

Дифференциальное уравнение переписывается как: d2wdz2={u2f(α,z)+g(z)}w\frac{d^2w}{dz^2} = \{u^2 f(α, z) + g(z)\} w

где: f(α,z)=(z+12α)2+1+αz2,g(z)=−14z2f(α, z) = \frac{(z + \frac{1}{2}α)^2 + 1 + α}{z^2}, \quad g(z) = -\frac{1}{4z^2}

3. Анализ точек поворота

Точки поворота расположены в: z1,2(α)=±iσ−12α,σ=1+αz_{1,2}(α) = ±iσ - \frac{1}{2}α, \quad σ = \sqrt{1 + α}

Через переменную Лиувилля ζζ и переменную LG ξξ: 23ζ3/2=ξ=∫z1(α)zf1/2(α,t)dt\frac{2}{3}ζ^{3/2} = ξ = \int_{z_1(α)}^z f^{1/2}(α, t) dt

Асимптотические разложения нулей

Нули удовлетворяют неявному уравнению: Z(u,a,u−1tm(u,a))=u−2/3amZ(u, a, u^{-1}t_m(u, a)) = u^{-2/3}a_m

где ama_m — mm-й отрицательный нуль функции Эйри Ai(x)Ai(x).

Форма разложения: tm(u,a)∼u∑s=0∞τm,s(α)u2s(u→∞)t_m(u, a) \sim u \sum_{s=0}^{\infty} \frac{τ_{m,s}(α)}{u^{2s}} \quad (u → ∞)

Вычисление основных коэффициентов

  • Ведущий член τm,0τ_{m,0}: определяется сложным неявным уравнением (3.16)
  • Последующие члены: вычисляются через рекуррентные соотношения (3.20)-(3.23):

τm,1=−Υm,1ζm,0′τ_{m,1} = -\frac{Υ_{m,1}}{ζ'_{m,0}}

τm,2=−12ζm,0′{τm,12ζm,0′′+2τm,1Υm,1′+2Υm,2}τ_{m,2} = -\frac{1}{2ζ'_{m,0}}\{τ_{m,1}^2 ζ''_{m,0} + 2τ_{m,1} Υ'_{m,1} + 2Υ_{m,2}\}

Проектирование численного алгоритма

1. Итерационная функция

Tn(a,z)=z−1Ωn(a,z)arctan⁡(Ωn(a,z)wn(0)(z;a)∂wn(0)(z;a)/∂z)T_n(a, z) = z - \frac{1}{\sqrt{Ω_n(a, z)}} \arctan\left(\frac{\sqrt{Ω_n(a, z)}w_n^{(0)}(z; a)}{\partial w_n^{(0)}(z; a)/\partial z}\right)

2. Функция шага

Hn(a,z)=z+πΩn(a,z)H_n(a, z) = z + \frac{π}{\sqrt{Ω_n(a, z)}}

3. Разложение в ряд Тейлора

Использование рекуррентных соотношений (4.7) для вычисления производных высокого порядка, избегая прямого вычисления значений функции.

Экспериментальная установка

Параметры тестирования

  • Диапазон степеней: n=15,30,50,200,500,1000,2000n = 15, 30, 50, 200, 500, 1000, 2000
  • Значения параметров: a=1.01,20.2a = 1.01, 20.2 и т.д.
  • Индексы нулей: m=1,3,10,15,25m = 1, 3, 10, 15, 25 и т.д.

Показатели оценки

  • Относительная ошибка: сравнение с высокоточной реализацией в Maple
  • Время вычисления: время выполнения на CPU
  • Сходимость: точность сходимости итерационного алгоритма ε=10−12ε = 10^{-12}

Детали реализации

  • Асимптотические разложения используют первые 5 коэффициентов
  • Реализация в Maple использует 60-битную точность
  • Алгоритм в Matlab тестировался на Dell Latitude 7410

Экспериментальные результаты

Основные результаты

1. Высокоточные приближения

Относительные ошибки, показанные в таблицах 1 и 2:

  • Для a=1.01a = 1.01: относительная ошибка порядка 10−1810^{-18} до 10−1510^{-15}
  • Для a=20.2a = 20.2: относительная ошибка порядка 10−1810^{-18} до 10−1510^{-15}
  • Относительная ошибка во всех тестовых случаях близка или превосходит 10−1510^{-15}

2. Анализ зависимости от параметров

На рисунке 1 показано, что с увеличением параметра aa относительная ошибка уменьшается, максимальная относительная ошибка составляет менее 4×10−154×10^{-15}.

3. Эффективность алгоритма

Результаты времени вычисления из таблицы 3:

  • n=30n = 30 (15 нулей): 2.8×10−32.8×10^{-3} сек
  • n=2000n = 2000 (1000 нулей): 3.7×10−23.7×10^{-2} сек
  • При увеличении количества нулей в 67 раз время вычисления увеличивается только в 13 раз, что демонстрирует высокую эффективность алгоритма

Анализ сходимости

На рисунке 4 показано, что даже при вычислении большого количества нулей относительная ошибка остаётся хорошо контролируемой, алгоритм обладает отличной численной стабильностью.

Проверка на примерах

На рисунке 2 представлена визуализация решения уравнения F(w)=0F(w) = 0, метод Ньютона точно локализует нули.

Связанные работы

Историческое развитие

  1. Carpenter (1992) 2: ранний асимптотический анализ нулей обобщённых полиномов Бесселя
  2. Pasquini (2000) 13: численные методы точного вычисления нулей обобщённых полиномов Бесселя
  3. Segura (2013) 14: универсальные методы вычисления комплексных нулей специальных функций
  4. de Bruin и др. (1981) 3,4: исследование области распределения нулей обобщённых полиномов Бесселя

Преимущества данной работы

  • Равномерность: разложения равномерно эффективны для всех нулей
  • Точность: достигается уровень машинной точности
  • Эффективность: сложность алгоритма значительно превосходит существующие методы
  • Применимость: охватывает более широкий диапазон параметров

Заключение и обсуждение

Основные выводы

  1. Теоретический вклад: установлена полная асимптотическая теория нулей обратного обобщённого полинома Бесселя
  2. Численные достижения: разработан высокоточный и эффективный алгоритм вычисления нулей
  3. Практическая ценность: предоставлены надёжные вычислительные инструменты для связанных приложений

Ограничения

  1. Ограничения параметров: метод применим для диапазона параметров −Δ1n+32≤a≤Δ2n-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n
  2. Сложность: вывод асимптотических разложений математически довольно сложен
  3. Требования реализации: требуется поддержка высокоточной арифметики

Направления будущих исследований

  1. Расширение на более общие диапазоны параметров
  2. Исследование аналогичных методов для нулей других специальных функций
  3. Оптимизация вычислительной эффективности алгоритма
  4. Изучение приложений в практических инженерных задачах

Глубокая оценка

Достоинства

  1. Теоретическая строгость: строгий математический вывод, основанный на теории точек поворота и методе WKB
  2. Методологическая инновация: искусное объединение асимптотических разложений и численных алгоритмов
  3. Выдающаяся точность: достижение относительной ошибки порядка 10−1510^{-15}
  4. Высокая эффективность: сложность алгоритма превосходит традиционные методы
  5. Сильная практичность: предоставлены полная реализация и открытый исходный код

Недостатки

  1. Теоретическая сложность: математическая база метода довольно глубока, высокий порог понимания
  2. Зависимость от параметров: эффективность метода зависит от предположений о конкретном диапазоне параметров
  3. Сложность реализации: требуется вычисление нескольких сложных коэффициентов и производных

Влияние

  1. Академическая ценность: предоставляет важные теоретические инструменты для теории специальных функций
  2. Перспективы приложений: потенциальные приложения в обработке сигналов, теории управления и других областях
  3. Методологическое значение: служит примером для исследования аналогичных проблем

Применимые сценарии

  1. Научные вычисления: случаи, требующие высокоточного вычисления нулей полиномов Бесселя
  2. Инженерные приложения: проектирование фильтров, анализ систем управления и т.д.
  3. Теоретические исследования: дальнейшие исследования специальных функций и асимптотического анализа

Список литературы

Данная работа цитирует 14 важных источников, охватывающих классические и передовые работы в областях асимптотического анализа, специальных функций и численных алгоритмов, обеспечивая прочную теоретическую базу для исследования.

Общая оценка: Это высококачественная статья по математическому анализу, достигшая высокого уровня как в теоретическом выводе, так и в численной реализации. Инновационность метода, точность результатов и эффективность алгоритма заслуживают похвалы и предоставляют ценные инструменты для исследований и приложений в соответствующих областях.