4-laboratoriya mashg‘uloti: hisoblash fizikasi va sonli usullar aniqligi

WolframAlpha va AI (Code Interpreter) muhitida kompyuter modellarining hisoblash aniqligini (integratsiya qadamini) sinab ko‘rish va xatolik grafiklarini qurish

Kompyuter modellarida fizik qonuniyatlar uzluksiz differensial shakldan diskret sonli qadamlarga ajratiladi. Mazkur laboratoriya mashg‘ulotida oddiy mayatnik tebranishi misolida 1-tartibli Eyler usuli bilan aniq analitik yechim taqqoslanadi. Talaba turli integratsiya qadamlari ($\Delta t = 0.2\text{s}, 0.1\text{s}, 0.05\text{s}, 0.01\text{s}$) uchun xatoliklar konvergentsiyasini sinovdan o‘tkazadi, fazaviy energiya portlashini tahlil qiladi hamda WolframAlpha va ChatGPT (Code Interpreter) orqali qiyosiy xatolik grafiklarini quradi.

📋 4-laboratoriya mashg‘uloti rejasi va didaktik konveyeri:

0. Diagnostika va individual reja:
Sonli integratsiya, Taylor qatori va xatoliklar tartibi bo‘yicha kognitiv diagnostika va shaxsiy o‘quv yo‘nalishi.

1. AI integratsiyasi (Claude va Wolfram):
Claude yordamida algoritmik kod yozish va WolframAlpha orqali analitik qadamlarni audit qilish.

2. Sonli integratsiya algoritmi:
Eyler ($\mathcal{O}(\Delta t)$), Eyler–Kromer (Simplektik) va RK4 ($\mathcal{O}(\Delta t^4)$) usullarini matematik qiyoslash.

3. Analitik hisoblash va grafiklar:
Code Interpreter (Python) va WolframAlpha orqali logarifmik xatolik konvergentsiyasini vizuallashtirish.

1. AI integratsiyasi: Claude va WolframAlpha hisoblash bosqichlarini tekshirish va xatolikni tahlil qilishda yordamchi

Sinergik tandem: neyrotarmoq va simvolli matematika

Fizika modellarini kompyuterda sinashda ikkita tubdan farq qiluvchi texnologiya sinergiyasi talab etiladi: Generativ Katta Til Modeli (Claude) hamda Deterministik Simvolli Hisoblash Dvigateli (WolframAlpha). Har ikki tizimning rollari quyidagicha taqsimlanadi:

Generativ LLM va kod sintezi claude.ai ↗

Claude — mantiqiy arxitektor

Differensial tenglamani birinchi tartibli sistemaga keltirish algoritmini tushuntiradi.
Sonli integratsiya (Eyler, Eyler–Kromer, RK4) uchun toza, optimallashgan Python va JavaScript kodlarini sintezlaydi.
Natijalarni jadvallashtirish, xatoliklar dreyfini izohlash va metodik xulosalar tuzishda ko‘maklashadi.

💡 Asosiy vazifasi: Algoritmik kod yozish, gipoteza ilgari surish va tahliliy hisobot yaratish.

Simvolli dvigatel va qat’iy matematika wolframalpha.com ↗

WolframAlpha — aniq hisoblash etaloni

Mayatnik tenglamasining analitik yechimini topadi: kichik burchaklarda garmonik yechim, katta burchaklarda esa davr to‘liq elliptik integral $T = 4\sqrt{L/g}\,K(\sin(\theta_0/2))$ orqali ifodalanadi.
Taylor qatoriga yoyish orqali integratsiya xatoligining nazariy qoldiq hadini ($\mathcal{O}(\Delta t^2)$) hisoblaydi.
Simvolli hisoblashda oraliq natijalar yaxlitlanmaydi; biroq tabiiy tildagi so‘rovni noto‘g‘ri talqin qilishi mumkin, shuning uchun so‘rov sintaksisi va natija talabaning o‘zi tomonidan tekshiriladi.

💡 Asosiy vazifasi: LLM javoblarini mustaqil tekshirish uchun deterministik simvolli hisoblash vositasi.

1.3. WolframAlpha uchun tayyor simvolli so‘rovlar (sintaksis):

Quyidagi so‘rovlarni WolframAlpha qidiruv maydoniga kiritib, analitik hisob-kitoblar etalonini oling:

1-So‘rov: Mayatnikning Aniq Differensial Yechimi (Kichik burchaklar)

