🔢 mathematics

Implementation of Milstein Schemes for Stochastic Delay-Differential Equations with Arbitrary Fixed Delays

यह योगदान मनमाने निश्चित विलंब (arbitrary fixed delays) वाले स्टोकेस्टिक डिले डिफरेंशियल इक्वेशन्स पर लागू यूलर-मारियौमा (Euler–Maruyama) और मिल्स्टीन (Milstein) योजनाओं के लिए संख्यात्मक कार्यान्वयन तकनीकों को प्रस्तुत करता है, जो 1/2 का अभिसरण क्रम (convergence order) प्राप्त करने के लिए रैखिक इंटरपोलेशन (linear interpolation) का उपयोग करता है और 1 का अभिसरण क्रम प्राप्त करने के लिए परिवर्तनशील स्टेप साइज़ के साथ एक विस्तारित टाइम ग्रिड का उपयोग करता है, जिससे कम्युन्सरेट विलंब (commensurate delays) तक सीमित पूर्ववर्ती विधियों की सीमाओं को दूर किया जा सके।

मूल लेखक: Mitchell T. Griggs, Kevin Burrage, Pamela M. Burrage

प्रकाशित 2026-05-05

📖 6 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: Mitchell T. Griggs, Kevin Burrage, Pamela M. Burrage

मूल पेपर CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

कल्पना कीजिए कि आप एक नदी में बहती हुई एक नाव के मार्ग की भविष्यवाणी करने की कोशिश कर रहे हैं। नदी शांत नहीं है; यह लहरों से उथल-पुथल भरी हुई है (यह "स्टोकेस्टिक" या यादृच्छिक भाग है)। लेकिन यहाँ एक मोड़ है: नाव का स्टीयरिंग केवल उन लहरों पर प्रतिक्रिया नहीं देता जो इस समय उससे टकरा रही हैं। बल्कि, यह उन लहरों पर भी प्रतिक्रिया देता है जो पाँच मिनट पहले या यहाँ तक कि दस मिनट पहले उससे टकराई थीं।

वैज्ञानिक इसे स्टोकेस्टिक डिले डिफरेंशियल इक्वेशन (SDDE) कहते हैं। यह एक गणितीय समस्या है जिसका उपयोग स्टॉक की कीमतों, जैविक आबादी या रासायनिक प्रतिक्रियाओं जैसे मॉडलों को दर्शाने के लिए किया जाता है, जहाँ अतीत वर्तमान को एक यादृच्छिक तरीके से प्रभावित करता है।

यह लेख एक बहुत ही विशिष्ट, पेचीदा समस्या से संबंधित है: क्या होता है जब ये "अतीत के समय" आपकी घड़ी के साथ सटीक रूप से मेल नहीं खाते?

समस्या: "गलत तालमेल वाली घड़ी" (The Misaligned Clock)

इन समीकरणों को कंप्यूटर पर हल करने के लिए, वैज्ञानिक आमतौर पर समय को छोटे, समान खंडों में विभाजित करते हैं, जैसे कि हर मिलीमीटर पर निशान वाला एक पैमाना। इसे यूनिफॉर्म टाइम ग्रिड (समान समय ग्रिड) कहा जाता है।

सरल मामला (विभाजित विलंब/Divisible Delays): कल्पना करें कि विलंब (delay) ठीक 5 सेकंड का है और आपकी घड़ी हर सेकंड टिक करती है। यदि कंप्यूटर को यह जानने की आवश्यकता है कि 5 सेकंड पहले क्या हुआ था, तो वह बस 5 टिक पीछे के निशान को देख लेता है। यह बहुत आसान है।
कठिन मामला (अविभाज्य विलंब/Indivisible Delays): अब कल्पना करें कि विलंब $\pi$ सेकंड (लगभग 3.14159...) है, लेकिन आपकी घड़ी अभी भी हर सेकंड टिक करती है। यदि कंप्यूटर को $\pi$ सेकंड पहले की स्थिति जानने की आवश्यकता है, तो वहाँ कोई टिक मार्क नहीं होगा! वह 3-सेकंड के निशान और 4-सेकंड के निशान के बीच फँस जाएगा।

