我正在编写一个MATLAB代码来快速解决以下问题:设X是一个根据P(x)分布的随机变量,取X的两个独立副本,称之为X1和X2,找到Y = f的分布(X1, X2)其中f(,)是已知函数.
为了解决上述问题,我从两个向量x和p开始,使得p(i)= P(x(i)).假设它们都包含n个元素.我可以容易地计算n×n矩阵y,使得y(i,j)= f(x(i),x(j)).此外,我可以计算n×n矩阵p_out,使得p_out(i,j)= p(i)* p(j).这意味着P(Y = y(i,j))= p(i,j).
现在,如果y的所有元素都是不同的,我们差不多完成了.它仍然只是将矩阵转换为向量,并且可能对它们进行排序以获得良好的输出.假设我们也通过设置来做到这一点
y = y(:);
p_out = p_out(:);
[y, idx] = sort(y);
p_out = p_out(idx);
然而,问题是y的元素通常不是唯一的.因此,我必须合并y的相同元素如下:如果y(i)= y(j)(记住现在y被转换为向量)则删除y(j)并设置p(i)= p( i)p(j).执行此操作的一种肮脏方式是使用for循环(因为y现在已经排序,我们只需要将每个元素与其后续元素进行比较).但是,我想知道是否存在更好的方式.
我知道唯一会删除矢量的重复元素(因此,如果我们只需要y就足够了).我也知道它返回两个索引向量,以某种方式指示重复元素的位置.但是,我想不出任何好的方法来使用它的输出来适当地合并p的元素.