DSolve[{theta''[t] + (9.8/1.0)*theta[t] == 0, theta[0] == 0.2, theta'[0] == 0}, theta[t], t]

Natija: $\theta(t) = 0.2 \cos(\sqrt{9.8} t) \approx 0.2 \cos(3.1305 t)$. Bu bizning nazariy etalonimiz hisoblanadi.

2-So‘rov: Taylor Qatoriga Yoyish va Eyler Usulining Qoldiq Hadi

Series[theta0*Cos[omega*t], {t, 0, 3}]

Natija: $\theta(t) = \theta_0 - \frac{1}{2}\theta_0 \omega^2 t^2 + \mathcal{O}(t^4)$. Eyler usuli esa faqat $\theta_0 + 0 \cdot t$ hadini olgani sababli, birinchi qadamdanoq xatolik $\Delta = \frac{1}{2}\theta_0 \omega^2 (\Delta t)^2$ ga teng bo‘ladi!

🌐 WWW Ilmiy-Uslubiy Manba va Simvolli AI (2-Reja):

Wolfram Research. (n.d.). Wolfram|Alpha: Computational Intelligence [Hisoblash tizimi]. https://www.wolframalpha.com
Anthropic. (2024). Claude Model Card Addendum. https://www-cdn.anthropic.com/fed9cc193a14b84131812372d8d5857f8f304c52/Model_Card_Claude_3_Addendum.pdf
Manba mazmuni: WolframAlpha — Wolfram Language asosidagi simvolli hisoblash tizimi (DSolve, Series kabi buyruqlarni bajaradi). Model kartasi Claude’ning kod yozish va mulohaza qilish imkoniyatlari bilan birga uning xato qilishi mumkinligini ham qayd etadi — shu sababli LLM yozgan kod va formulalar WolframAlpha hamda analitik yechim bilan tekshiriladi.

2. Sonli integratsiya algoritmi: Claude va WolframAlpha

Matematik model va xatoliklar dekompozitsiyasi

Matematik mayatnik — uzunligi $L$ bo‘lgan vaznsiz ipga osilgan massasi $m$ bo‘lgan moddiy nuqtadir. Tizimning harakat tenglamasi Nyutonning 2-qonunidan kelib chiqadi:

$$\frac{d^2\theta}{dt^2} + \frac{g}{L} \sin\theta = 0$$

Kichik burchaklar yaqinlashishida ($\sin\theta \approx \theta$) xos tebranish chastotasi: $\omega_0 = \sqrt{\frac{g}{L}}$.

2.1. 1-tartibli differensial tenglamalar tizimiga ajratish:

Kompyuterda sonli yechish uchun 2-tartibli tenglama burchak siljishi $\theta$ va burchak tezligi $\omega = \frac{d\theta}{dt}$ orqali ikkita 1-tartibli tenglamaga keltiriladi:

1-tenglama (burchak o‘zgarishi): $$\frac{d\theta}{dt} = \omega$$

2-tenglama (burchak tezlanishi): $$\frac{d\omega}{dt} = -\frac{g}{L} \sin\theta$$

2.2. Uchta asosiy sonli integratsiya usuli tahlili:

Usul nomi	Rekurrent iteratsiya formulalari	Global xatolik	Energiya saqlanishi
Oddiy Eyler (Forward Euler)	$$\theta_{n+1} = \theta_n + \omega_n \Delta t$$ $$\omega_{n+1} = \omega_n - \frac{g}{L}\sin(\theta_n) \Delta t$$	$\mathcal{O}(\Delta t)$ (1-tartib)	Saqlanmaydi (energiya dreyfi eksponensial o‘sadi)
Eyler–Kromer (Simplektik Eyler)	$$\omega_{n+1} = \omega_n - \frac{g}{L}\sin(\theta_n) \Delta t$$ $$\theta_{n+1} = \theta_n + \omega_{n+1} \Delta t$$	$\mathcal{O}(\Delta t)$ (1-tartib)	Barqaror (fazaviy hajm saqlanadi, energiya tebranadi)
Runge–Kutta 4 (RK4)	$$k_1, k_2, k_3, k_4 \text{ 4 ta oraliq qiyaliklar yig‘indisi}$$ $$y_{n+1} = y_n + \frac{\Delta t}{6}(k_1 + 2k_2 + 2k_3 + k_4)$$	$\mathcal{O}(\Delta t^4)$ (4-tartib)	Yuqori aniqlik (chiziqli mayatnikda $\Delta t=0.1$ s, 5 s oralig‘ida xatolik $\sim 2\cdot10^{-4}$ rad); lekin simplektik emas — juda uzoq vaqtda energiya sekin kamayadi

Claude uchun algoritmik so‘rov (prompt):

Ushbu yo‘riqnomani Claude oynasiga yuboring va algoritmlarning JavaScript / Python realizatsiyasini oling:

Mavzu: Oddiy matematik mayatnik tebranishining sonli integratsiya xatoligi tahlili.
Parametrlar: L = 1.0 m, g = 9.8 m/s^2, theta0 = 0.2 rad, omega0 = 0 rad/s.
Vazifa:
1. Klassik Forward Euler hamda Simplektik Euler-Cromer algoritmlarini Python tilida (funksiya ko‘rinishida) yozib bering.
2. dt = 0.05 s qadam bilan t = 0 dan t = 10 s oralig‘ida mexanik energiyaning qanday o‘zgarishini hisoblab bering: E(t) = 0.5 * m * (L*omega)^2 + m*g*L*(1 - cos(theta)).
3. Eyler usulida har bir qadamda fazaviy hajm qanday o‘sishini (Jakobian determinanti |det(J)| = 1 + (g/L)*dt^2) analitik tushuntirib bering.

🌐 WWW Ilmiy-Uslubiy Manba va Sonli Differensiallash (3-Reja):

Hairer, E., Nørsett, S. P., va Wanner, G. (1993). Solving Ordinary Differential Equations I: Nonstiff Problems (2nd rev. ed.). Springer. https://doi.org/10.1007/978-3-540-78862-1
Butcher, J. C. (2016). Numerical Methods for Ordinary Differential Equations (3rd ed.). John Wiley va Sons. https://doi.org/10.1002/9781119121534
Cromer, A. (1981). Stable Solutions Using the Euler Approximation. American Journal of Physics, 49(5), 455–459. https://doi.org/10.1119/1.12478
Manba mazmuni: Hairer va hammualliflar hamda Butcher Eyler va Runge–Kutta usullarining lokal va global xatolik tartiblarini nazariy asoslaydi. Cromer (1981) Eyler usulidagi kichik o‘zgartirish tebranma masalalarda barqaror yechim berishini ko‘rsatgan. Gamilton tizimlarida simplektik integratorlar (Eyler–Kromer, Verlet) uzoq muddatli hisoblashlarda energiyani nosimplektik usullarga (standart RK4) qaraganda yaxshiroq saqlaydi — lekin qisqa vaqtdagi aniqlikda RK4 ustun.

3. Analitik hisoblash va qiyosiy grafiklar: WolframAlpha yoki ChatGPT (Code Interpreter)

Avtomatlashtirilgan Python tahlili va logarifmik konvergentsiya

Sonli usulning aniqlik tartibini isbotlashning eng qat’iy ilmiy yo‘li — xatoliklar konvergentsiyasi grafigini logarifmik masshtabda ($\log(\text{Xatolik})$ ga nisbatan $\log(\Delta t)$) qurishdir. Nazariyaga ko‘ra:

$$\text{Xatolik} \approx C \cdot (\Delta t)^p \implies \ln(\text{Xatolik}) \approx p \cdot \ln(\Delta t) + \ln(C)$$

Bu yerda to‘g‘ri chiziq qiyaligi $p$ sonli usulning tartibini ko‘rsatadi: Eyler va Eyler–Kromer uchun $p \approx 1$, RK4 uchun $p \approx 4$.
Muhim: analitik etalon $\theta_0\cos(\omega_0 t)$ faqat chiziqli tenglama ($\sin\theta \approx \theta$) yechimi. Agar sonli usulda $\sin\theta$ ishlatilsa, xatolik qadam kichrayganda ham model xatoligi darajasida (masalan, $\theta_0 = 0.2$ rad da $\sim 7\cdot10^{-3}$ rad) to‘xtab qoladi va qiyalik noto‘g‘ri chiqadi. Shu sababli quyidagi skriptda chiziqli tenglama integrallanadi.

ChatGPT (Code Interpreter / Advanced Data Analysis) uchun to‘liq Python skripti:

Foydalanish tartibi: Ushbu kodni nusxalab, ChatGPT (yoki Google Colab / mahalliy Python) oynasiga joylang va "Enter" tugmasini bosing:

import numpy as np
import matplotlib.pyplot as plt

# Mayatnik parametrlari
L = 1.0       # ip uzunligi, m
g = 9.8       # erkin tushish tezlanishi, m/s^2
theta0 = 0.2  # boshlang‘ich burchak, rad (11.5 gradus)
omega0_nat = np.sqrt(g / L)
t_final = 5.0 # modellashtirish vaqti, s

# Sinovdan o‘tkaziladigan integratsiya qadamlari
dt_list = [0.2, 0.1, 0.05, 0.01, 0.005]
max_errors_euler = []
max_errors_cromer = []
max_errors_rk4 = []

for dt in dt_list:
    N = int(t_final / dt)
    t_vals = np.linspace(0, t_final, N + 1)
    
    # Aniq analitik yechim
    theta_exact = theta0 * np.cos(omega0_nat * t_vals)
    
    # Eslatma: analitik etalon chiziqli tenglama (sin(theta) ~ theta) yechimi,
    # shuning uchun konvergentsiya tartibini to‘g‘ri o‘lchash uchun sonli
    # usullar ham aynan shu chiziqli tenglamani integrallaydi.
    accel = lambda th: -(g / L) * th

    # 1. Forward Euler
    th_e = np.zeros(N + 1)
    om_e = np.zeros(N + 1)
    th_e[0] = theta0
    for i in range(N):
        th_e[i+1] = th_e[i] + om_e[i] * dt
        om_e[i+1] = om_e[i] + accel(th_e[i]) * dt
    max_errors_euler.append(np.max(np.abs(th_e - theta_exact)))

    # 2. Euler-Cromer
    th_c = np.zeros(N + 1)
    om_c = np.zeros(N + 1)
    th_c[0] = theta0
    for i in range(N):
        om_c[i+1] = om_c[i] + accel(th_c[i]) * dt
        th_c[i+1] = th_c[i] + om_c[i+1] * dt
    max_errors_cromer.append(np.max(np.abs(th_c - theta_exact)))

    # 3. Runge–Kutta 4
    th_r = np.zeros(N + 1)
    om_r = np.zeros(N + 1)
    th_r[0] = theta0
    for i in range(N):
        th, om = th_r[i], om_r[i]
        k1t, k1o = om, accel(th)
        k2t, k2o = om + 0.5*dt*k1o, accel(th + 0.5*dt*k1t)
        k3t, k3o = om + 0.5*dt*k2o, accel(th + 0.5*dt*k2t)
        k4t, k4o = om + dt*k3o, accel(th + dt*k3t)
        th_r[i+1] = th + dt/6 * (k1t + 2*k2t + 2*k3t + k4t)
        om_r[i+1] = om + dt/6 * (k1o + 2*k2o + 2*k3o + k4o)
    max_errors_rk4.append(np.max(np.abs(th_r - theta_exact)))

# Log-Log konvergentsiya grafigini chizish
plt.figure(figsize=(9, 5))
plt.loglog(dt_list, max_errors_euler, 'o-', color='red', label='Eyler usuli (O(dt^1))')
plt.loglog(dt_list, max_errors_cromer, 's-', color='green', label='Eyler–Kromer usuli (O(dt^1))')
plt.loglog(dt_list, max_errors_rk4, '^-', color='purple', label='RK4 usuli (O(dt^4))')
plt.grid(True, which="both", ls="--", alpha=0.5)
plt.xlabel("Integratsiya qadami Δt (soniya)")
plt.ylabel("Maksimal xatolik max|θ_sonli - θ_aniq| (rad)")
plt.title("Oddiy Mayatnik: Integratsiya Qadami va Xatolik Bog‘liqligi")
plt.legend()
plt.show()

🌐 WWW Ilmiy-Uslubiy Manba va Code Interpreter (4-Reja):

OpenAI. (2023). GPT-4 Technical Report (arXiv:2303.08774). arXiv. https://doi.org/10.48550/arXiv.2303.08774
Press, W. H., Teukolsky, S. A., Vetterling, W. T., va Flannery, B. P. (2007). Numerical Recipes: The Art of Scientific Computing (3rd ed.). Cambridge University Press.
Manba mazmuni: Numerical Recipes ODE integratorlari (Runge–Kutta, qadamni boshqarish) va xatolikni baholashning amaliy qo‘llanmasi. ChatGPT texnik hisoboti model imkoniyatlari bilan birga uning ishonchsiz bo‘lishi va "gallyutsinatsiya" qilishi mumkinligini ham ta’kidlaydi — shu sababli Code Interpreter natijalari (kod, grafik, qiyalik) talaba tomonidan analitik baho bilan solishtirilib tekshiriladi.

4. Amaliy tadqiqot: sonli integratsiya xatoligini hisoblash va interaktiv dastgoh

Jonli mayatnik, fazaviy portret va dinamik xatolik monitori

Ushbu interaktiv laboratoriya stendida siz integratsiya qadami $\Delta t$ ni o‘zgartirib, Eyler usulining etalon yechimdan qanday uzoqlashib ketishini (amplitudaning sun’iy eksponensial o‘sishini) hamda Eyler–Kromer usulining fazaviy barqarorligini real vaqtda kuzatishingiz mumkin. Dastgoh katta burchaklarda ham to‘g‘ri ishlashi uchun nochiziqli tenglama ($\sin\theta$) yechiladi; etalon sifatida juda kichik qadamli ($\Delta t = 10^{-4}$ s) RK4 yechimi olinadi.

⏱ Integratsiya qadami (Δt): 0,050 s

0,005 s (aniq) 0,25 s (dag‘al)

📐 Boshlang‘ich burchak (θ₀): 0,20 rad (11,5°)

0,05 rad (3°) 0,80 rad (46°)

📏 Ip uzunligi (L): 1,00 m

0,5 m 2,0 m

⚙️ Sonli integratsiya usuli:

Vaqt: 0,00 s

Etalon θ: 0,2000 rad

Sonli θ: 0,2000 rad

Xatolik |Δθ|: 0 rad

Energiya E: 0,195 J

🪐 Mayatnik tebranishi va fazaviy portret ($\theta$ ga nisbatan $\omega$):

📈 Burchak $\theta(t)$ va absolyut xatolik $|\theta_{\text{sonli}} - \theta_{\text{etalon}}|$ grafigi:

📊 4 xil qadam uchun avtomatik qiyosiy xatoliklar jadvali ($t = 5{,}0\text{ s}$ holatida):

Ushbu jadval joriy parametrlar ($L$, $\theta_0$) asosida $t = 5{,}0\text{ s}$ vaqt nuqtasida $\Delta t \in \{0{,}2;\ 0{,}1;\ 0{,}05;\ 0{,}01\}\text{ s}$ uchun sonli yechimlarning etalon yechimdan og‘ishini dinamik hisoblab beradi:

Qadam ($\Delta t$)	Etalon ($\theta_{\text{etalon}}$, RK4, $\Delta t=10^{-4}$)	Oddiy Eyler ($\theta_E$) va xatolik	Eyler–Kromer ($\theta_C$) va xatolik	RK4 ($\theta_{RK4}$) va xatolik

5. Baholash rubrikasi (100 ball) va ilmiy manbalar (APA 7)

Akademik baholash va bibliografiya

4-laboratoriya mashg‘uloti talabaning differensial tenglamalarni sonli integrallash nazariyasini tushunishi, WolframAlpha va Code Interpreter orqali xatolik tahlilini to‘g‘ri o‘tkaza olishini 100 ballik shkalada baholaydi:

📊 Interaktiv rubrika kalkulyatori (4-laboratoriya):

1. Nazariy Tushuncha va Diagnostika Testi
Differensial sistemaga ajratish, Eyler usulida energiya dreyfi sabablarini to‘g‘ri tushunish.

Max: 15b

2. Claude va WolframAlpha Simvolli Tahlili
WolframAlpha’da analitik yechim va Taylor qatorini chiqarish, Claude yordamida algoritmni diskretlash.

Max: 25b

3. Turli Δt Qadamlari Uchun Xatolik Jadvalini Tuzish
0.2s, 0.1s, 0.05s, 0.01s qadamlarda Eyler va Eyler–Kromer natijalarini etalon bilan solishtirish.

Max: 25b

4. ChatGPT Code Interpreter (Python) Logarifmik Grafigi
Python skriptini yurgazib, log(xatolik) vs log(dt) konvergentsiya to‘g‘ri chizig‘ini ilmiy asoslash.

Max: 20b

5. Xulosa va Ilmiy Himoya
Simplektik integratorlar ahamiyati va Eyler usulining cheklovlarini asosli bayon etish.