Метод послідовних наближень для розв’язання інтегральних рівнянь актуарної математики

Специфіка системи інтегральних рівнянь для ймовірностей нерозорення на нескінченному інтервалі часу для процесу ризику у випадковому марковському середовищі. Характеристика та особливості класичного актуарного інтегрального рівняння типу Вольтерра.

Рубрика Математика
Вид автореферат
Язык украинский
Дата добавления 28.12.2015
Размер файла 79,6 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

Національна академія наук України

Інститут кібернетики імені В.М. Глушкова

Автореферат

дисертації на здобуття наукового ступеня

кандидата фізико-математичних наук

Метод послідовних наближень для розв'язання інтегральних рівнянь актуарної математики

01.05.01 - теоретичні основи інформатики та кібернетики

Норкін Богдан Володимирович

Київ - 2006

Дисертацією є рукопис.

Робота виконана в Інституті кібернетики ім. В.М. Глушкова НАН України.

Науковий керівник: доктор фізико-математичних наук, професор,

академік НАН України

Єрмольєв Юрій Михайлович,

Інститут кібернетики ім. В.М. Глушкова НАН України,

головний науковий співробітник.

Офіційні опоненти: доктор фізико-математичних наук, професор

Єлейко Ярослав Іванович,

Львівський національний університет імені Івана Франка,

механіко-математичний факультет,

завідувач кафедри теоретичної та прикладної статистики,

кандидат фізико-математичних наук

Журбенко Микола Георгійович,

Інститут кібернетики ім. В.М. Глушкова НАН України,

старший науковий співробітник.

Провідна установа: Дніпропетровський національний університет, факультет прикладної математики, кафедра методів оптимізації.

ЗАГАЛЬНА ХАРАКТЕРИСТИКА РОБОТИ

Актуальність теми. Наразі в Україні бурхливо розвивається страховий ринок. Страхування майна, життя, здоров'я, цивільної відповідальності автомо-білістів посіли чільне місце у нашому бутті. Про медичне і пенсійне страхування варто згадати окремо - розвиток цих видів страхування складає зараз чи не найважливішу частину державної соціальної програми. Сфера бізнесу вимагає розвитку страхування інвестицій, кредитів, банківських вкладів, перевезень та інших видів страхування. Тому вдосконалення методів актуарних розрахунків, методик вибору параметрів управління страховою компанією, страхових тарифів, параметрів страхового договору, ціни перестраховки, більш точна оцінка ймовірності банкрутства компанії дозволило б підняти якість і здешевити послуги страхових компаній для бізнес-структур та населення.

Головною рисою страхової діяльності є те, що бізнес тут ґрунтується на випадковості, тому відповідні моделі є стохастичними. Діяльність страхової компанії спрямована на захист клієнта від ризику великих випадкових втрат за рахунок малих, але детермінованих втрат, які називаються страховими преміями або страховими внесками. Здатність страхової компанії вести такий бізнес, по суті, заснована на законі великих чисел: працюючи з великою кількістю незалежних страхових договорів компанія забезпечує майже передбачену динаміку її капіталу.

Математичні моделі страхової діяльності розвивалися в роботах зарубіжних учених Ф. Лундберга, Г. Крамера, О. Лундберга, К. Борха, С. Андерсена, С. Асмусена, Х.У. Гербера, Н. Прабху та ін., російських вчених А.Н. Ширяєва, І.В. Євстігнеєва, В.Е. Бенінга, В.В. Калашнікова, О.В. Мельникова, Г.І. Фаліна та ін., українських вчених Ю.М. Єрмольєва, І.М. Коваленка, В.С. Королюка, Б.В. Бондарєва, М.С. Братійчука, Д.В. Гусака, Я.І. Єлейко, Ю.С. Мішури, О.М. Наконечного, О.І. Ястремського та ін.

Сучасна теорія страхування побудована на теорії випадкових процесів. Найпершою класичною моделлю такого роду була модель Ф. Лундберга (1903 р.) (класичний процес ризику або складний пуассонівський процес ризику), яка з'явилася приблизно в той же час, що і стохастична модель зміни цін акцій Башельє (модель броунівського руху). На відміну від моделі, що трактує динаміку ціноутворення акцій як броунівський рух, модель Лундберга розглядає стрибкоподібний випадковий процес зміни капіталу: за відсутності страхових вимог капітал монотонно зростає, а в момент надходження вимоги він миттєво зменшується на випадкову величину. Інтервали часу між приходом вимог випадкові і розподілені згідно з експоненційним законом. Ці інтервали і величини вимог незалежні й однаково розподілені. Дану модель у подальшому розробляли Г. Крамер, В. Феллер, С. Андерсен, К. Борх, Х.У. Гербер та ін.

