在R中绘制矩阵方程

Question

我是R的新手，我需要绘制二次矩阵方程：

x^T A x + b^T x + c = 0

在R ^ 2中，A为2x2，b为2x1，c为常数。该等式用于定义点类的边界。我需要为x0 = -6...6，x1 = -4...6绘制边界。我的第一个想法是生成一堆点并看到它们在零的位置，但它取决于数字之间的增量（很可能我不会猜测哪些点为零）。

有没有一种更好的方法，而不仅仅是生成一堆点，看看它是零还是乘以它？任何帮助将不胜感激，

谢谢。

Answer 1

假设你有一个对称矩阵A，

例如

 # A = | a    b/2 |
 #     | b/2  c   |

并且您的等式代表圆锥曲线，您可以使用conics package

您需要的是表示

的系数c(a,b,c,d,e,f)的向量

a.x^2 + b*x*y + c*y^2 + d*x + e*y + f

在你的情况下，说你有

 A <- matrix(c(2,1,1,2))

 B <- c(-20,-28)
 C <- 10


# create the vector
v <- append(c(diag(A),B,C),A[lower.tri(A)]*2), 1)



 conicPlot(v)

您可以轻松地将multiplication out打包成一个简单的函数

# note this does no checking for symmetry or validity of arguments

expand.conic <- function(A, B, C){
 append(c(diag(A),B,C),A[lower.tri(A)]*2), 1)
}