Сдвинутая синусоидальная интегральная функция
ssinint(
возвращает сдвинутую синусоидальную интегральную функцию X
)ssinint(X) = sinint(X) — pi/2
.
В зависимости от его аргументов, ssinint
возвращает результаты с плавающей точкой или точные символьные результаты.
Вычислите сдвинутую синусоидальную интегральную функцию для этих чисел. Поскольку эти числа не являются символическими объектами, ssinint
возвращает результаты с плавающей точкой.
A = ssinint([- pi, 0, pi/2, pi, 1])
A = -3.4227 -1.5708 -0.2000 0.2811 -0.6247
Вычислите сдвинутую синусоидальную интегральную функцию для чисел, преобразованных в символические объекты. Для многих символических (точных) чисел, ssinint
возвращает неразрешенные символические вызовы.
symA = ssinint(sym([- pi, 0, pi/2, pi, 1]))
symA = [ - pi - ssinint(pi), -pi/2, ssinint(pi/2), ssinint(pi), ssinint(1)]
Использовать vpa
для аппроксимации символьных результатов с числами с плавающей запятой:
vpa(symA)
ans = [ -3.4227333787773627895923750617977,... -1.5707963267948966192313216916398,... -0.20003415864040813916164340325818,... 0.28114072518756955112973167851824,... -0.62471325642771360428996837781657]
Постройте график сдвинутой интегральной функции синуса на интервале от -4*pi
на 4*pi
.
syms x fplot(ssinint(x),[-4*pi 4*pi]) grid on
Многие функции, такие как diff
, int
, и taylor
, может обрабатывать выражения, содержащие ssinint
.
Найдите первую и вторую производные сдвинутой синусоидальной интегральной функции:
syms x diff(ssinint(x), x) diff(ssinint(x), x, x)
ans = sin(x)/x ans = cos(x)/x - sin(x)/x^2
Найдите неопределенный интеграл сдвинутой синусоидальной интегральной функции:
int(ssinint(x), x)
ans = cos(x) + x*ssinint(x)
Найдите расширение ssinint(x)
серии Тейлора:
taylor(ssinint(x), x)
ans = x^5/600 - x^3/18 + x - pi/2
[1] Gautschi, W. and W. F. Cahill. «Экспоненциальные интегральные и связанные функции». Руководство по математическим функциям с формулами, графиками и математическими таблицами. (М. Абрамовиц и И. А. Штегун, эд.). Нью-Йорк: Дувр, 1972.