Математична теорія страхування включає такі розділи як моделювання розподілів вимог, принципи обчислення страхових премій, теорію перестраховки, оцінки ймовірності розорення, управління портфелем страхових договорів та інші.

Одним з важливих показників стабільності бізнесу є ймовірність неро-зорення страхової компанії, тому важливою компонентою цих досліджень було вивчення ймовірності (не) розорення страхової компанії як функції її початкового капіталу та інших параметрів. Це нетривіальна проблема, оскільки слід знайти ймовірність попадання реалізації випадкового процесу в негативну область на нескінченному інтервалі часу. Зазначену ймовірність можна обчислити загальним методом Монте-Карло шляхом прямого моделювання випадкового процесу. Складність полягає в тому, що за малої ймовірності розорення (а це, безумовно, найважливіший та найцікавіший випадок) загальний метод Монте-Карло стає неефективним, тобто він не дозволяє оцінити малу ймовірність з необхідною точністю за прийнятний час. Тому набули розвитку різні точні та наближені аналітичні методи, числові методи, практичні наближені оцінки ймовірності розорення. Вони побудовані на тому, що для класичного процесу ризику ймовірність розорення задовольняє інтегральному рівнянню типу Вольтерра. Виявляється, що для більш загальних моделей страхової діяльності ймовірність розорення як функція початкового капіталу також задовольняє деяким складнішим інтегральним рівнянням, для яких не відомі ні аналітичні, ні чисельні методи розв'язання. Дисертація продовжує і розвиває цю лінію досліджень: виводяться нові більш загальні інтегральні рівняння і системи рівнянь для ймовірності (не) розорення страхової компанії, що функціонує за різних умов, встановлюються умови існування і єдиності їх розв'язків, обґрунтовується загальний метод послідовних наближень для знаходження розв'язків.

Зв'язок роботи з науковими програмами, планами, темами. Робота виконувалася у відповідності з планами наукових досліджень відділу №130 Інституту кібернетики ім. В.М. Глушкова НАН України в рамках наступних науково-дослідних тем:

ВФК 130.09 "Розробка математичних методів та обчислювальних алгоритмів для аналізу, оптимізації та керування ризиком", виконувалася за постановою Президії НАН України, протокол бюро ВІ № 8 від 12.01.2002;

ВФ 130.07 "Розробити методи оцінки ризику та їх застосувань у економіці, фінансовій та страховій математиці, теорії надійності", виконувалася за постановою Президії НАН України, протокол бюро ВІ № 6-А від 10.02.2000 р.

Мета і задачі дослідження. Основною метою роботи є розробка методів оцінки ризику банкрутства страхових компаній. Ця мета досягається шляхом побудови математичних моделей страхової діяльності у вигляді стохастичних процесів ризику, виведення інтегральних рівнянь (або систем інтегральних рівнянь) для ймовірності банкрутства як функції початкових умов процесу, дослідження властивостей цих рівнянь, розробки методів їх розв'язання.

Наукова новизна отриманих результатів. Отримані результати є новими і узагальнюють деякі результати, що відомі в літературі лише для класичного процесу ризику. А саме отримані наступні нові результати:

розглянуто узагальнення класичного процесу ризику, що описує еволюцію капіталу страхової компанії, зокрема процеси зі змінними детермінованими преміями, з випадковими преміями, з непуассонівськими потоками премій і вимог, процеси ризику у випадковому марковському середовищі;

узагальнено інтегральні рівняння для ймовірності (не) розорення на нескінченному інтервалі часу як функції початкового капіталу компанії для розглянутих узагальнень класичного процесу ризику;

вперше виведено системи інтегральних рівнянь для ймовірностей нерозорення на нескінченному інтервалі часу для процесу ризику у випадковому марковському середовищі;

вперше встановлені необхідні і достатні умови існування та загальні достатні умови існування і єдиності розв'язків інтегральних рівнянь страхової математики для розглянутих процесів ризику;

вперше теоретично і практично обґрунтовано метод послідовних наближень для чисельного вирішення загальних інтегральних рівнянь актуарної математики, який був застосовний раніше лише для класичного актуарного інтегрального рівняння типу Вольтерра;

вперше отримані умови рівномірної збіжності та оцінки швидкості такої збіжності запропонованого методу послідовних наближень при розв'язанні інтегральних рівнянь страхової математики;

