Разлика между версии на „Теория на вероятностите“

От Администрация и управление
Направо към навигацията Направо към търсенето
Ред 1: Ред 1:
'''Методът Монте Карло''' е метод за '''разглеждане поведението на система''', на която всеки етап от [[промени]]те се моделира по известен [[закон]] с помощта на някакъв източник на [[случайност|случайни]] числа.
+
'''Теорията на вероятностите се занимава''' с [[изследване]] [[закономерност]]ите при случайните събития.  
  
Според тази предпоставка се определят съществените [[белег|белези]] на метода, неговото приложение и [[прогнозиране]] на [[риск]]овете, както и взаимодействието му с останалите техниики и закони в процеса на работа.
+
Всяко [[явление]] протича при известни условия. Обратно, при определени условия протича някакво явление, но ние не сме в състояние да опишем всички условия, при които протича дадено явление.  
  
==История==
+
Ние определяме само някои условия за протичането на дадено явление, което ще наричаме събитие. При реализирането на тези условия даденото [[събитие]] може да се сбъдне, а може и да не се сбъдне. Такова събитие, което при даден [[комплекс]] от условия може да се сбъдне, а може и да не се сбъдне, наричаме случайно.
  
  
  
Смята се че методите за [[статистически изследван]]ия - метода на изчислителни и [[приложна математика]], въз основа на [[симулация]]та на случайни величини и изграждане на статистически оценки за неизвестните величини, са същите като медота Монте Карло. През 1930 година Енрико Ферми използва методът „Монте Карло” за изследване на [[неутрон|неутронния]] поток. По-късно той разработва механично устройство, което се използва при изчисленията по ядрена физика. В момента връзката на подобни методи и началото на ерата на компютъните технологии дават възможност за извършване на редица компютърни [[експеримент]]и, включително и тези които водят до получаването на случайни числа.
+
==ВЪВЕДЕНИЕ В ТЕОРИЯ НА ВЕРОЯТНОСТИТЕ==
  
За основатели на метода „Монте Карло” се считат амрекинаските математици [[Стенли Уламовите]], [[Джон фон Нойман]] и Никола „Метрополис”. През четиридесетте години на ХХ век, Джон фон Нойман полага основите на метода - създаването на математическата основа на [[функция|функциите]], [[вероятността плътност]], интегрирани функции, обратност на разпределение и [[генератор]]и на случайни числа. Проучванията показват, че в тясно сътрудничество със Стенли Уламовите, той е първия който реализира идеята за необходимостта от [[компютър]], за да се изпълняват изчисленията на метод Монте Карло много по-лесно.
 
  
[[Произход]]ът на името на метода е свързан с град в Княжество Монако, където се намира и едно от най-известните казина в света – факт е, че случайните числа са в основата на „Монте Карло”, а рулетката е едно от най-лесните [[устройство|устройства]] за генериране на случайни числа.
 
  
Друг интересен [[факт]] е свързан с това, че някои случайни методи за изчисление и експерименти, са проектирани и изпълнени и в [[праисторическата ера]] на [[компютърни технологии|компютърните технологии]]. " Основната разлика на „Монте Карло” с по-ранни изследвания в [[област]]та на [[статистическо моделиране]] е:Симулацията „Монте Карло” се е превърнала в стандарт за това това как да решим проблема с помощта на инструментите от [[теорията на вероятностите]] и математическата статистика. Което значи че това е предпоставка за необходимостта да се за занимаваме с [[детерминация]]та на проблема, след което да използваме [[симулация]]та, за да проверим [[изчисления]]та, направени по-рано.
+
===Основни понятия на теория на вероятностите===
  
==„Монте Карло” – анализ, предимства и недостатъци==
 
  
  
 +
Най-често се използва [[понятието]] опит ([[експеримент]]). [[Опитът]] се състои в осъществяване на определен комплекс от условия, в резултат на което се случва нещо, което не може да бъде предсказано предварително със сигурност. При многократното повторение на опита могат да бъдат получени различни [[резултат]]и (изходи), които се представят като поява на определено събитие. Събитията се разграничават по видове, в съответствие с особеностите на тяхното проявление.
  
За провеждане на симулация със [[стохастични модели]] най-често се използва методът „Монте Карло”. Това е [[универсален]] метод за симулация, който намира приложение в различни области на научните изследвания и в практиката
 
  
Обобщената дефиниция на симулацията включва:
 
  
 +
'''Видове събития'''
  
  
* Разработване ([[конструиране]]) на [[модел]] на една [[система]] – обикновено математически и логически по природа; системата може да бъде действителна или теоретична (виртуална);
 
  
 +
'''Сигурни(достоверни) събития'''
  
  
* Описване на реалната система с [[термин]]и, приемливи за компютърните системи;
 
  
 +
Ако при многократно повтаряне на даден [[комплекс]] от условия се реализира едно и също събитие, то се нарича сигурно или достоверно събитие.
  
 +
Ето няколко примера. Ако подхвърлим с ръце еди тежък предмет нагоре, този предмет ще падне. Падането на предмета е сигурно събитие.
  
*  Използване (обикновено) компютър за да се изпълни симулация;
+
Пример 2: при падането на монета, върху която различаваме лицева и гербова страна, върху хоризонтална равнина тя не застава върху равнината отвестно, а пада върху една от страните си. Падането на монетата върху една от страните й е сигурно събитие.
  
 +
Пример 3: при подхвърляне върху хоризонтална равнина едно кубче застава върху една от стените си. Следователно сигурно събитие е падането на куба върху една стена. За конкретност вместо куб ще говорим за [[зар]](кубче което се използува при [[хазарт]]ни игри), стените на които са номерирани от 1 до 6 посредством съответен брой точки. Събитието падане на зара върху стена с брой на точките, по-малък от 7 е [[сигурно събитие]]. Сигурните събитя ще означаваме с буквата U.
  
  
* [[имитация|Имитира]] се действието на реалната система/процес;
 
  
 +
'''[[Невъзможни събития]]'''
  
  
* Симулацията е експериментиране;
 
  
 +
Ако при реализиране на определен комплекс от условия дадено събитие не се сбъдва, то се нарича невъзможно събитие.
  
 +
Пример за невъзможно събитие е при падане върху хоризонтална равнина една монета да застане отвесно. Друг такъв пример е заставането на зар върху ръб или даже върху връх. Невъзможно събитие е падането на зара върху стена със 7 точки. Невъзможните събития ще означаваме с буквата V.
  
* [[Цел]]та на системата е да намерим нещо, отнасящо се до реалната система.
 
  
  
 +
'''[[Противоположни събития]]'''
  
За да се реши проблема чрез метода, първо се изгражда възможен модел, който представлява желана [[стойност]] като [[многомерно]] неразделно под формата на очакването ( с други думи това е случаен [[процес]], който след това се симулира на компютър ). Естествено съществуват и трудности по пътя към прилагането на метода. Но не е необходима тяхната [[дефиниция]], а по скоро надеждната оценка на неизвестните величини. Постигането на това не е толкова лесно колкото изглежда и голяма роля играе изграждането на [[вероятностен модел]].
 
  
Друг важен [[компонент]] това е симулацията на случайни величини с дадена [[дистрибуция]]. Обикновенно такива модели се осъществяват чрез [[конвертиране]] на една или повече независими стойности на случйни числа и  равномерното им разпределние в [[интервал]]а. Последователността на случйните стойности обикновено се извършва от компютър, като се използват теоритичните [[алгоритъм|алгоритми]], сред които най-често се използва „[[Методът на остатъците]]”. Тези числа се наричат „псевдо” и се проверяват от статистически тестове. От тук идва и основната роля на качественото генериране на случайни числа.
 
  
Основната причина, поради която използваме метода „Монте Карло” е защото чрез него много по-бързо можем да получим отговор на даден въпрос, за който по принцип е необходимо доста време, за да бъде разрешен. Например метода „Монте Карло” много добре пасва за намирането на някои много сложни [[интеграл]]и, което е задача с която други методи не биха били толкова полезни. Друго предимиство на метода е лесното симулиране на сложни физически системи в областта на [[физика]]та, [[инженерство]]то и [[химия]]та.
+
Две събития ще наричаме противоположни, ако се реализират при еди и същи комплекс от условия и ако винаги, когато едното събитие се реализира, другото не се реализира.  
  
 +
Пример за противоположни събития са разгледаните по-горе сигурни и невъзможни събития. Падането на една монета върху лицевата й страна и падането й върху гербовата й срана, разгледани като две събития, са също противоположни, защото монетата не може да е паднала върху лицевата си страна, ако е паднала върху гербовата си страна и обратно.
  
 +
Падането на зара върху стена с [[четен брой]] точки и върху стена с [[нечетен брой]] е също невъзможно събитие. Не е трудно да се разбере, че всяко събитие има свое противоположно. [[противоположност|Противоположното]] на дадено събитие можем да дефинираме като събитие, което се реализира, ако даденото не се реализира, или обратно, което остава нереализирано в резултат на реализирането на даденото събитие.
  
