besseli

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

Описание

пример

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

пример

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

Примеры

свернуть все

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

z = 0:0.01:5;

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

I = zeros(5,501);
for nu = 0:4
    I(nu+1,:) = besseli(nu,z);
end

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

plot(z,I)
axis([0 5 0 8])
grid on
legend('I_0','I_1','I_2','I_3','I_4','Location','NorthWest')
title('Modified Bessel Functions of the First Kind for $\nu \in [0,4]$','interpreter','latex')
xlabel('z','interpreter','latex')
ylabel('$I_\nu(z)$','interpreter','latex')

Figure contains an axes object. The axes object with title Modified Bessel Functions of the First Kind for nu in bracketleft 0 , 4 bracketright contains 5 objects of type line. These objects represent I_0, I_1, I_2, I_3, I_4.

Вычислите масштабированную модифицированную функцию Бесселя первого рода Iν(z)e-|Re[Z]| для значений z в интервале [0,20] и для порядков ν между 0 и 3.

z = linspace(0,20);
scale = 1;
Is = zeros(4,100);
for nu = 0:3
  Is(nu+1,:) = besseli(nu,z,scale);
end

Постройте все функции на том же рисунке. Для больших значений z, масштабированные функции не переполняют пределов двойной точности, расширяя их область значений исчисляемости по сравнению с немасштабированными функциями.

plot(z,Is)
legend('I_0','I_1','I_2','I_3')
title('Scaled Mod. Bessel Functions of the First Kind for $\nu \in \left[0, 3 \right]$','interpreter','latex')
xlabel('z','interpreter','latex')
ylabel('$e^{-|{z}|} \cdot I_\nu(z)$','interpreter','latex')

Figure contains an axes object. The axes object with title Scaled Mod. Bessel Functions of the First Kind for nu in bracketleft 0 , 3 bracketright contains 4 objects of type line. These objects represent I_0, I_1, I_2, I_3.

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

свернуть все

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

Пример: besseli(3,Z)

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

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

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

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

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

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

  • 1 — Масштабируйте выход besseli exp(-abs(real(Z)))

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

Пример: besseli(nu,Z,1)

Больше о

свернуть все

Модифицированные функции Бесселя

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

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

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

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

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

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

Kν(z)=(π2)Iν(z)Iν(z)sin(νπ).

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

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

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

| | | |

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