вперше запропонована методика оцінки точності наближених розв'язків розглянутих рівнянь шляхом побудови наближень до розв'язку згори та знизу.

Практичне значення отриманих результатів. Метод послідовних наближень для розв'язання інтегральних рівнянь страхової математики, що розроблений в дисертації, дозволяє підвищити точність актуарних розрахунків для складних моделей страхової діяльності, зокрема, у рамках обраної моделі, дозволяє обчислити ймовірність банкрутства страхової компанії із будь-якою точніостю.

Особистий внесок здобувача. Всі публікації за темою дисертації виконані без співавторів.

Апробація результатів дисертації. Результати досліджень у рамках даної дисертації доповідалися на наукових семінарах в Інституті кібернетики імені В.М. Глушкова НАН України, Львівському національному університеті імені Івана Франка та на наступних конференціях:

міжнародні конференції "Problems of decision making under uncertainties (PDMU-2003, 2004, 2005)", мм. Алушта, Тернопіль, Бердянськ;

Десята міжнародна конференція імені академіка М. Кравчука, 13-15 травня 2004 р., м. Київ;

VI Всеукраїнська конференція молодих наукових співробітників "Інформаційні технології в освіті, науці і техніці (ІТОНТ-2004)", 28-30 квітня 2004 р., м. Черкаси;

міжнародна конференція “Functional methods in approximation theory, operator theory, stochastic analysis and statistics II” пам'яті А.Я. Дороговцева (1935-2004 рр.), 1-5 жовтня 2004 р., м. Київ;

міжнародна конференція "Modern problems and new trends in probability theory", 19-26 червня 2005 р., м. Чернівці.

Публікації. Результати роботи опубліковано в 5 статтях у фахових наукових виданнях і в 7 тезах конференцій.

Структура та обсяг роботи. Дисертаційна робота складається із вступу, семи розділів, висновків, списку використаних джерел (123 найменування). Загальний обсяг роботи 137 сторінок, основний текст роботи викладено на 125 сторінках.

ЗМІСТ РОБОТИ

У вступі обґрунтовується актуальність вибраної тематики, формулюються мета і задачі дослідження, наводиться зміст одержаних результатів.

Перший розділ присвячено огляду літератури за темою дисертаційної роботи та вибору напрямків досліджень.

У другому розділі розглядається класична модель страхової діяльності (Ф. Лундберга), яка ґрунтується на моделюванні динаміки капіталу страхової компанії як складного пуассонівського процесу ризику, та різні її узагальнення. В класичній моделі страхові премії постійні, тому капітал страхової компанії з часом лінійно монотонно зростає між моментами надходження вимог, а в момент надходження страхової вимоги він миттєво зменшується на випадкову величину. Інтервали часу між приходом вимог випадково розподілені згідно з експоненційним законом. Ці інтервали та розмір вимог незалежні й однаково розподілені. Розглянуто також моделі із змінюваним темпом надходження премій (залежним від поточного капіталу), з непуассонівським потоком вимог, моделі з випадковими преміями, що надходять у випадкові моменти часу, моделі функціонування страхової компанії у випадковому (марковському) середовищі. Для кожного типу процесу ризику наведено відповідні інтегральні рівняння для ймовірності небанкрутства як монотонної функції початкового капіталу страхової компанії.

Інтегральні рівняння страхової математики для ймовірності небанкротства мають загальний вигляд

, (1)

де - деякий лінійний інтегральний оператор, причому розв'язок інтегрального рівняння має бути монотонним, задовольняти нерівності та граничній умові

. (2)

Ця умова природна, вона означає, що компанія не розорюється маючи нескінченно великий початковий капітал. В дисертації пропонується розв'язувати всі інтегральні рівняння актуарної математики єдиним методом - методом послідовних наближень:

, (3)

де - номер ітерації. При цьому виявляється, що оператор не є стискаючим на усій своїй області визначення, отже обґрунтування методу послідовних наближень є неабиякою проблемою. У подальших розділах досліджуються властивості оператора для процесів ризику конкретних видів, умови існування і єдиності розв'язку задачі (1), (2), характер послідовних наближень та умови збіжності послідовності до розв'язку задачі (1), (2). Основним елементом дослідження задачі (1), (2) і методу (3) є перевірка існування монотонних функцій , , таких, що

(4)

