symulacje

Idź do spisu treści

Menu główne

3.2

Kody rozszerzone > Rozdział 3

##########################################################################################
###                                                                                      
Generator liniowy                                                                                     ###
##########################################################################################


### Generowanie 35 losowych wartości i prezentacja wykresu widmowego

x=7; a=11; c=0; m=31
### parametry początkowe
y=c(x/m,(a*x+c) %% m)
### znak %% oznacza modulo (reszta z dzielenia)
plot(x[1],x[2],xlim=c(0,1),ylim=c(0,1))
for (i in 1:35)
{

x=a*x+c

  #%% - operator modulo

x=x %% m
y=c(y,x/m)
Sys.sleep(.5)
### spowolnienie wykonania programu, możemy obserwować kolejne punkty na wykresie
points(y[i],y[i+1])

}

### Generowanie 10000 losowych wartości i prezentacja wykresu widmowego
x=
111; a=71; c=7; m=1371 ### parametry początkowe
y=c(x/m,(a*x+c) %% m)
### znak %% oznacza modulo (reszta z dzielenia)
plot(x[1],x[2],xlim=c(0,1),ylim=c(0,1))
for (i in 1:35)
{
    x=a*x+c
    x=x %% m
    y=c(y,x/m)
    points(y[i],y[i+1])
}

 
Wróć do spisu treści | Wróć do menu głównego