2017-03-25 88 views
0

我的函數f採用形式爲(1,x,y)的向量,輸出是一個實數。如何在python中繪製隱式定義函數的圖形?

我想繪製一條線,使x(y,y)平面中的f(1,x,y)== 0。

我嘗試:

delta = 0.025 
xrange = arange(0, 12, delta) 
yrange = arange(0, 12, delta) 
p, q = meshgrid(xrange, yrange) 

mesh_point = [] 
for i in range(len(p[0])): 
    for j in range (len(q[0])): 
     mesh_point.append([p[0][i], q[0][j]]) 

for i in range (0, len(p[0])): 
plt.contour(p, q, (y1([1, p[0][i], q[0][i]])-y2([1, p[0][i], q[0][i]])) , [0]) 
plt.show() 

這裏F = Y1(...) - Y2(...)

任何人都可以幫助嗎?

回答

0

通過將contour的左側提供並選擇單個級別(0),可以在f(1,x,y) = 0上繪製輪廓線。

import matplotlib.pyplot as plt 
import numpy as np 

delta = 0.025 
x = np.arange(0, 12, delta) 
y = np.arange(0, 12, delta) 
p, q = np.meshgrid(x, y) 
# define some function f(n,x,y) 
f = lambda n, x, y: np.sin(np.sqrt(np.abs(x-3-0.5*y)**1.5))+np.exp(-(y-6)**2)-1 
z=f(1, p,q) 

# plot contour line of f(1,x,y)==0 
plt.contour(p, q, z , [0], colors=["k"]) 

#make legend 
proxy, = plt.plot([], color="k") 
plt.legend(handles=[proxy], labels=["f(1,x,y) = 0"]) 
plt.show() 

enter image description here