Виявляється, що як завжди можна взяти , а як - узагальнену межу Крамера - Лундберга, залежну від оператора . В дисертації показано, що для розглянутих операторів умова (4) є необхідною і достатньою для існування розв'язку задачі (1), (2), а послідовні наближення, що стартують з , монотонно збігаються згори, а ті, що стартують з , монотонно збігаються знизу до розв'язку задачі (1), (2). Єдиність розв'язку задачі (1), (2) доводиться для кожного рівняння окремо з використанням властивостей відповідного оператора .

У третьому розділі досліджено метод послідовних наближень Пікара для розв'язання інтегрального рівняння відновлення (типу Вольтерра), якому задовольняє ймовірність (не) банкрутства класичного процесу ризику.

Класичний процес ризику, що описує еволюцію капіталу страхової компанії, задається співвідношенням

,(5)

де - час; - початковий капітал страхової компанії; - інтенсивність надходження премій; - агреговані виплати вимог до моменту , ; - незалежні однаково розподілені випадкові величини (вимоги) з функцією розподілу і середнім значенням ; - число виплат до моменту (пуассонівський процес з інтенсивністю ). Як відомо, функція ймовірності не банкрутства за початкового капіталу задо-вольняє інтегральному рівнянню (відновлення)

. (6)

Як відомо, рівняння Вольтерра має єдиний розв'язок, який можна знайти за допомогою методу послідовних наближень Пікара:

(7)

де - деяка початкова функція; - номер ітерації.

Метод послідовних наближень дозволяє також знаходити аналітичні наближені розв'язки для ймовірності як функції та інших параметрів процесу (5). Як відомо, метод Пікара для інтегральних рівнянь Вольтерра збігається на кожному скінченному інтервалі значень згідно з принципом стискуючих відображень. Проте, в загальному випадку, метод не є монотонним і стиснення має місце тільки після деякого числа ітерацій, тому гарантована збіжність методу за великих може бути повільною.

У розділах 4, 5 узагальнено результати розділу 3 для числового розв'язання загальних (вже не вольтерровських) інтегральних рівнянь теорії процесів ризику. Розглянемо процес ризику вигляду

, (8)

де інтенсивність надходження премій залежить від поточного капіталу компанії і є кусково-неперервною функцією; - агреговані виплати вимог до моменту , ; - незалежні однаково розподілені випадкові величини (вимоги) з функцією розподілу і середнім значенням ; - число виплат до моменту . Позначимо функцію розподілу інтервалів часу між приходом послідовних вимог та часу надходження першої вимоги. Випадок пуассонівського потоку вимог з розглянуто у розділі 4, а загальний випадок - у розділі 5. Визначимо функцію як розв'язок задачі Коші: . Для ймовірності небанкрутства страхової компанії з початковим резервом на нескінченному інтервалі часу справедливо наступне інтегральне рівняння (для див.):

,(9)

де інтеграли розуміються у сенсі Лебега - Стільтьєса. У випадку та лінійної функції це рівняння з урахуванням вимоги (2) зводиться до (6), проте у загальному нелінійному або непуассонівському випадку воно не зводиться до рівняння Вольтерра.

У розділі 6 досліджується процес ризику, в якому стохастичними є як вимоги, так і премії. Він більш адекватно, ніж класична модель Крамера -Лундберга, описує стохастичний процес надходження клієнтів у страхову компанію. В роботі А.В. Бойкова для подібного процесу з незалежними пуассонівськими потоками випадкових премій і вимог було виведено інтегральне рівняння для ймовірності нерозорення і знайдені деякі аналітичні розв'язки цього рівняння. В даному розділі виведено більш загальні інтегральні рівняння для ймовірності нерозорення у разі, можливо, непуассонівських потоків випадкових премій і вимог та за наявності детермінованих платежів, а також обґрунтовано метод послідовних наближень для розв'язання цих рівнянь. Встановлено необхідні та достатні умови існування і загальні достатні умови єдиності розв'язку рівнянь, доведено збіжність методу послідовних наближень та оцінено швидкість збіжності.

У сьомому розділі розглянуто загальний процес ризику, що описує стохастичну еволюцію капіталу страхової компанії в марковському випадковому середовищі за нелінійного надходження премій. Показано, що ймовірність розорення процесу (компанії) як функція початкового стану в загальному випадку

задовольняє деякій системі інтегральних рівнянь з граничними умовами на нескінченності. Встановлено достатні умови існування розв'язку цієї системи і обґрунтовано метод послідовних наближень для обчислення ймовірності розорення. При старті з одиничного початкового наближення ітерації збігаються до розв'язку задачі згори, а при старті з узагальненої межі Крамера - Лундберга - знизу. Таким чином, завжди є можливість оцінити точність одержаного наближеного розв'язку. Працездатність методу перевірена на числовому прикладі процесу ризику на марковському ланцюгу з двома станами.

