import matplotlib.pyplot as plt

import numpy as np

# داده‌ها

روزها = ['روز 0', 'روز 7', 'روز 14', 'روز 28', 'روز 56']

روزها_عددی = np.array([0, 7, 14, 28, 56])

# ویژگی‌های تخمیری

pH_WCC = [6.20, 4.38, 4.28, 4.22, 4.18]

pH_C_TMR = [5.78, 4.45, 4.36, 4.30, 4.26]

لاکتیک_WCC = [4.2, 62.4, 71.5, 75.8, 78.3]

لاکتیک_C_TMR = [6.8, 54.7, 62.9, 66.3, 69.0]

استیک_WCC = [2.6, 18.5, 20.9, 21.7, 22.4]

استیک_C_TMR = [3.9, 14.2, 17.5, 19.1, 20.2]

پروپیونیک_WCC = [0.5, 1.2, 1.4, 1.5, 1.5]

پروپیونیک_C_TMR = [0.7, 0.9, 1.2, 1.3, 1.4]

بوتیریک_WCC = [0.2, 0.5, 0.4, 0.3, 0.3]

بوتیریک_C_TMR = [0.1, 0.3, 0.2, 0.2, 0.2]

آمونیاک_WCC = [44.6, 63.2, 61.7, 59.8, 57.9]

آمونیاک_C_TMR = [32.3, 48.7, 46.9, 45.2, 44.0]

# جمعیت میکروبی

باکتری_لاکتیک_WCC = [5.88, 8.71, 8.58, 8.46, 8.35]

باکتری_لاکتیک_C_TMR = [6.23, 8.42, 8.30, 8.18, 8.08]

باکتری_هوازی_WCC = [6.14, 4.02, 3.65, 3.40, 3.20]

باکتری_هوازی_C_TMR = [5.73, 3.77, 3.44, 3.25, 3.11]

مخمر_WCC = [5.20, 3.40, 3.00, 2.75, 2.60]

مخمر_C_TMR = [4.80, 3.20, 2.90, 2.65, 2.50]

# ترسیم نمودارها

fig, axs = plt.subplots(2, 1, figsize=(12, 10), sharex=True)

# نمودار ویژگی‌های تخمیری

axs[0].plot(روزها_عددی, pH_WCC, marker='o', label='pH - WCC')

axs[0].plot(روزها_عددی, pH_C_TMR, marker='o', label='pH - C-TMR')

axs[0].plot(روزها_عددی, لاکتیک_WCC, marker='s', label='اسید لاکتیک - WCC')

axs[0].plot(روزها_عددی, لاکتیک_C_TMR, marker='s', label='اسید لاکتیک - C-TMR')

axs[0].plot(روزها_عددی, استیک_WCC, marker='^', label='اسید استیک - WCC')

axs[0].plot(روزها_عددی, استیک_C_TMR, marker='^', label='اسید استیک - C-TMR')

axs[0].plot(روزها_عددی, پروپیونیک_WCC, marker='v', label='اسید پروپیونیک - WCC')

axs[0].plot(روزها_عددی, پروپیونیک_C_TMR, marker='v', label='اسید پروپیونیک - C-TMR')

axs[0].plot(روزها_عددی, بوتیریک_WCC, marker='d', label='اسید بوتیریک - WCC')

axs[0].plot(روزها_عددی, بوتیریک_C_TMR, marker='d', label='اسید بوتیریک - C-TMR')

axs[0].plot(روزها_عددی, آمونیاک_WCC, marker='x', label='آمونیاک-N - WCC')

axs[0].plot(روزها_عددی, آمونیاک_C_TMR, marker='x', label='آمونیاک-N - C-TMR')

axs[0].set_title('ویژگی‌های تخمیری در طول زمان')

axs[0].set_ylabel('مقدار (g/kg DM یا pH)')

axs[0].grid(True)

axs[0].legend(loc='upper right', fontsize=8)

# نمودار جمعیت میکروبی

axs[1].plot(روزها_عددی, باکتری_لاکتیک_WCC, marker='o', label='باکتری لاکتیک - WCC')

axs[1].plot(روزها_عددی, باکتری_لاکتیک_C_TMR, marker='o', label='باکتری لاکتیک - C-TMR')

axs[1].plot(روزها_عددی, باکتری_هوازی_WCC, marker='s', label='باکتری هوازی - WCC')

axs[1].plot(روزها_عددی, باکتری_هوازی_C_TMR, marker='s', label='باکتری هوازی - C-TMR')

axs[1].plot(روزها_عددی, مخمر_WCC, marker='^', label='مخمر - WCC')

axs[1].plot(روزها_عددی, مخمر_C_TMR, marker='^', label='مخمر - C-TMR')

axs[1].set_title('جمعیت میکروبی در طول زمان')

axs[1].set_xlabel('روز')

axs[1].set_ylabel('log₁₀ cfu/g FM')

axs[1].grid(True)

axs[1].legend(loc='upper right', fontsize=8)

plt.tight_layout()

plt.show()