4-ma’ruza Hisoblash aniqligi Sonli barqarorlik va simplektik tizimlar

Zamonaviy kompyuter modeliga
qo‘yiladigan talablar va hisoblash aniqligi

Fizik kompyuter modellarining uch fundamental mezoni — Aniqlik (Accuracy), Barqarorlik (Numerical Stability) va Hisoblash Tezligi (Efficiency) o‘rganiladi. Forward Euler, Euler-Cromer, Heun (RK2) va Klassik Runge-Kutta 4 (RK4) usullarining xatolik darajalari, Teylor qatori kesish xatosi (LTE va GTE) hamda integratsiya qadami ($\Delta t$) ning animatsiya sifatiga ta’siri matematik va interaktiv tarzda tahlil qilinadi.

📋 Ma’ruza rejasi va AI integratsiyasi bosqichlari:

1. AI diagnostika va individual reja:

MagicSchool AI orqali talabaning differensial modellar va sonli usullar bo‘yicha tayyorgarligini testlash.

2. Sonli integratsiya hisobi va tahlil:

Claude yordamida Forward Euler, Heun va RK4 diskretlash qadamlarini annotatsiya bilan keltirib chiqarish.

3. Analitik va sonli yechimlar qiyosi:

WolframAlpha yordamida Koshi masalasining aniq analitik yechimini topish va $L_2$ xatolik normasini o‘rnatish.

4. Xatolik grafiklarini qurish:

ChatGPT Advanced Data Analysis orqali $\log(E)$ vs $\log(\Delta t)$ Log-Log konvergensiya egri chiziqlarini qurish.

5. Integratsiya qadamini ($\Delta t$) vizual sinash:

Claude Artifacts asosidagi 4-kamerali interaktiv fizik simulyator va telemetriya sandboxingi.

6. Hisoblash barqarorligi nazariyasi:

Gemini orqali Fon Neyman (von Neumann) spektral tahlili va simplektik fazaviy hajm saqlanishini isbotlash.

🔬 4-ma’ruza bo‘yicha ilmiy AI-pipeline konveyeri:

1-bosqich: Diagnostika va individual reja

MagicSchool AI bilan talaba tayyorgarligini baholash

Fizik modellashtirishda sonli integrallash usullarini tanlashdan oldin, talabaning differensial tenglamalar, Koshi masalasi va diskretlash qadami tushunchalari bo‘yicha bazaviy bilimlari diagnostika qilinadi.

🤖 MagicSchool AI diagnostic prompt namunasi:

Rol: Fizika va hisoblash matematikasi bo‘yicha oliy toifali dotsent.
Vazifa: 3-kurs fizika talabalari uchun "Oddiy differensial tenglamalarni sonli integrallash va hisoblash xatolari" mavzusida 3 darajali diagnostik test va individual o‘quv trayektoriyasi tuzib ber.
Format:
1. Bazaviy daraja: Teylor qatori va 1-tartibli Eyler usuli diskretlash qadami.
2. O‘rta daraja: Konservativ tebranishlarda energiyaning sun’iy dreyfi va faza xatosi.
3. Ilg‘or daraja: Runge-Kutta 4 usulining oraliq qiyaliklari (k1..k4) va simplektik integratsiya mohiyati.
Chiqish formati: Har bir xato javob uchun individual metodik yo‘llanma berilsin.

🎯 Diagnostika maqsadlari:

Talabaning analitik integrallash bilan sonli diskretlash farqini bilishi;
Lokal kesish xatosi ($\text{LTE} \propto \Delta t^2$) va global xato ($\text{GTE} \propto \Delta t$) ni ajrata olishi;
Kompyuter animatsiyasida $\Delta t$ ning 60 FPS kadrlar chastotasi bilan bog‘liqligini tushunishi.

📈 Individual rejalashtirish:

A-trayektoriya (bazaviy): Eyler usuli va xatoliklar vizualizatsiyasidan boshlash;
B-trayektoriya (standart): Euler-Cromer va Heun (RK2) usullarini solishtirish;
C-trayektoriya (ilg‘or): RK4, Verlet va adaptiv qadamli (RK45) tizimlarni dasturlash.

2-bosqich: Sonli integratsiya hisobi

Differensial tenglamalarni sonli diskretlash nazariyasi