ВИСНОВКИ

В дисертації отримано наступні основні результати.

1. Виведені інтегральні рівняння для ймовірності (не) розорення різних узагальнень класичного процесу ризику, зокрема процесу з детермінованими преміями, залежними від поточного капіталу, з випадковими преміями, з непуассонівськими потоками премій і вимог та для процесу ризику у випадковому марковському середовищі.

2. Встановлені загальні необхідні та достатні, а також конкретні достатні умови існування і єдиності розв'язків розглянутих інтегральних рівнянь і систем рівнянь страхової математики.

3. Теоретично і практично обґрунтовано метод послідовних наближень для числового розв'язання розглянутих інтегральних рівнянь страхової математики. інтегральний вольтерр марковський

4. Розроблена методика оцінки точності наближеного розв'язку інтегральних рівнянь страхової математики шляхом побудови наближень до точного розв'язку знизу і згори.

5. Запропонований метод послідовних наближень для розв'язання інтегральних рівнянь страхової математики апробовано на ряді числових прикладів, проведено його порівняння з методом Монте-Карло та з відомими наближеними оцінками ймовірності банкрутства.

Розроблений метод послідовних наближень для розв'язання інтегральних рівнянь страхової математики дозволяє:

· підвищити точність актуарних розрахунків;

· обчислити (у рамках обраної моделі) ймовірність розорення страхової компанії з будь-якою наперед заданою точністю;

· провести перевірку точності і застосовності різних відомих наближених формул для ймовірності розорення;

· за необхідності поліпшити точність емпіричних апроксимацій ітеративним шляхом;

· оцінити точність та скорегувати параметри загального Монте - Карло при обчисленні ймовірності розорення за допомогою імітаційного моделювання.

ОСНОВНІ РЕЗУЛЬТАТИ ДИСЕРТАЦІЇ ОПУБЛІКОВАНО В ТАКИХ ПРАЦЯХ

1. Норкин Б.В. Система интегро-дифференциальных уравнений для вероятности банкротства процесса риска в Марковской среде // Теорія оптимальних рішень. Вип.1. - К.: Ін-т кібернетики ім. В.М. Глушкова НАН України, 2002. - С. 21-29.

2. Норкин Б.В. О методе последовательных приближений для вычисления вероятности банкротства классического процесса риска // Теорія оптимальних рішень. Вип.2 - К.: Ін-т кібернетики ім. В.М. Глушкова НАН України, 2003. - С. 10-18.

3. Норкин Б.В. Метод последовательных приближений для решения интегральных уравнений теории процессов риска // Кибернетика и системный анализ. - 2004. - № 4. - С.10-18.

4. Норкин Б.В. Метод последовательных приближений для вычисления вероятности банкротства процесса риска в марковской среде // Кибернетика и системний анализ - 2004. - №6. - C. 149-161.

5. Норкин Б.В. О вычислении вероятности банкротства непуассоновского процесса риска методом последовательных приближений // Проблемы управления и информатики. - 2005. - №2. - C. 133-144.

6. Норкін Б.В. Інтегральні рівняння для дивідендів і ймовірності банкрутства страхової компанії // Матеріали десятої міжнар. конф. ім. акад. М. Кравчука, 13-15 травня 2004 р., Київ. - К.: ТОВ “Задруга”, 2004. - С. 193.

7. Норкин Б.В. Математическое моделирование выплат дивидендов и риска разорения страховой компании // Матеріали IV Всеукр. конф. молодих науковців “Інформаційні технології в освіті, науці і техніці (ІТОНТ 2004)”, 28-30 квітня 2004 р. - Черкаси: Черкаський національний ун-т ім. Богдана Хмельницького, 2004. Ч. 1, - С. 194.

8. Norkin B.V. On calculation of ruin probability // Abstracts of intern. conf. “Problems of decision making under uncertainties (PDMU-2003), September 8-12, 2003, Alushta (Crimea). - Kyiv: Taras Shevchenko Kyiv National University, 2003. - P. 44-45.

9. Norkin B.V. Calculation of ruin probabilities for a risk process in Markovian environment // Abstracts of intern. workshop “Prediction and decision making under uncertainties” (PDMU-2004, 25-30 May 2004), Ternopil, Ukraine. - Kyiv: Taras Shevchenko Kyiv National University, 2004. -P. 42-43.

