2-बिंदु गॉसियन और समलंब चतुर्भुज के बीच तुलना। नीला वक्र उस फलन को दर्शाता है जिसका अंतराल पर निश्चित समाकल है [−1, 1] की गणना (एकीकृत) की जानी है। समलम्बाकार नियम फ़ंक्शन को रैखिक फ़ंक्शन के साथ अनुमानित करता है जो अंतराल के अंत बिंदुओं पर समाकलित के साथ मेल खाता है और नारंगी धराशायी रेखा द्वारा दर्शाया जाता है। सन्निकटन स्पष्ट रूप से अच्छा नहीं है, इसलिए त्रुटि बड़ी है (समलम्बाकार नियम समाकल के बराबर सन्निकटन देता है y(–1) + y(1) = –10, जबकि सही मान 2⁄3 है )। अधिक त्रुटिहीन परिणाम प्राप्त करने के लिए, अंतराल को कई उप-अंतरालों में विभाजित किया जाना चाहिए और फिर समग्र समलम्बाकार नियम का उपयोग किया जाना चाहिए, जिसके लिए बहुत अधिक गणनाओं की आवश्यकता होती है। गॉसियन चतुर्भुज इसके अतिरिक्त अधिक उपयुक्त बिंदु चुनता है, इसलिए रैखिक फ़ंक्शन भी फ़ंक्शन को श्रेष्ठतर (काली धराशायी रेखा) अनुमानित करता है। जैसा कि समाकलित घात 3 का बहुपद है (y(x) = 7x3 – 8x2 – 3x + 3), 2-बिंदु गॉसियन चतुर्भुज नियम भी त्रुटिहीन परिणाम देता है।
संख्यात्मक विश्लेषण में, चतुर्भुज नियम फ़ंक्शन (गणित) के समाकल का अनुमान है, जिसे सामान्यतः समाकलन के डोमेन के भीतर निर्दिष्ट बिंदुओं पर फ़ंक्शन मानों के भारित योग के रूप में कहा जाता है। (चतुर्भुज (गणित) नियमों पर अधिक जानकारी के लिए संख्यात्मक समाकलन देखें।) An n-बिंदु गॉसियन चतुर्भुज नियम, कार्ल फ्रेडरिक गॉस के नाम पर,[1] घात के बहुपद के लिए त्रुटिहीन परिणाम प्राप्त करने के लिए निर्मित चतुर्भुज नियम है 2n − 1 या उससे कम नोड्स के उपयुक्त विकल्प द्वारा xi और भार wi के लिए i = 1, …, n. 1826 में कार्ल गुस्ताव जैकोबी द्वारा ऑर्थोगोनल बहुपद का उपयोग करते हुए आधुनिक सूत्रीकरण विकसित किया गया था।[2] इस प्रकार के नियम के लिए समाकलन का सबसे साधारण डोमेन लिया जाता है [−1, 1], इसलिए नियम के रूप में कहा गया है
जो 2n − 1 या उससे कम घात वाले बहुपदों के लिए है । इस त्रुटिहीन नियम को गॉस-लेजेंड्रे चतुर्भुज नियम के रूप में जाना जाता है। चतुष्कोण नियम उपरोक्त समाकलन के लिए केवल त्रुटिहीन सन्निकटन होगा यदि f (x) को [−1, 1] पर घात 2n − 1 या उससे कम के बहुपद द्वारा अच्छी प्रकार से अनुमानित किया गया है।
एड्रियन मैरी लीजेंड्रे क्वाडरेचर नियम सामान्यतः समापन बिंदु विलक्षणता (गणित) के साथ पूर्णांक कार्यों के लिए उपयोग नहीं किया जाता है। इसके अतिरिक्त, यदि समाकलित को इस रूप में लिखा जा सकता है
जहाँ g(x) कम-घात बहुपद, फिर वैकल्पिक नोड्स द्वारा अच्छी प्रकार से अनुमानित है xi' और भार wi' सामान्यतः अधिक त्रुटिहीन चतुर्भुज नियम देगा। इन्हें गॉस-जैकोबी चतुष्कोण नियम के रूप में जाना जाता है, अर्थात,
यह दिखाया जा सकता है (देखें प्रेस, एट अल।, या स्टोअर और बुलिरश) कि चतुर्भुज नोड्स xi ओर्थोगोनल बहुपदों के वर्ग से संबंधित बहुपद के फलन के मूल हैं (भारित आंतरिक-उत्पाद के संबंध में वर्ग ऑर्थोगोनल)। गॉस क्वाडरेचर नोड्स और वजन की गणना के लिए यह महत्वपूर्ण अवलोकन है।
ऊपर बताई गई सरलतम समाकलन समस्या के लिए, अर्थात, f(x) पर बहुपदों द्वारा अच्छी प्रकार से अनुमानित है , संबंधित ऑर्थोगोनल बहुपद लीजेंड्रे बहुपद हैं, जिन्हें द्वारा दर्शाया गया है Pn(x). साथ n-वाँ बहुपद देने के लिए सामान्यीकृत Pn(1) = 1, द i-वां गॉस नोड, xi, है i-की जड़ Pn और भार सूत्र द्वारा दिए गए हैं[3]
कुछ निम्न-क्रम द्विघात नियम नीचे सारणीबद्ध हैं (अंतराल पर [−1, 1], अन्य अंतरालों के लिए नीचे दिया गया अनुभाग देखें)।
बिंदुओं की संख्या, n
अंक, xi
भार, wi
1
0
2
2
±0.57735...
1
3
0
0.888889...
±0.774597...
0.555556...
4
±0.339981...
0.652145...
±0.861136...
0.347855...
5
0
0.568889...
±0.538469...
0.478629...
±0.90618...
0.236927...
मध्यान्तर का परिवर्तन
[a, b] समाकल पर ऊपर को [−1, 1] समाकल ऊपर में बदला जाना चाहिए , गाऊसी चतुर्भुज नियम लागू करने से पहले। मध्यान्तर का यह परिवर्तन निम्न प्रकार से किया जा सकता है:
साथ लागू करना बिंदु गाऊसी चतुर्भुज नियम तो निम्नलिखित सन्निकटन में परिणाम:
दो-बिंदु गॉस चतुर्भुज नियम का उदाहरण
किसी रॉकेट द्वारा स्थापित की गई दूरी को मीटर में अनुमानित करने के लिए दो-बिंदु गॉस चतुर्भुज नियम का उपयोग करें को जैसा दिया गया है
सीमाएं बदलें जिससे कि तालिका 1 में दिए गए भार और भुज का उपयोग किया जा सके। साथ ही, पूर्ण सापेक्ष सत्य त्रुटि भी पाएं। सही मान 11061.34 मीटर के रूप में दिया गया है।
समाधान
सबसे पहले, समाकलन की सीमाओं को बदलना को देता है
अगला, दो-बिंदु नियम के लिए भार कारक और फ़ंक्शन तर्क मान तालिका 1 से प्राप्त करें,
अब हम गॉस चतुर्भुज सूत्र का उपयोग कर सकते हैं
तब से
यह देखते हुए कि सही मान 11061.34 मीटर है, पूर्ण सापेक्ष सत्य त्रुटि है
अन्य रूप
सकारात्मक भार फलन प्रारंभ करके समाकलन समस्या को थोड़ा और सामान्य विधियों से व्यक्त किया जा सकता है ω समाकलित में, और इसके अतिरिक्त अंतराल की अनुमति देता है [−1, 1]. अर्थात समस्या गणना करने की है
कुछ विकल्पों के लिए a, b, और ω के लिए a = −1, b = 1, और ω(x) = 1, समस्या वही है जो ऊपर मानी गई है। अन्य विकल्प अन्य समाकलन नियमों की ओर ले जाते हैं। इनमें से कुछ नीचे सारणीबद्ध हैं। अब्रामोवित्ज़ और स्टेगुन (A & S) के लिए समीकरण संख्याएँ दी गई हैं।
होने देना pn घात का गैर-तुच्छ बहुपद हो n ऐसा है कि
ध्यान दें कि यह उपरोक्त सभी ऑर्थोगोनल बहुपदों के लिए सत्य होगा, क्योंकि प्रत्येक pn का निर्माण अन्य बहुपदों के लिए ओर्थोगोनल होने के लिए किया गया है pj के लिए j<n, और xk उस समुच्चय की अवधि में है।
यदि हम चुनते हैं n नोड्स xi का शून्य होना pn, तो वहाँ उपस्तिथ हैं n भार wi जो सभी बहुपदों के लिए गॉस-चतुर्भुज परिकलित समाकल को त्रुटिहीन बनाता है h(x) घात 2n − 1 या कम। इसके अतिरिक्त , ये सभी नोड्स xi खुले अंतराल में होगा (a, b).[4]इस प्रमाण के पहले भागों को सिद्ध करने के लिए होने देना h(x) कोटि का कोई भी बहुपद हो 2n − 1 या कम। प्राप्त करने के लिए इसे ओर्थोगोनल बहुपद pn से विभाजित करें
जहाँ q(x) भागफल है, घात का n − 1 या उससे कम (क्योंकि इसकी घात का योग और विभाजक का pn लाभांश के बराबर होना चाहिए), और r(x) शेष है, घात का भी n − 1 या उससे कम (क्योंकि शेष की घात सदैव भाजक की घात से कम होती है)। तब से pn से कम घात के सभी एकपदीयों के लिए ऑर्थोगोनल है n, यह भागफल के लिए ओर्थोगोनल q(x) होना चाहिए । इसलिए
शेष के बाद से r(x) घात का है n − 1 या उससे कम, हम इसका उपयोग करके इसे प्रक्षेपित कर सकते हैं nलैग्रेंज बहुपद के साथ प्रक्षेप बिंदु li(x), जहाँ
अपने पास
तब इसका समाकल बराबर होगा
जहाँ wi, नोड से xi जुड़ा भार , li(x) के भारित समाकल के बराबर परिभाषित किया गया है (भार के लिए अन्य सूत्रों के लिए नीचे देखें)। किन्तु सभी xi की जड़ें हैं pn, तो उपरोक्त विभाजन सूत्र हमें बताता है
सभी के लिए i. इस प्रकार अंत में हमारे पास है
यह सिद्ध करता है कि किसी भी बहुपद के लिए h(x) घात 2n − 1 या उससे कम, इसका समाकल बिल्कुल गाऊसी चतुर्भुज योग द्वारा दिया जाता है।
प्रमाण के दूसरे भाग को सिद्ध करने के लिए, pn बहुपद के गुणनखंडित रूप पर विचार करें । कोई भी जटिल संयुग्मी जड़ें द्विघात कारक उत्पन्न करेंगी जो पूरी वास्तविक रेखा पर या तो सख्ती से सकारात्मक या सख्ती से नकारात्मक है। a को b से अंतराल के बाहर जड़ों के लिए कोई कारक उस अंतराल पर चिह्न नहीं बदलेगा। अंत में, जड़ों से संबंधित कारकों के लिए xi अंतराल के अंदर से a को b जो विषम बहुलता के हैं, गुणा करें pn और गुणक द्वारा नया बहुपद बनाने के लिए
यह बहुपद a से b तक के अंतराल पर चिह्न नहीं बदल सकता, क्योंकि इसकी सारी जड़ें अब भी बहुलता की हैं। तो समाकल
भार फलन के बाद से ω(x) सदैव गैर-ऋणात्मक होता है। किन्तु pn घात के सभी बहुपदों के लिए ऑर्थोगोनल है n-1 या उससे कम, तो उत्पाद की डिग्री
कम से कम n होना चाहिए . इसलिए pn के n भिन्न मूल हैं, सभी वास्तविक हैं, a से b के अंतराल में ।
भार के लिए सामान्य सूत्र
भार के रूप में व्यक्त किया जा सकता है
(1)
जहाँ का गुणांक में है। यह सिद्ध करना करने के लिए, ध्यान दें कि लैग्रेंज प्रक्षेप का उपयोग करके r(x) व्यक्त किया जा सकता है के अनुसार, जैसा
क्योंकि r(x) से कम घात है n और इस प्रकार इसे प्राप्त होने वाले मूल्यों द्वारा स्थापित किया जाता है n विभिन्न बिंदु। दोनों पक्षों को से गुणा करना ω(x) और से समाकलित a को b प्रस्तुतीकरण
भार wi इस प्रकार दिए गए हैं
के लिए यह समाकल अभिव्यक्ति ऑर्थोगोनल बहुपदों और के संदर्भ में निम्नानुसार व्यक्त किया जा सकता है।
हम लिख सकते हैं
जहाँ का गुणांक है में . की सीमा ले रहा है x को ल'हॉपिटल नियम का उपयोग करके प्रस्तुतीकरण
हम इस प्रकार भार के लिए समाकल अभिव्यक्ति लिख सकते हैं
(2)
समाकलित में, लेखन
प्रस्तुतीकरण
परंतु , क्योंकि
घात का बहुपद है k − 1 जो तब ओर्थोगोनल है . तो यदि q(x) हमारे पास अधिकतम n वें घात का बहुपद है
हम के लिए दाहिने ओर समाकल का मूल्यांकन इस प्रकार कर सकते हैं। क्योंकि घात का बहुपद है n − 1, अपने पास
जहाँ s(x) घात का बहुपद है । तब से s(x) ओर्थोगोनल है अपने पास
हम तब लिख सकते हैं
कोष्ठक में शब्द घात का बहुपद है , जो इसलिए ओर्थोगोनल है । समाकलन को इस प्रकार लिखा जा सकता है
समीकरण के अनुसार (2), भार को इससे विभाजित करके प्राप्त किया जाता है और वह समीकरण में (1)अभिव्यक्ति देता है।
ऑर्थोगोनल बहुपदों और अब के संदर्भ में भी व्यक्त किया जा सकता है। 3-टर्म पुनरावृत्ति संबंध में के साथ शब्द गायब हो जाता है, इसलिए Eq में (1) द्वारा प्रतिस्थापित किया जा सकता है ।
प्रमाण है कि भार सकारात्मक हैं
घात के निम्नलिखित बहुपद पर विचार करें
जहां, ऊपर के रूप में, xj बहुपद के मूल हैं । स्पष्ट रूप से . चूंकि की घात इससे कम है , गाऊसी चतुर्भुज सूत्र से प्राप्त भार और नोड्स सम्मलित हैं लागू होता है। तब से j के लिए i के बराबर नहीं, हमारे पास है
चूंकि दोनों और गैर-नकारात्मक कार्य हैं, यह इस प्रकार है .
गौसियन चतुर्भुज नियमों की गणना
नोड्स की गणना के लिए कई एल्गोरिदम हैं xi और भार wi गाऊसी चतुर्भुज नियम। सबसे लोकप्रिय गोलूब-वेल्श एल्गोरिदम की आवश्यकता है O(n2) संचालन, हल करने के लिए न्यूटन की विधि ओर्थोगोनल बहुपदों का उपयोग करना हैं। पुनरावृत्ति संबंध मूल्यांकन के लिए तीन-अवधि की पुनरावृत्ति की आवश्यकता होती है O(n2) संचालन और बड़े n आवश्यकता के लिए स्पर्शोन्मुख सूत्र O(n) संचालन।
पुनरावृत्ति संबंध
ऑर्थोगोनल बहुपद साथ के लिए अदिश उत्पाद के लिए , घात और प्रमुख गुणांक, अर्थात मोनिक बहुपद ऑर्थोगोनल बहुपद पुनरावृत्ति संबंध को संतुष्ट करते हैं
और अदिश उत्पाद परिभाषित
के लिए, जहाँ n अधिकतम घात है जिसे अनंत माना जा सकता है, और जहाँ . सबसे पहले, से प्रारंभ होने वाले पुनरावृत्ति संबंध द्वारा परिभाषित बहुपद अग्रणी गुणांक और सही घात है। द्वारा प्रारंभिक बिंदु दिया गया , की रूढ़िवादिता प्रवर्तन द्वारा दिखाया जा सकता है। के लिए, किसी के पास
अब यदि ओर्थोगोनल हैं, फिर भी , क्योंकि
पहले और को छोड़कर सभी अदिश उत्पाद गायब हो जाते हैं समान लंबकोणीय बहुपद को पूरा करता है। इसलिए,
चूंकि, यदि अदिश उत्पाद संतुष्ट करता है जो गौसियन चतुर्भुज के स्थितियों में है, पुनरावृत्ति संबंध तीन-अवधि के पुनरावृत्ति संबंध में कम हो जाता है। के लिए से कम या बराबर घात का बहुपदr − 1 है। वहीं दूसरी ओर, से कम या बराबर घात के हर बहुपद के लिए r − 1 ओर्थोगोनल है । इसलिए,s < r − 1 के लिए और है। पुनरावृत्ति संबंध तब सरल हो जाता है
या
(सम्मेलन के साथ ) जहाँ
(आखिरी के कारण , तब से से मतभेद होना घात से कम है r).
गोलूब-वेल्श एल्गोरिथम
तीन-पद पुनरावृत्ति संबंध को आव्यूह के रूप में लिखा जा सकता है, जहाँ , है मानक आधार सदिश, अर्थात, , और J तथाकथित जैकोबी आव्यूह है:
शून्य घात तक बहुपदों का n, जो गाऊसी चतुर्भुज के लिए नोड्स के रूप में उपयोग किया जाता है, इस त्रिकोणीय आव्यूह के आइजन मूल्य की गणना करके पाया जा सकता है। इस प्रक्रिया को गोलूब-वेल्श एल्गोरिथम के रूप में जाना जाता है।
भार और नोड्स की गणना के लिए, सममित त्रिभुज आव्यूह पर विचार करना श्रेष्ठतर होता है तत्वों के साथ
J और समान आव्यूह हैं और इसलिए समान आइजन मूल्य( नोड्स) हैं। वज़न की गणना संबंधित आइजन वैक्टर से की जा सकती है, यदि सामान्यीकृत आइजन वैक्टर है (अर्थात , यूक्लिडियन मानदंड के बराबर आइजन वैक्टर ) आइजन मूल्य से जुड़ा हुआ है xj, इस आइजन वैक्टर के पहले घटक से संबंधित भार की गणना की जा सकती है, अर्थात्:
गॉसियन चतुर्भुज नियम की त्रुटि निम्नानुसार बताई जा सकती है।[5] समाकलित के लिए जिसके पास 2n निरंतर व्युत्पन्न है,
(a, b) कुछ के लिए ξ में , जहाँ pn मोनिक है (अर्थात अग्रणी गुणांक है 1) n घात का ऑर्थोगोनल बहुपद और जहाँ
ω(x) = 1 के महत्वपूर्ण विशेष स्थितियों में , हमारे पास त्रुटि अनुमान है[6]
स्टॉयर और बुलिरश टिप्पणी करते हैं कि यह त्रुटि अनुमान अभ्यास में असुविधाजनक है, क्योंकि आदेश का अनुमान लगाना कठिन हो सकता है 2n व्युत्पन्न और इसके अतिरिक्त वास्तविक त्रुटि व्युत्पन्न द्वारा स्थापित सीमा से बहुत कम हो सकती है। अन्य दृष्टिकोण विभिन्न आदेशों के दो गॉसियन चतुर्भुज नियमों का उपयोग करना है, और दो परिणामों के बीच अंतर के रूप में त्रुटि का अनुमान लगाना है। इस प्रयोजन के लिए, गॉस-क्रोनरोड चतुर्भुज नियम उपयोगी हो सकते हैं।
यदि अंतराल [a, b] उप-विभाजित है, नए उप-अंतरालों के गॉस मूल्यांकन बिंदु पिछले मूल्यांकन बिंदुओं (विषम संख्याओं के लिए शून्य को छोड़कर) के साथ कभी मेल नहीं खाते हैं और इस प्रकार प्रत्येक बिंदु पर पूर्णांक का मूल्यांकन किया जाना चाहिए। गॉस-क्रोनरोड नियम गौस चतुर्भुज नियमों के विस्तार हैं जो n + 1 जोड़कर उत्पन्न होते हैं n की ओर संकेत करता है 2n + 1 -बिंदु नियम इस प्रकार है कि परिणामी नियम क्रम का हो। यह निम्न-क्रम अनुमान के फ़ंक्शन मानों का पुन: उपयोग करते हुए उच्च-क्रम अनुमानों की गणना करने की अनुमति देता है। गॉस चतुर्भुज नियम और इसके क्रोनरोड प्रसार के बीच का अंतर अधिकांशतः सन्निकटन त्रुटि के अनुमान के रूप में उपयोग किया जाता है।
गॉस-लोबेटो नियम
लोबेटो चतुर्भुज के रूप में भी जाना जाता है,[7] डच गणितज्ञ रेहुएल लोबेटो के नाम पर। यह निम्नलिखित अंतरों के साथ गॉसियन चतुर्भुज के समान है:
समाकलन बिंदुओं में समाकलन अंतराल के अंतिम बिंदु सम्मलित होते हैं।
यह घात तक के बहुपदों के लिए त्रुटिहीन है 2n – 3, जहाँ n समाकलन बिंदुओं की संख्या है।[8]
कार्य का लोबेटो चतुर्भुज f(x) अंतराल पर [−1, 1]:
भुज: xi है गोसाई शून्य ,यहाँ m-वें घात के मानक लेजेंड्रे बहुपद को दर्शाता है और डैश निरूपित को दर्शाता है।
वजन,
शेष,
कुछ भार हैं,
बिंदुओं की संख्या, n
अंक, xi
भार, wi
2 आंतरिक नोड्स के साथ इस एल्गोरिथम का अनुकूली संस्करण[9]जीएनयू सप्तक और मतलब में quadl और integrate पाया जाता है। [10][11]
Gautschi, Walter (1970). "On the construction of Gaussian quadrature rules from modified moments". Math. Comp. Vol. 24. pp. 245–260. doi:10.1090/S0025-5718-1970-0285117-6. MR0285177.
Piessens, R. (1971). "Gaussian quadrature formulas for the numerical integration of Bromwich's integral and the inversion of the laplace transform". J. Eng. Math. Vol. 5, no. 1. pp. 1–9. Bibcode:1971JEnMa...5....1P. doi:10.1007/BF01535429.
Danloy, Bernard (1973). "Numerical construction of Gaussian quadrature formulas for and ". Math. Comp. Vol. 27, no. 124. pp. 861–869. doi:10.1090/S0025-5718-1973-0331730-X. MR0331730.
Riener, Cordian; Schweighofer, Markus (2018). "Optimization approaches to quadrature: New characterizations of Gaussian quadrature on the line and quadrature with few nodes on plane algebraic curves, on the plane and in higher dimensions". Journal of Complexity. 45: 22–54. arXiv:1607.08404. doi:10.1016/j.jco.2017.10.002.
Stoer, Josef; Bulirsch, Roland (2002), Introduction to Numerical Analysis (3rd ed.), Springer, ISBN978-0-387-95452-3.
Walter Gautschi: "A Software Repository for Gaussian Quadratures and Christoffel Functions", SIAM, ISBN978-1-611976-34-2 (2020).
Teresa Laudadio, Nicola Mastronardi & Paul Van Dooren: "Computing Gaussian quadrature rules with high relative accuracy", Numerical Algorithms, vol.92(2023), pp.767–793.
Specific
↑Methodus nova integralium valores per approximationem inveniendi. In: Comm. Soc. Sci. Göttingen Math. Band 3, 1815, S. 29–76, Gallica, datiert 1814, auch in Werke, Band 3, 1876, S. 163–196.
↑C. G. J. Jacobi: Ueber Gauß' neue Methode, die Werthe der Integrale näherungsweise zu finden. In: Journal für Reine und Angewandte Mathematik. Band 1, 1826, S. 301–308, (online), und Werke, Band 6.