Что вычисляет этот калькулятор
Инструмент находит s-й положительный нуль функций Бесселя первого рода \(J_v(x)\) и второго рода \(Y_v(x)\) (её также называют функцией Неймана или Вебера) для вещественного порядка \(v\). Эти нули, обозначаемые \(j_{v,s}\) и \(y_{v,s}\), встречаются повсюду в физике и инженерии: колебания круглой мембраны (барабана), теплопроводность в цилиндрах, электромагнитные волноводы, ряды Фурье — Бесселя. Это чисто математический инструмент для работы со спецфункциями — он универсален и не зависит ни от региона, ни от системы единиц.
Как пользоваться
Введите порядок \(v\) (вещественное число, обычно от 0 до 200) и индекс нуля \(s\) (целое положительное число 1, 2, 3, …). Калькулятор вернёт \(j_{v,s}\) — s-й положительный корень \(J_v(x)\), и \(y_{v,s}\) — s-й положительный корень \(Y_v(x)\). Например, при \(v = 0\), \(s = 1\) вы получите основную моду колебаний барабана.
Формула и метод
Обе функции являются решениями уравнения Бесселя \(x^2 y'' + x y' + (x^2 - v^2) y = 0\). Значения \(J_v\) вычисляются по сходящемуся степенному ряду; для \(Y_v\) используется выражение $$Y_v = \frac{J_v \cos(v\pi) - J_{-v}}{\sin(v\pi)}$$ причём случай целого порядка обрабатывается как численный предел. Чтобы найти s-й нуль, мы берём начальное приближение по асимптотической формуле Мак-Магона, локализуем корень в интервале, а затем уточняем его методом деления отрезка пополам (бисекции) до достижения сходимости. Поскольку \(Y_v(x) \to -\infty\) при \(x = 0\), поиск начинается с небольшого положительного \(x\), чтобы обойти логарифмическую особенность.
Разбор примера
При \(v = 0\), \(s = 1\): первый положительный нуль \(J_0(x)\) равен \(2{,}4048255577\), а первый положительный нуль \(Y_0(x)\) равен \(0{,}8935769663\). При \(v = 1\), \(s = 1\): \(j_{1,1} = 3{,}8317059702\) и \(y_{1,1} = 2{,}1971413260\).
Частые вопросы
Почему первый нуль \(Y_v\) меньше единицы? \(Y_0\) уходит в \(-\infty\) в начале координат и пересекает нуль вблизи \(x \approx 0{,}894\), ещё до своего первого максимума, поэтому её первый нуль гораздо меньше, чем у \(J_0\).
Может ли \(v\) быть нецелым? Да. Определяющая формула для \(Y_v\) напрямую справедлива при любом нецелом \(v\), а целые порядки вычисляются как гладкий предел.
Насколько точны результаты? Расчёты выполняются с двойной точностью, что даёт примерно 10 значащих цифр при умеренных значениях \(v\) и \(s\).
Что означают нули в приложениях
Нули функций Бесселя — это не просто абстрактное любопытство: это дискретные собственные значения, которые появляются всякий раз, когда волновая, тепловая или потенциальная задача ставится на круговой или цилиндрической области. Граничное условие заставляет радиальную часть решения обращаться в нуль (или её производную обращаться в нуль) на границе, и это условие удовлетворяется только при нулях \(j_{v,s}\).
Колеблющаяся круглая мембрана (барабан)
Для идеального круглого барабана радиусом \(a\) с зажатым краем смещение разделяется на моды, помеченные \((v,s)\), где \(v\) подсчитывает угловые узловые диаметры, а \(s\) подсчитывает радиальные узловые окружности. Допустимые собственные частоты имеют вид \(f_{v,s}=\frac{c}{2\pi a}\,j_{v,s}\), где \(c\) — скорость волны. Основной тон использует \(j_{0,1}=2.4048256\); больший \(s\) и больший \(v\) оба повышают высоту звука, и поскольку \(j_{v,s}\) не являются целыми кратными друг друга, обертоны барабана не гармоничны.
Теплопроводность в цилиндре
При решении уравнения теплопроводности в длинном цилиндре с фиксированной температурой на поверхности радиальные собственные функции имеют вид \(J_0(\lambda_s r/a)\) с \(\lambda_s=j_{0,s}\). Каждая мода затухает во времени как \(\exp\!\left(-\alpha (j_{0,s}/a)^2 t\right)\), поэтому наименьший нуль \(j_{0,1}\) определяет самый медленно затухающий, наиболее долгоживущий тепловой профиль. Больший \(s\) даёт большие собственные значения и, следовательно, более быстрое затухание.
Частоты отсечки волновода
В полом круглом металлическом волноводе радиусом \(a\) поперечно-магнитные (ТМ) моды отсекаются на частотах, заданных \(j_{v,s}\), а поперечно-электрические (ТЕ) моды — на нулях производной \(J_v'\). Для ТМ мод отсечка происходит при \(f_{c}=\frac{c\,j_{v,s}}{2\pi a}\); только выше этой частоты мода распространяется. Самая низкая ТМ мода (ТМ\(_{01}\)) снова использует \(j_{0,1}\).
Ряды Фурье–Бесселя
Любая разумная функция на диске может быть разложена как \(f(r)=\sum_{s=1}^{\infty} c_s\,J_v\!\left(j_{v,s}\,r/a\right)\). Масштабированные нули \(j_{v,s}/a\) действуют точно так же, как волновые числа \(n\pi/L\) обычного синус-ряда Фурье, и ортогональность \(J_v(j_{v,s}r/a)\) на диске (с весом \(r\)) позволяет вычислять коэффициенты путём интегрирования, \(c_s=\frac{2}{a^2 J_{v+1}^2(j_{v,s})}\int_0^a f(r)\,J_v\!\left(j_{v,s}r/a\right) r\,dr\).
Как сдвигаются нули
Для фиксированного порядка увеличение индекса \(s\) увеличивает \(j_{v,s}\) шагами, приближающимися к \(\pi\) (осцилляция становится почти периодической вдали от начала координат). Для фиксированного индекса увеличение порядка \(v\) смещает первый нуль наружу примерно как \(j_{v,1}\approx v + 1.8557\,v^{1/3}+\cdots\) при большом \(v\), отражая то, как центрифугальный член \(v^2/x^2\) в уравнении Бесселя задерживает начало осцилляции. Практический вывод: большая угловая сложность (\(v\)) и больше радиальных узлов (\(s\)) оба соответствуют большим собственным значениям и, следовательно, более высоким частотам или более быстрому затуханию.
Для конкретного примера барабана с \(c=100\ \text{м/с}\) и \(a=0.30\ \text{м}\), основная частота равна \(f_{0,1}=\frac{100}{2\pi(0.30)}\,(2.4048256)\approx 127.6\ \text{Гц}\). Первый обертон использует \(j_{1,1}=3.8317060\), давая \(\approx 203\ \text{Гц}\), поэтому отношение обертона к основному тону равно \(j_{1,1}/j_{0,1}\approx 1.593\) — слышно не октава и не квинта, поэтому барабаны звучат не транспонированно по сравнению со струнами.