10. Norkin B.V. On solution of the basic actuarial integral equation // Abstracts of сonf. “Functional methods in approximation theory, operator theory, stochastic analysis and statistics II” dedicated to the memory of A.Ya. Dorogovtsev (1935-2004), October 1-5, 2004, Kiev, Ukraine. - Kyiv: Taras Shevchenko Kyiv National University, 2004. - P. 90.

11. Norkin B.V. A successive approximation method for solution of actuarial integral equations // Abstracts II of Intern. conf. "Modern problems and new trends in probability theory", June 19-26, 2005, Chernivtsi, Ukraine. - Kyiv: Institute of mathematics of the National Academy of Sciences of Ukraine, 2005. - P. 51.

12. Norkin B.V. On the theory of actuarial integral equations // Abstracts of intern. workshop “Problems of decision making under uncertainties” (PDMU-2005, 12-17 October 2005), Berdyansk, Ukraine, 2005. - Kyiv: Taras Shevchenko Kyiv National University, 2005. - P. 40-42.

АНОТАЦІЯ

Норкін Б.В. Метод послідовних наближень для розв'язання інтегральних рівнянь актуарної математики. - Рукопис.

Дисертація на здобуття наукового ступеня кандидата фізико-математичних наук за спеціальністю 01.05.01 - теоретичні основи інформатики та кібернетики. - Інститут кібернетики ім. В.М. Глушкова НАН України, Київ, 2006.

В дисертації одержано наступні основні результати. Розглянуто різні узагальнення класичного процесу ризику (модель Крамера - Лундберга), що описує стохастичну еволюцію капіталу страхової компанії. Зокрема, розглянуті процеси ризику з непостійними детермінованими преміями, з випадковими преміями, з непуассонівськими потоками премій і вимог, процеси ризику у випадковому марковському середовищі. Виведено інтегральні рівняння для ймовірності (не) розорення як функції початкового капіталу компанії для різних узагальнень класичного процесу ризику. Для процесу ризику у випадковому марковському середовищі одержані системи інтегральних рівнянь для набору ймовірності нерозорення з різних початкових станів середовища як функцій початкового капіталу компанії. Встановлені загальні необхідні і достатні, а також конкретні достатні умови існування та єдиності рішень розглянутих інтегральних рівнянь і систем рівнянь страхової математики. Теоретично і практично обґрунтований метод послідовних наближень для чисельного або аналітичного розв'язання розглянутих інтегральних рівнянь математики страховки, зокрема, доведена його рівномірна збіжність і одержані оцінки швидкості збіжності. Розроблена методика оцінки точності наближеного розв'язку інтегральних рівнянь математики страховки шляхом побудови наближень до точного розв'язку зверху та знизу. Запропонований метод послідовних наближень апробований на низці числових прикладів, проведено його порівняння з методом Монте-Карло, з відомими наближеними оцінками розв'язку. Розроблений метод послідовних наближень для розв'язання інтегральних рівнянь математики страховки дозволяє підвищити точність актуарних розрахунків, а саме - обчислити в рамках обраної моделі ймовірність розорення страхової компанії з будь-якою наперед заданою точністю, провести перевірку точності і застосовності різних відомих наближених формул для ймовірності розорення. За необхідності поліпшити точність емпіричних апроксимацій ітеративним шляхом, оцінити точність і скоректувати параметри Монте-Карло при обчисленні ймовірності розорення за допомогою імітаційного моделювання.

Ключові слова: актуарна математика, страхова математика, процес ризику, ймовірність банкрутства, ймовірність розорення, інтегральні рівняння, існування та єдиність розв'язків, метод послідовних наближень.

Норкин Б.В. Метод последовательных приближений для решения интегральных уравнений актуарной математики. - Рукопись.

Диссертация на соискание ученой степени кандидата физико-математических наук по специальности 01.05.01 - теоретические основы информатики и кибернетики. - Институт кибернетики им. В.М. Глушкова НАН Украины, Киев, 2006.

В диссертации рассмотрены различные обобщения классического процесса риска (модель Крамера - Лундберга), описывающего стохастическую эволюцию капитала страховой компании. В частности, рассмотрены процессы риска с непостоянными детерминированными премиями (которые зависят от текущего капитала компании), процессы со случайными премиями (приходящими в случайные или детерминированные моменты времени), процессы с непуассоновскими моментами прихода премий и требований, процессы риска в случайной марковской среде.

