पुनरावर्ती अनुक्रम मॉडल
पुनरावर्ती नियम हर पद को उससे पहले के पद और एक शुरुआती मान से बनाता है। दो सरल नियम:
- समांतर: u(n+1) = u(n) + r। हर कदम पर r जुड़ता है। तब u(n) = u(0) + n·r।
- गुणोत्तर: u(n+1) = q·u(n)। हर कदम पर q से गुणा। तब u(n) = u(0)·qn।
बराबर सालाना बचत के लिए समांतर, ब्याज और जनसंख्या वृद्धि के लिए गुणोत्तर काम आता है।
समांतर-गुणोत्तर अनुक्रम
बहुत-से असली मॉडल दोनों करते हैं: u(n+1) = a·u(n) + b। अगर a ≠ 1 हो तो स्थिर बिंदु L ढूँढो, ऐसा मान जो बदलता नहीं: L = a·L + b, इसलिए L = b ÷ (1 − a)।
अब v(n) = u(n) − L लिखो। तब v(n+1) = a·v(n): यह गुणोत्तर है! इसलिए v(n) = an·v(0) और
u(n) = L + an·(u(0) − L)
अगर |a| < 1 तो an शून्य की ओर जाता है, इसलिए u(n) → L। अगर a > 1 तो यह बिना सीमा बढ़ता है (जब तक u(0) = L न हो)।
समीकरण y' = ay + b
अब समय सतत मानो। y के बदलने की दर y' = a·y + b है। हल है
y(x) = C·eax − b/a (a ≠ 0), जहाँ C शुरुआती मान से मिलता है।
स्थिर हल y = −b/a साम्य है, स्थिर बिंदु जैसा। a < 0 हो तो हर हल उसकी ओर जाता है; a > 0 हो तो उससे दूर।
असतत और सतत मॉडल
असतत मॉडल कदमों में बदलता है (साल में एक बार)। सतत मॉडल हर पल बदलता है। अगर u(n+1) = a·u(n) + b को y' = k·y + c से मिलाएँ जहाँ a = ek, तो वक्र हर बिंदु से ठीक गुज़रता है। निश्चित समय पर नापते हों (सालाना जनगणना) तो असतत लो; बदलाव सहज हो (ठंडा होना, क्षय) तो सतत।
लॉजिस्टिक मॉडल
सादी वृद्धि हमेशा नहीं चल सकती। सीमित जगह में वृद्धि सीमा K के पास धीमी हो जाती है। एक सरल लॉजिस्टिक नियम है
u(n+1) = u(n) + r·u(n)·(1 − u(n)/K)
u छोटा हो तो कोष्ठक 1 के पास है, इसलिए वृद्धि लगभग गुणोत्तर और तेज़। u, K के पास पहुँचे तो कोष्ठक 0 के पास है और वृद्धि रुक जाती है। ग्राफ़ S-आकार का वक्र है। सतत रूप y' = r·y·(1 − y/K) है।
यूलर विधि
ज़्यादातर समीकरणों का सुंदर सूत्र नहीं होता। यूलर विधि छोटे कदम h में ढलान के साथ चलती है:
y(n+1) = y(n) + h·f(y(n)), जहाँ y' = f(y)।
छोटा h बेहतर उत्तर देता है पर ज़्यादा कदम चाहिए। 3D के लॉजिस्टिक बार ठीक h = 1 वाले यूलर कदम हैं।
करके देखो: पहले अनुमान, फिर जाँच
3D के फ़्री-प्ले में a = 0.8, b = 4, शुरुआत = 2 रखो। b ÷ (1 − a) से सीमा का अनुमान लगाओ। हरी रेखा से जाँचो। अब a = 0.5, b = 4 लो: पहले अनुमान, फिर जाँच। आख़िर में a = 1 करो। यह कैसा अनुक्रम है?
मुख्य सूत्र और परिभाषाएँ
- समांतर: u(n) = u(0) + n·r
- गुणोत्तर: u(n) = u(0)·q<sup>n</sup>
- स्थिर बिंदु: L = b ÷ (1 − a), a ≠ 1
- u(n) = L + a<sup>n</sup>·(u(0) − L)
- y' = ay + b का हल y = C·e<sup>ax</sup> − b/a
- यूलर: y(n+1) = y(n) + h·f(y(n))
- लॉजिस्टिक: u(n+1) = u(n) + r·u(n)·(1 − u(n)/K)
हल किए गए उदाहरण
1. u(0) = 2 और u(n+1) = u(n) + 3। u(5) निकालो।
समांतर, r = 3। u(5) = 2 + 5 × 3 = 17।
2. u(0) = 5 और u(n+1) = 2·u(n)। u(4) निकालो।
गुणोत्तर, q = 2। u(4) = 5 × 2⁴ = 5 × 16 = 80।
3. u(0) = 2 और u(n+1) = 0.8·u(n) + 4। स्थिर बिंदु और u(3) निकालो।
L = 4 ÷ (1 − 0.8) = 20। u(n) = 20 − 18 × 0.8ⁿ। u(3) = 20 − 18 × 0.512 = 20 − 9.216 = 10.784। कदमों से जाँच: u1 = 5.6, u2 = 8.48, u3 = 10.784।
4. आप 1000 रुपये जमा करते हैं, साल में 5% ब्याज मिलता है और हर साल 100 रुपये और जोड़ते हैं। नियम लिखो और 2 साल बाद रक़म निकालो।
u(n+1) = 1.05·u(n) + 100। u1 = 1050 + 100 = 1150। u2 = 1.05 × 1150 + 100 = 1207.5 + 100 = 1307.5। यानी 1307.50 रुपये।
5. y' = 2y + 4 को y(0) = 1 के साथ हल करो और y(1) निकालो।
y = C·e^(2x) − 4/2 = C·e^(2x) − 2। y(0) = 1 से C − 2 = 1, इसलिए C = 3। y(1) = 3e² − 2 ≈ 3 × 7.389 − 2 ≈ 20.17।
6. y' = y, y(0) = 1 के लिए h = 0.5 से यूलर विधि लगाकर y(1) का अनुमान लगाओ। सही मान e ≈ 2.718 से तुलना करो।
y1 = 1 + 0.5 × 1 = 1.5। y2 = 1.5 + 0.5 × 1.5 = 2.25। अनुमान 2.25 है; सही मान 2.718 है, यानी यूलर लगभग 0.47 कम बताता है। छोटा h सही मान के और पास देता है।
आम गलतियाँ
- a = 1 होने पर भी L = b ÷ (1 − a) लगाना। तब स्थिर बिंदु नहीं होता (जब तक b = 0 न हो); अनुक्रम समांतर है।
- सोचना कि हर u(n+1) = a·u(n) + b की कोई सीमा होती है। सीमा तभी जब |a| < 1।
- y' = ay + b हल करते समय ऋण चिह्न भूल जाना: स्थिर हल −b/a है, b/a नहीं।
- यूलर विधि में बड़ा h लेकर नतीजे पर भरोसा करना। छोटे h से जाँच लो।