реклама на сайте
подробности

 
 
 
Reply to this topicStart new topic
> псевдослуч. послед-ть огромной длительности, посоветуйте
shf_05
сообщение Nov 25 2008, 11:48
Сообщение #1


Профессионал
*****

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



алгоритм формирования ПСП такой:
//------------
//algoritm: iran= 1664525L*iran+1013904223L;
//для генерации числа в диапазоне 0 - 2^32-1 достаточно простого умножения
//на мультипликатор и сложения с инкрементом.
//Деление по модулю будет произведено автоматически при переполнении.
//значения мультипликатора и инкремента для этого случая получены в исследованиях
//D. Knuth и H.W. Lewis.
//значения (0...MAX/2) и (MAX/2...MAX) равновероятны и можно интерпретировать число как знаковое в доп коде
//как -MAX/2...MAX/2
//------------

для 16 разр. проца беру старшие 16бит
//generate random signal from -32768 to 32767
//OUTPUT in sr0

интересующие меня вопросы:
1, последовательность повторяется через 2^32 или 2^64?
2, более важный вопрос:
для запуска генератора применены следующие НУ
#define RAND0_HI 0xa4ce //---inital conditions //HI
#define RAND0_LO 0x39a7 //LO
скажем досчитав до 2^32 или 2^64 поменяю нач. уловия (например подставлю значение входной переменной из АЦП оно случайно), получу ли новый набор чисел?
имхо- да, но есть люди, говорящие что это чушь и я запущу ту же посл-ть лишь с другого "места".
Go to the top of the page
 
+Quote Post
Artem_Petrik
сообщение Nov 25 2008, 11:56
Сообщение #2


Местный
***

Группа: Свой
Сообщений: 443
Регистрация: 22-07-06
Из: Украина, г. Харьков
Пользователь №: 19 006



Цитата(shf_05 @ Nov 25 2008, 15:48) *
интересующие меня вопросы:
1, последовательность повторяется через 2^32 или 2^64?
2, более важный вопрос:
для запуска генератора применены следующие НУ
#define RAND0_HI 0xa4ce //---inital conditions //HI
#define RAND0_LO 0x39a7 //LO
скажем досчитав до 2^32 или 2^64 поменяю нач. уловия (например подставлю значение входной переменной из АЦП оно случайно), получу ли новый набор чисел?
имхо- да, но есть люди, говорящие что это чушь и я запущу ту же посл-ть лишь с другого "места".


1. Если оперируете 32-х разрядными числами, период повторения больше чем 2^32 быть не может в принципе.
2. Люди полностью правы.

Чтоб получить другую последовательность, нужно менять те константы, которые у вас 1664525L и
1013904223L. Но нельзя брать "от фонаря", потому как только немногие комбинации дадут максимальный период повторения.
Go to the top of the page
 
+Quote Post
shf_05
сообщение Nov 25 2008, 12:19
Сообщение #3


Профессионал
*****

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



1. ага 2^32
2. по уму то вроде правы, я сам так думал, но смутило, что, смутно помню- моделили такой алгоритм, для разных НУ получались разные последовательности, были даже случаи, когда алгоритм давал "плохой" шум... надо повспоминать.

Сообщение отредактировал shf_05 - Nov 25 2008, 12:22
Go to the top of the page
 
+Quote Post
shf_05
сообщение Nov 25 2008, 14:22
Сообщение #4


Профессионал
*****

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



все верно, там изменялись коэффициенты, которые и были этими НУ, изменение "стартового числа" ни к чему хорошему не приводит

Сообщение отредактировал shf_05 - Nov 25 2008, 14:24
Go to the top of the page
 
+Quote Post
Михаил_K
сообщение Nov 25 2008, 14:51
Сообщение #5


Знающий
****

Группа: Свой
Сообщений: 552
Регистрация: 29-02-08
Пользователь №: 35 481



Цитата(shf_05 @ Nov 25 2008, 14:48) *
1, последовательность повторяется через 2^32 или 2^64?