Естествено методът си има и своите недостатъци. Все още не е общоприето дали „Монте Карло” може да се използва за симулиране на системи, които не са в [[равновесие]] (т.е. преходно състояние). Друг недостатък е, че за да се използва „Монте Карло” симулация трябва да се генерират голям брои [[проба|проби]], а това може да отнеме доста време до постигане на желаните резултати, понеже не може да се генерира само един [[образец]] и веднага да се използва в симулация, а е необходимо да се направят много проби и да се получи средната им стойност.
+
Ако с А означим даденото събитие, то с А ще означаваме неговото противоположно. [[логика|Логически]] казано, А е не А и А е не А
  
Резултатите от „Монте Карло” симулацията са само приблизителни на действителната стойност и съвсем не са толкова точни.
+
Съгласно въведените по-горе означения за сигурно и невъзможно събитие имаме равенствата
  
==Приложение==
+
U = V, V = U
  
  
 +
===Несъвместими и съвместими събития===
  
Методът Монте Карло не се използва за определяне на [[двумерни площи]]. За това се прилагат по-точни методи. Но при определяне на площи и [[обем]]и на тела в многомерното пространство Монте Карло методът се оказва единствената възможност да се реши задачата.
 
  
  
 +
Събития, които не могат да се реализират (едновременно) при един и същ комплекс от условия, се наричат несъвместими. Несъвместими събития са например противоположните. Но докато при противоположните се реализира едно от тях, щом другото не се е реализирало, при несъвместимите събития това не е задължително. Две несъвместими събития могат и едновременно да не се реализират. Тяхното едновременно реализиране обаче е отречено по определение. При хвърлянето на един зар несъвместими са например следните две събития: падането на зара върху стена с една точка и падането му върху стена с четен брой точки. Съвместимите събития ще определим като събития, които не са несъвместими.Пример на съвместими събития е падането на зар върху стена с една точка и падането на зар върху стена с нечетен брой точки
  
[[Image:snimka1.jpg|thumb|left|alt=A scheme.|Модел на Метода Монте Карло.]]
 
  
Монте Карло Методът има две важни особености:
 
  
1. Проста [[структура]] на изчислителния алгоритъм;
+
===Елементарни и сложни събития===
  
2. Грешката при изчисленията се отбелязва с формулата, където D=cost, a N – броят на изпитанията.
 
  
[[File:dn.png]]
 
  
 +
При даден [[комплекс]] от условия се реализират различни събития. Някои от тях могат да бъдат определени като основни, неизразени посредством други събития, но посредством които се изграждат останалите събития. Такива основни събития ще наричаме елементарни а останалите събития – сложни. Ще отбележим, че както елементарните, така и сложните събития се обуславят единствено от комплекса от условия. Не е изключено едно и също събитие да може да реализира при два различни комплекса от условия, но при единия комплекс то да е [[елементарно]], а при другия [[сложно]]. По принцип смятаме, че при различни комплекси от условия се реализират различни събития.
  
 +
Примери за елементарни и сложни събития:
  
От тук се вижда, че за да се намали грешката с един порядък (10 пъти) трябва да се увеличи броят на опитите N – 100 пъти. Методът Монте Карло позволява да се моделира произволен процес на чието протичане влияят случайни фактори. За много задачи несвързани със случайни фактори може изкуствено да се въведат вероятностни модели, както беше направено в примера с площите. Поради това за методът Монте Карло може да се говори като за универсален метод за решаване на задачи. Методът е много подходящ при оценката и [[прогнозиране на рискове]] за [[безопасност]] и при оценка на [[надеждност]]та на сложни системи. Решаването на задачи от този тип е свързано с отчитане на множество случайни фактори. Ако анализираме основната формула за изчисляване на [[риск]]а R.
+
При хвърлянето на една монета елементарните събития са две: падане върху лизевата страна и падане върху гербовата страна. Също така при хвърляне на зар имаме 6 [[елементарни събития]], за 6-те различни стени. Падането на зара върху стена с четен брой точки е сложно събитие, което се състои в реализирането на едно от елементарните събития за падане върху стена с 2, 4 или 6 точки. Ако хвърлим две монети едновременно, то елементарните събития са двойка елементарни събития за една монета. Ще означаваме с Л и Г съответно събитията падане на монета върху лицевата и падане върху гербовата й страна. За падане на 2 монети имаме следните четири елементарни събития: ЛЛ, ЛГ, ГЛ, ГГ.  
  
''R ''=''P''.''T''
+
Сложно събитие е падане на двете монети върху различни стени, тогава, когато едната монета е паднала върху Л страна, а другата върху Г или обратно(ЛГ или ГЛ). Така това сложно събитие се реализира само когато се реализира едно от двете елементарни събития(ЛГ или ГЛ), но то не е [[тъждествено]] с тях, а само е изградено посредством тези елементарни събития. И така при даден [[комплекс от условия]] всяко събитие може да се разглежда като изградено посредством елементарни събития. Ние ще разгледаме всяко събитие като [[множество]], чиито елементи са елементарни събития.  
  
където:
 
  
- P е вероятността да се случи [[инцидент]];
 
  
- Т е тежестта на последиците от инцидента.
+
===Действия със събития. Поле от събития===
  
Следва да отбележим, че и последиците от инцидентите, както и тяхната тежест, също са с вероятностен характер. Вероятността P да се случи инцидент на практика винаги е сложна вероятност, която описва възможността да се реализира [[каскада]] от събития.
 
  
Метода Монте Карло е изключително полезен при проектиране на [[предприятия]], [[технологии]], [[планиране на промени]] организация. Вместо скъпо струващи експерименти в [[натура]] може да се експериментира на компютър различни варианти на организационната работа и използване на [[оборудване]]то. В редица случаи компютърното моделиране е единствения начин да се получат сведения за поведението на проектираната система. Разбира се, разчетите не идват даром. Няма универсално решение, което да се избере от менюто на “магически” [[софтуер]], което да върне резултата веднага след като попълниш своето задание. Трябва да се създаде добър модел на изследваната система, за да се получат резултати с практическа [[ценност]]. Необходимо е да се изследват щателно действителните потоци от заявки, да се проведат [[хронометражни измервания]] за работата на отделните възли и т.н. Когато системата е в проектна фаза и все още не съществува в действителност се използват проектните параметри и технологичните граници. Въобще, следва да се знаят вероятностните закони на функциониране на отделните части на системата, с което се моделира начина, по който работата им се колебае около проектните стойности. Тогава методът Монте Карло позволява да се изчислят вероятностните закони на работата на цялата система, отчитайки [[взаимосвързаност]]та между отделните и части, колкото и сложна да е тя. По такъв начин се решават задачи за [[разчет]] на надеждност на изделие, [[агрегат]], система, състояща се от множество елементи. Задачи за ефективната поддръжка и др.
 
  
Ето още един пример за полезността на метода Монте Карло при една на пръв поглед проста задача за определяне на [[качество]]то и [[надеждност]]та на изделие състоящо се от даден брой елементи.
+
Известни са действията с множества. Тъй като ние разглеждаме събитията като множества, то имаме възможност лесно да извършваме [[действие|действия]] и със [[събитие|събития]].  
  
 +
Преди всичко, ако А е дадено събитие, изградено посредством елементарните събития а1, а2, а3,... аn,то смятаме, че събитието А се сбъднало, ако се е сбъднало поне едно от [[изграждащ]]ите го елементарни събития а1, а2, а3,... аn.
  
 +
Еквивалентни събития
  
==Пример за определяне на качество и надеждност:==
+
Казваме, че събитието А следва от събитието В или че А се съдържа в В, или В съдържа А, което записваме така; А ?(само че знака е легнал) В или В ?(само че знака е легналА). Занапред, когато говорим за различни събития , ще смятаме, че имаме един и същи комплекс от условия, ако не е казано противното. Ако А и В са две събития, такива че А следва от В и В следва от А, т.е. А?В и В?А Ще казваме че А и B са [[еквивалент]]ни(равносилни) или [[прост]]о равни и това ще записваме така А = B
  
Дадено е изделие състоящо се от множество елементи. Да допуснем, че качеството на изделието се определя от стойността на един изходен параметър - U , който може да се изчисли, знаейки стойностите на входните параметри на изделието R,C и др. Т.е. известна е функционалната зависимост:
 
  
''U = f(R1, R2...; C1, C2...;) където: U е [[напрежение]], R - [[съпротивление]], C – [[капацитет]]''
 
  
Съществен момент при моделирането е, че параметрите не се подбират с фиксирани стойности а се третират като случайни променливи разпределени в някакъв интервал. С метода Монте Карло може да се получи разпределението на изходния параметър U при различни комбинации от стойности на входните параметри. Изходният параметър U , като функция на случайни величини също се явява случайна променлива. Търси именно вида и интервала на разпределението му. Не винаги е ясно при кои набор от параметри ще се получи най-лош резултат. Освен това е малко вероятно всички елементи да се окажат с лоши параметри. За това е [[целесъобразно]] да се отчитат всички параметри, като случайни величини. А оценката на качеството и надеждността да се оцени статистически. Така в резултат Монте Карло симулация ще се получи средната стойност на изходния параметър и отклонението от него, което при достатъчно голям брой повторения на компютърния експеримент ще представляват на практика математическото очакване MU и дисперсията DU на изходният параметър при всички възможни комбинации от стойности на входните параметри.
+
==Вероятност. Методи за определяне на вероятности==
  
