Отправлено SM 06 июля 2005 г. 13:43 В ответ на: Получение "хорошего" псевдослучайного числа (+) отправлено DASM 06 июля 2005 г. 13:36

проверяется обычно рассчетом АКФ, которая должна иметь пик в нулевой точке и значение -1 в остальных (если нули рассматривать как минус единицы). Причем для битовой последовательности. LFSR никогда не давал хороших значений при таком применении как в данном случае.

P.S. А что значит "накапливаю" - складываешь что-ли что-то с чем-то?

А вообще ищи видимо другие варианты генераторов псевдослучайных чисел. их есть много. с другим распределением.

• накапливаю - 32 бита с выхода беру и побитно в свой аккумулятор кладу. АКФ красивая, без всплесков на выборке из пары тысяч, вблизи нуля ( нормирую СЧ -0.5 < x < 0.5) — DASM   (06.07.2005 13:59 217.118.78.1, пустое)
  • Так сделал бы LFSR сразу нужной тебе длины и брал бы прямо из него. А так (+) — SM   (06.07.2005 14:05 213.141.159.26, 280 байт)
    • ай, сам ступил. Регистр то у меня 39 бит, а я с него выдвигаю только 32. Оставшиеся 7 явно дадут вполне предсказуемое влияние на след выборку. Насчет LFSR нужной длины - не понял. Тогда соседние выборки будут отличатся одним сдвигом и старшим битом - тут уж точно зависимость будет жуткая — DASM   (06.07.2005 14:09 217.118.78.1, пустое)
      • да уж...я тоже не понял, ты так уверенно объяснял что делаешь что взять 32 из 39 это уже само собой разумеющееся...убаюкал:)) — ы   (06.07.2005 14:14 80.92.98.211, пустое)
      • Старший бит это очень много :) Соседние выборки могут отличаться ой-ой на сколько. — SM   (06.07.2005 14:13 213.141.159.26, пустое)
        • ты предлагаешь делать снимок регистра на каждом сдвиге ? Ну совсем неслучайно выйдет. На двумерном графике - две наклонные линии выходят, хотя автокореляцция прекрасно около нуля сидит — DASM   (06.07.2005 14:18 217.118.78.1, пустое)
          • Не упускай, что я говорил про поле Галуа, там выход регистра не один бит, а сразу сколько надо. — SM   (06.07.2005 14:19 213.141.159.26, пустое)
            • Вот сейчас беру снимок этого 39 битного регистра, но не после каждого сдвига, а после 39-ти. То есть как бы совсем новый. И все равно гадость такая. С Галуа еще разбираться ночь надо, ибо не в курсах — DASM   (06.07.2005 14:32 217.118.78.1, пустое, картинка)
              • Ну правильно - ты берешь не все значения, часть выкидываешь, отсюда и белые полосы. В идеале все должно почернеть. — SM   (06.07.2005 14:34 213.141.159.26, пустое)
                • ну так если брать все значения то вообще две линии выходят — DASM   (06.07.2005 14:35 217.118.78.1, пустое)
                  • Так я наверное не понял... Ты как делаешь? Я бы так проверял (+) — SM   (06.07.2005 14:38 213.141.159.26, 70 байт)
                    • не, x = XMAX * (myrand() / RAND_MAX) . В double все ессно конвертировано — DASM   (06.07.2005 14:40 217.118.78.1, пустое)
                      • А Y? — SM   (06.07.2005 14:43 213.141.159.26, пустое)
                        • аналогично. По сути рисуется зависимость следущего от предыдущего — DASM   (06.07.2005 14:48 217.118.78.1, пустое)
                          • вот посмотрел сурцы rand из NET. Это какой класс алгоритмов ? (+) — DASM   (06.07.2005 14:58 217.118.78.1, 349 байт)
                            • Это линейный конгруэнтный генератор. — SM   (06.07.2005 15:07 213.141.159.26, пустое)
                              • а что у него с периодом ? Распределение вроде на первой взгляд ничего — DASM   (06.07.2005 15:10 217.118.78.1, пустое)
                                • Максимум это (+) — SM   (06.07.2005 15:29 213.141.159.26, 196 байт)
                                  • ок, всем спасибо и тебе в частности — DASM   (06.07.2005 15:32 217.118.78.1, пустое)
                  • и вообще, для "хорошей" псевдослучайности последовательности не все ли равно, подряд брать или с прорежевинием ? — DASM   (06.07.2005 14:36 217.118.78.1, пустое)
                    • Конечно не все равно. При прореживании часть чисел просто никогда не появится на выходе прореживателя. Соотв. изменится распределение. — SM   (06.07.2005 14:43 213.141.159.26, пустое)
• Причем АКФ при накапливании расползается и получаются узоры. — argus98   (06.07.2005 13:48 81.22.205.230, пустое)
  • Ну эт естессно, накопление это есть фильтрация... и спектр становится совсем другим. — SM   (06.07.2005 13:51 213.141.159.26, пустое)
    • да не накапливаю, просто беру выборку 32 бита, получаю случайное число, где все биты - случайны — DASM   (06.07.2005 14:00 217.118.78.1, пустое)