1

我有这个：

    install.packages("mnormt")
    library(mnormt)


    ma = c(1,5)
    Sa = matrix(c(2,0,0,0.5),2,2)
    A = rmnorm(100, ma, Sa)


    mb = c(5,7)
    Sb = matrix(c(1,0.5,0.5,0.5),2,2)
    B = rmnorm(200, mb, Sb)


    plot(rbind(A,B), type = "n")
    points(A, col = "red")
    points(B, col = "blue")


    points(ma[1], ma[2], pch = 20)
    points(mb[1], mb[2], pch = 20)


    lines(ma, mb)
    centroreta = c((ma[1] + mb[1])/2,(ma[2] + mb[2])/2)
    points(centroreta[1], centroreta[2], pch = 20, col = "green")


    cangular = (mb[2]-ma[2]) / (mb[1]-ma[1])
    cangular = -1/cangular

    clinear = centroreta[2]-cangular*centroreta[1]
    abline(clinear,cangular, col = "green")

这将创建一条线，将 2D 中的两个点云分开。

我也想做同样的事情，但是是 3D 的。但我无法生成垂直于线的平面

    cloud1 <- rmnorm(100,mean=c(1,1,1),varcov=diag(.25,3)) 
    cloud2 <- rmnorm(75, mean=c(3,3,3),varcov=diag(.25,3))

    install.packages("rgl")
    library(rgl)
    plot3d(cloud1,box=F)
    points3d(cloud2,col="red")

    a <- matrix(c(1,1,1),1,3)
    b <- matrix(c(3,3,3),1,3)
    c <- (a+b)/2

    points3d(a,col="blue")
    points3d(b,col="blue")
    points3d(c,col="blue")

    segments3d(rbind(a,b))

我走到这一步了。现在我无法生成飞机。

谢谢

马塞尔

4

1 回答 1

2

飞机的法线是a-b ，我们想c在上面。

normal <- a-b
offset <- - sum( normal * c )
planes3d( normal[1], normal[2], normal[3], offset, alpha = .2 )

飞机

于 2013-10-04T06:40:44.427 回答