लंबे समय तक, कंप्यूटर केवल "सरल मामले" को ही संभाल सकते थे। यदि विलंब $\pi$ या $\sqrt{2}$ जैसे जटिल नंबर थे जो घड़ी के साथ फिट नहीं बैठते थे, तो वैज्ञानिकों को या तो:

अनुमान लगाना पड़ता था: खाली जगह को भरने के लिए 3-सेकंड और 4-सेकंड के निशानों के बीच एक सीधी रेखा खींचकर मान का अनुमान लगाना (लीनियर इंटरपोलेशन)। यह काम करता है, लेकिन यह धीमा और गलत हो सकता है।
हार माननी पड़ती थी: केवल उन्हीं समस्याओं का अध्ययन करना जहाँ विलंब "अच्छे" या सरल नंबरों में होता था।

समाधान: एक कस्टम रूलर (पैमाना) बनाना

लेखक मिचेल ग्रिग्स, केविन बुराग और पामेला बुराग कहते हैं: "विलंब के अनुरूप घड़ी को बदलने के बजाय, आइए एक कस्टम रूलर बनाते हैं जो विलंब के अनुकूल हो।"

उन्होंने एक विस्तारित समय ग्रिड (extended time grid) बनाने की विधि विकसित की है। इसे ऐसे समझें जैसे आप अपने मानक पैमाने में ठीक वहीं पर अतिरिक्त, कस्टम-मेड निशान जोड़ रहे हैं जहाँ "अतीत" के क्षण स्थित हैं, चाहे वे नंबर कितने भी जटिल क्यों न हों।

यहाँ उनका समाधान सरल शब्दों में बताया गया है:

1. "स्मार्ट" रूलर (विस्तारित ग्रिड)

कंप्यूटर को हर सेकंड चेक करने के लिए मजबूर करने के बजाय, कंप्यूटर उन हर उस क्षण की सूची बनाता है जिसे उसे चेक करने की आवश्यकता है।

यह वर्तमान को चेक करता है।
यह वर्तमान में से विलंब घटाकर चेक करता है।
यह वर्तमान में से दो बार विलंब घटाकर चेक करता है (क्योंकि कभी-कभी अतीत, अतीत के अतीत पर निर्भर करता है)।
यह इन सभी समयों को एक एकल, क्रमबद्ध सूची में व्यवस्थित करता है।

अब, भले ही विलंब $\pi$ हो, कंप्यूटर के पास ठीक $\pi$ सेकंड पीछे के लिए एक विशिष्ट टिक मार्क होगा। अब उसे अनुमान लगाने की आवश्यकता नहीं है।

2. "हाई-डेफिनिशन" कैमरा (मिलस्टीन स्कीम)

लेख एक विशिष्ट गणितीय उपकरण पर ध्यान केंद्रित करता है जिसे मिलस्टीन स्कीम (Milstein Scheme) कहा जाता है।

लो-रेज़ कैमरा (यूलर-मारौमा स्कीम/Euler-Maruyama): यह नाव की एक धुंधली फोटो की तरह है। यह तेज़ है, लेकिन बहुत सटीक नहीं है। यह "अनुमान लगाने" वाली विधि के साथ अच्छा काम करता है, लेकिन यह नाव की गति के सूक्ष्म विवरणों को नहीं पकड़ सकता।
हाई-रेज़ कैमरा (मिलस्टीन): यह एक सुपर-शार्प कैमरा है जो नाव की गति को पूरी तरह से कैप्चर करता है। हालाँकि, इसके लिए बहुत सटीक डेटा की आवश्यकता होती है। यदि आप इसे "अनुमानित" मान (धुंधली फोटो विधि से) देते हैं, तो छवि खराब हो जाती है और हाई-रेज़ कैमरा अपना लाभ खो देता है।

