秦时明星 发表于 2012-5-18 23:18

回复 7 # xiezhh 的帖子

谢老师,您好!看了您的《MATLAB统计分析与应用39个案例分析》后,我在编程中碰到一个问题:
在第四章的生成任意一元分布随机数的调用命令中,在comand命令中输入如下命令:
k=-0.345119445;
a=11.67962773;
e=68.61830219;
y='(-(log(1-('num2str(k)')*(x-'num2str(e)')/('num2str(a)')))/('num2str(k)'))';
pdffun = ['1/('num2str(a)')/exp((1-('num2str(k)'))*y)*(1+exp(-y))^2'];
x = crnd(pdffun, , 100, 1);
为什么输入后,一直显示 y='........'; 这个表达式有错呢?
希望谢老师您能帮我解答看看怎么解决这个问题呢!
先谢谢 谢老师您了!

秦时明星 发表于 2012-5-18 23:23

回复 13 # xiezhh 的帖子

刚刚的那个程序中,密度函数是:
f(x)=1/(a*exp((1-k)*y)*(1+exp(-y))^2)
其中,y=-1/k*ln(1-k*(x-e)/a)

秦时明星 发表于 2012-5-19 08:40

回复 7 # xiezhh 的帖子

谢老师,您好!我注意到,' num2str(k) ' 在引号的左右两侧是有空格的,修改后,产生的数怎么会是 复数呢?结果居然为 54.567865+0.5655785 i 之类的数!我表示很不解!
希望谢老师能帮忙解答一下!谢谢了!

fmlv2007 发表于 2012-11-4 17:40

您好谢老师:
    我想问一下如果想生成服从二元函数的随机数,能否用您给出的函数,最简单的如联合分布为二元正态分布。
页: 1 [2]
查看完整版本: 已知概率密度函数的随机数?