Отправлено SM 06 июля 2005 г. 14:05 В ответ на: накапливаю - 32 бита с выхода беру и побитно в свой аккумулятор кладу. АКФ красивая, без всплесков на выборке из пары тысяч, вблизи нуля ( нормирую СЧ -0.5 < x < 0.5) отправлено DASM 06 июля 2005 г. 13:59

ты делаешь выборку определенных значений из последовательности, тем самым получив свое хитрое распределение, не попортив спектральные хар-ки. Ну или можно извратиться, сделав LFSR в поле Галуа, там будеш на каждом "такте" получать новое значение. Однако там умножать потрахацца...

• ай, сам ступил. Регистр то у меня 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, пустое)