черчение матричного уравнения в R

Question

Я новичок в R, и мне нужно, чтобы построить квадратное уравнение матрицы:

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

в R^2, причем А 2х2, 2х1 ба, и постоянная ча , Уравнение предназначено для границы, определяющей классы точек. Мне нужно построить эту границу для 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) 
}