Garmonik tebranishning 2-tartibli differensial tenglamasi $\ddot{x} + \omega_0^2 x = 0$ ni 1-tartibli tenglamalar sistemasiga keltiramiz: $$\begin{cases} \frac{dx}{dt} = v \\ \frac{dv}{dt} = -\omega_0^2 x \end{cases}$$

1. Forward Euler (klassik Eyler)

1-tartibli oshkor usul ($\mathcal{O}(\Delta t)$):

                $x_{n+1} = x_n + v_n \Delta t$

                $v_{n+1} = v_n - \omega_0^2 x_n \Delta t$

⚠️ Xato: Faza maydoni har qadamda $(1 + \omega_0^2 \Delta t^2)$ koeffitsiyentga kattalashadi va energiya eksponensial portlaydi!

2. Euler-Cromer (simplektik Eyler)

1-tartibli simplektik usul ($\mathcal{O}(\Delta t)$):

                $v_{n+1} = v_n - \omega_0^2 x_n \Delta t$

                $x_{n+1} = x_n + v_{n+1} \Delta t$

✅ Afzallik: Yangilangan tezlikdan foydalanadi. Faza fazosi hajmi saqlanadi ($\det(J)=1$), tebranish uzoq vaqt barqaror!

3. Runge-Kutta 4 (klassik RK4)

4-tartibli yuqori aniqlik ($\mathcal{O}(\Delta t^4)$):

                $k_1, k_2, k_3, k_4$ qiyaliklari

                $x_{n+1} = x_n + \frac{\Delta t}{6}(k_1^x + 2k_2^x + 2k_3^x + k_4^x)$

⭐ Standart: Xatolik o‘ta kichik ($\sim 10^{-10}$), professional simulyatsiyalarning tayanch dvigateli.

🧠 Claude: RK4 oraliq qiyaliklarining fizik ma’nosi

$k_1 = f(t_n, y_n)$: Interval boshidagi boshlang‘ich qiyalik (Klassik Eyler vektori).
$k_2 = f(t_n + \frac{\Delta t}{2}, y_n + \frac{\Delta t}{2} k_1)$: $k_1$ yordamida oraliq nuqtaga ($\Delta t/2$) borib, u yerdagi yangilangan qiyalikni baholash.
$k_3 = f(t_n + \frac{\Delta t}{2}, y_n + \frac{\Delta t}{2} k_2)$: $k_2$ qiyaligi orqali oraliq nuqtadagi qiyalikni ikkinchi bor aniqlashtirish.
$k_4 = f(t_n + \Delta t, y_n + \Delta t k_3)$: Oraliq qiyalik $k_3$ asosida to‘liq qadam ($\Delta t$) oxiridagi yakuniy qiyalikni bashorat qilish.
Yakuniy natija Simpson qoidasi bo‘yicha o‘rtacha qiyalik sifatida olinadi: $k_{avg} = \frac{k_1 + 2k_2 + 2k_3 + k_4}{6}$.

3-bosqich: Analitik va sonli yechimlar qiyosi

WolframAlpha bilan Koshi masalasini yechish va xatolik normasi