Отбелязваме, че ''MU ''¹ ''f (MR1,MR2,…;MC1, MC2,…;)  ''където MR и MC са математическите очаквания съответно на R и C. Аналитичното определяне на разпределението на изходният параметър U при сложна функция е на трудно и в много практически случаи – невъзможно. Използвайки метода Монте Карло можем лесно да получим търсеното разпределение чрез компютърно моделиране.
 
  
Ролята на случайните разпределения в метода Монте Карло:
 
  
Ето някой от най-използваните теоретични разпределения на случайни величини.
+
Когато [[експериментът]] е проведен n пъти и събитие А се е появило m пъти, отношението
  
[[Image:snimka2.jpg|thumb|left|alt=A scheme of some methods.|Методи на разпределение.]]
+
[[File:teoriqnaveroqtnostite1.jpg]]
  
Нормалното разпределение описва много природни феномени като грешки при измервания, ръст на хората, [[коефициент]]и на [[интелигентност]], ръст на [[инфлация]] и много други. То е непрекъснато разпределение с два параметъра – средно и стандартно отклонение. Има три условия на които отговаря нормалното разпределение:
 
  
* От стойностите на случайната променлива най-често срещани са средните стойности;
 
* Случайната променлива дава стойности еднакво разположени около средната – разпределението е [[симетрия|симетрично]];
 
* Случайната променлива е по-вероятно да даде стойност близо до средната отколкото далече от нея. От всички стойностите на нормалното разпределение приблизително 68% са разположени на едно стандартно отклонение от средната.
 
  
Триъгълното разпределение показва броят на успехите на даден опит, когато са известни минимума, максимума и най-вероятната стойност. На пример може да опише броят коли продадени на седмица, когато се знаят трите параметъра от минали продажби. То е непрекъснато вероятностно разпределение с три параметъра - [[минимум]], [[максимум]]и най-вероятната стойност. Има три условия, на които отговаря триъгълното разпределение:
+
показва какъв дял от общия брой случаи се пада на събитие А. Тази [[величина]] се нарича [[относителна честота]] или честост, а числото m - [[честота]]. Най-често относителната честота се записва с
  
* Минималната стойност е фиксирана;
+
[[File:teoriqnaveroqtnostite2.jpg]]
* Максималната стойност е фиксирана;
 
* Най вероятната стойност лежи между минимума и максимума, показвайки че стойностите около долните ъгли на триъгълника са по-малко вероятни от тези около върха – най-вероятната стойност
 
  
[[разпределение на Пуасон|Разпределението на Пуасон]] описва броят пъти за които едно събитие се случва в даден [[интервал]], като брой телефонни разговори за час или брой на грешки на страница от документ или брой на пострадалите от инциденти по време на работа за месец. То е дискретно вероятностно разпределение с един параметър - честота. Има три условия на които отговаря пуасоновото разпределение:
 
  
* Броят на възможните случаи във всяка единица не е [[лимит|лимитиран]];
 
* Случайте са независими. Броят на случаите в една единица от проведеното измерване не влияе върху броят на случаите в другите единици;
 
* Средният брой случай трябва да остава същият от единица към единица
 
  
[[биномно разпределение|Биномното разпределение]] описва броят на успехите при фиксиран брой опити, като броят на падания на ези при 10 хвърляния на монета или броят на дефектните изделия от 50 тествани от склада с готова продукция. То е дискретно вероятностно разпределение с два параметъра – вероятност и брой опити. Има три условия на които отговаря Биномното разпределение:
+
Относителната честота има следните свойства:
  
* При всеки опит има два възможни изхода – успех и неуспех;
+
1) 0W(A)1, т.е. тя е [[неотрицателна величина]].
* Опитите са независими. Това което се получава в резултат от първия опит не засяга вторият и т.н.;
 
* Вероятността за успех остава същата от опит на опит.
 
  
[[логнормално разпределение|Логнормалното разпределение]] широко се използва в ситуации, при които стойностите са положително асиметрични (когато повечето от стойностите се получават близо до минималната стойност). Ситуация забелязана при финансовите анализи при оценката за сигурност или при недвижимите имоти при оценката на собствеността. Финансовите анализатори са забелязали, че цените на фондовата борса обикновено са положително свити а не симетрично разпределени, както при нормалното разпределение. [[ценни книжа|Ценните книжа]] проявяват тази тенденция, защото цените им не могат да паднат под нулевата долна граница, но могат да нараснат произволно. Подобно на тях цените на недвижимите имоти проявяват същата тенденция, защото не могат да станат отрицателни. Логнормалното разпределение е непрекъснато разпределение с два параметъра – средно и стандартно отклонение. Има три условия, на които отговаря Логнормалното разпределение:
+
2) Ако събитие К е достоверно (сигурно) събитие, W(K) = 1, т.е относителната честота на достоверно събитие, в случая събитие К, е равна на единица. Важи и обратното, че ако дадено събитие има [[относителна]] честота равна на единица, то е достоверно събитие.
  
* Случайната променлива може да нараства безгранично;
+
3) За [[редица]]та от събития A1, А2, ... , Ak
* Стойностите на случайната променлива имат положително свито разпределение;
 
* [[натурален логаритъм|Натурален логаритъм]] от случайната променлива води до нормалната крива.
 
  
[[еднородно разпределение|Еднородно (униформено) разпределение]], всички стойности между минимума и максимума могат да се получат с еднаква вероятност. Това е непрекъснато разпределение с два параметъра – минимум и максимум. Има три условия, на които отговаря Еднородното разпределение:
+
[[File:teoriqnaveroqtnostite3.jpg]]
  
* Минималната стойност е фиксирана;
 
* Максималната стойност е фиксирана;
 
* Всички стойности между минимума и максимума са еднакво вероятни.
 
  
[[експоненциално разпределение|Експоненциалното разпределение]] описва периода от време, което минава между събития, като времето между два отказа в електронно оборудване, или времето между посещенията на обслужващо гише. То е непрекъснато разпределение с един параметър – честота. Има две условия, на които отговаря Експоненциалното разпределение:
 
  
* Експоненциалното разпределение допълва Пуасоновото разпределение.
+
За събитията А и В с честоти m(А) и m(В) и честости W(A) и W(В) важи 
* Пуасоновото разпределение описва броят на събитията за единица време, Експоненциалното – периода от време между събитията. Експоненциалното разпределение не се влияе от предишни събития
 
  
[[геометрично разпределение|Геометричното разпределение]] описва броят на опитите до първият успех, такива като броят [[сондаж]]и, които трябва да се пробият докато избие нефтения фонтан. Това е дискретно разпределение с един параметър - [[вероятност]]. Има три условия, на които отговаря Геометричното разпределение:
+
[[File:teoriqnaveroqtnostite4.jpg]]
  
* Броят на опитите не е фиксиран;
 
* Опитите продължават до първи успех;
 
* Вероятността за успех остава същата от опит към опит;
 
  
[[разпределение на Вейбул|Разпределението на Вейбул]] е по същество семейство разпределения, което може да приеме свойствата на няколко други разпределения. Например в зависимост от дефинираната форма разпределението на Вейбул може да [[апроксимира]] Експоненциалното разпределение и това на Рейлай ([[Rayleigh]]). То е непрекъснато разпределение. С три параметъра – [[площ]], [[мащаб]] и [[форма]]. Разпределението на Вейбул намира приложение при оценката на надеждност.
 
  
[[Бета разпределение]]то е много гъвкаво разпределение подходящо за моделиране на вероятности, базирани на статистиката на [[Бейз]]. Статистиката на Бейз и бейзовата гледна точка са в основата на съвременната статистическата теория на решенията и оценката на риска. Също така се използва за описание на [[емпирични данни]] и прогнозиране на случайно поведение на количества и части от цялото чрез представянето им като проценти и дроби. То е непрекъснато вероятностно разпределение с три параметъра – алфа, бета и мащаб. Две са условията, на които отговаря Бета разпределението:
+
Относителната честота е [[характеристика]], която се получава след като е проведен експеримент или са извършени съответни [[наблюдение|наблюдения]]. В редица случаи са необходими характеристики за настъпването на разглежданите събития, които се получават преди да е извършен експеримент и описват техните [[честотни]] закономерности. Такива характеристики се наричат вероятности. Ако се разглежда събитие А, принадлежащо на [[множеството]] F (А F), вероятностната мера за събитие А се бележи с Р(А).
  
* Неизвестната променлива е случайна стойност между 0 и мащаба;
+
За вероятностите важи следното:
* Формата на разпределението се задава с две положителни стойности - a ,b
 
  
[[Хипергеометричното разпределение]] е подобно на Биномното; и двете описват броят на успехите при фиксиран брой опити. Обаче при Биномното разпределение опитите са независими, докато при хипергеометричното – вероятността се променя при следващ опит и са опити без заместване. Това е дискретно разпределение с три параметъра – вероятност, опити и размер. Параметъра вероятност определя началната вероятност. Три са условията, на които отговаря хипергеометричното разпределение:
+
#) За всяко събитие А F 0 Р(А) 1.
  
* Пълният брой елементи (размера) е фиксиран – крайна [[популация]];
+
#) Вероятността за поява на достоверно събитие е равна на [[единица]], т.е. Р() = 1.
* Примерния размер (броя на опитите) представя част от популацията;
 