लेखकों की बड़ी सफलता यह दिखाना है कि अपने कस्टम रूलर (विस्तारित ग्रिड) का उपयोग करके, वे हाई-रेज़ कैमरा (मिलस्टीन स्कीम) को वह सटीक डेटा दे सकते हैं जिसकी उसे आवश्यकता है, भले ही विलंब के नंबर जटिल क्यों न हों। यह कंप्यूटर को उन समस्याओं के लिए भी "हाई-डेफिनिशन" परिणाम प्राप्त करने की अनुमति देता है जिन्हें पहले सटीक रूप से हल करना असंभव था।

ट्रेड-ऑफ: अधिक निशान, अधिक काम

एक समस्या है। इस कस्टम रूलर को बनाने के लिए अधिक काम की आवश्यकता होती है।

यदि आपके पास एक सरल विलंब है, तो आपके रूलर में 100 निशान होंगे।
यदि आपके पास एक जटिल विलंब है, तो आपके रूलर को उन सभी "बीच के" क्षणों को कवर करने के लिए 1,000 निशानों की आवश्यकता हो सकती है।

लेख बताता है कि हालांकि इसके लिए अधिक कंप्यूटर मेमोरी और समय की आवश्यकता होती है, लेकिन परिणाम इसके लायक है। "हाई-रेज़" विधि (ऑर्डर 1) "लो-रेज़" विधि (ऑर्डर 1/2) की तुलना में काफी अधिक सटीक हो जाती है, खासकर जब आपको सटीक उत्तरों की आवश्यकता होती है।

लेख में वास्तविक उदाहरण

लेखकों ने इसे वित्तीय विकल्पों (Financial Options) (विशेष रूप से एक "डिजिटल कॉल ऑप्शन") के मॉडल पर परखा।

कल्पना कीजिए कि आप इस बात पर दांव लगा रहे हैं कि भविष्य में स्टॉक की कीमत एक निश्चित संख्या से ऊपर होगी या नहीं।
स्टॉक की कीमत उसके अतीत की अस्थिरता (volatility) पर निर्भर करती है।
यदि वह अतीत का समय एक जटिल नंबर है (जैसे $\pi$ सेकंड पहले), तो पुराने तरीके एक धुंधला और गलत दांव देंगे।
नया तरीका कस्टम रूलर बनाता है, सटीक अतीत की अस्थिरता की गणना करता है, और दांव के मूल्य का बहुत अधिक स्पष्ट और विश्वसनीय पूर्वानुमान प्रदान करता है।

सारांश

यह लेख अनिवार्य रूप से इस बारे में एक मार्गदर्शिका है कि जटिल समय विलंब (time delays) को साफ, समान समय अंतराल में कैसे फिट न किया जाए। इसके बजाय, यह हमें एक लचीला, कस्टम टाइमलाइन बनाने का तरीका सिखाता है जिसमें गणित द्वारा आवश्यक हर एक क्षण शामिल हो। यह हमें उन समस्याओं के लिए अपने सबसे शक्तिशाली, उच्च-परिशुद्धता वाले गणना उपकरणों को लागू करने की अनुमति देता है जो पहले सटीक रूप से हल करने के लिए बहुत जटिल थीं।

तकनीकी सारांश: मनमाने निश्चित विलंब (Arbitrary Fixed Delays) वाले स्टोकेस्टिक डिले डिफरेंशियल इक्वेशंस के लिए मिलस्टीन स्कीम्स का कार्यान्वयन

1. समस्या विवरण

यह शोध पत्र कई निश्चित विलंबों ( $\tau_1, \dots, \tau_K$ ) वाले इटो स्टोकेस्टिक डिले डिफरेंशियल इक्वेशंस (SDDEs) के संख्यात्मक सिमुलेशन को संबोधित करता है, जो एक समान समय ग्रिड (Uniform Time Grid) के साथ संरेखित नहीं होते हैं। हालांकि ऐसे समीकरणों के लिए $1/2$ (यूलर-मारौयामा) और $1$ (मिलस्टीन) की स्ट्रोंग कन्वर्जेंस ऑर्डर (Strong Convergence Order) के लिए सैद्धांतिक ढांचे मौजूद हैं, लेकिन उच्च-क्रम (High-order) स्कीम्स का व्यावहारिक कार्यान्वयन काफी हद तक "विभाज्य विलंबों" (Divisible Delays) तक ही सीमित है—अर्थात, उन मामलों में जहाँ सभी विलंब एक सामान्य स्टेप साइज $h$ के पूर्णांक गुणक होते हैं।