Sonli usullarning xatosini baholash uchun etalon (aniq analitik) yechim kerak bo‘ladi. WolframAlpha orqali garmonik ossilyator tenglamasini aniq integrallaymiz:

              🔍 WolframAlpha Query: DSolve[{x''[t] + omega^2 * x[t] == 0, x[0] == x0, x'[0] == 0}, x[t], t]
            

$$x_{exact}(t) = x_0 \cos(\omega_0 t), \quad v_{exact}(t) = -x_0 \omega_0 \sin(\omega_0 t)$$

📊 Asosiy xatolik metrikalari:

1. Nuqtaviy mutlaq xato:
$e(t_n) = |x_{num}(t_n) - x_{exact}(t_n)|$

2. $L_2$ integrallangan normasi:
$\|e\|_2 = \sqrt{\frac{1}{N} \sum_{n=1}^N \left(x_{num}(t_n) - x_{exact}(t_n)\right)^2}$

3. Faza xatosi ($\Delta \phi$):
$\Delta \phi(t) = \omega_{num} t - \omega_0 t$ (tebranish davrining sonli cho‘zilishi yoki qisqarishi)

4. Nisbiy energiya dreyfi:
$\frac{\Delta E(t)}{E_0} = \frac{|E_{num}(t) - E_0|}{E_0} \times 100\%$

📈 Vaqt bo‘yicha siljish grafigi:

$\Delta t = 0.25\text{ s}$ qadamda Forward Euler (qizil) amplitudasining eksponensial o‘sishi va RK4 (moviy) ning aniq yechim bilan mos kelishi.

4-bosqich: Xatolik grafiklarini qurish

ChatGPT bilan log-log konvergensiya tahlili

Sonli usulning nazariy tartibini ($p$) isbotlash uchun xatolikning qadamga bog‘liqligi logarifmik masshtabda o‘rganiladi: $$E(\Delta t) = C \cdot \Delta t^p \implies \log E(\Delta t) = p \cdot \log(\Delta t) + \log C$$ Log-Log grafigining qiyalik burchagi $\tan\alpha = p$ integratsiyalash tartibini to‘g‘ridan-to‘g‘ri ko‘rsatadi.

📉 Qiyalik burchaklari tahlili:

Usul	Qiyalik ($p$)	$\Delta t = 0.1\text{ s}$ xato	$\Delta t = 0.01\text{ s}$ xato
Forward Euler	$p \approx 1.0$	$\sim 2.5 \times 10^{-1}$	$\sim 2.5 \times 10^{-2}$
Heun / RK2	$p \approx 2.0$	$\sim 4.1 \times 10^{-3}$	$\sim 4.1 \times 10^{-5}$
Klassik RK4	$p \approx 4.0$	$\sim 1.2 \times 10^{-7}$	$\sim 1.2 \times 10^{-11}$

💡 Xulosa: Qadam $\Delta t$ ni 10 barobarga qisqartirganda: Forward Euler xatosi 10 barobar, Heun xatosi 100 barobar, RK4 xatosi esa $10^4 = 10\,000$ barobarga kamayadi!

5-bosqich: integrallash qadamini vizual sinash

Interaktiv xatolik va barqarorlik sinov maydoni

Vaqt qadami $\Delta t$, mayatnik uzunligi $L$ va boshlang‘ich burchak $\theta_0$ ni o‘zgartiring. To‘rt usulning bir vaqtda harakatini, to‘lqin shaklini, faza portretini va energiya dreyfini jonli kuzating:

Vaqt qadami ($\Delta t$): 0.150 s

Mayatnik uzunligi ($L$): 1.00 m

Boshlang‘ich burchak ($\theta_0$): 0.50 rad (28.6°)

Aniq yechim Oshkor Eyler Eyler–Kromer Heun (RK2) Runge–Kutta 4

🎥 Mayatnik harakati animatsiyasi Vaqt: 0,00 s

📈 Siljish grafigi $\theta(t)$, rad Qadamlar: 0

🌀 Faza portreti $\omega(\theta)$

⚡ To‘liq mexanik energiya dreyfi $E(t)/E_0$

📊 Real vaqtdagi hisoblash xatoliklari va energiya telemetriyasi:

Sonli usul	Nazariy tartib	Joriy burchak $\theta(t)$	Mutlaq xato $\|\theta_{\text{son}} - \theta_{\text{aniq}}\|$	Energiya o‘zgarishi $\Delta E/E_0$

6-bosqich: Hisoblash barqarorligi nazariyasi

Gemini: Fon Neyman tahlili va simplektik geometriya

Nima sababdan ba’zi sonli usullar vaqt o‘tishi bilan fizik qonuniyatlarni buzib, energiya portlashini keltirib chiqaradi? Buning sababi Fon Neyman (von Neumann) spektral barqarorligi va Liuvill teoremasi bilan tushuntiriladi.

🔬 Fon Neyman spektral barqarorlik sharti:

Chiziqli differensial operator uchun bitta qadam matritsasi $A$ tuziladi: $\mathbf{u}_{n+1} = A \mathbf{u}_n$. Barqarorlik sharti uning xos qiymatlari $\lambda$ moduli bo‘yicha belgilanadi: $$|\lambda| \le 1$$

Forward Euler: $\det(A) = 1 + \omega_0^2 \Delta t^2 > 1$. Barcha xos qiymatlar moduli birlik doiradan tashqarida ($|\lambda| > 1$). Har qanday garmonik tizimda muqarrar beqaror!
Euler-Cromer: $\det(A) = 1$. Matritsa determinanti aynan 1 ga teng, faza maydoni buzilmaydi ($|\lambda| = 1$).
RK4: Xos qiymatlar kompleks tekislikda keng barqarorlik doirasiga ega ($\omega_0 \Delta t < 2.83$).

🛡️ Simplektik integratorlar va Liuvill teoremasi:

Gamilton mexanikasida fazaviy hajm invariantdir: $d\Gamma = dq \wedge dp = \text{const}$. Simplektik integratorlar (Euler-Cromer, Verlet, Yoshidaning yuqori tartibli sxemalari) aynan shu geometrik simplektik 2-formani aynan saqlaydi.

💡 Veb animatsiyada oltin qoida (Sub-stepping):

Brauzerda ekranni yangilash kadrlar chastotasi o‘zgaruvchan (30 FPS dan 144 FPS gacha) bo‘ladi. Hech qachon brauzerning $\Delta t_{frame}$ ini to‘g‘ridan-to‘g‘ri fizik tenglamaga bermang! Doimo Fixed Timestep (masalan $dt = 0.005\text{ s}$) va kadr ichida bir necha kichik hisoblash qadami tashlaydigan accumulator mexanizmini qo‘llang.

O‘z-o‘zini tekshirish

4-ma’ruza diagnostika testi (3 ta ilmiy savol)

Hisoblash aniqligi, Teylor qatori xatosi va simplektik barqarorlik bo‘yicha bilimlaringizni sinab ko‘ring:

1-savol: Klassik forward Euler usulida garmonik ossilyator modellashtirilganda nima sababdan to‘liq mexanik energiya eksponensial o‘sib ketadi?

⚪ A) Gravitatsiya tezlanishi noto‘g‘ri o‘zgarmas deb olingani uchun

⚪ B) Forward Euler simplektik emas va uning o‘tish matritsasi determinanti det(A) = 1 + ω²Δt² > 1 bo‘lib, faza maydonini har qadamda sun’iy kengaytiradi