* Известната начална вероятност за успех се променя плавно след всеки опит.
 
  
Със смесеното разпределение, може да се опише серия от индивидуални стойности, дискретни интервали или непрекъснати интервали за уникална ситуация, която не може да се опише с друг тип разпределение. Параметрите на това разпределение са начална граница, крайна граница, вероятност и стъпка.
+
#) Вероятността на сума от множество две по две несъвместими събития Aj F, j = 1,2, ... , е равна на сумата от техните вероятности
  
Вероятностните разпределения са изведени като теоретично описание на явления и събития наблюдавани в природата и в живота през мирно и военно време и обобщават в себе си вековен опит. [[теоретична вероятност|Теоретичните вероятностни]] разпределения със съвременните компютърни технологии могат да бъдат превърнати в генератори на данни в случаите, когато трябва да решим даден проблем, но разполагаме само с основни познания, експертни мнения и опита от аналогии с други обекти. Тяхното използване лежи в основата на съвременното компютърно моделиране и реализация на Монте Карло метода.
+
[[File:teoriqnaveroqtnostite5.jpg|100px]]
  
С посочените примери илюстрирахме метода Монте Карло и видяхме как се определят математическото очакване и [[дисперсия]]та на интересуващи ни параметри за да изследваме качеството и надеждността на системи и изделия състоящи се от множество елементи или да оценим [[риск]]а от индустриален инцидент. Приведените примери показват, че методиката за разчет е достатъчно проста по идея. На практика почти във всички случаи се знаят или могат да се определят вероятностните характеристики на отделните елементи на интересуващата ни система. Случайността на параметрите се отчита в общото цяло на системата по пътя на моделирането. Използвайки съвременните технологии за [[компютърно имитационно моделиране]] чрез симулиране на множество варианти на работа на дадена система може да се получи ценна практическа информация за нейното поведение при стотици и хиляди комбинации на входните й параметри и вътрешни състояния.
 
  
==Методът „Монте Карло” и съвременните технологии==
 
  
 +
Вероятността за противоположното събитие (A) е
  
 +
[[File:teoriqnaveroqtnostite6.jpg]]
  
Изводите от направения литературен преглед относно възможностите на компютърното симулиране са следните:
 
  
Компютърните симулационни модели:
 
  
* Позволяват да се извлекат изводи за нова система без да се налага да се построи реално, или да се направят промени в съществуваща без да се разстройва работата й;
+
като Р(А) + Р(A) = 1.
* Позволяват на [[ръководител]]я да визуализира операциите на нова или съществуваща система при различни условия;
+
 
* Позволяват да се види начина, по който различни компоненти си взаимодействат и как това влияе върху общата системна производителност;
+
Съвкупността от множеството , определеното на него множество F и вероятностната [[мера]] Р(А) за поява на събитие А, принадлежащо на F, образува вероятностно пространство и е част от [[математически модел|математическия модел]] на експеримента. За определяне на вероятностите се използват т.нар. [[лапласов метод|класически (лапласов) метод]], [[статистически метод]] и метод на експертните оценки (субективните вероятности).
* Позволяват общо вникване в същността на процеса;
+
 
* Позволяват разпознаване на специфични проблеми и проблемни области в изследваната система;
+
 
* Подпомагат разработването на специфични политики и планове за [[процес]]и;
+
 
* Подобряват системната ефективност.
+
===Класически (лапласов) метод===
 +
 
 +
 
 +
 
 +
Разглежда се вероятностен експеримент Е с пространство от елементарни събития = {i / i = 1, 2, ... , n}. Събитията i са несъвместими помежду си и са равно възможни, т. е. Р(1) = Р(2) = ... = Р(n). В резултат от експеримента се появява едно от тях. Вероятността за поява на което и да е елементарно събитие е
 +
 
 +
[[File:teoriqnaveroqtnostite7.jpg]]
 +
 
 +
 
 +
 
 +
При този метод на задаване на вероятностите, вероятността на което и да е събитие А, състоящо се от появата на едно от елементарните събития 1, 1, ... , m , m < n, т.е. събитието А = {1, 2, , m} се изчислява по формулата
 +
 
 +
[[File:teoriqnaveroqtnostite8.jpg]]
 +
 
 +
 
 +
 
 +
Събитията 1, 1, ... , m се наричат благоприятни за събитие А. [[Формула]]та показва, че ако всички елементарни събития j са равно възможни, вероятността за настъпване на което и да е събитие А F е равна на отношението на броя на елементарните събития, благоприятни за събитие А, към общия брой на елементарните събития.
 +
 
 +
===Статистически метод===
 +
 
 +
 
 +
 
 +
При класическия метод се предполага, че е известен броят на елементарните събития, и че те са равно възможни. Това не винаги е изпълнено. В такива случаи се прилага [[статистика|статистическият]] метод. Нека е произволно множество от елементарни събития на експеримента Е. Експериментът се повтаря n пъти при еднакви условия, представени като комплекс от условия S. Относителната честота за поява на събитие А е
 +
 
 +
[[File:teoriqnaveroqtnostite9.jpg]]
 +
 
 +
 
 +
 
 +
При този метод относителната честота на събитие А се приема като оценка на вероятността за поява на това събитие с точност , т.е. Р(А) = W(A) + , където е грешката, която се допуска. В основата на метода е заложена [[хипотеза]]та, че съществува такава константа Р(А), около която се групират относителните честоти
 +
 
 +
[[File:teoriqnaveroqtnostite10.jpg]]
 +
 
 +
 
 +
 
 +
изчислени за многократни повторения на експеримента при еднакви условия S. С нарастването на n, относителните честоти се приближават към вероятността Р(А).
 +
 
 +
 
 +
 
 +
 
 +
 
 +
===Метод на експертните оценки===
 +
 
 +
 
 +
 
 +
Прилага се когато повторението на експеримента при еднакви условия е невъзможно, причина за което могат да бъдат естеството на [[процес]]а, [[разход]]ите за експеримента, [[ограничеността]] на времето и др. В тези случаи се прибягва към [[експертна оценка|експертни оценки]] на вероятностите за поява на даденото [[събитие]].Нека експертизата обхваща n [[експерт]]и и m от тях са изказали мнение, че събитие А ще се появи. Тогава вероятността за поява на събитието се определя като
 +
 
 +
[[File:teoriqnaveroqtnostite11.jpg]]
 +
 
 +
 
 +
 
 +
В този случай мненията на [[експерт]]ите се приемат за еднакво [[правдоподобни]]. Основание за това е предположението за еднаква степен на квалификация, информираност и неангажираност. В конкретни случаи могат да се дадат различни [[тегло|тегла]] на експертите. При метода на експертните оценки се ползва понятието [[субективна вероятност]]. По предложение на [[Кейнс]] вероятността за поява на събитие А, оценявана от експерти, разполагащи с комплекс информация I, се определя като P(A,I). Ако информацията се промени, например на I', тогава вероятността се определя като P(A,I').
 +
 
 +
==Действия с вероятности==
 +
 
 +
 
 +
 
 +
===Събиране на вероятности (вероятност на сума от събития)===
 +
 
 +
 
 +
 
 +
Събиране на вероятности се извършва като се има предвид дали събитията са съвместими или не са.
 +
 
 +
а) При несъвместими събития
 +
 
 +
Нека събитията A1, А2, ... , Ak са несъвместими помежду си. Сумата на вероятностите на тези събития представлява вероятността за поява на кое да е от тях, която се означава Р(А1 или А2 или ... или Ak). Тя се определя като
 +
 
 +
[[File:teoriqnaveroqtnostite12.jpg]]
 +
 
 +
 
 +
 
 +
б) При съвместими събития
 +
 
 +
Ако две събития, събитие А и събитие В, са съвместими, което означава, че тяхното произведение не е празно множество, вероятността за поява на което и да е от тях е равна на сумата от техните вероятности минус вероятността за съвместното им появяване, т.е.
 +
 
 +
P(AB) = P(A) + P(B) P(AB) (4.2.)
 +
 
 +
За три произволни съвместими събития А, В и С вероятността за сумата на тези събития се определя като
 +
 
 +
P(ABC) = P(A) + P(B) + P(C) P(AB) - P(AC) - P(BC) + P(ABC)
 +
 
 +
===Условна вероятност===
 +
 
 +
 
 +
 
 +
В определени случаи появата на дадено събитие зависи от това дали друго събитие, с което то е в определена вероятностна връзка, сее осъществило. Появата на второто събитие изменя условията и съответно вероятността за поява на разглежданото събитие. В общия случай нека е дадено вероятностното пространство (, F, Р) и А и В са произволни събития, като вероятността за поява на събитие В е положителна величина, т. е. Р(В)>О. Условната вероятност за настъпване на събитие А при [[условие]], че е настъпило събитие В, се означава обикновено с P(A/В). Тя се определя като отношение на произведението (сечението) на двете събития към вероятността за настъпване на събитие В
 +
 
 +
[[File:teoriqnaveroqtnostite13.jpg]]
 +
 
 +
 
 +
 
 +
По подобен начин се определя условната вероятност за настъпване на събитие В при условие, че е настъпило събитие А
 +
 
 +
[[File:teoriqnaveroqtnostite14.jpg]]
 +
 
 +
 
 +
 
 +
Събитията А и В са независими тогава, когато P(AB) = P(A)P(B). Това означава, че P(A/B) = P(A) и P(B/А) = P(B).
 +
 
 +
