2013-07-10 17 views
5

matplotlibを使用して線形不等式関数で囲まれた領域をどのように描画できますか?例えばPython Matplotlib:線形不等式の描画

、我々は3つの機能がある場合: Y < = -2 + 4X、Y> = 2 + 0.5×、Iはタングステンとしてsimmilarのsomethinを描画したいY < = 7 -0.3x

をアルファ:http://www3.wolframalpha.com/Calculate/MSP/MSP43251aca1dfd6ebcd862000067b9fd36a79h3igf?MSPStoreType=image/gif&s=39&w=200.&h=210.&cdf=Coordinates&cdf=Tooltips

+2

、あなたは何を試してみましたか?あなたはギャラリーを見ましたか?私は 'fill_between'があなたが望むことをさせると思うhttp://matplotlib.org/examples/pylab_examples/fill_between_demo.html – tacaswell

答えて

10

非常に簡単な例を書きましたが、問題にのみ有効ですが、簡単に拡張して一般化することができます。唯一のやり方は、所望のポリゴンを構築するためのルーツを見つける問題を単純化するためにsimpyを使用することです。だから、(http://docs.sympy.org/dev/modules/solvers/solvers.htmlからの引用)

import numpy as np 
import matplotlib.pyplot as plt 
from sympy.solvers import solve 
from sympy import Symbol 

def f1(x): 
    return 4.0*x-2.0 
def f2(x): 
    return 0.5*x+2.0 
def f3(x): 
    return -0.3*x+7.0 

x = Symbol('x') 
x1, = solve(f1(x)-f2(x)) 
x2, = solve(f1(x)-f3(x)) 
x3, = solve(f2(x)-f3(x)) 

y1 = f1(x1) 
y2 = f1(x2) 
y3 = f2(x3) 

plt.plot(x1,f1(x1),'go',markersize=10) 
plt.plot(x2,f1(x2),'go',markersize=10) 
plt.plot(x3,f2(x3),'go',markersize=10) 

plt.fill([x1,x2,x3,x1],[y1,y2,y3,y1],'red',alpha=0.5) 

xr = np.linspace(0.5,7.5,100) 
y1r = f1(xr) 
y2r = f2(xr) 
y3r = f3(xr) 

plt.plot(xr,y1r,'k--') 
plt.plot(xr,y2r,'k--') 
plt.plot(xr,y3r,'k--') 

plt.xlim(0.5,7) 
plt.ylim(2,8) 

plt.show() 

enter image description here

よろしく

+0

ありがとう、これはまさに私が必要なものです:) – N10

関連する問題