⚪ C) Vaqt qadami Δt juda kichik tanlangani tufayli

2-savol: Runge-Kutta 4 (RK4) usulida integratsiya qadami Δt 2 barobarga qisqartirilsa, global kesish xatosi (GTE) necha barobarga kamayadi?

⚪ A) 2 barobarga

⚪ B) 8 barobarga

⚪ C) 16 barobarga (chunki RK4 usulida xatolik tartibi O(Δt⁴) bo‘lib, 2⁴ = 16)

3-savol: Nega real vaqt o‘yinlarida va uzoq davom etuvchi planetar simulyatsiyalarda RK4 o‘rniga ko‘pincha Euler-Cromer yoki Verlet usullari afzal ko‘riladi?

⚪ A) Chunki Euler-Cromer va Verlet simplektik bo‘lib, faza fazosi hajmini aynan saqlaydi va 1 FLOP hisoblash sarfi bilan uzoq vaqt energetik barqarorlikni kafolatlaydi

⚪ B) Chunki ular 5-tartibli aniqlik beradi

⚪ C) Chunki ularda vaqt qadami Δt ni umuman belgilash shart emas

📚 4-ma’ruza bo‘yicha ilmiy-uslubiy adabiyotlar ro‘yxati:

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.
Goldberg, D. (1991). What Every Computer Scientist Should Know About Floating-Point Arithmetic. ACM Computing Surveys, 23(1).
Hairer, E., Nørsett, S. P., va Wanner, G. (1993). Solving Ordinary Differential Equations I: Nonstiff Problems (2nd ed.). Springer.
IEEE. (2019). IEEE Standard for Floating-Point Arithmetic (IEEE 754-2019). IEEE Standards Association.
Landau, R. H., Páez, M. J., va Bordeianu, C. C. (2015). Computational Physics: Problem Solving with Python (3rd ed.). Wiley-VCH.
Chapra, S. C., va Canale, R. P. (2021). Numerical Methods for Engineers (8th ed.). McGraw-Hill.
Sargent, R. G. (2013). Verification and Validation of Simulation Models. Journal of Simulation, 7(1).
Anderson, L. W., va Krathwohl, D. R. (2001). A Taxonomy for Learning, Teaching, and Assessing. Longman.

Zamonaviy kompyuter modeliga qo‘yiladigan talablar va hisoblash aniqligi