===Вероятност на сложно събитие (вероятност на произведение на събития)===
 +
 
 +
 
 +
 
 +
[[Сложно събитие]] е събитието, което се състои в съвместното осъществяване, едновременно или последователно, на две или повече събития. Ако това са събитията А и В, сложно събитие е едновременното появяване и на А, и на В, означено с (А и В) или, което е същото, с (АВ). Вероятността за сложното събитие (А и В) се определя като произведение на вероятността за настъпване на събитие А по условната вероятност за настъпване на събитие В при условие, че е настъпило събитие А
 +
 
 +
Р(АВ) = Р(А)Р(В/А) (4.5)
 +
 
 +
или като произведение на вероятността за настъпване на събитие В по условната вероятност за настъпване на събитие А при условие, че е настъпило събитие В
 +
 
 +
Р(АВ) = P(В)P(A/В).
 +
 
 +
===Теорема за умножение на вероятности.===
 +
 
 +
 
 +
 
 +
Вероятността на произведението на две събития е равна на произведението на вероятността на едното от тези събития по условната вероятност на второто събитие, изчислена при условие, че първото е настъпило:
 +
 
 +
Р(А В) = Р(А)Р(В/А) = Р(В)Р(А/В)
 +
 
 +
- Теорема за пълната вероятност. Нека събитие К да може да се появи с едно от несъвместимите събития Ai, i = 1,2, ... ,n, обединението на които съвпада с пространството от елементарни събития , т.е.
 +
 
 +
= {А1 U А2 U U Аn}, Ai U Aj и Ai Aj = , ij.
 +
 
 +
Събитията Ai обикновено се наричат хипотези.
 +
 
 +
Теоремата за пълната вероятност гласи, че вероятността на събитие К, което може да се появи съвместно с която и да е от хипотезите Аi, i = 1,2,..., n, е равна на сумата от произведенията, две по две, на вероятностите на всяка от тези хипотези по съответстващите им условни вероятности за настъпване на събитие К
 +
 
 +
[[File:teoriqnaveroqtnostite15.jpg]]
 +
 
 +
 
 +
 
 +
===Формула на Бейс===
 +
 
 +
 
 +
 
 +
Нека събитията А1, А2, ..., Аn представляват всички възможни изходи от предприето действие и събитие К е резултат, който може да настъпи с настъпването на което и да е от тези събития. Съгласно формулата за умножение на вероятности, вероятността за настъпване на събитие К и което и да е от събитията Аi , i = 1, 2, ... , n, е равна на
 +
 
 +
P(Ai K) = P(Ai)P(K/Ai) = P(K)P(Ai / K). (4.8)
 +
 
 +
 
 +
 
 +
От тук се определя
 +
 
 +
[[File:teoriqnaveroqtnostite16.jpg]]
 +
 
 +
 
 +
 
 +
Съгласно теоремата за пълната вероятност
 +
 
 +
[[File:teoriqnaveroqtnostite17.jpg]]
 +
 
 +
 
 +
 
 +
Като се замести се получава:
 +
 
 +
[[File:teoriqnaveroqtnostite18.jpg]]
 +
 
 +
 
 +
 
 +
==Елементи на теория на вероятностите==
 +
 
 +
 
 +
 
 +
===Случайни величини===
 +
 
 +
 
 +
 
 +
Под [[случайна]] величина ще разбираме величината приемаща в резултат на опита някакво числово значение. Пример: Числото, падащо при хвърляне на зар. Тази случайна величина може да приеме едно от числата 1, 2, 3, 4, 5 и 6. Случайна величина, приемаща краен брой значения се нарича дискретна случайна величина. Случайна величина, приемаща значенията си в някакъв интервал се нарича непрекъсната. За да се характеризира случайна величина е необходимо да се укажат възможните нейни значения (по-нататък ще разглеждаме дискретни случайни величини). Но това не е достатъчно: необходимо е да знаем колко често приема [[значенията]] си. Това най-добре се характеризира с вероятността за отделните нейни значения т.е. за случайната величина X трябва да укажем значенията и x1, x2, ., но и вероятностите за събитията X=xi.
 +
 
 +
pi=P(X=xi)
 +
 
 +
състоящи се в това, че X приема значение xi. Съответствието xi pi се нарича закон за разпределение на случайната величана. [[Математическо очакване]] и дисперсия на случайната величина. Математическо очакване на сл. величина X се нарича числото M(X) определено с равенството
 +
 
 +
[[File:teoriqnaveroqtnostite20.jpg]]
 +
 
 +
 
 +
 
 +
Задача. Компютърът може да се намира в едно следните състояния:
 +
 
 +
S1 изправен работи; S2 неизправен, спрян, търси се повреда; S3 неизправността е открита, тече ремонт; S4 ремонтът е завършен, тече подготовка за пуска на компютъра. Средното време за безотказна работа на компютъра е половин ден. За ремонта отиват средно 6 ч. Търсенето на повредата продължава средно половин час (0,5 ч.). След ремонта подготовката за пускане е средно 1 ч. Да се намерят пределните вероятности на състоянията на компютрите.
 +
 
 +
 
 +
 
 +
Решение.
 +
 
 +
Граф на състоянията на интензивностите 
 +
 
 +
[[File:teoriqnaveroqtnostite21.jpg]]
 +
 
 +
 
 +
 
 +
Времената t1,t2,t3 и t4 ще ги представим в дни
 +
 
 +
[[File:teoriqnaveroqtnostite22.jpg]]
 +
 
 +
 
 +
 
 +
(Време за безотказна работа)
 +
 
 +
[[File:teoriqnaveroqtnostite23.jpg]]
 +
 
 +
 
 +
 
 +
(Време за търсене на повреда)
 +
 
 +
[[File:teoriqnaveroqtnostite24.jpg]]
 +
 
 +
 
 +
 
 +
(Време за ремонт)
 +
 
 +
[[File:teoriqnaveroqtnostite25.jpg]]
 +
 
 +
 
 +
 
 +
(Време за подготовка за пускане)
 +
 
 +
От системата намираме:
 +
 
 +
[[File:teoriqnaveroqtnostite30.jpg]]
 +
 
 +
 
 +
 
 +
[[File:teoriqnaveroqtnostite31.jpg]]
 +
 
 +
 
 +
 
 +
[[File:teoriqnaveroqtnostite32.jpg]]
 +
 
 +
 
 +
 
 +
След което
 +
 
 +
[[File:teoriqnaveroqtnostite33.jpg]]
 +
 
 +
 
 +
 
 +
p1 = 0.615, p2 =0.026, p3 = 0,308, p4 = 0,0051
  
Компютърните симулационни модели не могат:
 
  
* Да [[оптимизация|оптимизират]] ... те могат само да генерират резултати от “КАКВО-АКО” запитвания;
 
* Да получат коректни резултати от неточни данни;
 
* Да опишат системни характеристики, които не са били заложени в модела;
 
* Да решават проблеми ... те могат само да предоставят [[информация]], която подпомага процеса на изработване на решение;
 
* Да дават лесни отговори на комплексни проблеми.
 
  
 
==Вижте още==
 
==Вижте още==
Ред 209: Ред 404:
  
  
* [[Монте Карло методи в областта на финансите]]
+
* [[алтернатива]]
* [[Алгоритмично решение]]
+
* [[вероятност]]
* [[Вероятност]]
+
* [[ефективност]]
* [[Експеримент]]
+
* [[хипотеза]]
* [[Ефкетивност]]
+
* [[цел]]
* [[Избор (Теория на решенията)]]
+
* [[теорема на Бейс]]
* [[Задача за решаване]]
+
* [[случайно число]]
* [[Информация]]
+
* [[следствие]]
* [[Модел]]
+
* [[решение]]
* [[Неопределеност]]
+
* [[риск]]
* [[Научен подход]]
+
* [[оптимизация]]
* [[Оптимизация]]
+
* [[планиране]]
* [[Предвиждане]]
+
* [[информация]]
* [[Риск]]
+
* [[критерий]]
* [[Решение]]
+
* [[задача за решаване]]
* [[Симулация]]
+
* [[избор]]
  
==Книги==
+
==Източници==
  
  
  
