Вычисление функции интеграла вероятности и ей обратной
В статье предложены алгоритмы для эффективного вычисления функции интеграла вероятностей и ей обратной, основанные на разложении этих функций в ряд невязок по полиномам Эрмита. Для вычисления квантилей интеграла вероятностей получено разложение функций в ряд невязок на основе обращения разложения фу...
Saved in:
| Date: | 2004 |
|---|---|
| Main Authors: | , |
| Format: | Article |
| Language: | Russian |
| Published: |
Інститут проблем математичних машин і систем НАН України
2004
|
| Subjects: | |
| Online Access: | https://nasplib.isofts.kiev.ua/handle/123456789/83909 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Journal Title: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Cite this: | Вычисление функции интеграла вероятности и ей обратной / Г.С. Теслер, Зы Хак Зунг // Мат. машини і системи. — 2004. — № 3. — С. 31-40. — Бібліогр.: 26 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1859904872864808960 |
|---|---|
| author | Теслер, Г.С. Зы Хак Зунг |
| author_facet | Теслер, Г.С. Зы Хак Зунг |
| citation_txt | Вычисление функции интеграла вероятности и ей обратной / Г.С. Теслер, Зы Хак Зунг // Мат. машини і системи. — 2004. — № 3. — С. 31-40. — Бібліогр.: 26 назв. — рос. |
| collection | DSpace DC |
| description | В статье предложены алгоритмы для эффективного вычисления функции интеграла вероятностей и ей обратной, основанные на разложении этих функций в ряд невязок по полиномам Эрмита. Для вычисления квантилей интеграла вероятностей получено разложение функций в ряд невязок на основе обращения разложения функции интеграла вероятностей, что требует значительно меньших затрат, чем разложение в ряд Тейлора. На основе этого разложения получены итерационные формулы произвольного порядка сходимости. Особое внимание уделено начальным приближениям на основе сегментной аппроксимации, вопросам табулирования функций и вычисления взаимно-обратных функций. Представлены материалы моделирования предложенных алгоритмов на компьютере.
В статті запропоновані алгоритми для ефективного обчислення функції інтеграла імовірностей та її оберненої, які базуються на розкладі цих функцій у ряд за нев’язками по поліномах Ерміта. Для обчислення квантилей інтеграла імовірностей отримано розклад функцій у ряд нев’язок на основі обернення розкладу функції інтеграла імовірностей, що потребує значно менших витрат, ніж розклад у ряд Тейлора. На базі цього розкладу отримані ітераційні формули довільного порядку збіжності. Особливу увагу присвячено початковим наближенням на основі сегментної апроксимації, питанням табулювання функцій та обчислення взаємно-обернених функцій. Представлені матеріали моделювання запропонованих алгоритмів на комп’ютері.
In the article therehase been proposed the algorithms for effective calculating probabilities integral function and the reverse bit founded on decomposing these functions in row of closure errors on polynomials Erpith. For calculating probabilities integral quantiles there has been received the functions decomposing in row of closure errors on the basis of converting the probabilities integral functions uhat requires considerably lesser expenses than Jailor’s decomposing in raw/ On the basis of this decomposing there have been received the iterative formulae of cowergency arbitrary arrangement. Special attention has been spared to initial estimates on the basis of segment approximation? To the problems of tabulating the functions and computing mutually-inverse functions. Shere have been represented the materials of modelling the proposed algorithms on a computer.
|
| first_indexed | 2025-12-07T15:59:18Z |
| format | Article |
| fulltext |
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
31
УДК 518.5
Г.С. ТЕСЛЕР, ЗЫ ХАК ЗУНГ
ВЫЧИСЛЕНИЕ ФУНКЦИИ ИНТЕГРАЛА ВЕРОЯТНОСТИ И ЕЙ ОБРАТНОЙ
Abstract: In the article therehase been proposed the algorithms for effective calculating probabilities integral function
and the reverse bit founded on decomposing these functions in row of closure errors on polynomials Erpith. For
calculating probabilities integral quantiles there has been received the functions decomposing in row of closure errors on
the basis of converting the probabilities integral functions uhat requires considerably lesser expenses than Jailor’s
decomposing in raw/ On the basis of this decomposing there have been received the iterative formulae of cowergency
arbitrary arrangement. Special attention has been spared to initial estimates on the basis of segment approximation? To
the problems of tabulating the functions and computing mutually-inverse functions. Shere have been represented the
materials of modelling the proposed algorithms on a computer.
Key words: probabilities integral, quantile, approximation, modelling, error.
Анотація: В статті запропоновані алгоритми для ефективного обчислення функції інтеграла імовірностей
та її оберненої, які базуються на розкладі цих функцій у ряд за нев’язками по поліномах Ерміта. Для обчислення
квантилей інтеграла імовірностей отримано розклад функцій у ряд нев’язок на основі обернення розкладу
функції інтеграла імовірностей, що потребує значно менших витрат, ніж розклад у ряд Тейлора. На базі цього
розкладу отримані ітераційні формули довільного порядку збіжності. Особливу увагу присвячено початковим
наближенням на основі сегментної апроксимації, питанням табулювання функцій та обчислення взаємно-
обернених функцій. Представлені матеріали моделювання запропонованих алгоритмів на комп’ютері.
Ключові слова: інтеграл імовірностей, квантиль, апроксимація, моделювання, похибки.
Аннотация: В статье предложены алгоритмы для эффективного вычисления функции интеграла
вероятностей и ей обратной, основанной на разложении этих функций в ряд невязок по полиномам Эрмита.
Для вычисления квантилей интеграла вероятностей получено разложение функций в ряд невязок на основе
обращения разложения функции интеграла вероятностей, что требует значительно меньших затрат, чем
разложение в ряд Тейлора. На основе этого разложения получены итерационные формулы произвольного
порядка сходимости. Особое внимание уделено начальным приближениям на основе сегментной аппроксимации,
вопросам табулирования функций и вычисления взаимно-обратных функций. Представлены материалы
моделирования предложенных алгоритмов на компьютере.
Ключевые слова: интеграл вероятностей, квантиль, аппроксимация, моделирование, погрешность.
1. Введение
Интеграл вероятности (функция Крампа, интеграл ошибок) и близкие к ним функции являются одними
из наиболее используемых в теории и практике специальных функций.
Это связано с тем, что данные функции непосредственно связаны с нормальным (гауссовским)
распределением случайной величины, которые во многих случаях являются адекватной моделью
многих процессов, происходящих в живой и искусственной природе. Они широко используются в
теории вероятности, математической статистике, теоретической и математической физике,
описывающей процессы теплопроводности, диффузии, броуновского движения, кинетической теории
газов и т.д., в теории связи (при исследовании влияния шума на полезный сигнал), при решении
дифференциальных и интегральных уравнений методами Монте Карло, в теории чисел
(распределение значений “коротких” чисел Вейля и других разделах теоретической и прикладной
математики, а также в теории надёжности и многих других приложениях).
Описание функции интеграла вероятностей либо их таблиц включено в большинство
известных справочников и монографий, посвящённых описанию специальных математических
функций [1 – 6]. Имеется обширная литература с подробными таблицами данных функций [7, 8]. В [8]
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
32
приведена обширная библиография по этой тематике. Имеются также отдельные монографии,
посвящённые функции интеграла вероятности [9, 26], и множество таблиц.
В работах [10 – 13] приведены различные методы и алгоритмы аппроксимации функции
вероятности, а также обширная библиография.
Особо хочется отметить тот факт, что в шестидесятые годы прошлого столетия в работах [14 –
15] ввели в рассмотрение функцию ( )xinverfy = , являющуюся обратной по отношению функции
интеграла вероятности. Это не означает, что ранее в неявном виде данные функции использовались
под названием квантиля распределения и процентных точек распределения.
Для эффективного вычисления как прямых, так и обратных функций на ЭВМ, Г.С. Теслером
было предложено использовать адаптивные аппроксимации, основанные на разложении функций по
невязкам (см. справочник [13] и приведенную в нём библиографию).
В настоящей статье предлагается дальнейшее развитие этих исследований для функций
( ) ∫
−=
x
t dtexerf
0
2
π
2
, ∫
−=Φ
x
t dtex
0
2
0
2
2
1
)(
π
и им обратным. Отметим, что ( ) )2(2 0 xxerf Φ= .
В свою очередь функция
)(
2
1
)]2/(1[
2
1
2
1
)( 0
2/2
xxerfdtex
x
t Φ+=−==Φ ∫
∞−
−
π
.
2. Вычисление ( )xerf и )(0 xΦ
Разложение ( )yxerf + в ряд Тейлора имеет вид
∑∑
∞
=
+
=+
+
=
∞
= +
+==+
0
1
1
1
0 )!1(
.)()(
!
.)()(
k
k
xyk
kk
xy
k
k
k
k
y
yerf
dy
d
xerf
k
y
yerf
dy
d
yxerf .
Учтём, что производная интеграла вероятностей
;
2)(
)(
2xe
dx
xderf
x −==
π
ϕ
...;2,1,0),().()1().(
2
)1()()(
2
1
1
)( =−=−== −
+
+
nxxHexHxerf
dx
d
x n
nx
n
n
n
n
n ϕ
π
ϕ
),()(,0)()()( )0(11 xxxnxxx nnn ϕϕϕϕϕ ==++ −+
)(xH n – многочлены Эрмита.
Отметим, что
.)1()(
22 x
n
n
xn
n e
dx
d
exH −−=
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
33
Произведя несложные преобразования и учтя, что невязка xxxyinverfz −=−= 000 )( и 0y –
приближённое значение ( )xerfy = на заданном интервале, получим
.
)!1(
).(2
)!1(
)).((
.
2
))(()(
0
1
00
0
0
1
00))((
000
2
0
2
0
∑
∑
∞
=
+
−
∞
=
+
−
+
−=
=
+
−=−=
k
k
kx
k
k
kyinverf
k
zxH
ey
k
zyinverfH
eyzyinverferfxerf
π
π
(1)
Многочлены Эрмита могут быть вычислены на основе рекуррентных соотношений
)(2)(2)( 11 xnHxxHxH nnn −+ −= , где .2)(,1)( 10 xxHxH ==
Для уменьшения погрешности за счёт резкого возрастания порядков чисел ряда (1) перейдём к
другому рекуррентному соотношению:
)),()((
1
2
)( 11 xxx
k
x kkk −+ −
+
= γγγ (2)
где .
!
)(
)(,...,2)(,1)( 10 k
xH
xxxx k
k === γγγ
В результате формула (1) может вычисляться на основе следующего рекуррентного соотношения:
.
1
).(.
2 1
0
01
0
+
−=
+
−
+ k
z
xeyy
k
k
x
kk γ
π
(3)
Для вычисления погрешности при оставлении в ряде (1) m членов следует подставить в
невязку величину ,00 ∆+= yy где 0∆ – абсолютная погрешность приближения на заданном
интервале. В результате получим [13]:
1
0
101)(
0 )(
)!1(
)(
.
2 2
0
2
0 +++− ∆
2+
≈∆ mmymxm e
m
xH
e
π
π
.
Начальные приближения целесообразно выбирать на основе сегментной (разновидность
сплайн) аппроксимации. При использовании постоянных приближений на соответствующем
подынтервале с учётом соответствующих норм погрешностей для ],[ 1+∈ ii xxx будет [13 ,16]:
)
2
(
;2))()((
1
0
10
+
+
+
=
+=
ii
i
iii
xx
erfy
xerfxerfy
либо ).(xthyi =
При программной реализации алгоритмов на ЭВМ необходимая точность вычислений задается
в виде ,10.5,0 αε −= где степень α определяет количество верных значащих цифр результата. Тогда
выбираем условие прекращения:
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
34
ε≤−+−
+
+−
1
11
k
kk
k
kk
y
yy
y
yy
или .
1
).(.
2
1
1
1
0
2
0 εγ
π
≤
+∑
−=
+
+
−
k
ki
i
i
i
x y
i
z
xe (4)
Для расчета по формуле (3) задается начальное приближение ),( 00 yx . В частности, можно задать
).0,0( 00 == yx Для стабилизации времени счёта при любых значениях аргумента вводится система
опорных точек. Диапазон изменения аргумента x , начиная с нуля, разбивается на интервалы с шагом
h . В качестве опорной точки в каждом интервале берем его середину и вычисляем )( ii xerfy = ,
причем ).
2
1
(
2
1
)1( −=+−= ihhihxi Из формул (3) и (4) видно, что количество шагов итераций
зависит от h .
Для примера приведем две системы опорных точек ( 5,0=h и 25,0=h ) и их соответствующие
количества шагов итераций при программной реализации.
Таблица 1. Система опорных Таблица 2. Система опорных точек
точек при 5,0=h при 25,0=h
0x 0y
0x 0y 0x 0y
0,25 0,276326390168236933 0,125 0,140316204801333839 2,625 0,999794624263858738
0,75 0,711155633653515132 0,375 0,404116909434822325 2,875 0,999952145160256182
1,25 0,922900128256458230 0,625 0,623240882188417924 3,125 0,999990103270201730
1,75 0,986671671219182444 0,875 0,784075061059859613 3,375 0,999998184723399535
2,25 0,998537283413318848 1,125 0,888388231701707819 3,625 0,999999704859807446
2,75 0,999899378077880363 1,375 0,948170072782090378 3,875 0,999999957486055875
3,25 0,999995697225363250 1,625 0,978443733239983643 4,125 0,999999994576599199
3,75 0,999999886272743430 1,875 0,991990057670119949 4,375 0,999999999387516704
4,25 0,999999998149425820 2,125 0,997345970640517631 4,625 0,999999999938783840
4,75 0,999999999981514951 2,375 0,999217061782108854 4,875 0,999999999994586539
… … …
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
35
Рис. 1. График количества шагов итерации (линия соединения точек) при 5,0=h и 15 ,10 ,7=α
Рис. 2. График количества шагов итерации (линия соединения точек) при 25,0=h и 15 ,10 ,7=α
3. Вычисление )(xinverf и )(0 xinvΦ
В работах [13, 16] приведены методы получения итерационных формул произвольного порядка
сходимости на основе соответствующих функциональных уравнений и методов дифференцирования,
что для получения итерационных формул высоких порядков представляется весьма трудоёмкой
задачей. Итерационная формула второго порядка:
.).)((
2
1
iy
iii exyerfyy −
2
−=+
π
Итерационная формула третьего порядка:
..))((.).)((
22 22
1
ii y
ii
y
iii exyerfyexyerfyy −
4
+−
2
−=+
ππ
(5)
Однако более просто получить эти формулы за счёт обращения ряда (1).
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
36
В основе обоснования такого типа обращений в соответствии с [17] рассмотрим решение
уравнений рядами. Предположим, что функция )(xfy = представлена в неявном виде:
0),( =yxF (6)
и ),( yxF в окрестности точки ),( 00 yx разлагается в ряд по степеням 0xx − и 0yy − , причём
постоянный член в нём равен нулю, а коэффициент при 0yy − отличен от нуля. Тогда и функция
)(xfy = , определяемся уравнением (6) в окрестности указанной точки, также разлагается в ряд по
степеням 0xx − вблизи 0xx = , т.е. если функция ),( yxF является аналитической в точке ),( 00 yx ,
то и функция )(xfy = , определяемая уравнением (6), будет аналитической в точке 0x .
Отсюда следует, что функция )(xfy = в некоторой окрестности точки 0xx − представляется
рядом:
...)(...)()( 0
2
02010 +−++−+−=− n
n xxaxxaxxayy , (7)
где 0y – свободный член, то при 01 ≠a в окрестности 0yy = , из (7) величина х определяется как
функция от y , разлагающаяся, в свою очередь, в ряд по степеням 0yy − . Таким образом, если y
является аналитической функцией от х в точке 0x , то в соответствующей точке 0y (при выполнении
указанного условия) обратная функция будет аналитической).
Так как )(xerf является функцией аналитической, то можно обратить степенной ряд (1).
Техника такого обращения хорошо освещена в работе [13]. В современных условиях для применения
таких преобразований целесообразно использовать системы компьютерной алгебры [19 – 24].
В работах [23, 25] рассмотрен несколько иной подход, чем было указано выше. Так, вместо
невязки 00 xxz −= , удовлетворяющей функциональному уравнению 0)(),( =−= yfxyxF ,
используется невязка )('))(( 000 yfyfxz −= , удовлетворяющая функциональному уравнению
0)('))((),( =−= yfyfxyxF , так что в конечном счёте итерационные формулы для обоих видов
невязок совпадают с точностью перегруппировки членов. Однако для получения подобных формул уже
используются прямые методы разложения по невязкам. В этом случае для )(xinverf невязка имеет
вид .2)).((
2
0
00
yeyerfxxz −=
Очевидно, такой тип невязок наиболее подходит при вычислении интегральных функций.
Используя первый вид невязок, т.е. )( 00 yerfxz −= на основе обращения ряда невязок (1)
(таблица) и подставляя значения полиномов Эрмита и невязки получим
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
37
....
!3
).)(14(2
!2
).(2
.))(()(
3
00
2
00
0000
2
0
2
0
2
0
+2
+
−2+
+
2
−=−==
zeyzey
zeyzyerfinverfxinverfy
yy
y
ππ
π
Для невязки
2
0)).((
2 00
yeyerfx
x
z −= получим [23]
....)..
3
4
3
1
(. 3
0
2
0
2
0000 +++++= zyzyzyy
Погрешности при оставлении двух, трёх и т.д. членов составляют [13, 23]
.).
15
14
5
38
5
24
(.).712(
6
1
;.).54(
2
1
).
3
1
3
4
(
;).
3
2
3
2
(.
524
4
2
4
4232
3
322
2
∆−−+∆+−≈∆
∆++∆+≈∆
∆−−∆−≈∆
yyyy
yyy
yy
Отметим, что погрешности для первого и второго подходов выбора невязок совпадают.
Для сравнения полученных результатов с известными приведём результаты алгоритма
H.E. Fetts [13] для вычисления )(xinverfcy = по формулам
( ) ( )
),
1
2/3
1
2/2
11
(),(
;)),(ln(
));)ln((ln(
222
21
1
21
0
tttt
xxyF
yxFy
xy
nn
π
π
=
−=
−−=
+
где .1 yt =
Необходимое количество членов непрерывной дроби i и итераций k при точности 12
десятичных знаков:
x 610−
410− 210− 0,05 0,1 0,2 0,3
i 29 33 51 74 98 150 220
k 7 8 11 14 26 21 27
Отметим, что для значительно большей разрядности, достигающей 30 десятичных разрядов,
при счёте по предложенным итерационным формулам количество итераций колеблется в диапазоне 4
– 7 для всего интервала изменения аргумента ]5,0[∈x , что ведёт к значительно меньшему времени
счёта по сравнению с вышеприведенным алгоритмом.
+ + + +…
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
38
При программной реализации итерационной формулы (5) сходимость вычислений тесно
связана с начальным приближением. Авторы предлагают следующее: для маленьких значений
аргумента, примерно 2,0≤x , в качестве начального приближения следует выбирать
логарифмическую формулу )1ln( 2
0 xy −−−= , а для остальных значений аргумента, воспользуемся
линейной сегментной аппроксимацией: DtCy +−= )1(0 , где коэффициенты C и D определяются в
зависимости от интервала изменения аргумента x .
Таблица 3. Коэффициенты линейной сегментной аппроксимации функции ( )xinverf
Описанный в статье подход для вычисления интеграла вероятностей нормального
распределения и соответствующих квантилей может быть использован не только для вычисления
отдельных значений этих характеристик случайного процесса с фиксированной разрядностью, но и
для вычисления их с произвольной разрядностью и табулирования функций. В последнем случае
время вычисления отдельных значений значительно сокращается, а начальные приближения при
вычислении интеграла вероятности могут быть использованы в качестве опорных значений, что
гарантирует невозрастание погрешностей округления.
Помимо этого, описанный в статье подход для вычисления интеграла вероятностей и
квантилей позволяет, когда это необходимо, порождать соответствующие алгоритмы, адаптивные к
внутренним и внешним условиям применения, что обеспечит повышение уровня их эффективности.
Напомним, что квантилем t называют значение аргумента x функции распределения,
которому с заданной вероятностью Р соответствует условие tx ≤ , т.е. выполняется равенство
∫
∞−
=
t
Pdxxf ,)(
где )(xf – функция плотности распределения.
Интервал изменения
аргумента x
C D
[ )2101;0 −− 1,8398 1,5559
[ )42 101;101 −− −− 93,902 2,5665
[ )64 101;101 −− −− 7149,4 3,3124
[ )86 101;101 −− −− 599310 3,9270
[ )108 101;101 −− −− 5,2592. 710 4,4620
[ )1210 101;101 −− −− 4,7147. 910 4,9399
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
39
В нашем случае квантиль t адекватен функциям )(xinverf либо ).(0 xinvΦ
4. Проверка согласованности результатов вычислений
Для проверки согласованности вычисляем величины:
xxinverferfпр −=∆ ))(( и xxerfinverfобр −=∆ ))(( при различных значениях аргумента x и
точности вычислений α . Напомним, что .10.5,0 αε −= Результаты такой проверки приведены на рис. 3
и 4
а) при 0 и 9 ≈= xα б) при 1 и 9 ≈= xα в) при 1 и 12 ≈= xα
Рис. 3. Результаты вычислений величины xxinverferfпр −=∆ ))((
а) при 0 и 12 ≈= xα б) при 1 и 12 ≈= xα в) при 1 и 15 ≈= xα
Рис. 4. Результаты вычислений величину xxerfinverfобр −=∆ ))((
5. Выводы
В статье приведены алгоритмы вычисления функции интеграла вероятностей для нормального
распределения и квантилей на основе разложения этих функций в ряд невязок и приведены
результаты их программной реализации.
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
40
Для получения алгоритмов при вычислении функции, обратной интегралу вероятностей,
использован эффективный метод обращения функциональных рядов, который является значительно
менее затратным, чем прямые методы, требующие использования обратных производных.
Показана перспективность предложенного подхода для табулирования функции интеграла
вероятностей и ей обратной, а также для их вычисления с произвольной разрядностью и порождения
адаптивных алгоритмов к внутренним и внешним условиям применения.
СПИСОК ЛИТЕРАТУРЫ
1. Справочник по специальным функциям с формулами, графиками и математическими таблицами / Под ред. М.
Абрамовица и И. Стиган: Пер. с англ. – М.: Наука, 1979. – 832 с.
2. Градштейн И.С. , Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. – М.: Физматгиз, 1963. – 1100с.
3. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. – М.: Физматгиз, 1978. – 832с.
4. Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся втузов. – М.: Наука,
1981. – 720 с.
5. Янке Е., Эмде Ф., Лёш Ф. Специальные функции. Формулы, графики, таблицы. – М.: Наука, 1964. – 344с.
6. Прудников А.П., Брычков Ю.А., Марычев О.И. Интегралы и ряды. Специальные функции. – М.: Наука, 1983. –
752 с.
7. Таблицы нормального интеграла вероятностей, нормальной плотности и её нормированных производных / Под
ред. чл. – корр. АН СССР Н.В. Смирнова. – М.: Издательство АН СССР, 1960. – 136 с.
8. Большев Л.Н., Смирнов Н.В. Таблицы математической статистики. – М.: Наука, 1983. – 416 с.
9. Хаджи П.И. Функция вероятности. Интегралы, ряды и некоторые обобщения. – Кишинев: Издательство АН
Молдавской ССР, 1971. – 398 с.
10. Fike C.T. Computer evaluation of mathematical function. – New Jersey: Prentice – Hall, 1968. – 228 p.
11. Computer approximations. – Ed. by J.F. Hart. New York: Wiley, 1968. – 344 p.
12. Люк Ю. Специальные математические функции и их аппроксимации. – М.: Мир. – 1980. – 608 с.
13. Попов Б.А. , Теслер Г.С. Вычисление функций на ЭВМ: Справочник. – Киев: Наукова думка, 1984. – 600 с.
14. Philip J.R. The function inverfθ . “ Austral. J. Phys.” 13. – 1960. – P. 13 – 20.
15. Strecok A.J. On the calculation of the inverse of the error function. “ Math. Comput.” 22. – 1968. – P. 144 – 158.
16. Попов Б.А., Теслер Г.С. Приближение функций для технических приближений. – Киев: Наукова думка, 1980. –
356 с.
17. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления. – М.: Физматлит, 1962. – 807 с. (С.
501).
18. Фильчаков П.Ф. Численные и графические методы прикладной математики: Справочник. – Киев: Наукова
думка, 1970. – 800 с.
19. Heck A. Introduction to Maple. – N.Y. : Springer – Verlag, 1993. – 498 p.
20. Maeder R.E. Programming in Mathematica. – Redwood City, California: Addison – Weslagerg, 1990. – 267 p.
21. Hearn A.C. REDUCE User Manual. – Santa Monica: The Rand Corporation, 1995. – 480 p.
22. Misniewski M., Zero A. MathCAD Plus 6.0. – Warszawa: EXIT, 1998. – 322 p.
23. Попов Б.О. Розв’язування математичних задач у системах комп’ютерної алгебри. – MAPLEV. – Київ: V i P, 2001.
– 312 с.
24. Морозов А.А., Клименко В.П., Фишман Ю.С. АНАЛИТИК – 2000 // Математические машины и системы. – 2001.
№ 1. – C. 66 – 99.
25. Попов Б.О., Лушник О.І. Побудова оптимального алгоритму для обчислення функції, оберненої до інтегралу
імовірності // Волинський мат. вісник. – 1999. – Вип. 6. – C. 107 – 112.
26. Митропольский А.К. Интеграл вероятностей. – Л.: Изд-во Ленинградского университета, 1972. – 86 с.
|
| id | nasplib_isofts_kiev_ua-123456789-83909 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 1028-9763 |
| language | Russian |
| last_indexed | 2025-12-07T15:59:18Z |
| publishDate | 2004 |
| publisher | Інститут проблем математичних машин і систем НАН України |
| record_format | dspace |
| spelling | Теслер, Г.С. Зы Хак Зунг 2015-06-28T09:42:31Z 2015-06-28T09:42:31Z 2004 Вычисление функции интеграла вероятности и ей обратной / Г.С. Теслер, Зы Хак Зунг // Мат. машини і системи. — 2004. — № 3. — С. 31-40. — Бібліогр.: 26 назв. — рос. 1028-9763 https://nasplib.isofts.kiev.ua/handle/123456789/83909 518.5 В статье предложены алгоритмы для эффективного вычисления функции интеграла вероятностей и ей обратной, основанные на разложении этих функций в ряд невязок по полиномам Эрмита. Для вычисления квантилей интеграла вероятностей получено разложение функций в ряд невязок на основе обращения разложения функции интеграла вероятностей, что требует значительно меньших затрат, чем разложение в ряд Тейлора. На основе этого разложения получены итерационные формулы произвольного порядка сходимости. Особое внимание уделено начальным приближениям на основе сегментной аппроксимации, вопросам табулирования функций и вычисления взаимно-обратных функций. Представлены материалы моделирования предложенных алгоритмов на компьютере. В статті запропоновані алгоритми для ефективного обчислення функції інтеграла імовірностей та її оберненої, які базуються на розкладі цих функцій у ряд за нев’язками по поліномах Ерміта. Для обчислення квантилей інтеграла імовірностей отримано розклад функцій у ряд нев’язок на основі обернення розкладу функції інтеграла імовірностей, що потребує значно менших витрат, ніж розклад у ряд Тейлора. На базі цього розкладу отримані ітераційні формули довільного порядку збіжності. Особливу увагу присвячено початковим наближенням на основі сегментної апроксимації, питанням табулювання функцій та обчислення взаємно-обернених функцій. Представлені матеріали моделювання запропонованих алгоритмів на комп’ютері. In the article therehase been proposed the algorithms for effective calculating probabilities integral function and the reverse bit founded on decomposing these functions in row of closure errors on polynomials Erpith. For calculating probabilities integral quantiles there has been received the functions decomposing in row of closure errors on the basis of converting the probabilities integral functions uhat requires considerably lesser expenses than Jailor’s decomposing in raw/ On the basis of this decomposing there have been received the iterative formulae of cowergency arbitrary arrangement. Special attention has been spared to initial estimates on the basis of segment approximation? To the problems of tabulating the functions and computing mutually-inverse functions. Shere have been represented the materials of modelling the proposed algorithms on a computer. ru Інститут проблем математичних машин і систем НАН України Обчислювальні системи Вычисление функции интеграла вероятности и ей обратной Обчислення функції інтеграла імовірностей та її оберненої Computing probability integral function and reverse Article published earlier |
| spellingShingle | Вычисление функции интеграла вероятности и ей обратной Теслер, Г.С. Зы Хак Зунг Обчислювальні системи |
| title | Вычисление функции интеграла вероятности и ей обратной |
| title_alt | Обчислення функції інтеграла імовірностей та її оберненої Computing probability integral function and reverse |
| title_full | Вычисление функции интеграла вероятности и ей обратной |
| title_fullStr | Вычисление функции интеграла вероятности и ей обратной |
| title_full_unstemmed | Вычисление функции интеграла вероятности и ей обратной |
| title_short | Вычисление функции интеграла вероятности и ей обратной |
| title_sort | вычисление функции интеграла вероятности и ей обратной |
| topic | Обчислювальні системи |
| topic_facet | Обчислювальні системи |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/83909 |
| work_keys_str_mv | AT teslergs vyčisleniefunkciiintegralaveroâtnostiieiobratnoi AT zyhakzung vyčisleniefunkciiintegralaveroâtnostiieiobratnoi AT teslergs občislennâfunkcíííntegralaímovírnosteitaííobernenoí AT zyhakzung občislennâfunkcíííntegralaímovírnosteitaííobernenoí AT teslergs computingprobabilityintegralfunctionandreverse AT zyhakzung computingprobabilityintegralfunctionandreverse |