bessely

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

Описание

пример

Y = bessely(nu,Z) вычисляет Функцию Бесселя второго доброго Y ν (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')

Вычислите немасштабированное (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')

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')

Входные параметры

свернуть все

Порядок уравнения в виде скаляра, вектора, матрицы или многомерного массива. 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)

Больше о

свернуть все

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

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

z2d2ydz2+zdydz+(z2ν2)y=0.

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

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

Jν(z)=(z2)ν(k=0)(z24)kk!Γ(ν+k+1).

Можно вычислить использование функций Бесселя первого рода besselj.

Функции Бесселя второго доброго, обозначенного Y ν (z), сформируйте второе решение уравнения функции Бесселя, которое линейно независимо от J ν (z). Y ν (z) задан

Yν(z)=Jν(z)cos(νπ)Jν(z)sin(νπ).

Советы

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

Hν(1)(z)=Jν(z)+iYν(z)Hν(2)(z)=Jν(z)iYν(z).

Hν(K)(z) besselh, J ν (z) является besselj, и Y ν (z) является bessely. Функции Ганкеля также формируют основной набор решений уравнения функции Бесселя (см. besselh).

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

Смотрите также

| | |

Представлено до R2006a
Для просмотра документации необходимо авторизоваться на сайте