* B. Kaye, Chaos & Comlexity, VCH, Weinheim, New York 1993
+
* Pierre Simon de Laplace (1812). Analytical Theory of Probability.
* И .М. Соболь, Метод Монте Карло. Популярньiе лекции по математике вьiпуск 46, Москва “Наука”, Главная редакция физико-математической лератуьi, 1978.
+
* Andrei Nikolajevich Kolmogorov (1950). Foundations of the Theory of Probability.
* С.М. Ермаков, Г.А.Михайлов, Курс статистического моделирования, Москва, “Наука”, 1976.
+
* Patrick Billingsley (1979). Probability and Measure. New York, Toronto, London: John Wiley and Sons.
* С. М. Ермаков, Метод Монте Карло и смежньiе вопросьi, Москва, “Наука” 1971.
+
* Olav Kallenberg; Foundations of Modern Probability, 2nd ed. Springer Series in Statistics. (2002). 650 pp.  
* Н. П. Бусленко и др. (по редакцией Ю. А. Шрейдера), Метод статистических изпитаний (метод Монте Карло), СМБ, Москва, Физматгиз, 1962.
+
* Henk Tijms (2004). Understanding Probability. Cambridge Univ. Press.
* D. D. McCracken, “The Monte Carlo Method”, Scientific American, May 1955, pp. 162-165, reprinted in Computers and Computations, W. H. Freeman, San Farncisco 1971.
+
* Olav Kallenberg; Probabilistic Symmetries and Invariance Principles. Springer -Verlag, New York (2005). 510 pp.  
* G. Gamow, “One, Two, Three … Infinity, Viking Press, New York, 1974.
+
* Gut, Allan (2005). Probability: A Graduate Course. Springer-Verlag.
* Sturgul J.R. Mine Design - Examples Using Simulation, SME Inc., 2000
+
* ET. Jaynes. ''Probability Theory: The Logic of Science'' Cambridge University Press, (2003).
* Russell A.C. Sizing Up Simulation, Jurnal of Mining Engineering, No:26, August 2001.
+
* Stephen M. Stigler (1986) The history of statistics. Harvard University press
* Smith R.D. Enciclopedia Article. Simulation Article, 4th Edition, N.Y. July 2000
+
* Richard T. Cox, Algebra of Probable Inference, The Johns Hopkins University Press
* MINExpo 2000 – Web Site.
+
* Bishop, CM., Pattern Recognition and Machine Learning. Springer, 2007
* Стефанов Т.П., В.Т. Томов, И.С. Велчев. Анализ Вениляционньiх систем с помощью ЭБМ 16. Международная конференция НИИ по безопасности в горной промишлености, Вашингтон, США, 1975.
+
* Stephen M. Stigler (1986) The history of statistics. Harvard University press
* Власева Е., Т. Стефанов. Изграждане на умение за анализ на поведението на рудничната вентилационна система чрез компютърно симулиране. Научно-техническа конференция по охрана на труда в подземните и откритите рудници и кариерите, Варна, 8-11 юни, 1998.
+
* Wolpert, RL. (2004) A conversation with James O. Berger, Statistical science
* Elena Vlaseva. Control volumes technique applied to gas dynamical problems in underground mines. Proceeding of Finite Volumes for Complex Applications II. Hermes Science Publications, Paris, 1999
+
* Berger, James O (1985). Statistical Decision Theory and Bayesian Analysis. Springer Series in Statistics (Second ed.). Springer-Verlag
 +
* de Finetti, Bruno. "Probabilism: A Critical Essay on the Theory of Probability and on the Value of Science," (translation of 1931 article) in Erkenntnis, volume 31, September 1989
 +
* Hacking, I (1988) "Slightly More Realistic Personal Probability". 1967 article partly reprinted in: Gärdenfors, Peter and Sahlin, Nils-Eric. (1988) Decision, Probability, and Utility: Selected Readings. 1988. Cambridge University Press
 +
* Jaynes E.T. (2003) Probability Theory: The Logic of Science, CUP
 +
* Morgenstern, Oskar (1978). "Some Reflections on Utility". In Andrew Schotter. Selected Economic Writings of Oskar Morgenstern. New York University Press
 +
* Pfanzagl, J (1967). "Subjective Probability Derived from the Morgenstern-von Neumann Utility Theory". In Martin Shubik. Essays in Mathematical Economics In Honor of Oskar Morgenstern. Princeton University Press
 +
 
 +
 
  
 
==Външни препратки==
 
