R 语言隐函数作图（二元二次函数）

Posted 2020-10-09

在进行论文写作时，需要绘制椭圆图像，隐函数表示如下：

a*(x^2)+b*(y^2)+c*(x*y)-d^2 = 0

其中abcd均为已知值

在matlab中用ezplot函数可以很方便的绘制出该图像：

 ezplot(‘a*x^2+c*x*y+b*y^2=d‘,[-min max ],[-min max])

R语言中没有找到类似相对应的函数，只能使用描点法，假设有两个变元x和y，你先通过分析方程获得x的变化范围，然后通过方程解出y，再画出y对x的图像。

https://stackoverflow.com/questions/12702248/plotting-implicit-function

原文中代码如下：

x<-seq(-1.1,1.1,length=1000)
y<-seq(-1,1,length=1000)
z<-outer(x,y,function(x,y) 4*y^2+x^2+x*y -1 )
contour(x,y,z,levels=0)

根据所需参数替换具体代码就可以绘制成功。