Выведены интегральные уравнения для вероятности (не) разорения как функции начального капитала компании для различных обобщений классического процесса риска. В случае классического процесса риска соответствующее уравнение сводится к уравнению Вольтерра и известно много подходов к его решению. Этому случаю посвящено огромное число работ в страховой математике. Однако для рассмотренных обобщений классического процесса риска оказывается, что соответствующие актуарные интегральные уравнения не сводятся к уравнениям Вольтерра. Все эти уравнения имеют следующий общий вид, искомая функция равняется значению некоторого оператора от этой же функции. Это линейные однородные интегральные уравнения, не принадлежащие к известным типам интегральных уравнений. Они содержат интегралы по неограниченным областям. Операторы, входящие в уравнения, не являются сжимающими на своих областях определения. Данные уравнения имеют много решений, но ищется решение, удовлетворяющее граничному условию на бесконечности, а именно, решение, стремящееся к единице при росте капитала к бесконечности. Такое условие означает, что компания не разоряется из бесконечного начального капитала.

Для процесса риска в случайной Марковской среде получена система линейных интегральных уравнений для набора вероятностей неразорения как функций начального капитала компании при различных начальных состояниях среды. Математически задача состоит в том, чтобы найти решение уравнения или системы, удовлетворяющее указанному граничному условию на бесконечности. В диссертации установлены общие необходимые и достаточные, а также конкретные достаточные условия существования и единственности решений рассмотренных интегральных уравнений и систем интегральных уравнений страховой математики. Эти условия сводятся к существованию некоторой монотонной функции (обобщенной нижней границы Крамера - Лундберга) со значениями между нулем и единицей такой, что интегральные операторы уравнений переводят эту функцию в некоторую большую функцию.

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

В силу общности и сложности рассматриваемых интегральных уравнений страховой математики получить их точное аналитическое решение невозможно, за исключением очень редких специальных случаев. Поэтому в диссертации предложен, теоретически и практически обоснован метод последовательных приближений для численного или приближенного аналитического решения рассмотренных интегральных уравнений страховой математики, в частности, доказана его равномерная сходимость к решению и получены оценки скорости сходимости. Показано, что метод равномерно сходится в любую окрестность решения со скоростью геометрической прогрессии (знаменатель которой зависит от размера окрестности). Предложена методика оценки точности приближенных решений интегральных уравнений страховой математики путем организации приближений к точному решению сверху и снизу. Показано, что запущенный с единичного начального приближения метод порождает монотонно убывающую последовательность приближений, а запущенный с обобщенной нижней границы Крамера - Лундберга - монотонно возрастающую последовательность приближений, сходящуюся к решению. Все эти результаты справедливы и для решения систем интегральных уравнений, описывающих вероятности неразорения процесса риска в марковской случайной среде.

Предложенный метод последовательных приближений апробирован на ряде численных примеров, проведено его сравнение с методом Монте - Карло, с известными приближенными оценками вероятности разорения.

Разработанный метод последовательных приближений для решения интегральных уравнений страховой математики позволяет повысить точность актуарных расчетов, а именно, вычислить в рамках выбранной модели вероятность разорения страховой компании с любой наперед заданной точностью; провести проверку точности и применимости различных известных приближенных формул для вероятности разорения; при необходимости улучшить точность эмпирических аппроксимаций итеративным путем; оценить точность и скорректировать параметры общего метода Монте - Карло при вычислении вероятности разорения с помощью имитационного моделирования.

Ключевые слова: актуарная математика, страховая математика, процесс риска, вероятность банкротства, вероятность разорения, интегральные уравнения, существование и единственность решения, метод последовательных прибли-жений.

Norkin B.V. A successive approximation method for solution of actuarial integral equations. - Manuscript.

Dissertation for a scientific candidate degree of physical and mathematical sciences by speciality 01.05.01 - theoretical bases of computer science and cybernetics. - V.M. Glushkova Institute of Cybernetics of the National Academy of Sciences of Ukraine, Kiev, 2006.

In the dissertation various generalizations of the classical risk process describing a stochastic evolution of the capital of an insurance company (Kramer - Lundberg model) are considered. In particular risk processes with variable deterministic premiums, with random premiums, with non Poisson flows of premiums and claims, risk processes in a stochastic markovian environment are studied. Integral equations for the probability of nonruin as a function of the initial capital of the company for various generalizations of the classical risk process are deduced. For a risk process in a stochastic markovian environment system of integral equations for a set of nonruin probabilities from various initial states of the process are obtained. A general necessary and sufficient, and also concrete sufficient conditions of the existence and uniqueness of solutions of the considered integral equations and systems of equations are established. A successive approximation method for a numerical or analytical solution of the considered integral equations of insurance mathematics is theoretically and practically validated, in particular its uniform convergence and rate of such convergence are established. A technique of estimation of the accuracy of approximate solutions of integral equations of insurance mathematics is developed by construction of approximations to the exact solution from above and from below. The suggested method of consecutive approximations was tested on a number of numerical examples, its comparison with Monte Carlo method, with known approximations of solutions was carried out.

