r8mat_normal.m

用于生成一个正态分布的伪随机数数值序列

function [ r, seed ] = r8mat_normal ( m, n, a, b, seed )

%% R8MAT_NORMAL returns a scaled pseudonormal R8MAT.
%
%  Modified:
%
%    17 July 2006
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Paul Bratley, Bennett Fox, Linus Schrage,
%    A Guide to Simulation,
%    Springer Verlag, pages 201-202, 1983.
%
%    Bennett Fox,
%    Algorithm 647:
%    Implementation and Relative Efficiency of Quasirandom
%    Sequence Generators,
%    ACM Transactions on Mathematical Software,
%    Volume 12, Number 4, pages 362-376, 1986.
%
%    Peter Lewis, Allen Goodman, James Miller,
%    A Pseudo-Random Number Generator for the System/360,
%    IBM Systems Journal,
%    Volume 8, pages 136-143, 1969.
%
%  Parameters:
%
%    Input, integer M, N, the number of rows and columns in the array.
%
%    Input, real A, B, the mean and standard deviation.
%
%    Input/output, integer SEED, the "seed" value, which should NOT be 0.
%    On output, SEED has been updated.
%
%    Output, real R(M,N), the array of pseudonormal values.
%
  for j = 1 : n

    [ v, seed ] = r8vec_normal ( m, a, b, seed );
    r(1:m,j), seed ] = v(1:m)';

  end