我想在下面的問題中填寫3行之間的內容。這是代碼：import numpy as npimport matplotlib.pyplot as plt%matplotlib inline# Construct lines# x > 0x = np.linspace(0, 20, 2000)# C1y1 = (36-2*x)/6# C2y2 = (30-5*x)/3# C3y3 = (40-8*x)/2# C4# y4 = 0*x# Make plotplt.plot(x, y1, label=r'$2 x_{1} + 6 x_{2}\leq 36$')plt.plot(x, y2, label=r'$x_{1} + 3 x_{2}\leq 30$')plt.plot(x, y3, label=r'$x_{1} + 2 x_{2}\leq 40$')# plt.plot(x, y4, label=r'$x_{1}, x_{2}\geq 0$')plt.xlim((0, 16))plt.ylim((0, 11))plt.xlabel(r'$x_1$')plt.ylabel(r'$x_2$')# Fill feasible regiony5 = np.minimum(0, 0)y6 = np.maximum(y2, y3)plt.fill_between(x, y1, y2, color='grey', alpha=0.5, interpolate=True)plt.legend(bbox_to_anchor=(1.05, 1), loc=2, borderaxespad=0.)我想填充下圖中紅色陰影部分（在 y1、y2、y3 和零之間）
2 回答
守候你守候我
TA貢獻1802條經驗 獲得超10個贊
您可以在零和三條曲線的最小值之間進行填充：
import numpy as np
import matplotlib.pyplot as plt
x = np.linspace(0, 20, 2000)
y1 = (36 - 2 * x) / 6
y2 = (30 - 5 * x) / 3
y3 = (40 - 8 * x) / 2
plt.plot(x, y1, label=r'$2 x_{1} + 6 x_{2}\leq 36$')
plt.plot(x, y2, label=r'$x_{1} + 3 x_{2}\leq 30$')
plt.plot(x, y3, label=r'$x_{1} + 2 x_{2}\leq 40$')
plt.xlim((0, 16))
plt.ylim((0, 11))
plt.xlabel(r'$x_1$')
plt.ylabel(r'$x_2$')
plt.fill_between(x, y1, y2, color='grey', alpha=0.5,
interpolate=True)
plt.fill_between(x, 0, np.min([y1, y2, y3], axis=0), color='red', alpha=0.5, hatch='//',
interpolate=True, label='$intersection$')
plt.legend(bbox_to_anchor=(1.05, 1), loc=2, borderaxespad=0.)
plt.tight_layout()
plt.show()
犯罪嫌疑人X
TA貢獻2080條經驗 獲得超4個贊
注意：
np.vstack([y1, y2, y3])從 3 個“y”數組創建一個 3 行數組。….min(0)計算每列的最小值，因此它實際上是 3 個源數組中的最小值（對于較高的x也具有負值）。….clip(min=0)將上述負數元素轉換為0。
因此，添加到您的代碼中：
plt.fill_between(x, 0, np.vstack([y1, y2, y3]).min(0).clip(min=0), color='yellow', alpha=0.5, interpolate=True)
例如在你的第一個fill_ Between之后。
對于您的數據以及添加了上述指令的代碼，我得到：
如果需要，請將填充顏色更改為適合您需要的顏色。
添加回答
舉報
0/150
提交
取消