如何在matplotlib中的线之间填充区域

Question

我想在matplotlib中绘制后根据以下等式填充最大面积 尝试了所有可能性，但无法填充所需的区域。

import numpy as np
import matplotlib.pyplot as plt

A = np.linspace(0, 100, 2000)

# 3A+4B≤30

y1 = (30 - A * 3 ) /4
# 5A+6B≤60
y2 = (60 - A * 5)/6
# 1.5A+3B≤21
y3 = (21 - A * 1.5)/3.0

plt.plot(A, y1, label=r'$3A+4B\leq30$')
plt.plot(A, y2, label=r'$5A+6B\leq60$')
plt.plot(A, y3, label=r'$1.5A+3B\leq21$')


plt.xlim((0, 20))
plt.ylim((0, 15))
plt.xlabel(r'$x values$')
plt.ylabel(r'$y values$')

plt.fill_between(A, y3, where = y2<y3,color='grey', alpha=0.5)
plt.legend(bbox_to_anchor=(.80, 1), loc=2, borderaxespad=0.1)
plt.show()

想要填充x = 2.0和y = 6.0的格言区域

Answer 1

这是基于this链接的一种解决方案。与链接解决方案的唯一区别在于，对于您的情况，我不得不使用fill_betweenx覆盖曲线共有的整个x轴，并切换x和Y的顺序。想法是首先找到在一定公差范围内的交点，然后从左一条曲线到该点的曲线和另一条曲线在该交点的曲线获取值。我还必须在[0]中添加一个额外的ind才能使其正常工作

import numpy as np
import matplotlib.pyplot as plt

A = np.linspace(0, 100, 2000)

y1 = (30 - A * 3 ) /4
y2 = (60 - A * 5)/6
y3 = (21 - A * 1.5)/3.0

plt.plot(A, y1, label=r'$3A+4B\leq30$')
plt.plot(A, y2, label=r'$5A+6B\leq60$')
plt.plot(A, y3, label=r'$1.5A+3B\leq21$')

plt.xlim((0, 20))
plt.ylim((0, 12))
plt.xlabel(r'$x values$')
plt.ylabel(r'$y values$')

plt.legend(bbox_to_anchor=(.65, 0.95), loc=2, borderaxespad=0.1)

def fill_below_intersection(x, S, Z):
    """
    fill the region below the intersection of S and Z
    """
    #find the intersection point
    ind = np.nonzero( np.absolute(S-Z)==min(np.absolute(S-Z)))[0][0]
    # compute a new curve which we will fill below
    Y = np.zeros(S.shape)
    Y[:ind] = S[:ind]  # Y is S up to the intersection
    Y[ind:] = Z[ind:]  # and Z beyond it
    plt.fill_betweenx(Y, x, facecolor='gray', alpha=0.5) # <--- Important line

fill_below_intersection(A, y3, y1)

Answer 2

我假设您想填充y1和y3之间的区域，直到它们彼此相交，因为您指定了（2，6）作为点？然后使用：

plt.fill_between(A, y1, y3, where = y1<y3)

如果要表示另一条曲线，则将y3替换为y2。正如@gmds已经评论的那样，“最大面积”有点误导。