मुख्य कठिनाई इसलिए उत्पन्न होती है क्योंकि उच्च-क्रम स्कीम्स (विशेष रूप से मिलस्टीन स्कीम) को विलंबित समय बिंदुओं ( $t_n - \tau_k$ ) और विलंबित पुनरावृत्त स्टोकेस्टिक इंटीग्रल्स ( $I^{\tau_k}_{ij}$ ) के मूल्यांकन की आवश्यकता होती है। जब विलंब "अविभाज्य" (Indivisible) होते हैं (अर्थात, विलंब और स्टेप साइज के बीच कोई सामान्य विभाजक मौजूद नहीं होता), तो आवश्यक समय बिंदु सिमुलेशन ग्रिड के साथ मेल नहीं खाते। फलस्वरूप, मानक फिक्स्ड-स्टेप कार्यान्वयन इन पदों का सटीक मूल्यांकन करने में विफल रहते हैं, जिससे अनुमान लगाने के लिए मजबूर होना पड़ता है जो कन्वर्जेंस ऑर्डर को घटाकर $1/2$ कर देता है।

2. कार्यप्रणाली (Methodology)

लेखक मनमाने निश्चित विलंबों वाले SDDEs के लिए स्ट्रोंग कन्वर्जेंस ऑर्डर (OoC) 1 प्राप्त करने के लिए एक व्यवस्थित दृष्टिकोण प्रस्तावित करते हैं। कार्यप्रणाली को दो प्राथमिक रणनीतियों में विभाजित किया गया है:

A. OoC 1/2 के लिए लीनियर इंटरपोलेशन (LI)

एक गणनात्मक रूप से कुशल बेंचमार्क प्रदान करने के लिए, लेखक एक निश्चित स्टेप साइज $h$ के साथ स्कीम्स लागू करते हैं। जब विलंबित समय बिंदु ( $t_n - \tau_k$ ) ग्रिड बिंदुओं के बीच स्थित होते हैं, तो विलंबित समाधान मानों को निकटतम सिम्युलेटेड मानों के बीच लीनियर इंटरपोलेशन (LI) के माध्यम से अनुमानित किया जाता है।

स्कीम्स: यूलर-मारौयामा (EM), मिलस्टीन, मैग्नस-यूलर-मारौयामा (MEM), और मैग्नस-मिलस्टीन (MM) को LI का उपयोग करके अनुकूलित किया गया है।
परिणाम: यह दृष्टिकोण सरल और कुशल है लेकिन रिफाइंड इंटीग्रल एप्रोक्सिमेशन का उपयोग करने के बावजूद कन्वर्जेंस ऑर्डर को $1/2$ तक सीमित कर देता है, क्योंकि इंटरपोलेशन एरर ग्लोबल एरर पर हावी हो जाता है।

B. OoC 1 के लिए विस्तारित समय ग्रिड (Extended Time Grid)

सैद्धांतिक OoC 1 प्राप्त करने के लिए, लेखक एक विस्तारित समय ग्रिड का निर्माण करते हैं जो मिलस्टीन स्कीम के मूल्यांकन के लिए आवश्यक सभी समय बिंदुओं को स्पष्ट रूप से शामिल करता है।

