假设已知总体密度是 p(x1,x2)=0.1*N(mu1,sigma1)+0.7N(mu2,sigma2)+0.2N(mu3,sigma3), 其中N均是二维Gaussian密度函数。问题是:如何产生服从p(x1,x2)的随机数。
我的思路是:
1.产生多项分布随机数1000个,n=100; p=[0.1 0.7 0.2]; m=mnrnd(n,p,1000);
2。如果模拟结果m=20 32 48, 48为最大数,则由N(mu3,sigma3)产生随机数;
如果模拟结果m=50 32 18, 50为最大数,则由N(mu1,sigma1)产生随机数;
依次产生其他随机数。
请问:这个思路有问题吗?