随机数的产生
常用函数
A=rand(m,n); %生成m*n随机数矩阵,每个元素的选取服从[0,1]上的平均分布
B=randn(m,n); %生成m*n随机数矩阵,每个元素的选取服从标准正态分布
C=randi(nmax,m,n); %生成m*n随机数矩阵,每个元素是1和nmax之间的随机选取的正整数
D=randi([a,b],m,n); %m*n随机数矩阵,每个元素是[a,b]上的随机选取的正整数
一般随机数的产生
y=random('name',A,m,n) %单参数分布的随机数数组产生
y=random('name',A,B,m,n) %双参数分布的随机数数组产生
例1:生成4行5列的标准正态分布的随机数数组
A=random('Norm',0,1,4,5)
例2:生成1-10的4*5离散均匀分布数组
B=random('unid', 10, 4,5)
利用蒙特卡洛方法计算积分
考虑二重积分