ПАЙТОН БОДЛОГО №1 (Сорьцлолтын төлөөлөх дээжний хамгийн бага жин)

Энэхүү python программ нь GY-ийн томьёогоор төлөөлөх дээжний хамгийн бага жинг олно.

#  Gy-ийн томьёоны Python функц
# Author: G. Davaatseren - Python for Mineral Processing Engineers

def gy_min_sample_mass(C: float, d_mm: float, rho: float, f: float) -> float:
     # мм-ийг метр болгох
    d_m = d_mm / 1000.0
    
    # Gy-ийн томьёо
    m_tonnes = (C * (d_m ** 3) * rho) / (f ** 2)
    return m_tonnes
    
#  Жишээ бодлогын өгөгдөл

C = 0.5      # Gy constant
d_mm = 50.0  # top size, мм
rho = 2.7    # t/m3
f = 0.05     # 5% relative error

m_min_t = gy_min_sample_mass(C, d_mm, rho, f)
m_min_kg = m_min_t * 1000  # кг руу хөрвүүлэх

print(f"Хамгийн бага сорьцын масс ≈ {m_min_t:.4f} тн")
print(f"Эсвэл ≈ {m_min_kg:.1f} кг")

#  Ялгаатай ширхэглэл (d) болон алдааны түвшин (f) -ийн нөлөөг харуулах

import numpy as np
import matplotlib.pyplot as plt

C = 0.5
rho = 2.7

# 20 мм-ээс 100 мм хүртэл
d_values = np.linspace(20, 100, 9)  # 20, 30, ..., 100 мм
f_values = [0.1, 0.05, 0.03]        # 10%, 5%, 3%

plt.figure(figsize=(7,5))

for f in f_values:
    m_list = []
    for d_mm in d_values:
        m_list.append(gy_min_sample_mass(C, d_mm, rho, f) * 1000)  # кг
    plt.plot(d_values, m_list, marker="o", label=f"f = {f*100:.0f} %")

plt.xlabel("Хамгийн том ширхэглэл, d (мм)")
plt.ylabel("Хамгийн бага сорьцын масс, m (кг)")
plt.title("Gy-ийн томьёогоор   өөрчлөгдөхөд сорьцын масс")
plt.grid(True)
plt.legend()
plt.show()