==Външни препратки==
Ред 249: Ред 451:
  
  
* [http://translate.googleusercontent.com/translate_c?hl=bg&langpair=en|bg&u=http://www.puc-rio.br/marco.ind/monte-carlo.html&rurl=translate.google.bg&usg=ALkJrhjNWRXn6HoxFL0zvPbmzylJtzVUHg Реални модели и опции за прилагането на метода Монте Карло]
+
* [http://www.mgu.bg/drugi/ebooks/nikolina/ Случйни величини]
* [http://translate.googleusercontent.com/translate_c?hl=bg&langpair=en|bg&u=http://www.palisade.com/risk/monte_carlo_simulation.asp&rurl=translate.google.bg&usg=ALkJrhgXqOfUOQgD6QY1IFhOAl6yT68w8w  Монте Карло симулация]
+
* [http://www.math10.com/forumbg/viewforum.php?f=33 Математическа статистика]
* [http://translate.googleusercontent.com/translate_c?hl=bg&langpair=en|bg&u=http://www.global-derivatives.com/index.php?option=com_content&task=view&id=21&rurl=translate.google.bg&usg=ALkJrhjJmrfNpnLo1NgA05R2uE2A33lTrQ Финансова математика (речник)]
+
* [http://books.balkanatolia.com/c/sl-b/cid-27/p-p/id-8858/teoria-na-veroiatnostite-i-statistika-krasimir-kalinov.html Теория на вероятностите и статистика]
* [http://www1.ecs.uni-ruse.bg/fbm/bsa/bsa-L4.pdf Бизнес симулация]
+
* [http://www.kaminata.net/3-t44485.html Висша математика]
* [http://www.kaminata.net/excel-t45794.html Имитационно моделиране]
+
* [http://www.math10.com/forumbg/viewtopic.php?t=2779 Теория на числата]
* [http://parallel.bas.bg/CE_SuperCA/REP_Dec09_WP4_F.pdf Монте Карло методи за анализ на чувствителността на големи математически модели]
+
* [http://arsstat.com/MitakaStatLabVirtual/VTULectures/ralitsa/RALITSASLAVOVAANDDIMITERTSVETKOVANNUALNVU2010.htm Надеждност и вероятност]
 +
* [http://docs.google.com/viewer?a=v&q=cache:L4d3IFV1yGgJ:www.mgu.bg/drugi/ebooks/nikolina/lessons/content.pdf+%D1%82%D0%B5%D0%BE%D1%80%D0%B8%D1%8F+%D0%BD%D0%B0+%D0%B2%D0%B5%D1%80%D0%BE%D1%8F%D1%82%D0%BD%D0%BE%D1%81%D1%82%D1%82%D0%B0&hl=bg&gl=bg&pid=bl&srcid=ADGEESgTUzm4cR-6oZ2DiMJyBl5YBqusscWc5AX8-Q0O_aJ2kAeIPIHnjqyFfJXDM_G4K8a03AYTNoFV9k8wLsW-scMVILfhTKEOZxsqqbXt70daXOTxFsBZS0i0MQwzhzmqQPoUHTXb&sig=AHIEtbSFKDDHcMc8OER5KzXfV2zwQCk4gw Вероятности и статистика]
 +
* [http://tuj.asenevtsi.com/Informatica2/I013.htm Аспекти на информацията]
 +
 
 
[[category:Управленски решения и риск]]
 
[[category:Управленски решения и риск]]

Версия от 16:57, 24 януари 2011

Теорията на вероятностите се занимава с изследване закономерностите при случайните събития.

Всяко явление протича при известни условия. Обратно, при определени условия протича някакво явление, но ние не сме в състояние да опишем всички условия, при които протича дадено явление.

Ние определяме само някои условия за протичането на дадено явление, което ще наричаме събитие. При реализирането на тези условия даденото събитие може да се сбъдне, а може и да не се сбъдне. Такова събитие, което при даден комплекс от условия може да се сбъдне, а може и да не се сбъдне, наричаме случайно.


ВЪВЕДЕНИЕ В ТЕОРИЯ НА ВЕРОЯТНОСТИТЕ

Основни понятия на теория на вероятностите

Най-често се използва понятието опит (експеримент). Опитът се състои в осъществяване на определен комплекс от условия, в резултат на което се случва нещо, което не може да бъде предсказано предварително със сигурност. При многократното повторение на опита могат да бъдат получени различни резултати (изходи), които се представят като поява на определено събитие. Събитията се разграничават по видове, в съответствие с особеностите на тяхното проявление.


Видове събития


Сигурни(достоверни) събития


Ако при многократно повтаряне на даден комплекс от условия се реализира едно и също събитие, то се нарича сигурно или достоверно събитие.

Ето няколко примера. Ако подхвърлим с ръце еди тежък предмет нагоре, този предмет ще падне. Падането на предмета е сигурно събитие.

Пример 2: при падането на монета, върху която различаваме лицева и гербова страна, върху хоризонтална равнина тя не застава върху равнината отвестно, а пада върху една от страните си. Падането на монетата върху една от страните й е сигурно събитие.

Пример 3: при подхвърляне върху хоризонтална равнина едно кубче застава върху една от стените си. Следователно сигурно събитие е падането на куба върху една стена. За конкретност вместо куб ще говорим за зар(кубче което се използува при хазартни игри), стените на които са номерирани от 1 до 6 посредством съответен брой точки. Събитието падане на зара върху стена с брой на точките, по-малък от 7 е сигурно събитие. Сигурните събитя ще означаваме с буквата U.


Невъзможни събития


Ако при реализиране на определен комплекс от условия дадено събитие не се сбъдва, то се нарича невъзможно събитие.

Пример за невъзможно събитие е при падане върху хоризонтална равнина една монета да застане отвесно. Друг такъв пример е заставането на зар върху ръб или даже върху връх. Невъзможно събитие е падането на зара върху стена със 7 точки. Невъзможните събития ще означаваме с буквата V.


Противоположни събития


Две събития ще наричаме противоположни, ако се реализират при еди и същи комплекс от условия и ако винаги, когато едното събитие се реализира, другото не се реализира.

Пример за противоположни събития са разгледаните по-горе сигурни и невъзможни събития. Падането на една монета върху лицевата й страна и падането й върху гербовата й срана, разгледани като две събития, са също противоположни, защото монетата не може да е паднала върху лицевата си страна, ако е паднала върху гербовата си страна и обратно.

Падането на зара върху стена с четен брой точки и върху стена с нечетен брой е също невъзможно събитие. Не е трудно да се разбере, че всяко събитие има свое противоположно. Противоположното на дадено събитие можем да дефинираме като събитие, което се реализира, ако даденото не се реализира, или обратно, което остава нереализирано в резултат на реализирането на даденото събитие.

Ако с А означим даденото събитие, то с А ще означаваме неговото противоположно. Логически казано, А е не А и А е не А

Съгласно въведените по-горе означения за сигурно и невъзможно събитие имаме равенствата

U = V, V = U


Несъвместими и съвместими събития

Събития, които не могат да се реализират (едновременно) при един и същ комплекс от условия, се наричат несъвместими. Несъвместими събития са например противоположните. Но докато при противоположните се реализира едно от тях, щом другото не се е реализирало, при несъвместимите събития това не е задължително. Две несъвместими събития могат и едновременно да не се реализират. Тяхното едновременно реализиране обаче е отречено по определение. При хвърлянето на един зар несъвместими са например следните две събития: падането на зара върху стена с една точка и падането му върху стена с четен брой точки. Съвместимите събития ще определим като събития, които не са несъвместими.Пример на съвместими събития е падането на зар върху стена с една точка и падането на зар върху стена с нечетен брой точки


Елементарни и сложни събития

При даден комплекс от условия се реализират различни събития. Някои от тях могат да бъдат определени като основни, неизразени посредством други събития, но посредством които се изграждат останалите събития. Такива основни събития ще наричаме елементарни а останалите събития – сложни. Ще отбележим, че както елементарните, така и сложните събития се обуславят единствено от комплекса от условия. Не е изключено едно и също събитие да може да реализира при два различни комплекса от условия, но при единия комплекс то да е елементарно, а при другия сложно. По принцип смятаме, че при различни комплекси от условия се реализират различни събития.

Примери за елементарни и сложни събития:

При хвърлянето на една монета елементарните събития са две: падане върху лизевата страна и падане върху гербовата страна. Също така при хвърляне на зар имаме 6 елементарни събития, за 6-те различни стени. Падането на зара върху стена с четен брой точки е сложно събитие, което се състои в реализирането на едно от елементарните събития за падане върху стена с 2, 4 или 6 точки. Ако хвърлим две монети едновременно, то елементарните събития са двойка елементарни събития за една монета. Ще означаваме с Л и Г съответно събитията падане на монета върху лицевата и падане върху гербовата й страна. За падане на 2 монети имаме следните четири елементарни събития: ЛЛ, ЛГ, ГЛ, ГГ.

Сложно събитие е падане на двете монети върху различни стени, тогава, когато едната монета е паднала върху Л страна, а другата върху Г или обратно(ЛГ или ГЛ). Така това сложно събитие се реализира само когато се реализира едно от двете елементарни събития(ЛГ или ГЛ), но то не е тъждествено с тях, а само е изградено посредством тези елементарни събития. И така при даден комплекс от условия всяко събитие може да се разглежда като изградено посредством елементарни събития. Ние ще разгледаме всяко събитие като множество, чиито елементи са елементарни събития.


Действия със събития. Поле от събития

Известни са действията с множества. Тъй като ние разглеждаме събитията като множества, то имаме възможност лесно да извършваме действия и със събития.

Преди всичко, ако А е дадено събитие, изградено посредством елементарните събития а1, а2, а3,... аn,то смятаме, че събитието А се сбъднало, ако се е сбъднало поне едно от изграждащите го елементарни събития а1, а2, а3,... аn.

Еквивалентни събития

Казваме, че събитието А следва от събитието В или че А се съдържа в В, или В съдържа А, което записваме така; А ?(само че знака е легнал) В или В ?(само че знака е легналА). Занапред, когато говорим за различни събития , ще смятаме, че имаме един и същи комплекс от условия, ако не е казано противното. Ако А и В са две събития, такива че А следва от В и В следва от А, т.е. А?В и В?А Ще казваме че А и B са еквивалентни(равносилни) или просто равни и това ще записваме така А = B


Вероятност. Методи за определяне на вероятности

Когато експериментът е проведен n пъти и събитие А се е появило m пъти, отношението

Teoriqnaveroqtnostite1.jpg


показва какъв дял от общия брой случаи се пада на събитие А. Тази величина се нарича относителна честота или честост, а числото m - честота. Най-често относителната честота се записва с

Teoriqnaveroqtnostite2.jpg


Относителната честота има следните свойства:

1) 0W(A)1, т.е. тя е неотрицателна величина.

2) Ако събитие К е достоверно (сигурно) събитие, W(K) = 1, т.е относителната честота на достоверно събитие, в случая събитие К, е равна на единица. Важи и обратното, че ако дадено събитие има относителна честота равна на единица, то е достоверно събитие.

3) За редицата от събития A1, А2, ... , Ak

Teoriqnaveroqtnostite3.jpg


За събитията А и В с честоти m(А) и m(В) и честости W(A) и W(В) важи

Teoriqnaveroqtnostite4.jpg


Относителната честота е характеристика, която се получава след като е проведен експеримент или са извършени съответни наблюдения. В редица случаи са необходими характеристики за настъпването на разглежданите събития, които се получават преди да е извършен експеримент и описват техните честотни закономерности. Такива характеристики се наричат вероятности. Ако се разглежда събитие А, принадлежащо на множеството F (А F), вероятностната мера за събитие А се бележи с Р(А).

За вероятностите важи следното:

  1. ) За всяко събитие А F 0 Р(А) 1.
  1. ) Вероятността за поява на достоверно събитие е равна на единица, т.е. Р() = 1.
  1. ) Вероятността на сума от множество две по две несъвместими събития Aj F, j = 1,2, ... , е равна на сумата от техните вероятности

Teoriqnaveroqtnostite5.jpg


Вероятността за противоположното събитие (A) е

Teoriqnaveroqtnostite6.jpg


като Р(А) + Р(A) = 1.

Съвкупността от множеството , определеното на него множество F и вероятностната мера Р(А) за поява на събитие А, принадлежащо на F, образува вероятностно пространство и е част от математическия модел на експеримента. За определяне на вероятностите се използват т.нар. класически (лапласов) метод, статистически метод и метод на експертните оценки (субективните вероятности).


Класически (лапласов) метод

Разглежда се вероятностен експеримент Е с пространство от елементарни събития = {i / i = 1, 2, ... , n}. Събитията i са несъвместими помежду си и са равно възможни, т. е. Р(1) = Р(2) = ... = Р(n). В резултат от експеримента се появява едно от тях. Вероятността за поява на което и да е елементарно събитие е

Teoriqnaveroqtnostite7.jpg


При този метод на задаване на вероятностите, вероятността на което и да е събитие А, състоящо се от появата на едно от елементарните събития 1, 1, ... , m , m < n, т.е. събитието А = {1, 2, , m} се изчислява по формулата

Teoriqnaveroqtnostite8.jpg


Събитията 1, 1, ... , m се наричат благоприятни за събитие А. Формулата показва, че ако всички елементарни събития j са равно възможни, вероятността за настъпване на което и да е събитие А F е равна на отношението на броя на елементарните събития, благоприятни за събитие А, към общия брой на елементарните събития.

Статистически метод

При класическия метод се предполага, че е известен броят на елементарните събития, и че те са равно възможни. Това не винаги е изпълнено. В такива случаи се прилага статистическият метод. Нека е произволно множество от елементарни събития на експеримента Е. Експериментът се повтаря n пъти при еднакви условия, представени като комплекс от условия S. Относителната честота за поява на събитие А е

Teoriqnaveroqtnostite9.jpg


При този метод относителната честота на събитие А се приема като оценка на вероятността за поява на това събитие с точност , т.е. Р(А) = W(A) + , където е грешката, която се допуска. В основата на метода е заложена хипотезата, че съществува такава константа Р(А), около която се групират относителните честоти

Teoriqnaveroqtnostite10.jpg


изчислени за многократни повторения на експеримента при еднакви условия S. С нарастването на n, относителните честоти се приближават към вероятността Р(А).



Метод на експертните оценки

Прилага се когато повторението на експеримента при еднакви условия е невъзможно, причина за което могат да бъдат естеството на процеса, разходите за експеримента, ограничеността на времето и др. В тези случаи се прибягва към експертни оценки на вероятностите за поява на даденото събитие.Нека експертизата обхваща n експерти и m от тях са изказали мнение, че събитие А ще се появи. Тогава вероятността за поява на събитието се определя като

Teoriqnaveroqtnostite11.jpg


В този случай мненията на експертите се приемат за еднакво правдоподобни. Основание за това е предположението за еднаква степен на квалификация, информираност и неангажираност. В конкретни случаи могат да се дадат различни тегла на експертите. При метода на експертните оценки се ползва понятието субективна вероятност. По предложение на Кейнс вероятността за поява на събитие А, оценявана от експерти, разполагащи с комплекс информация I, се определя като P(A,I). Ако информацията се промени, например на I', тогава вероятността се определя като P(A,I').

Действия с вероятности

Събиране на вероятности (вероятност на сума от събития)

Събиране на вероятности се извършва като се има предвид дали събитията са съвместими или не са.

а) При несъвместими събития

Нека събитията A1, А2, ... , Ak са несъвместими помежду си. Сумата на вероятностите на тези събития представлява вероятността за поява на кое да е от тях, която се означава Р(А1 или А2 или ... или Ak). Тя се определя като

Teoriqnaveroqtnostite12.jpg


б) При съвместими събития

Ако две събития, събитие А и събитие В, са съвместими, което означава, че тяхното произведение не е празно множество, вероятността за поява на което и да е от тях е равна на сумата от техните вероятности минус вероятността за съвместното им появяване, т.е.

P(AB) = P(A) + P(B) P(AB) (4.2.)

За три произволни съвместими събития А, В и С вероятността за сумата на тези събития се определя като

P(ABC) = P(A) + P(B) + P(C) P(AB) - P(AC) - P(BC) + P(ABC)

Условна вероятност

В определени случаи появата на дадено събитие зависи от това дали друго събитие, с което то е в определена вероятностна връзка, сее осъществило. Появата на второто събитие изменя условията и съответно вероятността за поява на разглежданото събитие. В общия случай нека е дадено вероятностното пространство (, F, Р) и А и В са произволни събития, като вероятността за поява на събитие В е положителна величина, т. е. Р(В)>О. Условната вероятност за настъпване на събитие А при условие, че е настъпило събитие В, се означава обикновено с P(A/В). Тя се определя като отношение на произведението (сечението) на двете събития към вероятността за настъпване на събитие В

Teoriqnaveroqtnostite13.jpg


По подобен начин се определя условната вероятност за настъпване на събитие В при условие, че е настъпило събитие А

Teoriqnaveroqtnostite14.jpg


Събитията А и В са независими тогава, когато P(AB) = P(A)P(B). Това означава, че P(A/B) = P(A) и P(B/А) = P(B).

Вероятност на сложно събитие (вероятност на произведение на събития)

Сложно събитие е събитието, което се състои в съвместното осъществяване, едновременно или последователно, на две или повече събития. Ако това са събитията А и В, сложно събитие е едновременното появяване и на А, и на В, означено с (А и В) или, което е същото, с (АВ). Вероятността за сложното събитие (А и В) се определя като произведение на вероятността за настъпване на събитие А по условната вероятност за настъпване на събитие В при условие, че е настъпило събитие А

Р(АВ) = Р(А)Р(В/А) (4.5)

или като произведение на вероятността за настъпване на събитие В по условната вероятност за настъпване на събитие А при условие, че е настъпило събитие В

Р(АВ) = P(В)P(A/В).

Теорема за умножение на вероятности.

Вероятността на произведението на две събития е равна на произведението на вероятността на едното от тези събития по условната вероятност на второто събитие, изчислена при условие, че първото е настъпило:

Р(А В) = Р(А)Р(В/А) = Р(В)Р(А/В)

- Теорема за пълната вероятност. Нека събитие К да може да се появи с едно от несъвместимите събития Ai, i = 1,2, ... ,n, обединението на които съвпада с пространството от елементарни събития , т.е.

= {А1 U А2 U U Аn}, Ai U Aj и Ai Aj = , ij.

Събитията Ai обикновено се наричат хипотези.

Теоремата за пълната вероятност гласи, че вероятността на събитие К, което може да се появи съвместно с която и да е от хипотезите Аi, i = 1,2,..., n, е равна на сумата от произведенията, две по две, на вероятностите на всяка от тези хипотези по съответстващите им условни вероятности за настъпване на събитие К

Teoriqnaveroqtnostite15.jpg


Формула на Бейс

Нека събитията А1, А2, ..., Аn представляват всички възможни изходи от предприето действие и събитие К е резултат, който може да настъпи с настъпването на което и да е от тези събития. Съгласно формулата за умножение на вероятности, вероятността за настъпване на събитие К и което и да е от събитията Аi , i = 1, 2, ... , n, е равна на

P(Ai K) = P(Ai)P(K/Ai) = P(K)P(Ai / K). (4.8)


От тук се определя

Teoriqnaveroqtnostite16.jpg


Съгласно теоремата за пълната вероятност

Teoriqnaveroqtnostite17.jpg


Като се замести се получава:

Teoriqnaveroqtnostite18.jpg


Елементи на теория на вероятностите

Случайни величини

Под случайна величина ще разбираме величината приемаща в резултат на опита някакво числово значение. Пример: Числото, падащо при хвърляне на зар. Тази случайна величина може да приеме едно от числата 1, 2, 3, 4, 5 и 6. Случайна величина, приемаща краен брой значения се нарича дискретна случайна величина. Случайна величина, приемаща значенията си в някакъв интервал се нарича непрекъсната. За да се характеризира случайна величина е необходимо да се укажат възможните нейни значения (по-нататък ще разглеждаме дискретни случайни величини). Но това не е достатъчно: необходимо е да знаем колко често приема значенията си. Това най-добре се характеризира с вероятността за отделните нейни значения т.е. за случайната величина X трябва да укажем значенията и x1, x2, ., но и вероятностите за събитията X=xi.

pi=P(X=xi)

състоящи се в това, че X приема значение xi. Съответствието xi pi се нарича закон за разпределение на случайната величана. Математическо очакване и дисперсия на случайната величина. Математическо очакване на сл. величина X се нарича числото M(X) определено с равенството

Teoriqnaveroqtnostite20.jpg


Задача. Компютърът може да се намира в едно следните състояния:

S1 изправен работи; S2 неизправен, спрян, търси се повреда; S3 неизправността е открита, тече ремонт; S4 ремонтът е завършен, тече подготовка за пуска на компютъра. Средното време за безотказна работа на компютъра е половин ден. За ремонта отиват средно 6 ч. Търсенето на повредата продължава средно половин час (0,5 ч.). След ремонта подготовката за пускане е средно 1 ч. Да се намерят пределните вероятности на състоянията на компютрите.


Решение.

Граф на състоянията на интензивностите

Teoriqnaveroqtnostite21.jpg


Времената t1,t2,t3 и t4 ще ги представим в дни

Teoriqnaveroqtnostite22.jpg


(Време за безотказна работа)

Teoriqnaveroqtnostite23.jpg


(Време за търсене на повреда)

Teoriqnaveroqtnostite24.jpg


(Време за ремонт)

Teoriqnaveroqtnostite25.jpg


(Време за подготовка за пускане)

От системата намираме:

Teoriqnaveroqtnostite30.jpg


Teoriqnaveroqtnostite31.jpg


Teoriqnaveroqtnostite32.jpg


След което

Teoriqnaveroqtnostite33.jpg


p1 = 0.615, p2 =0.026, p3 = 0,308, p4 = 0,0051


Вижте още

Източници

  • Pierre Simon de Laplace (1812). Analytical Theory of Probability.
  • Andrei Nikolajevich Kolmogorov (1950). Foundations of the Theory of Probability.
  • Patrick Billingsley (1979). Probability and Measure. New York, Toronto, London: John Wiley and Sons.
  • Olav Kallenberg; Foundations of Modern Probability, 2nd ed. Springer Series in Statistics. (2002). 650 pp.
  • Henk Tijms (2004). Understanding Probability. Cambridge Univ. Press.
  • Olav Kallenberg; Probabilistic Symmetries and Invariance Principles. Springer -Verlag, New York (2005). 510 pp.
  • Gut, Allan (2005). Probability: A Graduate Course. Springer-Verlag.
  • ET. Jaynes. Probability Theory: The Logic of Science Cambridge University Press, (2003).
  • Stephen M. Stigler (1986) The history of statistics. Harvard University press
  • Richard T. Cox, Algebra of Probable Inference, The Johns Hopkins University Press
  • Bishop, CM., Pattern Recognition and Machine Learning. Springer, 2007
  • Stephen M. Stigler (1986) The history of statistics. Harvard University press
  • Wolpert, RL. (2004) A conversation with James O. Berger, Statistical science
  • Berger, James O (1985). Statistical Decision Theory and Bayesian Analysis. Springer Series in Statistics (Second ed.). Springer-Verlag
  • de Finetti, Bruno. "Probabilism: A Critical Essay on the Theory of Probability and on the Value of Science," (translation of 1931 article) in Erkenntnis, volume 31, September 1989
  • Hacking, I (1988) "Slightly More Realistic Personal Probability". 1967 article partly reprinted in: Gärdenfors, Peter and Sahlin, Nils-Eric. (1988) Decision, Probability, and Utility: Selected Readings. 1988. Cambridge University Press
  • Jaynes E.T. (2003) Probability Theory: The Logic of Science, CUP
  • Morgenstern, Oskar (1978). "Some Reflections on Utility". In Andrew Schotter. Selected Economic Writings of Oskar Morgenstern. New York University Press
  • Pfanzagl, J (1967). "Subjective Probability Derived from the Morgenstern-von Neumann Utility Theory". In Martin Shubik. Essays in Mathematical Economics In Honor of Oskar Morgenstern. Princeton University Press


Външни препратки