निर्माण: ग्रिड का निर्माण एक प्रारंभिक अवलोकन समय बिंदुओं के सेट (जिसमें अंतिम समय $T$ और विलंब के गुणक शामिल हैं) से शुरू करके और अंतराल $[0, T]$ के भीतर विलंब ऑपरेशंस के तहत क्लोज्ड होने तक सभी आवश्यक विलंबित बिंदुओं ( $t_n - \tau_k$ , $t_n - \tau_l - \tau_k$ , आदि) को पुनरावर्ती रूप से जोड़कर किया जाता है।
वेरिएबल स्टेप साइज: इसके परिणामस्वरूप एक नॉन-यूनिफॉर्म ग्रिड प्राप्त होता है जिसमें वेरिएबल स्टेप साइज होता है।
इंटीग्रल सिमुलेशन: लेखक पुनरावृत्त स्टोकेस्टिक इंटीग्रल्स के लिए ट्रेपेज़ॉइडल एप्रोक्सिमेशन विधि (जो विभाज्य विलंबों के लिए मूल रूप से विकसित की गई थी) को इस विस्तारित ग्रिड के लिए विस्तारित करते हैं। वे एक एक्सटेंडेड रिफाइंड टाइम ग्रिड (ARTM) को एक महीन प्रारंभिक स्टेप साइज के साथ लागू करके परिभाषित करते हैं। यह सुनिश्चित करता है कि स्कीम ग्रिड में प्रत्येक अंतराल $[t_n, t_{n+1}]$ के लिए, संगत विलंबित अंतराल $[t_n - \tau_k, t_{n+1} - \tau_k]$ को ARTM द्वारा पूरी तरह से हल किया गया है, जिससे ट्रेपेज़ॉइडल नियम का उपयोग करके विलंबित पुनरावृत्त इंटीग्रल्स $I^{\tau_k}_{ij}$ का सटीक सिमुलेशन संभव हो पाता है।

3. मुख्य योगदान

विस्तारित ग्रिड के लिए ढांचा: यह पेपर एक रचनात्मक एल्गोरिदम (एल्गोरिदम 1) पेश करता है जो एक ऐसा समय ग्रिड तैयार करता है जो मनमाने निश्चित विलंबों को समायोजित करता है, जिससे बिना इंटरपोलेशन के विलंबित स्कीम मानों का सटीक मूल्यांकन संभव हो जाता है।
ट्रेपेज़ॉइडल इंटीग्रेशन का विस्तार: लेखक ARTM का उपयोग करके, अविभाज्य विलंबों के संदर्भ में पुनरावृत्त स्टोकेस्टिक इंटीग्रल्स के लिए ट्रेपेज़ॉइडल एप्रोक्सिमेशन को सामान्यीकृत करते हैं, जिससे यह सुनिश्चित होता है कि विलंबित एकीकरण अंतराल ग्रिड बिंदुओं के साथ संरेखित हों।
कन्वर्जेंस सत्यापन: यह पेपर एक सैद्धांतिक प्रमाण (थ्योरम 4) और संख्यात्मक साक्ष्य प्रदान करता है कि मिलस्टीन और मैग्नस-मिलस्टीन स्कीम्स, विस्तारित ग्रिड पर रिफाइंड इंटीग्रल्स के साथ लागू होने पर, अविभाज्य विलंबों के लिए भी अपने सैद्धांतिक स्ट्रोंग कन्वर्जेंस ऑर्डर 1 को प्राप्त करती हैं।
तुलनात्मक विश्लेषण: यह कार्य फिक्स्ड-स्टेप स्कीम्स (LI के साथ, OoC 1/2) और विस्तारित ग्रिड (OoC 1) पर आधारित स्कीम्स के प्रदर्शन की व्यवस्थित तुलना करता है, यह प्रदर्शित करते हुए कि जबकि विस्तारित ग्रिड अधिक गणनात्मक लागत वहन करता है, ऑर्डर-1 सटीकता प्राप्त करने के लिए यह आवश्यक है।

4. परिणाम

