R 语言隐函数作图(二元二次函数)


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

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的图像。
原文中代码如下:
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)

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


免责声明!

本站转载的文章为个人学习借鉴使用,本站对版权不负任何法律责任。如果侵犯了您的隐私权益,请联系本站邮箱yoyou2525@163.com删除。



 
粤ICP备18138465号  © 2018-2025 CODEPRJ.COM