через 2^32-1;
Go to the top of the page
 
+Quote Post
shf_05
сообщение Nov 25 2008, 19:23
Сообщение #6


Профессионал
*****

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



Цитата(Михаил_K @ Nov 25 2008, 19:51) *
через 2^32-1;

простите за неточность, под 2^32 или 2^64 я и имел ввиду 32-1.(63)
вобщем почесав репу, отдохнув немного во всем разобрался, тема закрыта, всем спасибо.
Go to the top of the page
 
+Quote Post
shf_05
сообщение Nov 26 2008, 05:19
Сообщение #7


Профессионал
*****

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



однако...

http://www.xakep.ru/post/25259/default.asp?print=true

"Получаемая последовательность зависит от выбора стартового числа X0 и при разных его значениях получаются различные последовательности случайных чисел. В то же время, многие свойства последовательности Xk определяются выбором коэффициентов в формуле и не зависят от выбора стартового числа. Ясно, что последовательность чисел, генерируемая таким алгоритмом, периодична с периодом, не превышающим m. "

они неправы или неправильно поняты? есть специалисты в этой теме?

PS чтобы не заморачиваться использовал ал-м Marsaglia-Multicarry.
Go to the top of the page
 
+Quote Post
Guest_TSerg_*
сообщение Nov 26 2008, 08:01
Сообщение #8





Guests






Цитата(shf_05 @ Nov 26 2008, 08:19) *
однако...

http://www.xakep.ru/post/25259/default.asp?print=true

"Получаемая последовательность зависит от выбора стартового числа X0 и при разных его значениях получаются различные последовательности случайных чисел. В то же время, многие свойства последовательности Xk определяются выбором коэффициентов в формуле и не зависят от выбора стартового числа. Ясно, что последовательность чисел, генерируемая таким алгоритмом, периодична с периодом, не превышающим m. "

они неправы или неправильно поняты? есть специалисты в этой теме?

PS чтобы не заморачиваться использовал ал-м Marsaglia-Multicarry.


"Да что Вы всякую щваль читаете" (С)

Хотя бы:
http://www.agner.org/
http://crypto.mat.sbg.ac.at/results/karl/server/node1.html
Go to the top of the page
 
+Quote Post
shf_05
сообщение Nov 26 2008, 08:18
Сообщение #9


Профессионал
*****

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



виноват, что нагуглил
пробовал зайти на [url=http://www.agner.org/]- не пустил пров,- сказал- сайт забанен- порно((, с др. места попробую

имхо: эта "шваль" передирает все с тех же сайтов, тока переводит- наверное неточно)

огромное спасибо за ссылки!

Сообщение отредактировал shf_05 - Nov 26 2008, 08:29
Go to the top of the page
 
+Quote Post
alex_os
сообщение Nov 26 2008, 08:38
Сообщение #10


Знающий
****

Группа: Свой
Сообщений: 521
Регистрация: 12-05-06
Пользователь №: 17 030



Цитата(shf_05 @ Nov 26 2008, 11:18) *
виноват, что нагуглил
пробовал зайти на [url=http://www.agner.org/]- не пустил пров,- сказал- сайт забанен- порно((, с др. места попробую

имхо: эта "шваль" передирает все с тех же сайтов, тока переводит- наверное неточно)

огромное спасибо за ссылки!


Да все правильно почему последовательность будет зависеть от стартового числа? Есть стартовое число (т.е. состояние генератора ПСП ) М, при любом стартовом состоянии X, мы придем к состоянию М за число итераций <=2^32-1. Т.е начальное состояние задает какбы "фазу" этой ПСП, видимо это имелось ввиду.


--------------------
ну не художники мы...
Go to the top of the page
 
+Quote Post

Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 18th June 2025 - 14:10
Рейтинг@Mail.ru


Страница сгенерированна за 0.01408 секунд с 7
ELECTRONIX ©2004-2016