Generacia nahodnych cisel

[Marek Peca]

> Zda existuje rovnomerny barevny sum ovsem nejsem si jist, cekam
> na pouceni.

Bohuzel nemam na to, abych se pokusil na to prijit nejak exaktne,
nicmene mi to nedalo a ucinil jsem pokus na zaklade selskeho rozumu.
Prohnal jsem barevny gaussovsky sum fci erf a podival se na
histogram, skutecne slo prakticky o rovnomerne rozlozeny sum, rovnez
spektrum bylo _podobne_ (ne stejne) barevne, jako v zdrojovem sumu.

Takze i rovnomerne rozlozeny barevny sum patrne existuje.

Pokud to ale cte nekdo, kdo tomu opravdu rozumi, prosil bych o
potvrzeni ci vyvraceni zaveru, neb vim, ze takovehle zjednodusene
uvahy nekdy vedou na zcesti.

% --- barevsum.m
% pokus o vytvoreni barevneho rovnomerne rozlozeneho sumu

% gaussovsky bily
u1=randn(10000,1);
[b,a]=butter(2,0.2);
% gaussovsky barevny
u2=filter(b,a,u1);
u2=u2/std(u2);

figure; hist([u1 u2],100);

% rovnomerny bily, barevny
y1=erf(u1/sqrt(2));
y2=erf(u2/sqrt(2));

figure; hist([y1 y2],100);

% overeni spekter
f=linspace(0,1,500)';
h=spectrum.welch;
Hu1=psd(h,u1);
Hu2=psd(h,u2);
Hy1=psd(h,y1);
Hy2=psd(h,y2);

figure; hold on;
plot(Hu1.Frequencies,Hu1.Data,'b');
plot(Hu2.Frequencies,Hu2.Data,'g');
plot(Hy1.Frequencies,Hy1.Data,'c');
plot(Hy2.Frequencies,Hy2.Data,'r');

% -- KOHEII, --
%
% pekny den preje Marek P.