Exponenta Banner
exponenta event 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')

Скрипт 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')

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

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

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

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

| | |

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

Документация MATLAB
Поддержка