#!/usr/bin/octave-cli --persist
% Zufallszahlen ergeben in einer Autokorrelation eine Dirac-Funktion
# Das kann man zum Testen der Qualität verwenden.
% ----------------------------------
pkg load signal;
addpath("~/octave");
y = importdata("random-numbers");
N = length(y)
sy = (sum(y)/N)^2
#
Zb = xcorr(y,y,'biased');
Zu = xcorr(y,y,'unbiased');
Zc = xcorr(y,y,'coeff');
#
fig = figure(1);
plot(y,"og")
title("Zufallszahlen von random.org",'FontSize',12)
grid on
#
fig = figure(2);
plot(Zb)
title("Autokorrelation von Zufallsreihen (biased)",'FontSize',12);
grid on
#
fig = figure(3);
plot(Zu)
title("Autokorrelation von Zufallsreihen (unbiaased)",'FontSize',12);
grid on
#
fig = figure(4);
plot(Zc)
title("Autokorrelation von Zufallsreihen (coeff)",'FontSize',12);
grid on
#
dtn=importdata("random-dt");
Nd=length(dtn)
[hdtn1,xdt1]=hist(dtn,Nd/100);
fig = figure(5);
plot(xdt1,hdtn1,"or") #,xdt2,hdtn2*10,"+b",xdt3,hdtn3*100,".y")
title("Histogramm der Differenzanalyse der Zufallszahlen",'FontSize',12);
#xlim([1 N*1.05]);
grid on
fig = figure(6);
plot(dtn,"or")
title("DifferenzZeitreihenanalyse der Zufallszahlen",'FontSize',12);
#xlim([1 N*1.05]);
grid on

