rndLCvm¶
Purpose¶
Computes von Mises pseudo-random numbers.
Note
This function is deprecated but remains for backward compatibility.
Format¶
-
{ x, newstate } =
rndLCvm
(r, c, m, k, state)¶ - Parameters:
r (scalar) – number of rows of resulting matrix.
c (scalar) – number of columns of resulting matrix.
m (matrix, vector or scalar) – means for von Mises distribution, scalar or ExE conformable matrix with r and c.
k (matrix, vector or scalar) – shape argument von Mises distribution, scalar or ExE conformable matrix with r and c.
state (scalar or vector) –
scalar case
state = starting seed value only. System default values are used for the additive and multiplicative constants.
The defaults are 1013904223, and 1664525, respectively. These may be changed with
rndcon
andrndmult
.If state = -1, GAUSS computes the starting seed based on the system clock.
3x1 vector case
[1]
the starting seed, uses the system clock if -1
[2]
the multiplicative constant
[3]
the additive constant
4x1 vector case
state = the state vector returned from a previous call to one of the
rndLC
random number generators.
- Returns:
x (RxC matrix) – von Mises distributed random numbers.
newstate (4x1 vector) –
[1]
the updated seed
[2]
the multiplicative constant
[3]
the additive constant
[4]
the original initialization seed
Remarks¶
r and c will be truncated to integers if necessary.
Technical Notes¶
This function uses a linear congruential method, discussed in Kennedy, W.J. Jr., and J.E. Gentle, Statistical Computing, Marcel Dekker, Inc. 1980, pp. 136-147. Each seed is generated from the preceding seed using the formula
where %
is the mod operator and where a is the multiplicative constant
and c is the additive constant.
Source¶
randlc.src