求解答

注意第二个参数应该是BETA的倒数（已知X的分布，1/X的分布可以自己推一下，任何一半概率的书上都有）
1/X~IG(ALPHA,1/BETA)

IG是BAYESIAN常用的一个分布，而且如果你用MATLAB，知道这个性质会很有帮助滴，嘿嘿……

非常感谢您的帮助！！
我按照您说的修改了程序，结果总算不一直是0了。我要用自助法(Bootstrap法)随机抽样2000次，2次的耗时是1.2秒，预计20分钟后结果就出来了。
我提这个问题是为了生成一个【服从逆伽马分布的随机数】，这个随机数是自助法中要用到的一个参数，但matlab中只提供了服从gamma分布的随机数生成函数：gamrnd。参数alpha和belta是已知的，我现在的做法是：
>>x=gamrnd(alpha,1/belta);
>>x=1/x;

不对。（1/X）不服从伽马分布

这是密度函数，注意，如果写成关于（1/X）的分布的话，那么中间那一项（1/X）的指数应该是alpha+1。
而如果1/X不服从伽马分布，中间应该是alpha-1的。（也就是说，中间那一项会是X^(-alpha+1)）

还想请教个问题~~就是逆伽玛分布和贝塔分布都是指数型分布族么？？