Точность переменной использования численного интегрирования
vpaintegral(___,
дополнительные опции использования заданы одним или несколькими Name,Value
)Name,Value
парные аргументы.
Численно интегрируйте символьное выражение x^2
от 1
к 2
.
syms x vpaintegral(x^2, 1, 2)
ans = 2.33333
Численно интегрируйте символьную функцию y (x) = x 2 от 1
к 2
.
syms y(x) y(x) = x^2; vpaintegral(y, 1, 2)
ans = 2.33333
vpaintegral
арифметика переменной точности использования, в то время как MATLAB® integral
функционируйте использует арифметику с двойной точностью. Используя значения по умолчанию допуска, vpaintegral
может обработать значения, которые вызывают MATLAB integral
функционируйте, чтобы переполниться или потерять значимость.
Интегрируйте besseli(5,25*u).*exp(-u*25)
при помощи обоих integral
и vpaintegral
. integral
функция возвращает NaN
и выдает предупреждение в то время как vpaintegral
возвращает правильный результат.
syms u x f = besseli(5,25*x).*exp(-x*25); fun = @(u)besseli(5,25*u).*exp(-u*25); usingIntegral = integral(fun, 0, 30) usingVpaintegral = vpaintegral(f, 0, 30)
Warning: Infinite or Not-a-Number value encountered. usingIntegral = NaN usingVpaintegral = 0.688424
digits
функция не влияет vpaintegral
. Вместо этого увеличьте точность vpainteral
путем уменьшения допусков интегрирования. С другой стороны увеличьте скорость численного интегрирования путем увеличения допусков. Управляйте допуском, используемым vpaintegral
путем изменения относительной погрешности RelTol
и абсолютная погрешность AbsTol
, которые влияют на интегрирование через условие
Численно интегрируйте besselj(0,x)
от 0
к pi
, к 32 значащим цифрам установкой RelTol
к 10^(-32)
. Выключите AbsTol
путем установки его на 0
.
syms x vpaintegral(besselj(0,x), [0 pi], 'RelTol', 1e-32, 'AbsTol', 0)
ans = 1.3475263146739901712314731279612
Используя более низкий допуск значения увеличивает точность за счет скорости.
Интегрируйте 1/(2*z-1)
по треугольному пути от 0
к 1+1i
к 1-1i
назад к 0
путем определения waypoints.
syms z vpaintegral(1/(2*z-1), [0 0], 'Waypoints', [1+1i 1-1i])
ans = - 8.67362e-19 - 3.14159i
Инвертирование направления интеграла, путем изменения порядка waypoints и обмена пределами, изменяет знак результата.
Выполните несколько интегрирование вложенными вызовами vpaintegral
. Объединяться
syms x y vpaintegral(vpaintegral(x*y, x, [1 3]), y, [-1 2])
ans = 6.0
Пределы интегрирования могут быть символьными выражениями или функциями. Интегрируйте по треугольной области 0 ≤ x ≤ 1 и |y | < x путем определения пределов интегрирования по y
в терминах x
.
vpaintegral(vpaintegral(sin(x-y)/(x-y), y, [-x x]), x, [0 1])
ans = 0.89734