求助一个程序怎么编

假设有两个随机变量X和N，其中N服从N（0，σ），X具体不知道。然后有Y=X+N，Y可以进行观察，有10000个样本。首先通过样本来估计出Y的密度函数，再通过X=Y-N来计算X的密度函数，再使用贝叶斯法来计算使得P（X|Y）最大的X。

也就是求\[\widehat{X}=\arg\max exp(-\frac{(Y-X)^2}{2\sigma^2}) \int_{-\infty}^{+\infty} f(X+N)*exp(-\frac{N^2}{2\sigma^2})dN\]



还有就是怎么用R读取.mat数据集。


求助。。。。先谢谢了

没人吗？是不是实现不了

这个如果估计出来y的密度函数的话，公式得到是一个向量。直接对公式结果求which.max，相应的y-N不就是你想要的X么