如何在R中找到函数之间的交点路径?

时间:2014-10-13 22:18:27

标签: r optimization intersection gaussian

我试图弄清楚如何在一个多变量函数等于另一个的二维场中找到点。更具体地说,高斯人:

g1 <- function(x){x1 = x[1]; x2 = x[2]; return (1.2 / (sqrt(4 * pi * pi * 81)) * exp(-1/18*((x1-12)**2 + (x2-3))) + .4 / (sqrt(4 * pi * pi * 81)) * exp(-1/18*((x1-3)**2 + (x2-5)**2)) )}

g2 <- function(x){x1 = x[1]; x2 = x[2]; return (.8 / (sqrt(4 * pi * pi * 81)) * exp(-1/18*((x1-4)**2 + (x2-12))) + .2 / (sqrt(4 * pi * pi * 81)) * exp(-1/18*((x1-3)**2 + (x2-5)**2)) )}

g3 <- function(x){x1 = x[1]; x2 = x[2]; return (1.2 / (sqrt(4 * pi * pi * 81)) * exp(-1/18*((x1-4)**2 + (x2-12))) + .4 / (sqrt(4 * pi * pi * 81)) * exp(-1/18*((x1-12)**2 + (x2-3)**2)) )}

在这里,我已经将矩阵Sigma分解为简单代数。我想找到代表g1(x)== g2(x),g1(x)== g3(x)和g2(x)== g3(x)的曲线。我是否需要简单地进行双循环并使用一些阈值(例如,abs(g1-g2)&lt; 1E-5)?我喜欢更健壮的东西。

提前致谢!

0 个答案:

没有答案