Exponenta Banner

bessely

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

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

Синтаксис

Y = bessely(nu,Z)
Y = bessely(nu,Z,scale)

Описание

пример

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

пример

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

Примеры

свернуть все

Скрипт Open Live Script

Задайте область.

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.

Скрипт Open Live Script

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

Больше о

свернуть все

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

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

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)+i Yν(z)Hν(2)(z)=Jν(z)−i Yν(z).

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

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

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

| | |

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

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

Поддержка