The developed method of successive approximations for the solution of integral equations of insurance mathematics allows to increase accuracy of actuarial calculations, namely, to calculate within a chosen framework a probability of ruin of an insurance company with any prescribed accuracy, to check applicability and accuracy of various empirical approximate formulas for the probability of ruin, if necessary to improve accuracy of empirical approximations in an iterative way, to estimate accuracy and to correct parameters of Monte Carlo simulation method for calculation of the probability of ruin.

Keywords: actuarial mathematics, insurance mathematics, risk process, ruin probability, probability of bankruptcy, integral equations, existence and uniqueness of solutions, successive approximation method.

Размещено на Allbest.ru


Подобные документы

  • Чисельні методи розв’язання систем нелінійних рівнянь: лінійні і нелінійні рівняння, метод простих ітерацій, метод Ньютона. Практичне використання методів та особливості розв’язання систем нелінійних рівнянь у пакеті Mathcad, Excel та на мові С++.

    курсовая работа [2,0 M], добавлен 30.11.2010

  • Розв’язання систем лінійних рівнянь методом Жордана-Гауса. Еквівалентні перетворення системи, їх виконання як елемент методів розв’язування системи рівнянь. Базисні та вільні змінні. Лінійна та фундаментальна комбінації розв’язків, таблиці коефіцієнтів.

    контрольная работа [170,2 K], добавлен 16.05.2010

  • Лінійні діофантові рівняння. Невизначені рівняння вищих порядків. Невизначене рівняння Ферма. Приклади розв’язання лінійних діофантових рівнянь та системи лінійних діофантових рівнянь. Алгоритми знаходження всіх цілочисельних розв’язків рівнянь.

    курсовая работа [1,7 M], добавлен 29.12.2010

  • Історія створення теорії алгебраїчних рівнянь. Сутність системи лінійних алгебраїчних рівнянь в лінійній алгебрі. Повна характеристика методів розв'язання рівнянь: точні, ітераційні та ймовірнісні. Особливості теорем Гауса-Жордана та Габріеля Крамера.

    реферат [543,7 K], добавлен 23.04.2015

  • Прийоми розв’язання задач в першому і другому степені на Далекому Сході та Греції. Досягнення арабських математиків в області алгебраїчних рівнянь. Розв'язання похідного кубічного рівняння. Найвидатніші теореми про радикали вищих степенів, їх розв’язання.

    курсовая работа [536,1 K], добавлен 23.02.2014

  • Ознайлення з базовими поняттями, фактами, методами та найпростішими застосуваннями рівняння Пфаффа. Виконання завдань щодо розв’язання рівнянь Пфаффа. Аналітичний запис задачі про відшукання інтегральних поверхонь максимально можливої вимірності.

    курсовая работа [489,2 K], добавлен 30.12.2013

  • Системи лінійних алгебраїчних рівнянь, головні означення. Коротка характеристика головних особливостей матричного способу, методу Жордано-Гаусса. Формули Крамера, теорема Кронекера-Капеллі. Практичний приклад розв’язання однорідної системи рівнянь.

    курсовая работа [690,9 K], добавлен 25.04.2013

  • Графічний спосіб розв'язку рівнянь. Комбінований метод пошуку та відокремлення коренів. Метод Ньютона (метод дотичних або лінеаризації). Процедура Ейткена прискорення збіжності. Метод половинного поділу та простих ітерацій уточнення коренів рівняння.

    лекция [1,9 M], добавлен 27.07.2013

  • Етапи розв'язування задачі дослідження певного фізичного явища чи процесу, зведення її до диференціального рівняння. Методика та схема складання диференціальних рівнянь. Приклади розв'язування прикладних задач за допомогою диференціального рівняння.

    контрольная работа [723,3 K], добавлен 07.01.2016

  • Визначення поняття "рівняння з параметрами", розгляд принципів рішення даних рівнянь на загальних випадках. Особливості методів розв'язання рівнянь із параметрами, зв'язаних із властивостями показовою, логарифмічною й тригонометричною функціями.

    реферат [68,3 K], добавлен 15.02.2011

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.