Plot Числовой диапазон матрицы в R - PullRequest
Новая
       31

Plot Числовой диапазон матрицы в R

0 голосов
/ 07 апреля 2020

Как я могу построить числовой диапазон эрмитовой матрицы в R? Например, у меня есть следующая матрица: 


a = matrix(c(8+5i,6,0,14+5i),2,2,byrow = T)
a
eigen(a)

Числовой диапазон представляет собой эллипс с фокусами собственных значений матрицы. Например. 

qplot(Re(c(14,8)), Im(c(5i,5i)), geom="point")

Но как я могу нарисовать эллипс в R. Я попробовал это, но не получилось: 


ellipseFun <- function(center = c(0,0),axis = c(0,0), npoints = 100){
  phi =  pi/3
  tt  =  seq(0,2*pi,length.out = npoints)
  xx  =  center[1] + axis[1]*cos(t)*cos(phi) - axis[2]*sin(t)*sin(phi)
  yy  =  center[2] + axis[1]*cos(t)*cos(phi) + axis[2]*sin(t)*cos(phi)
  return(data.frame(x = xx, y = yy))
}
dat <- ellipseFun(c(14,8),c(20,20),npoints = 100)
head(dat)
dat1 = data.frame(x=c(14,8),y=c(0,0));dat1
p = ggplot (dat1, aes (x = x, y = y))+
  geom_point()
p
p+geom_path(data=dat)
...