1. 程式人生 > >R 語言隱函數作圖(二元二次函數)

R 語言隱函數作圖(二元二次函數)

sta cit get out 替換 隱函數 blog matlab function

在進行論文寫作時,需要繪制橢圓圖像,隱函數表示如下:

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)

根據所需參數替換具體代碼就可以繪制成功。

技術分享

R 語言隱函數作圖(二元二次函數)