c8_normal_01.m

来自「用于生成一个正态分布的伪随机数数值序列」· M 代码 · 共 32 行

function [ z, seed ] = c8_normal_01 ( seed )

%% C8_NORMAL_01 returns a unit pseudonormal C8.
%
%  Discussion:
%
%    The value has mean 0 and standard deviation 1.
%
%  Modified:
%
%    17 July 2006
%
%  Author:
%
%    John Burkardt
%
%  Parameters:
%
%    Input, integer SEED, a seed for the random number generator.
%
%    Output, complex Z, a sample of the PDF.
%
%    Output, integer SEED, a seed for the random number generator.
%
  [ v1, seed ] = r8_uniform_01 ( seed );
  [ v2, seed ] = r8_uniform_01 ( seed );

  x_r = sqrt ( - 2.0 * log ( v1 ) ) * cos ( 2.0 * pi * v2 );
  x_c = sqrt ( - 2.0 * log ( v1 ) ) * sin ( 2.0 * pi * v2 );

  z = x_r + x_c * i;