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

Примеры

свернуть все

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

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)потому что(νπ)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

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