कन्वर्जेंस ऑर्डर्स: संख्यात्मक प्रयोग पुष्टि करते हैं कि फिक्स्ड-स्टेप स्कीम्स और विस्तारित ग्रिड पर आधारित स्कीम्स दोनों ही विभाज्य विलंबों के लिए अपेक्षित व्यवहार करती हैं। हालांकि, अविभाज्य विलंबों के लिए, लीनियर इंटरपोलेशन का उपयोग करने वाली फिक्स्ड-स्टेप स्कीम्स घटकर OoC 1/2 रह जाती हैं। इसके विपरीत, विस्तारित ग्रिड पर आधारित स्कीम्स OoC 1 बनाए रखती हैं।
ग्रिड साइज बनाम विलंब परिमाण: विस्तारित ग्रिड का आकार विलंब के गुणकों की संख्या के साथ द्विघात (quadratically) रूप से बढ़ता है। लेखक नोट करते हैं कि यह विधि तब सबसे व्यावहारिक है जब विलंब अंतिम समय $T$ के सापेक्ष बड़ा हो। छोटे विलंबों के लिए, ग्रिड का आकार अत्यधिक बड़ा हो सकता है।
कंप्यूटेशनल दक्षता: तालिका 3 दिखाती है कि जबकि विस्तारित ग्रिड को LI वाले फिक्स्ड-स्टेप मेथड्स की तुलना में अधिक गणना समय की आवश्यकता होती है, फिर भी उच्च सटीकता की आवश्यकता होने पर यह LI मेथड्स की तुलना में अधिक कुशल हो सकता है। उदाहरण के लिए, विशिष्ट विन्यास में अविभाज्य विलंबों के साथ एक निश्चित त्रुटि स्तर तक पहुँचने के लिए विस्तारित ग्रिड पर मिलस्टीन स्कीम, EM-LI स्कीम की तुलना में तेज़ सिद्ध हुई।
वित्तीय मॉडलों पर अनुप्रयोग: इस विधि को डिजिटल कॉल ऑप्शंस की प्राइसिंग के लिए एक डिलेड लोकल वोलैटिलिटी मॉडल (DLVM) पर लागू किया गया था। विस्तारित ग्रिड पर मिलस्टीन स्कीम ने निम्न-क्रम के तरीकों की तुलना में मोटे समय रिज़ॉल्यूशन पर भी ऑप्शन वैल्यू के स्थिर और सटीक अनुमान प्रदान किए।

5. महत्व और दावे

यह पेपर एक महत्वपूर्ण अंतर को पाटने का दावा करता है जो सामान्य विलंबों के लिए उच्च-क्रम SDDE सॉल्वर के सैद्धांतिक अस्तित्व और उनके व्यावहारिक कार्यान्वयन के बीच है।

व्यावहारिक कार्यान्वयन: लेखक इस बात पर जोर देते हैं कि जबकि निरंतर मिलस्टीन स्कीम को क्लोडेन और शार्डलो [14] द्वारा सैद्धांतिक रूप से स्थापित किया गया था, अविभाज्य विलंबों के लिए कोई मानक कार्यान्वयन रणनीति मौजूद नहीं थी। यह कार्य वही रणनीति प्रदान करता है।
मजबूती (Robustness): यह विधि विभिन्न विलंब कॉन्फ़िगरेशन, जिसमें कई विलंब और अपरिमेय (irrational) विलंब मान (जैसे $\pi/4$ ) शामिल हैं, में मजबूत साबित हुई है।
सीमाएं और भविष्य का कार्य: लेखक विनम्रतापूर्वक स्वीकार करते हैं कि विस्तारित ग्रिड का आकार छोटे विलंबों के लिए एक बाधा बन सकता है। वे सुझाव देते हैं कि ऐसे मामलों के लिए, LI विधि (OoC 1/2) एक व्यवहार्य विकल्प बनी रहती है। भविष्य के विस्तार में समय-निर्भर या वितरित विलंबों को संभालना, साथ ही गणनात्मक लागत को कम करने के लिए एडेप्टिव रिफाइनमेंट या पैरेलल इम्प्लीमेंटेशन की खोज करना शामिल है।

संक्षेप में, यह पेपर निष्कर्ष निकालता है कि एक समान समय ग्रिड की प्रतिबंधात्मकता को त्यागकर एक व्यवस्थित रूप से निर्मित विस्तारित ग्रिड को अपनाकर, मनमाने निश्चित विलंबों वाले SDDEs के लिए स्ट्रोंग कन्वर्जेंस ऑफ ऑर्डर 1 प्राप्त करना संभव हो जाता है, जिससे वित्त और पारिस्थितिकी जैसे क्षेत्रों में जटिल स्टोकेस्टिक मॉडलों के लिए अधिक सटीक सिमुलेशन सक्षम होता है।