如何找到一个代码下面的算法

2视图(30天)
可能
可能 2013年9月20日
我想找到z随机变量的分布假设 X1, X2,…, Xn n 相互独立的随机变量,让Z是它们的和:
Z = X1 + X2 +……+ Xn
的分布 Z 可以递归地使用两个随机变量的结果总结如下功能:
函数[答案]= Sum_Of_2_Random_Variables (z)
F = @ (y) normpdf (z-y)。* normpdf (y);
回答=积分(F,负无穷,正);
结束
该算法是:
算法一:
  1. 首先,定义:日元= X1 + X2和计算的分布日元;
  2. 然后,定义:Y2 = y₁+ X3和计算的分布Y2;
  3. 等等,直到Z的分布计算:Z = Yn = Y (n - 1) + Xn
现在我不知道如何写的代码算法a .任何帮助将不胜感激。谢谢你!
主要的问题我有如下:假设我想找的分布 日元= X1 + X2 使用下面的代码:
F = @ (y) Sum_Of_2_Random_Variables (z-y)。* normpdf (y);
回答=积分(F,负无穷,正);
这不会工作,因为输入的函数 Sum_Of_2_Random_Variables 不是双型。
  • 我知道正常随机变量之和是正常随机变量,随机变量的分布的可以是任何分布,来简化这个问题我选择了标准正态分布。(在我的问题分布是不正常)。
1评论
可能
可能 2013年9月20日
我只是想知道如何调用一个函数包括积分递归。

登录置评。

答案(1)

马特·J
马特·J 2013年9月21日
编辑:马特·J 2013年9月21日
看起来像一个不必要的痛苦和暴力的方式来做这件事。而不是使用n次旋转,似乎最好使用特征函数。即。,以fft算法的pdf (Xi)的特色功能,它们相乘,然后进行傅里叶反变换的结果。
5个评论
马特·J
马特·J 2013年9月22日
编辑:马特·J 2013年9月22日
没有明确的理由这样做你需要象征性地而不是数值。反正你已经尝试数值积分,用积分()命令。
样品你的pdf(或您的特点功能如果你知道)和使用fft()和传输线()而不是试图向前做分析,傅里叶变换。

登录置评。

s manbetx 845

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!