赞
踩
问题:今天在研究中需要以一定的概率产生随机数,要生成一个矩阵c,行标和列标一致时为零,其他情况以概率为p生成[1,10]的正整数,以概率为(1-p)生成0。
找到了一个很好的函数:
randsrc(100,1,[alphabet; prob])
alphabet表示需要产生的数字,prob表示对应的概率,注意:两者的维数要一致。
比如:要产生 1、4、 6这三个数。它们分别出现的概率为 0.1 、0.4、 0.5 如何设计程序使得按照这个概率产生100个随机数呢?
alphabet = [1 4 6]; prob = [0.1 0.4 0.5];
randsrc(100,1,[alphabet; prob])
对于我的问题,我的想法是当角标一样时,赋值为0,当不一样时,生成一个[1,10]的正整数a,与0组合成一个新的矩阵b,此时以一定的概率添加到c中。
当c的大小为5*5时,程序为:
clc
clear all
n=5;
p=30;
prob=[0.3,1-p/100];
c=[];
for i=1:n
for j=1:n
if i==j
c(i,j)=0;
else
a=round(rand(1,1)*9+1);
b=[a 0];
c(i,j)=randsrc(1,1,[b;prob]);
end
end
end
问题得到解决。
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。