matlab中随机数的问题

来源:百度文库 编辑:神马文学网 时间:2024/04/27 17:06:32
norm(a)  一个向量的长度,也就是二阶范数
randn(2,20)  生成均值为0,方差为1的正态分布
scatter(a(1,:),a(2,:))  画散点图
关键问题不是用法而是意义,rand是0-1的均匀分布,randn是均值为0方差为1的正态分布
用法自己help rand或者help randn一下就知道了
rand(n)或randn(n)生成n*n的随机数矩阵
rand(n,m)或randn(m,n)生成m*n的随机数矩阵
取模(mod)与取余(rem)的区别——Matlab学习笔记【09-11-15】
本帖来自:数学中国 作者: 木长春 日期: 2009-11-15 19:51 您是本帖第664个浏览者
昨天在学习Matlab的数学函数时,教程中提到取模(mod)与取余(rem)是不同的,今天在网上具体查了一下:
通常取模运算也叫取余运算,它们返回结果都是余数.rem和mod唯一的区别在于:
当x和y的正负号一样的时候,两个函数结果是等同的;当x和y的符号不同时,rem函数结果的符号和x的一样,而mod和y一样。
这是由于这两个函数的生成机制不同,rem函数采用fix函数,而mod函数采用了floor函数(这两个函数是用来取整的,fix函数向0方向舍入,floor函数向无穷小方向舍入)。
rem(x,y)命令返回的是x-n.*y,如果y不等于0,其中的n = fix(x./y),而mod(x,y)返回的是x-n.*y,当y不等于0时,n=floor(x./y)
两个异号整数取模取值规律            (当是小数时也是这个运算规律,这一点好像与C语言的不太一样)
先将两个整数看作是正数,再作除法运算
①能整除时,其值为0
②不能整除时,其值=除数×(整商+1)-被除数
例:mod(36,-10)=-4
即:36除以10的整数商为3,加1后为4;其与除数之积为40;再与被数之差为(40-36=4);取除数的符号。所以值为-4。
例:mod(9,1.2)=0.6
例:
>> mod(5,2)
ans =1                   %“除数”是正,“余数”就是正
>> mod(-5,2)
ans =1
>> mod(5,-2)
ans =-1                  %“除数”是负,“余数‘就是负
>> mod(-5,-2)
ans =-1                  %用rem时,不管“除数”是正是负,“余数”的符号与“被除数”的符号相同
>> rem(5,2)
ans =1                   %“被除数”是正,“余数”就是正
>> rem(5,-2)
ans =1
>> rem(-5,2)
ans =-1                 %“被除数”是负,“余数”就是负
>> rem(-5,-2)
ans =-1
慢慢体会,两者确实不一样
问一下各位,我用mvnrnd函数生成多维正态分布函数,
R = MVNRND(MU,SIGMA,CASES),其中,系统提示SIGMA must be a positive semi-definite matrix.而我的程序中途,出现了SIGMA不是正定矩阵的情况,导致无法运行,请问各位,该怎么解决呢?谢谢!
那就把SIGMA弄成正定矩阵就行了,例如
mu = [2 3];
SIGMA = [1 1.5; 1.5 3];
r = mvnrnd(mu,SIGMA,100);
plot(r(:,1),r(:,2),'+')
为什么要转化呢?可以加一个判断SIGMA是否为正定矩阵的语句,若SIGMA是正定矩阵,就生成多维正态分布随机数,若不是,就不执行。
判断SIGMA是否为正定矩阵可以用chol函数。
因为不正定的话,密度函数不存在,此时是一个退化的分布,你可以参考一下多元统计分析。