Използването на генератор на псевдо-случайни числа различно ли е от реалното произволно число във веригата на Марков Монте-Карло?

Искам да използвам matlab, за да направя MCMC приближение до вероятностно разпределение. Открих обаче, че има много шум в сравнение с теоретичните резултати.

Въпросът ми е възможно ли е шумът да идва от генераторите на псевдослучайни числа в самия matlab?


person user2173168    schedule 09.11.2015    source източник
comment
Определено трябва да сте по-конкретни. Какъв шум? Можете ли да дадете вашия код и вашите теоретични разсъждения? Случайните генератори са изключително добре проучени и тествани, за да бъдат безпристрастни и че генерираните последователности наистина следват равномерно разпределение. Така че е трудно да се повярва, че във вашия случай случайните числа не са достатъчно случайни.   -  person cel    schedule 09.11.2015


Отговори (1)


Като общо правило не се постига близко съответствие между която и да е форма на симулация Монте Карло и теоретично разпределение на вероятностите, освен ако не се използват много голям брой проби. Обикновено човек би очаквал, че „шумът“, който наблюдавате, ще намалее толкова бавно, колкото единица върху корен квадратен от броя на пробите (т.е. може да са необходими 100 пъти повече проби, за да се намали шумът с коефициент около 10). При много специфични обстоятелства може да се открие, че прост генератор на произволни числа (напр. линеен конгруентен) може да не е подходящ за равномерно генериране на точки в многомерно пространство, но вероятно си струва да помислите за ефектите от вземането на проби от Монте-Карло, преди да се тревожите за качество на самия генератор на произволни числа. Очевидно е трудно да бъдем по-конкретни без конкретни подробности за вашето приложение.

person rwp    schedule 23.02.2018