import numpy as np
import matplotlib.pyplot as plt

初始化参数

a = 1 - 0.2 * (1 / 12)
m = 1.109 * 10**5
w3 = 17.86
w4 = 22.99
X = []
Z = []

计算X和Z的值

for k in np.arange(0, 0.875, 0.001):
x1 = 1.221011 * (1 - 1 / (m * (1 - a - 0.42 * k)8 * (1 - a)24 * (0.5 + (1 - a - k)8 * (1 - a)4)))
x2 = (1 - a) * 12 * x1
x3 = (1 - a) * 12 * x2
x4 = (1 - a - 0.42 * k) * 8 * (1 - a)4 * x3
X.append([x1, x2, x3, x4])
z = (0.42 * k * w3 * (1 - (1 - a - 0.42 * k)8) / (a + 0.42 * k) * x3 + k * w4 * (1 - (1 - a - k)8) / (a + k) * x4)
Z.append(z)

找到最大生产量及其对应的k值

mz = max(Z)
ind = Z.index(mz)
print("最大生产量:", mz)
print('各年龄组鱼群数(在最大生产量对应的k值下):', X[ind])
print('对应的k值(用于索引计算，可能需要根据实际情况转换为实际k值):', ind * 0.001)

绘制总生产量Z随捕捞强度k的变化图

plt.plot(np.arange(0, 0.875, 0.001), Z)
plt.ylabel('Z')
plt.xlabel('k')
plt.title('总生产量Z随捕捞强度k的变化')
plt.grid(True)
plt.show()

print("学号后两位:08")
结果如下图所示：