Matlab随机数产生的问题.docVIP

Matlab随机数产生的问题2009-02-13 08:13这2天发现这个问题讨论较多，所以就搜索了一些资料。发现自己之前的理解有些很有问题，同时欢迎大家继续讨论。先澄清一下几个容易弄错的地方（也不一定全对）（1）用计算机产生的是“伪随机数”。用投色子计数的方法产生真正的随机数 , 但电脑若也这样做 , 将会占用大量内存 ; 用噪声发生器或放射性物质也可产生真正的随机数 , 但不可重复 . 而用数学方法产生最适合计算机 , 这就是周期有限 , 易重复的 ” 伪随机数 ”（2）随机数的产生需要有一个随机的种子，因为用计算机产生的随机数是通过递推的方法得来的，必须有一个初始值。（3）用同一台电脑，且在初始值和递推方法相同的情况下，可以产生相同的随机序列（由于以前每次使用randn或者rand得到都是不同值，所以曾经误以为相同的seed无法产生相同的序列）[size=2][color=red][b]一 matlab里产生随机数的方法[/b][/color][/size]matlab里和随机数有关的函数：（1） rand：产生均值为0.5、幅度在0~1之间的伪随机数（2） randn：产生均值为0、方差为1的高斯白噪声（3） randperm(n):产生1到n的均匀分布随机序列（4） normrnd(a,b,c,d)：产生均值为a、方差为b大小为cXd的随机矩阵还有很多的扩展函数，不再一一列出。不过他们都调用的是rand或者randn函数，由此可见在matlab里rand和randn是产生随机数的关键所在。看来只有看他们的源文件了function [varargout] = randn(varargin)%%%help 文档的内容略去%%%if nargout == 0?? builtin(randn, varargin{:});else?? [varargout{1:nargout}] = builtin(randn, varargin{:});end从这里也看不出到底是怎么产生的，就只看到builtin。而builtin函数的源文件是这样的：%BUILTIN?? Execute built-in function from overloaded method.% BUILTIN is used in methods that overload built-in functions to execute% the original built-in function. If F is a string containing the name% of a built-in function then BUILTIN(F,x1,...,xn) evaluates that% function at the given arguments.%% BUILTIN(...) is the same as FEVAL(...) except that it will call the% original built-in version of the function even if an overloaded one% exists (for this to work, you must never overload BUILTIN).%% [y1,..,yn] = BUILTIN(F,x1,...,xn) returns multiple output arguments.%% See also FEVAL.% Copyright 1984-2002 The MathWorks, Inc. % $Revision: 1.9 $?? $Date: 2002/04/15 04:16:04 $% Built-in function.后来发现matlab基本函数的源文件都是这么一个结构function [varargout] = functionname(varargin)%%%help 文档的内容%%%if nargout == 0?? builtin(functionname, varargin{:});else[varargout{1:nargout}] = builtin(functionname, varargin{:});end其中一些是.m文件一些是.bi文件，数学里最基本的运算就是加减乘除（对于计算机来说就只有加法器和乘法器）。那么最简单的sin函数如何用四则运算求值？反正我搞不清楚，sin.m内容也是一个builtin函数。扯了这么多，只是为了说明在使用matlab基本函数的时候，很多情况下看源文件并不能知道其中具体用的是什么数值方法。有了ra