资源描述：
《matlab随机数产生的问题》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、Matlab随机数产生的问题2009-02-1308:13这2天发现这个问题讨论较多，所以就搜索了一些资料。发现自己之前的理解有些很有问题，同时欢迎大家继续讨论。先澄清一下几个容易弄错的地方（也不一定全对）（1）用计算机产生的是“伪随机数”。用投色子计数的方法产生真正的随机数,但电脑若也这样做,将会占用大量内存;用噪声发生器或放射性物质也可产生真正的随机数,但不可重复.而用数学方法产生最适合计算机,这就是周期有限,易重复的”伪随机数”（2）随机数的产生需要有一个随机的种子，因为用计算机产生的随机数是通过递推的方法得来的，必须有一个初始值。（3）用同一台电脑，且在初始值和
2、递推方法相同的情况下，可以产生相同的随机序列（由于以前每次使用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的随机矩阵还有很多的扩展函数，不再一
3、一列出。不过他们都调用的是rand或者randn函数，由此可见在matlab里rand和randn是产生随机数的关键所在。看来只有看他们的源文件了function[varargout]=randn(varargin)%%%help文档的内容略去%%%ifnargout==0 builtin('randn',varargin{:});else [varargout{1:nargout}]=builtin('randn',varargin{:});end从这里也看不出到底是怎么产生的，就只看到builtin。而builtin函数的源文件是这样的：%BUILTIN Ex
4、ecutebuilt-infunctionfromoverloadedmethod.%BUILTINisusedinmethodsthatoverloadbuilt-infunctionstoexecute%theoriginalbuilt-infunction.IfFisastringcontainingthename%ofabuilt-infunctionthenBUILTIN(F,x1,...,xn)evaluatesthat%functionatthegivenarguments.%%BUILTIN(...)isthesameasFEVAL(...)except
5、thatitwillcallthe%originalbuilt-inversionofthefunctionevenifanoverloadedone%exists(forthistowork,youmustneveroverloadBUILTIN).%%[y1,..,yn]=BUILTIN(F,x1,...,xn)returnsmultipleoutputarguments.%%SeealsoFEVAL.%Copyright1984-2002TheMathWorks,Inc.%$Revision:1.9$ $Date:2002/04/1504:16:04$%Buil
6、t-infunction.后来发现matlab基本函数的源文件都是这么一个结构function[varargout]=functionname(varargin)%%%help文档的内容%%%ifnargout==0 builtin('functionname',varargin{:});else[varargout{1:nargout}]=builtin('functionname',varargin{:});end其中一些是.m文件一些是.bi文件，数学里最基本的运算就是加减乘除（对于计算机来说就只有加法器和乘法器）。那么最简单的sin函数如何用四则运算求值？反正
7、我搞不清楚，sin.m内容也是一个builtin函数。扯了这么多，只是为了说明在使用matlab基本函数的时候，很多情况下看源文件并不能知道其中具体用的是什么数值方法。有了rand和randn就可以产生轻松产生均匀分布和正态分布的随机数了（1）产生在[a,b]区间服从均匀分布随机序列的方法 (b-a)*rand(m,n)+a>>3*rand(2)+2ans=2.81662.04582.59644.2404（2）产生服从正态分布的随机数>>randn('state',2)>>a=normrnd(0,1,1,6)a=1.74910.132
