Exponenta Banner
exponenta event banner

bessely

Функция Бесселя второго рода

свернуть все на странице

Синтаксис

Y = бессильно (nu, Z)
Y = bessely (ню, Z, масштаб)

Описание

пример

Y = bessely(nu,Z) вычисляет функцию Бесселя второго рода Ystart( z) для каждого элемента в массиве Z.

пример

Y = bessely(nu,Z,scale) определяет необходимость экспоненциального масштабирования функции Бесселя второго рода во избежание переполнения или потери точности. Если scale является 1, то выход bessely масштабируется на коэффициент exp(-abs(imag(Z))).

Примеры

свернуть все

Открыть сценарий в реальном времени

Определите домен.

z = 0:0.1:20;

Вычислите первые пять функций Бесселя второго рода. Каждая строка Y содержит значения одного порядка функции, вычисленные в точках в z.

Y = zeros(5,201);
for i = 0:4
    Y(i+1,:) = bessely(i,z);
end

Постройте график всех функций на одном рисунке.

plot(z,Y)
axis([-0.1 20.2 -2 0.6])
grid on
legend('Y_0','Y_1','Y_2','Y_3','Y_4','Location','Best')
title('Bessel Functions of the Second Kind for $\nu \in [0, 4]$','interpreter','latex')
xlabel('z','interpreter','latex')
ylabel('$Y_\nu(z)$','interpreter','latex')

Figure contains an axes. The axes with title Bessel Functions of the Second Kind for $\nu \in [0, 4]$ contains 5 objects of type line. These objects represent Y_0, Y_1, Y_2, Y_3, Y_4.

Открыть сценарий в реальном времени

Вычислите немасштабированное значение (Y) и масштабированы (Ys) функция Бесселя второго рода Y2 (z) для комплексных значений z .

x = -10:0.35:10;
y = x';
z = x + 1i*y;
scale = 1;
Y = bessely(2,z);
Ys = bessely(2,z,scale);

Сравните графики мнимой части масштабированных и немасштабированных функций. Для больших значений abs(imag(z)), немасштабированная функция быстро перетекает пределы двойной точности и перестает быть вычисляемой. Масштабированная функция удаляет это доминирующее экспоненциальное поведение из вычисления и, таким образом, имеет больший диапазон вычисляемости по сравнению с немасштабированной функцией.

surf(x,y,imag(Y))
title('Bessel Function of the Second Kind','interpreter','latex')
xlabel('real(z)','interpreter','latex')
ylabel('imag(z)','interpreter','latex')

Figure contains an axes. The axes with title Bessel Function of the Second Kind contains an object of type surface.

surf(x,y,imag(Ys))
title('Scaled Bessel Function of the Second Kind','interpreter','latex')
xlabel('real(z)','interpreter','latex')
ylabel('imag(z)','interpreter','latex')

Figure contains an axes. The axes with title Scaled Bessel Function of the Second Kind contains an object of type surface.

Входные аргументы

свернуть все

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

Пример: bessely(3,0:5)

Типы данных: single | double

Функциональная область, заданная как скаляр, вектор, матрица или многомерный массив. bessely является действительным, где Z является положительным. nu и Z должен быть одинакового размера, или один из них может быть скалярным.

Пример: bessely(1,[1-1i 1+0i 1+1i])

Типы данных: single | double
Поддержка комплексного номера: Да

Переключение на функцию масштабирования, заданную как одно из следующих значений:

  • 0 (по умолчанию) - без масштабирования

  • 1 - Масштабировать выход bessely около exp(-abs(imag(Z)))

На комплексной плоскости величина bessely быстро растет по мере того, как значение abs(imag(Z)) увеличивается, поэтому экспоненциальное масштабирование выходного сигнала полезно для больших значений abs(imag(Z)) где результаты в противном случае быстро теряют точность или переполняют пределы двойной точности.

Пример: bessely(3,0:5,1)

Подробнее

свернуть все

Функции Бесселя

Это дифференциальное уравнение, где start- вещественная константа, называется уравнением Бесселя:

z2d2ydz2 + zdydz + (z2 − start2) y = 0.

Его решения известны как функции Бесселя.

Функции Бесселя первого рода, обозначенные Jstart( z) и J-start( z), образуют фундаментальный набор решений уравнения Бесселя для неинтегрирующих start. Jstart( z) определяется

Jstart( z) = (z2) ν∑ (k = 0) ∞ (− z24) kk! Γ (start+ k + 1).

Функции Бесселя первого рода можно вычислить с помощью besselj.

Функции Бесселя второго рода, обозначенные Ystart( z), образуют второе решение уравнения Бесселя, линейно не зависящее от Jstart( z). Ystart( z) определяется как

Y, (z) = J, (z) cos (startδ) − J −

Совет

Функции Бесселя связаны с функциями Ганкеля, также называемыми функциями Бесселя третьего рода:

H, (1) (z) = J, (z) + i Y, (z) H, (2), (z) = J, (z) − i, Y, (z).

Hstart( K) (z) =besselh, Jstart( z) являетсяbesselj, и Ystart( z) являетсяbessely. Функции Ганкеля также образуют фундаментальный набор решений уравнения Бесселя (см. besselh).

Расширенные возможности

..

См. также

| | |

Представлен до R2006a

Документация MATLAB

Поддержка