क्यूआर एल्गोरिदम: Difference between revisions

From Vigyanwiki
No edit summary
No edit summary
 
(6 intermediate revisions by 3 users not shown)
Line 1: Line 1:
{{Short description|Algorithm to calculate eigenvalues}}  
{{Short description|Algorithm to calculate eigenvalues}}  


[[संख्यात्मक रैखिक बीजगणित]] में, '''क्युआर एल्गोरिथ्म''' या क्युआर पुनरावृत्ति [[eigenvalue एल्गोरिथ्म|आइजेनवैल्यू एल्गोरिथ्म]] है: अर्थात, [[मैट्रिक्स (गणित)|आव्युह (गणित)]] के आइजेनवैल्यू ​​और आइजेनवेक्टर  की गणना करने की प्रक्रिया। क्यूआर एल्गोरिदम को 1950 के दशक के अंत में जॉन जी.एफ. फ्रांसिस और वेरा एन. कुब्लानोव्स्काया द्वारा स्वतंत्र रूप से काम करते हुए विकसित किया गया था।<ref>J.G.F. Francis, "The QR Transformation, I", ''[[The Computer Journal]]'', '''4'''(3), pages 265&ndash;271 (1961, received October 1959). [[doi:10.1093/comjnl/4.3.265]]</ref><ref>{{cite journal |first=J. G. F. |last=Francis |title=क्यूआर परिवर्तन, II|journal=The Computer Journal |volume=4 |issue=4 |pages=332–345 |year=1962 |doi=10.1093/comjnl/4.4.332 |doi-access=free }}</ref><ref>
[[संख्यात्मक रैखिक बीजगणित]] में, '''क्युआर एल्गोरिथ्म''' या क्युआर पुनरावृत्ति [[eigenvalue एल्गोरिथ्म|आइजेनवैल्यू एल्गोरिथ्म]] है: अर्थात, [[मैट्रिक्स (गणित)|आव्युह (गणित)]] के आइजेनवैल्यू ​​और आइजेनसदिश की गणना करने की प्रक्रिया होती है । तथा क्यूआर एल्गोरिदम को 1950 के दशक के अंत में जॉन जी.एफ. फ्रांसिस और वेरा एन. कुब्लानोव्स्काया द्वारा स्वतंत्र रूप से काम करते हुए विकसित किया गया था।<ref>J.G.F. Francis, "The QR Transformation, I", ''[[The Computer Journal]]'', '''4'''(3), pages 265&ndash;271 (1961, received October 1959). [[doi:10.1093/comjnl/4.3.265]]</ref><ref>{{cite journal |first=J. G. F. |last=Francis |title=क्यूआर परिवर्तन, II|journal=The Computer Journal |volume=4 |issue=4 |pages=332–345 |year=1962 |doi=10.1093/comjnl/4.4.332 |doi-access=free }}</ref><ref>
Vera N. Kublanovskaya, "On some algorithms for the solution of the complete eigenvalue problem," ''USSR Computational Mathematics and Mathematical Physics'', vol. 1, no. 3, pages 637–657 (1963, received Feb 1961). Also published in: ''Zhurnal Vychislitel'noi Matematiki i Matematicheskoi Fiziki'', vol.1, no. 4, pages 555–570 (1961). [[doi:10.1016/0041-5553(63)90168-X]]</ref> मूल विचार [[क्यूआर अपघटन]] करना है, आव्युह को [[ऑर्थोगोनल मैट्रिक्स|ऑर्थोगोनल आव्युह]] और ऊपरी [[त्रिकोणीय मैट्रिक्स|त्रिकोणीय आव्युह]] के उत्पाद के रूप में लिखना, कारकों को रिवर्स ऑर्डर में गुणा करना और पुनरावृत्त करना है।
Vera N. Kublanovskaya, "On some algorithms for the solution of the complete eigenvalue problem," ''USSR Computational Mathematics and Mathematical Physics'', vol. 1, no. 3, pages 637–657 (1963, received Feb 1961). Also published in: ''Zhurnal Vychislitel'noi Matematiki i Matematicheskoi Fiziki'', vol.1, no. 4, pages 555–570 (1961). [[doi:10.1016/0041-5553(63)90168-X]]</ref> मूल विचार [[क्यूआर अपघटन]] करना है, आव्युह को [[ऑर्थोगोनल मैट्रिक्स|ऑर्थोगोनल आव्युह]] और ऊपरी [[त्रिकोणीय मैट्रिक्स|त्रिकोणीय आव्युह]] के उत्पाद के रूप में लिखना, कारकों को रिवर्स ऑर्डर में गुणा करना और पुनरावृत्त करना है।  


==व्यावहारिक क्युआर एल्गोरिथ्म==
==व्यावहारिक क्युआर एल्गोरिथ्म                                                                                       ==
औपचारिक रूप से, मान लीजिए कि A वास्तविक आव्युह है जिसके आइजेनवैल्यू ​​​​की गणना हम करना चाहते हैं, और A को मान लीजिए<sub>0</sub>:=. K-वें चरण पर (k = 0 से प्रारंभ करके), हम क्युआर अपघटन A की गणना करते हैं<sub>''k''</sub>=प्र<sub>''k''</sub>R<sub>''k''</sub> कहां प्र<sub>''k''</sub> ऑर्थोगोनल आव्युह है (यानी, Q<sup>टी</sup> = क्यू<sup>−1</sup>) और आर<sub>''k''</sub> ऊपरी त्रिकोणीय आव्युह है. फिर हम A बनाते हैं<sub>''k''+1</sub> = आर<sub>''k''</sub>Q<sub>''k''</sub>. ध्यान दें कि
सामान्य रूप से, मान लीजिए कि A वास्तविक आव्युह है जिसके आइजेनवैल्यू ​की गणना हम करना चाहते हैं, और मान लीजिए ''A<sub>0</sub>:=A. K''-वें चरण पर (k = 0 से प्रारंभ करके), हम क्युआर अपघटन ''A<sub>k</sub>=Q<sub>k</sub>R<sub>k</sub>'' की गणना करते हैं जहाँ ''Q<sub>k</sub>'' ऑर्थोगोनल आव्युह है (अर्थात, ''Q<sup>T</sup> = Q<sup>−1</sup>'') और ''R<sub>k</sub>'' ऊपरी त्रिकोणीय आव्युह है. फिर हम ''A<sub>k+1</sub> = R<sub>k</sub>Q<sub>k</sub>''. बनाते हैं ध्यान दें कि  
<math display="block"> A_{k+1} = R_k Q_k = Q_k^{-1} Q_k R_k Q_k = Q_k^{-1} A_k Q_k = Q_k^{\mathsf{T}} A_k Q_k, </math>
<math display="block"> A_{k+1} = R_k Q_k = Q_k^{-1} Q_k R_k Q_k = Q_k^{-1} A_k Q_k = Q_k^{\mathsf{T}} A_k Q_k, </math>
तो सभी ए<sub>''k''</sub> [[समान मैट्रिक्स|समान आव्युह]]  हैं और इसलिए उनके आइजेनवैल्यू  ​​​​समान हैं। एल्गोरिथ्म [[संख्यात्मक स्थिरता]] है क्योंकि यह ऑर्थोगोनल समानता परिवर्तनों द्वारा आगे बढ़ता है।


खास शर्तों के अन्तर्गत,<ref name="golubvanloan">{{cite book |last1=Golub |first1=G. H. |last2=Van Loan |first2=C. F. |title=मैट्रिक्स संगणना|edition=3rd |publisher=Johns Hopkins University Press |location=Baltimore |year=1996 |isbn=0-8018-5414-8 }}</ref> आव्युह  ए<sub>''k''</sub> त्रिकोणीय आव्युह  में अभिसरण करें, ए का [[शूर रूप]]। त्रिकोणीय आव्युह  के आइजेनवैल्यू  ​​​​विकर्ण पर सूचीबद्ध हैं, और आइजेनवैल्यू  समस्या हल हो गई है। अभिसरण के परीक्षण में सटीक शून्य की आवश्यकता अव्यावहारिक है, लेकिन [[गेर्शगोरिन सर्कल प्रमेय]] त्रुटि पर सीमा प्रदान करता है।


इस कच्चे रूप में पुनरावृत्तियाँ अपेक्षाकृत महंगी हैं। इसे पहले आव्युह  ए को ऊपरी [[हेसेनबर्ग फॉर्म|हेसेनबर्ग रूप]]  में लाकर कम किया जा सकता है (जिसकी लागत है <math>\begin{matrix}\frac{10}{3}\end{matrix} n^3 + \mathcal{O}(n^2)</math> घरेलू परिवर्तन पर आधारित तकनीक का उपयोग करके अंकगणितीय संचालन), ऑर्थोगोनल समानता परिवर्तनों के सीमित अनुक्रम के साथ, कुछ हद तक दो-तरफा क्यूआर अपघटन की तरह।<ref name=Demmel>{{cite book |first=James W. |last=Demmel |author-link=James W. Demmel |title=अनुप्रयुक्त संख्यात्मक रैखिक बीजगणित|publisher=SIAM |year=1997 }}</ref><ref name=Trefethen>{{cite book |first1=Lloyd N. |last1=Trefethen |author-link=Lloyd N. Trefethen |first2=David |last2=Bau |title=संख्यात्मक रैखिक बीजगणित|publisher=SIAM |year=1997 }}</ref> (क्यूआर अपघटन के लिए, हाउसहोल्डर रिफ्लेक्टर को केवल बाईं ओर गुणा किया जाता है, लेकिन हेसेनबर्ग मामले के लिए उन्हें बाएं और दाएं दोनों पर गुणा किया जाता है।) ऊपरी हेसेनबर्ग आव्युह  लागत के क्यूआर अपघटन का निर्धारण <math>6 n^2 + \mathcal{O}(n)</math> अंकगणितीय आपरेशनस। इसके अलावा, क्योंकि हेसेनबर्ग रूप  पहले से ही लगभग ऊपरी-त्रिकोणीय है (इसमें प्रत्येक विकर्ण के नीचे केवल गैर-शून्य प्रविष्टि है), इसे शुरुआती बिंदु के रूप में उपयोग करने से क्यूआर एल्गोरिदम के अभिसरण के लिए आवश्यक चरणों की संख्या कम हो जाती है।


यदि मूल आव्युह  [[सममित मैट्रिक्स|सममित आव्युह]]  है, तो ऊपरी हेसेनबर्ग आव्युह  भी सममित है और इस प्रकार त्रिविकर्ण आव्युह  है, और सभी ए भी हैं<sub>''k''</sub>. इस प्रक्रिया में लागत आती है <math>\begin{matrix}\frac{4}{3}\end{matrix} n^3 + \mathcal{O}(n^2)</math> हाउसहोल्डर रिडक्शन पर आधारित तकनीक का उपयोग करके अंकगणितीय परिचालन।<ref name=Demmel/><ref name=Trefethen/>एक सममित त्रिविकर्ण आव्युह  लागत के क्यूआर अपघटन का निर्धारण <math>\mathcal{O}(n)</math> परिचालन.<ref>{{cite journal |first1=James M. |last1=Ortega |first2=Henry F. |last2=Kaiser |title=The ''LL<sup>T</sup>'' and ''QR'' methods for symmetric tridiagonal matrices |journal=The Computer Journal |volume=6 |issue=1 |pages=99–101 |year=1963 |doi=10.1093/comjnl/6.1.99 |doi-access=free }}</ref>
तो सभी ''A<sub>k</sub>'' [[समान मैट्रिक्स|समान आव्युह]] हैं और इसलिए उनके आइजेनवैल्यू ​​​​समान हैं। एल्गोरिथ्म [[संख्यात्मक स्थिरता]] है क्योंकि यह ऑर्थोगोनल समानता परिवर्तनों द्वारा आगे बढ़ता है।
अभिसरण की दर आइजेनवैल्यू  ​​​​के बीच अलगाव पर निर्भर करती है, इसलिए व्यावहारिक एल्गोरिदम अलगाव को बढ़ाने और अभिसरण में तेजी लाने के लिए स्पष्ट या अंतर्निहित बदलावों का उपयोग करेगा। विशिष्ट सममित क्यूआर एल्गोरिदम केवल या दो पुनरावृत्तियों के साथ प्रत्येक आइजेनवैल्यू को अलग करता है (फिर आव्युह  के आकार को कम करता है), जिससे यह कुशल और मजबूत हो जाता है।


== विज़ुअलाइज़ेशन ==
कुछ नियमों के अन्तर्गत,<ref name="golubvanloan">{{cite book |last1=Golub |first1=G. H. |last2=Van Loan |first2=C. F. |title=मैट्रिक्स संगणना|edition=3rd |publisher=Johns Hopkins University Press |location=Baltimore |year=1996 |isbn=0-8018-5414-8 }}</ref> आव्युह ''A<sub>k</sub>'' त्रिकोणीय आव्युह में परिवर्तित हो जाते हैं जहाँ A का [[शूर रूप|स्चुर रूप]] दर्शाता है । त्रिकोणीय आव्युह के आइजेनवैल्यू ​​​​विकर्ण पर सूचीबद्ध किया जाता हैं, और आइजेनवैल्यू समस्या हल हो जाती है। अभिसरण के परीक्षण में स्पष्ट शून्य की आवश्यकता अव्यावहारिक है, किन्तु [[गेर्शगोरिन सर्कल प्रमेय]] त्रुटि पर सीमा प्रदान करता है।
[[File:QR and LR visualization illustrating fixed points (corrected).gif|thumb|चित्र 1: क्यूआर या एलआर एल्गोरिथ्म के एकल पुनरावृत्ति का आउटपुट उसके इनपुट के साथ कैसे भिन्न होता है]]मूल क्यूआर एल्गोरिदम की कल्पना उस स्थिति में की जा सकती है जहां ए सकारात्मक-निश्चित सममित आव्युह है। उस स्थिति में, A को 2 आयामों में दीर्घवृत्त या उच्च आयामों में दीर्घवृत्त के रूप में दर्शाया जा सकता है। एल्गोरिथम के इनपुट और एकल पुनरावृत्ति के बीच संबंध को चित्र 1 (एनीमेशन देखने के लिए क्लिक करें) के रूप में दर्शाया जा सकता है। ध्यान दें कि एलआर एल्गोरिदम को क्यूआर एल्गोरिदम के साथ दर्शाया गया है।


एक एकल पुनरावृत्ति के कारण दीर्घवृत्त x-अक्ष की ओर झुक जाता है या गिर जाता है। ऐसी स्थिति में जहां दीर्घवृत्त का बड़ा [[अर्ध-प्रमुख और अर्ध-लघु अक्ष]] | अर्ध-अक्ष x-अक्ष के समानांतर है, क्युआर  का पुनरावृत्ति कुछ नहीं करता है। और स्थिति जहां एल्गोरिदम कुछ नहीं करता है वह यह है कि जब बड़ा अर्ध-अक्ष x-अक्ष के बजाय y-अक्ष के समानांतर होता है। उस घटना में, दीर्घवृत्त को किसी भी दिशा में गिरने में सक्षम हुए बिना अनिश्चित रूप से संतुलन बनाने के रूप में सोचा जा सकता है। दोनों स्थितियों में, आव्युह  विकर्ण है। ऐसी स्थिति जहां एल्गोरिथम की पुनरावृत्ति कुछ नहीं करती, उसे [[निश्चित बिंदु (गणित)]] कहा जाता है। एल्गोरिथम द्वारा नियोजित रणनीति [[निश्चित-बिंदु पुनरावृत्ति]]|एक निश्चित-बिंदु की ओर पुनरावृत्ति है। ध्यान दें कि निश्चित बिंदु स्थिर है जबकि दूसरा अस्थिर है। यदि दीर्घवृत्त को अस्थिर निश्चित बिंदु से बहुत कम मात्रा में झुकाया जाता है, तो क्यूआर के पुनरावृत्ति के कारण दीर्घवृत्त निश्चित बिंदु की ओर झुकने के बजाय दूर झुक जाएगा। चूँकि  अंततः, एल्गोरिदम अलग निश्चित बिंदु पर परिवर्तित हो जाएगा, लेकिन इसमें लंबा समय लगेगा।
इस क्रूड रूप में पुनरावृत्तियाँ अपेक्षाकृत मूल्यवान हैं। इसे पहले आव्युह A को ऊपरी [[हेसेनबर्ग फॉर्म|हेसेनबर्ग रूप]] में लाकर कम किया जा सकता है (जिसमें हाउसहोल्डर रिडक्शन पर आधारित विधि का उपयोग करके अंकगणितीय संचालन की निवेश <math>\begin{matrix}\frac{10}{3}\end{matrix} n^3 + \mathcal{O}(n^2)</math> होती है), तथा ऑर्थोगोनल समानता परिवर्तनों के सीमित अनुक्रम के साथ, कुछ सीमा तक दो-तरफा क्यूआर अपघटन की तरह।<ref name="Demmel">{{cite book |first=James W. |last=Demmel |author-link=James W. Demmel |title=अनुप्रयुक्त संख्यात्मक रैखिक बीजगणित|publisher=SIAM |year=1997 }}</ref><ref name="Trefethen">{{cite book |first1=Lloyd N. |last1=Trefethen |author-link=Lloyd N. Trefethen |first2=David |last2=Bau |title=संख्यात्मक रैखिक बीजगणित|publisher=SIAM |year=1997 }}</ref> (क्यूआर अपघटन के लिए, हाउसहोल्डर रिफ्लेक्टर को केवल बाईं ओर गुणा किया जाता है, किन्तु हेसेनबर्ग स्तिथियाँ के लिए उन्हें बाएं और दाएं दोनों पर गुणा किया जाता है।) ऊपरी हेसेनबर्ग आव्युह निवेश के क्यूआर अपघटन का निर्धारण करने में <math>6 n^2 + \mathcal{O}(n)</math> अंकगणितीय संचालन की निवेश होती है।। इसके अतिरिक्त , क्योंकि हेसेनबर्ग रूप पहले से ही लगभग ऊपरी-त्रिकोणीय है (इसमें प्रत्येक विकर्ण के नीचे केवल गैर-शून्य प्रविष्टि है), इसे प्रारम्भिक बिंदु के रूप में उपयोग करने से क्यूआर एल्गोरिदम के अभिसरण के लिए आवश्यक चरणों की संख्या कम हो जाती है।


=== आइजनवैल्यू ढूंढना बनाम आइजेनवेक्टर ढूंढना ===
यदि मूल आव्युह [[सममित मैट्रिक्स|सममित आव्युह]] है, तो ऊपरी हेसेनबर्ग आव्युह भी सममित है और इस प्रकार त्रिविकर्ण आव्युह है, और इसी तरह सभी A<sub>''k''</sub> भी हैं. इस प्रक्रिया में हाउसहोल्डर रिडक्शन पर आधारित तकनीक का उपयोग करके<math>\begin{matrix}\frac{4}{3}\end{matrix} n^3 + \mathcal{O}(n^2)</math> अंकगणितीय परिचालन की निवेश आती है।।<ref name="Demmel" /><ref name="Trefethen" /> सममित त्रिविकर्ण आव्युह निवेश के क्यूआर अपघटन का निर्धारण <math>\mathcal{O}(n)</math> संचालन की निवेश आती है।.<ref>{{cite journal |first1=James M. |last1=Ortega |first2=Henry F. |last2=Kaiser |title=The ''LL<sup>T</sup>'' and ''QR'' methods for symmetric tridiagonal matrices |journal=The Computer Journal |volume=6 |issue=1 |pages=99–101 |year=1963 |doi=10.1093/comjnl/6.1.99 |doi-access=free }}</ref>  
[[File:Qr lr eigenvalue clash.gif|thumb|चित्र 2: जब दो आइजेनवैल्यू  ​​​​एक दूसरे के पास आते हैं तो क्युआर  या LR के एकल पुनरावृत्ति का आउटपुट कैसे प्रभावित होता है]]यह इंगित करने योग्य है कि सममित आव्युह का भी आइजनवेक्टर ढूंढना गणना योग्य नहीं है ([[गणना योग्य विश्लेषण]] में परिभाषाओं के अनुसार सटीक वास्तविक अंकगणित में)<ref>{{Cite web|title=linear algebra - Why is uncomputability of the spectral decomposition not a problem?|url=https://mathoverflow.net/questions/369930/why-is-uncomputability-of-the-spectral-decomposition-not-a-problem|access-date=2021-08-09|website=MathOverflow}}</ref> यह कठिनाई तब मौजूद होती है जब किसी आव्युह  के आइजेनवैल्यू  ​​​​की बहुलताएं जानने योग्य नहीं होती हैं। दूसरी ओर, आइजेनवैल्यू  ​​​​खोजने के लिए वही समस्या मौजूद नहीं है। आव्युह  के आइजेनवैल्यू  ​​​​हमेशा गणना योग्य होते हैं।                                   


अब हम चर्चा करेंगे कि बुनियादी क्यूआर एल्गोरिदम में ये कठिनाइयाँ कैसे प्रकट होती हैं। इसे चित्र 2 में दर्शाया गया है। (थंबनेल पर क्लिक करना याद रखें)। याद रखें कि दीर्घवृत्त सकारात्मक-निश्चित सममित आव्युह  का प्रतिनिधित्व करते हैं। जैसे ही इनपुट आव्युह  के दो आइगेनवैल्यू एक-दूसरे के करीब आते हैं, इनपुट दीर्घवृत्त सर्कल में बदल जाता है। वृत्त पहचान आव्युह  के गुणज से मेल खाता है। निकट-वृत्त पहचान आव्युह  के निकट-गुणक से मेल खाता है जिसका आइजेनवैल्यू  ​​​​आव्युह  की विकर्ण प्रविष्टियों के लगभग बराबर है। इसलिए उस मामले में लगभग आइजेनवैल्यू  ​​​​खोजने की समस्या आसान दिखाई देती है। लेकिन ध्यान दें कि दीर्घवृत्त के अर्ध-अक्षों का क्या होता है। क्यूआर (या एलआर) की पुनरावृत्ति अर्ध-अक्षों को कम से कम झुकाती है क्योंकि इनपुट दीर्घवृत्त वृत्त होने के करीब पहुंचता है। आइजनवेक्टर केवल तभी ज्ञात हो सकते हैं जब अर्ध-अक्ष x-अक्ष और y-अक्ष के समानांतर हों। निकट-समानांतरता प्राप्त करने के लिए आवश्यक पुनरावृत्तियों की संख्या बिना किसी सीमा के बढ़ जाती है क्योंकि इनपुट दीर्घवृत्त अधिक गोलाकार हो जाता है।
अभिसरण की दर आइजेनवैल्यू ​​​​के बीच सेपरेशन पर निर्भर करती है, इसलिए व्यावहारिक एल्गोरिदम सेपरेशन को बढ़ाने और अभिसरण में तेजी लाने के लिए स्पष्ट या अंतर्निहित शिफ्टो का उपयोग करेगा। विशिष्ट सममित क्यूआर एल्गोरिदम केवल या दो पुनरावृत्तियों के साथ प्रत्येक आइजेनवैल्यू को भिन्न करता है (फिर आव्युह के आकार को कम करता है), जिससे यह कुशल और प्रबल हो जाता है।  


हालांकि मनमाना सममित आव्युह  के आव्युह  के ईजेंडेकंपोजिशन की गणना करना असंभव हो सकता है, आव्युह को इच्छानुसार  छोटी राशि से परेशान करना और परिणामी आव्युह  के ईजेंडेकंपोजीशन की गणना करना हमेशा संभव होता है। ऐसे मामले में जब आव्युह  को निकट-वृत्त के रूप में दर्शाया गया है, आव्युह  को उस आव्युह  से बदला जा सकता है जिसका चित्रण पूर्ण वृत्त है। उस स्थिति में, आव्युह  पहचान आव्युह  का गुणक है, और इसका आइजेनडीकम्पोजिसन तत्काल है। चूँकि  सावधान रहें कि परिणामी [[अपना आधार]] मूल आइजेनबासिस से काफी दूर हो सकता है।
== विज़ुअलाइज़ेशन                                                                          ==
[[File:QR and LR visualization illustrating fixed points (corrected).gif|thumb|चित्र 1: क्यूआर या एलआर एल्गोरिथ्म के एकल पुनरावृत्ति का आउटपुट उसके इनपुट के साथ कैसे भिन्न होता है]]मूल क्यूआर एल्गोरिदम की कल्पना उस स्थिति में की जा सकती है जहां A धनात्मक -निश्चित सममित आव्युह है। उस स्थिति में, A को 2 आयामों में दीर्घवृत्त या उच्च आयामों में दीर्घवृत्त के रूप में दर्शाया जा सकता है। एल्गोरिथम के इनपुट और एकल पुनरावृत्ति के बीच संबंध को चित्र 1 (एनीमेशन देखने के लिए क्लिक करें) के रूप में दर्शाया जा सकता है। ध्यान दें कि एलआर एल्गोरिदम को क्यूआर एल्गोरिदम के साथ दर्शाया गया है।


== अंतर्निहित क्यूआर एल्गोरिदम ==
इस प्रकार की एकल पुनरावृत्ति के कारण दीर्घवृत्त x-अक्ष की ओर झुक जाता है या गिर जाता है। तब ऐसी स्थिति में जहां दीर्घवृत्त का बड़ा [[अर्ध-प्रमुख और अर्ध-लघु अक्ष]] | अर्ध-अक्ष x-अक्ष के समानांतर होता है, क्युआर का पुनरावृत्ति कुछ नहीं करता है। और स्थिति जहां एल्गोरिदम कुछ नहीं करता है और इसका कारण यह है कि जब बड़ा अर्ध-अक्ष x-अक्ष के अतिरिक्त y-अक्ष के समानांतर होता है। उस घटना में, दीर्घवृत्त को किसी भी दिशा में गिरने में सक्षम हुए बिना अनिश्चित रूप से संतुलन बनाने के रूप में सोचा जा सकता है। तथा दोनों स्थितियों में, आव्युह विकर्ण है। ऐसी स्थिति जहां एल्गोरिथम की पुनरावृत्ति कुछ नहीं करती, उसे [[निश्चित बिंदु (गणित)]] कहा जाता है। और एल्गोरिथम द्वारा नियोजित रणनीति [[निश्चित-बिंदु पुनरावृत्ति]] या एक निश्चित-बिंदु की ओर पुनरावृत्ति है। ध्यान दें कि निश्चित बिंदु स्थिर है जबकि दूसरा अस्थिर है। यही कारण है कि दीर्घवृत्त को अस्थिर निश्चित बिंदु से बहुत कम मात्रा में झुकाया जाता है, तब क्यूआर के पुनरावृत्ति के कारण दीर्घवृत्त निश्चित बिंदु की ओर झुकने के अतिरिक्त दूर झुक जाता है । चूँकि अंततः एल्गोरिदम भिन्न दुसरे निश्चित बिंदु पर परिवर्तित हो तो जाएगा, किन्तु इसमें लंबा समय लग सकता है
आधुनिक कम्प्यूटेशनल अभ्यास में, क्यूआर एल्गोरिदम को अंतर्निहित संस्करण में निष्पादित किया जाता है जो कई बदलावों के उपयोग को प्रारंभ  करना आसान बनाता है।<ref name="golubvanloan" /> आव्युह  को पहले ऊपरी हेसेनबर्ग रूप  में लाया जाता है <math>A_0=QAQ^{\mathsf{T}}</math> जैसा कि स्पष्ट संस्करण में है; फिर, प्रत्येक चरण पर, का पहला कॉलम <math>A_k</math> के पहले कॉलम में छोटे आकार के घरेलू समानता परिवर्तन के माध्यम से रूपांतरित किया गया है <math>p(A_k)</math>  (या <math>p(A_k)e_1</math>), कहाँ <math>p(A_k)</math>, डिग्री का <math>r</math>, वह बहुपद है जो स्थानांतरण रणनीति को परिभाषित करता है (अक्सर <math>p(x)=(x-\lambda)(x-\bar\lambda)</math>, कहाँ <math>\lambda</math> और <math>\bar\lambda</math> अनुगामी के दो आइजेनवैल्यू  ​​हैं <math>2 \times 2</math> का प्रमुख सबआव्युह  <math>A_k</math>, तथाकथित अंतर्निहित डबल-शिफ्ट)। फिर आकार का क्रमिक गृहस्वामी परिवर्तन <math>r+1</math> कार्यशील आव्युह  को वापस करने के लिए किया जाता है <math>A_k</math> ऊपरी हेसेनबर्ग रूप में। एल्गोरिदम के चरणों के साथ आव्युह की गैर-शून्य प्रविष्टियों के अजीब आकार के कारण, इस ऑपरेशन को उभार पीछा के रूप में जाना जाता है। पहले संस्करण की तरह, उप-विकर्ण प्रविष्टियों में से के रूप में ही अपस्फीति का प्रदर्शन किया जाता है <math>A_k</math> पर्याप्त रूप से छोटा है.


===नाम बदलने का प्रस्ताव                                      ===
=== आइजनवैल्यू खोजना बनाम आइजेनसदिश खोजना ===
चूंकि प्रक्रिया के आधुनिक अंतर्निहित संस्करण में कोई क्यूआर अपघटन स्पष्ट रूप से नहीं किया जाता है, कुछ लेखक, उदाहरण के लिए वॉटकिंस,<ref>{{cite book |last=Watkins |first=David S. |title=The Matrix Eigenvalue Problem: GR and Krylov Subspace Methods |publisher=SIAM |location=Philadelphia, PA |year=2007 |isbn=978-0-89871-641-2 }}</ref> इसका नाम बदलकर फ्रांसिस एल्गोरिथम रखने का सुझाव दिया। जीन एच. गोलूब और चार्ल्स एफ. वैन लोन फ्रांसिस क्यूआर स्टेप शब्द का उपयोग करते हैं।
[[File:Qr lr eigenvalue clash.gif|thumb|चित्र 2: जब दो आइजेनवैल्यू ​​​​एक दूसरे के पास आते हैं तो क्युआर या LR के एकल पुनरावृत्ति का आउटपुट कैसे प्रभावित होता है]]यह संकेतिक करने योग्य है कि सममित आव्युह का भी आइजनसदिश खोजना गणना योग्य नहीं है अर्थात ([[गणना योग्य विश्लेषण]] में परिभाषाओं के अनुसार स्पष्ट वास्तविक अंकगणित में)।<ref>{{Cite web|title=linear algebra - Why is uncomputability of the spectral decomposition not a problem?|url=https://mathoverflow.net/questions/369930/why-is-uncomputability-of-the-spectral-decomposition-not-a-problem|access-date=2021-08-09|website=MathOverflow}}</ref> यह कठिनाई तब उपस्तिथ होती है जब किसी आव्युह के आइजेनवैल्यू की बहुलताएं जानने योग्य नहीं होती हैं। दूसरी ओर आइजेनवैल्यू ​​​​खोजने के लिए वही समस्या उपस्तिथ नहीं है। आव्युह के आइजेनवैल्यू ​​​​सदैव गणना योग्य होते हैं।                                    


== व्याख्या और अभिसरण                  ==
अब हम चर्चा करेंगे कि मूलभूत क्यूआर एल्गोरिदम में ये कठिनाइयाँ कैसे प्रकट होती हैं। इसे चित्र 2 में दर्शाया गया है। (थंबनेल पर क्लिक करना याद रखें)। याद रखें कि दीर्घवृत्त धनात्मक -निश्चित सममित आव्युह का प्रतिनिधित्व करते हैं। जैसे ही इनपुटआव्युह के दो आइगेनवैल्यू एक-दूसरे के समीप आते हैं, तब इनपुट दीर्घवृत्त सर्कल में बदल जाता है। वृत्त पहचान आव्युह के गुणज से मेल खाता है। निकट-वृत्त पहचान आव्युह के निकट-गुणक से मेल खाता है तथा जिसका आइजेनवैल्यू ​​​​आव्युह की विकर्ण प्रविष्टियों के लगभग समान्तर है। इसलिए उस स्तिथियाँ में लगभग आइजेनवैल्यू ​​​​खोजने की समस्या आसान दिखाई देती है। किन्तु ध्यान दें कि दीर्घवृत्त के अर्ध-अक्षों का क्या होता है। क्यूआर (या एलआर) की पुनरावृत्ति अर्ध-अक्षों को कम से कम झुकाती है क्योंकि इनपुट दीर्घवृत्त वृत्त होने के समीप पहुंचता है। आइजनसदिश केवल तभी ज्ञात हो सकते हैं जब अर्ध-अक्ष x-अक्ष और y-अक्ष के समानांतर हों। निकट-समानांतरता प्राप्त करने के लिए आवश्यक पुनरावृत्तियों की संख्या बिना किसी सीमा के बढ़ जाती है क्योंकि इनपुट दीर्घवृत्त अधिक गोलाकार हो जाता है।                                                                    
क्यूआर एल्गोरिदम को बुनियादी पावर पुनरावृत्ति के अधिक परिष्कृत बदलाव के रूप में देखा जा सकता है पावर आइजेनवैल्यू एल्गोरिदम। याद रखें कि पावर एल्गोरिदम बार-बार वेक्टर को ए से गुणा करता है, प्रत्येक पुनरावृत्ति के बाद सामान्य हो जाता है। वेक्टर सबसे बड़े आइजेनवैल्यू के आइजेनवेक्टर  में परिवर्तित हो जाता है। इसके बजाय, क्यूआर एल्गोरिदम वैक्टर के पूर्ण आधार के साथ काम करता है, क्यूआर अपघटन का उपयोग करके पुनर्सामान्यीकरण (और ऑर्थोगोनलाइज़) करता है। सममित आव्युह  A के लिए, अभिसरण पर, AQ = QΛ, जहां Λ आइजेनवैल्यू  ​​​​का विकर्ण आव्युह  है जिसमें A अभिसरण करता है, और जहां Q वहां पहुंचने के लिए आवश्यक सभी ऑर्थोगोनल समानता परिवर्तनों का संयोजन है। इस प्रकार Q के कॉलम आइजेनवेक्टर  हैं।


== इतिहास            ==
चूँकि इच्छा से सममित आव्युह के ईजेंडेकंपोजिशन की गणना करना असंभव हो सकता है, तथा आव्युह को इच्छानुसार छोटी राशि से परेशान करना और परिणामी आव्युह के ईजेंडेकंपोजीशन की गणना करना सदैव संभव होता है। ऐसे स्तिथियाँ में जब आव्युह को निकट-वृत्त के रूप में दर्शाया गया है, आव्युह को उस आव्युह से बदला जा सकता है जिसका चित्रण पूर्ण वृत्त है। उस स्थिति में, आव्युह पहचान आव्युह का गुणक है, और इसका ईजेंडेकंपोजिशन तत्काल है। चूँकि तब भी सावधान रहें कि परिणामी [[अपना आधार]] मूल आइजेनबासिस से अधिक दूर हो सकता है।
क्यूआर एल्गोरिदम एलआर एल्गोरिदम से पहले था, जो क्यूआर अपघटन के बजाय [[एलयू अपघटन]] का उपयोग करता है। क्यूआर एल्गोरिदम अधिक स्थिर है, इसलिए आजकल एलआर एल्गोरिदम का उपयोग शायद ही कभी किया जाता है। चूँकि , यह क्युआर  एल्गोरिदम के विकास में महत्वपूर्ण कदम का प्रतिनिधित्व करता है।


एलआर एल्गोरिथ्म को 1950 के दशक की शुरुआत में [[हेंज रूटीशौसर]] द्वारा विकसित किया गया था, जो उस समय [[ईटीएच ज्यूरिख]] में [[एडवर्ड बूट्स]] के अनुसंधान सहायक के रूप में काम करते थे। स्टिफ़ेल ने सुझाव दिया कि रूटीशौसर क्षणों y के अनुक्रम का उपयोग करें<sub>0</sub><sup>टी</sup><sup></sup>x<sub>0</sub>, k = 0, 1, … (जहाँ x<sub>0</sub> और य<sub>0</sub> इच्छानुसार  वेक्टर हैं) ए के आइजेनवैल्यू  ​​​​को खोजने के लिए। रुतिशौसर ने इस कार्य के लिए [[अलेक्जेंडर ऐटकेन]] का एल्गोरिदम लिया और इसे भागफल-अंतर एल्गोरिदम या क्यूडी एल्गोरिदम में विकसित किया। गणना को उपयुक्त आकार में व्यवस्थित करने के बाद, उन्होंने पाया कि क्यूडी एल्गोरिदम वास्तव में पुनरावृत्ति ए है<sub>''k''</sub> = एल<sub>''k''</sub>U<sub>''k''</sub> (एलयू अपघटन), ए<sub>''k''+1</sub> = यू<sub>''k''</sub>L<sub>''k''</sub>, त्रिविकर्ण आव्युह  पर लागू किया जाता है, जिसमें से एलआर एल्गोरिदम अनुसरण करता है।<ref>{{Citation | last1=Parlett | first1=Beresford N. | last2=Gutknecht | first2=Martin H. | title=From qd to LR, or, how were the qd and LR algorithms discovered? | doi=10.1093/imanum/drq003 | year=2011 | journal=IMA Journal of Numerical Analysis | issn=0272-4979 | volume=31 | issue=3 | pages=741–754| hdl=20.500.11850/159536 | url=http://doc.rero.ch/record/299139/files/drq003.pdf }}</ref>
== अंतर्निहित क्यूआर एल्गोरिदम                                    ==
आधुनिक कम्प्यूटेशनल अभ्यास में, क्यूआर एल्गोरिदम को अंतर्निहित संस्करण में निष्पादित किया जाता है जो कई शिफ्टो के उपयोग को प्रारंभ करना आसान बनाता है।<ref name="golubvanloan" /> आव्युह को पहले ऊपरी हेसेनबर्ग रूप <math>A_0=QAQ^{\mathsf{T}}                                                                                                                                    </math> में लाया जाता है जैसा कि स्पष्ट संस्करण में है; फिर, प्रत्येक चरण पर, <math>A_k</math> का पहला स्तम्भ के पहले स्तम्भ में छोटे आकार के घरेलू समानता परिवर्तन के माध्यम से <math>p(A_k)</math> (या <math>p(A_k)e_1                                                                                                                    </math>) रूपांतरित किया गया है, जहाँ <math>p(A_k)</math>, डिग्री का <math>r</math>, वह बहुपद है जो स्थानांतरण रणनीति <math>p(x)=(x-\lambda)(x-\bar\lambda)                                                                                                </math> को परिभाषित करता है (प्रायः, जहाँ <math>\lambda</math> और <math>\bar\lambda</math> ) अनुगामी 2 के <math>2 \times 2</math> आइजेनवैल्यू ​​हैं इसका <math>A_k</math> प्रमुख उपआव्यूह है , तथाकथित अंतर्निहित डबल-शिफ्ट) होती है । फिर आकार का क्रमिक गृहस्वामी परिवर्तन <math>r+1</math> कार्यशील आव्युह <math>A_k</math> को ऊपरी हेसेनबर्ग रूप में वापस करने के लिए किया जाता है । एल्गोरिदम के चरणों के साथ आव्युह की गैर-शून्य प्रविष्टियों के विचित्र आकार के कारण, इस ऑपरेशन को उभार पीछा के रूप में जाना जाता है। पहले संस्करण की तरह, उप-विकर्ण प्रविष्टियों में से के रूप में ही अपस्फीति का प्रदर्शन किया जाता है जहाँ <math>A_k</math> पर्याप्त रूप से छोटा है.


===नाम बदलने का प्रस्ताव                        ===
चूंकि प्रक्रिया के आधुनिक अंतर्निहित संस्करण में कोई क्यूआर अपघटन स्पष्ट रूप से नहीं किया जाता है, कुछ लेखक, उदाहरण के लिए वॉटकिंस,<ref>{{cite book |last=Watkins |first=David S. |title=The Matrix Eigenvalue Problem: GR and Krylov Subspace Methods |publisher=SIAM |location=Philadelphia, PA |year=2007 |isbn=978-0-89871-641-2 }}</ref> इसका नाम बदलकर फ्रांसिस एल्गोरिथम रखने का सुझाव दिया जाता है । जीन एच. गोलूब और चार्ल्स एफ. वैन लोन फ्रांसिस क्यूआर स्टेप शब्द का उपयोग करते हैं।
== व्याख्या और अभिसरण        ==
क्यूआर एल्गोरिदम को मूलभूत शक्ति पुनरावृत्ति के अधिक परिष्कृत परिवर्तन के रूप में देखा जा सकता है तथा शक्ति आइजेनवैल्यू एल्गोरिदम होती है। अर्थात यह याद रखें कि शक्ति एल्गोरिदम निरंतर सदिश को A से गुणा करता है, जब प्रत्येक पुनरावृत्ति के बाद सामान्य हो जाता है। तब सदिश सबसे बड़े आइजेनवैल्यू के आइजेनसदिश में परिवर्तित हो जाता है। इसके अतिरिक्त ,क्यूआर एल्गोरिदम सदिश के पूर्ण आधार के साथ कार्य करता है, क्यूआर अपघटन का उपयोग करके पुनर्सामान्यीकरण (और ऑर्थोगोनलाइज़) करता है। और सममित आव्युह A के लिए, अभिसरण पर, AQ = QΛ, जहां Λ आइजेनवैल्यू ​​​​का विकर्ण आव्युह है जिसमें A अभिसरण करता है, और जहां Q वहां पहुंचने के लिए आवश्यक सभी ऑर्थोगोनल समानता परिवर्तनों का संयोजन है। इस प्रकार Q के स्तम्भ आइजेनसदिश हैं।
== इतिहास                ==
क्यूआर एल्गोरिदम एलआर एल्गोरिदम से पहले था, जो क्यूआर अपघटन के अतिरिक्त [[एलयू अपघटन]] का उपयोग करता है। क्यूआर एल्गोरिदम अधिक स्थिर है, इसलिए आजकल एलआर एल्गोरिदम का उपयोग संभवता ही कभी किया जाता है। चूँकि यह क्युआर एल्गोरिदम के विकास में महत्वपूर्ण कदम का प्रतिनिधित्व करता है।             
एलआर एल्गोरिथ्म को 1950 के दशक की प्रारंभ में [[हेंज रूटीशौसर]] द्वारा विकसित किया गया था, जो उस समय [[ईटीएच ज्यूरिख]] में [[एडवर्ड बूट्स]] के अनुसंधान सहायक के रूप में कार्य करते थे। तथा स्टिफ़ेल ने सुझाव दिया कि रूटीशौसर A के आइजेनवैल्यू ​​​​को खोजने के लिए क्षणों के अनुक्रम ''y<sub>0</sub><sup>T</sup>A<sup>k</sup>x<sub>0</sub>, k = 0, 1'' , … (जहाँ x<sub>0</sub> और y<sub>0</sub> इच्छानुसार सदिश हैं) का उपयोग करते है। रुतिशौसर ने इस कार्य के लिए [[अलेक्जेंडर ऐटकेन]] का एल्गोरिदम लिया और इसे भागफल-अंतर एल्गोरिदम या क्यूडी एल्गोरिदम में विकसित किया। गणना को उपयुक्त आकार में व्यवस्थित करने के बाद, उन्होंने पाया कि क्यूडी एल्गोरिदम वास्तव में पुनरावृत्ति A<sub>''k''</sub> = L<sub>''k''</sub>U<sub>''k''</sub> (एलयू अपघटन), A<sub>''k''+1</sub> = U<sub>''k''</sub>L<sub>''k''</sub>, द्वारा होती है त्रिविकर्ण आव्युह पर प्रयुक्त किया जाता है, जिसमें से एलआर एल्गोरिदम अनुसरण करता है।<ref>{{Citation | last1=Parlett | first1=Beresford N. | last2=Gutknecht | first2=Martin H. | title=From qd to LR, or, how were the qd and LR algorithms discovered? | doi=10.1093/imanum/drq003 | year=2011 | journal=IMA Journal of Numerical Analysis | issn=0272-4979 | volume=31 | issue=3 | pages=741–754| hdl=20.500.11850/159536 | url=http://doc.rero.ch/record/299139/files/drq003.pdf }}</ref>
== अन्य प्रकार          ==
क्यूआर एल्गोरिथ्म का प्रकार, गोलूब-कहान-रीन्स्च एल्गोरिथ्म सामान्य आव्युह को द्विभुजीय आव्युह में कम करने के साथ प्रारंभ होता है।<ref>Bochkanov Sergey Anatolyevich. ALGLIB User Guide - General Matrix operations - Singular value decomposition . ALGLIB Project. 2010-12-11. URL:[http://www.alglib.net/matrixops/general/svd.php] Accessed: 2010-12-11. (Archived by WebCite at https://www.webcitation.org/5utO4iSnR?url=http://www.alglib.net/matrixops/general/svd.php</ref> [[एकवचन मान|एकवचन मानों]] की गणना के लिए क्युआर एल्गोरिदम के इस संस्करण का वर्णन सबसे पहले किसके द्वारा किया गया था? अर्थात {{harvtxt|गोलुब |कहन|1965}}. [[LAPACK|लापैक]] सबरूटीन [http://www.netlib.org/lapack/double/dbdsqr.f डीबीडीएसयूआर] के द्वारा किया गया था | तब उस स्तिथियाँ को कवर करने के लिए कुछ संशोधनों के साथ इस पुनरावृत्त विधि को कार्यान्वित करता है जहां एकवचन मान बहुत छोटे होते हैं {{harv|डेमेल|कहन|1990}}. हाउसहोल्डर रिफ्लेक्शन और, यदि उपयुक्त हो, क्यूआर अपघटन का उपयोग करते हुए पहले चरण के साथ, यह एकवचन मूल्य अपघटन की गणना के लिए [http://www.netlib.org/lapack/double/dgesvd.f डीजीईएसवीडी] रूटीन बनाता है। क्यूआर एल्गोरिदम को संबंधित अभिसरण परिणामों के साथ अनंत आयामों में भी प्रयुक्त किया जा सकता है।<ref>{{cite journal |first1=Percy |last1=Deift |first2=Luenchau C. |last2=Li |first3=Carlos |last3=Tomei|title=टोडा अनंत अनेक चरों के साथ बहती है|journal=Journal of Functional Analysis |volume=64 | number=3| pages=358–402 |year=1985 |doi=10.1016/0022-1236(85)90065-5|doi-access=free }}</ref><ref>{{cite journal |first1=Matthew J. |last1=Colbrook |first2=Anders C.|last2=Hansen|title=अनंत-आयामी क्यूआर एल्गोरिदम पर|journal=Numerische Mathematik |volume=143 | number=1| pages=17–83 |year=2019 |doi=10.1007/s00211-019-01047-5|arxiv=2011.08172|doi-access=free }}</ref>


== अन्य प्रकार                                                                  ==
क्यूआर एल्गोरिथ्म का प्रकार, गोलूब-कहान-रीन्स्च एल्गोरिथ्म सामान्य आव्युह  को द्विभुजीय आव्युह  में कम करने के साथ प्रारंभ  होता है।<ref>Bochkanov Sergey Anatolyevich. ALGLIB User Guide - General Matrix operations - Singular value decomposition . ALGLIB Project. 2010-12-11. URL:[http://www.alglib.net/matrixops/general/svd.php] Accessed: 2010-12-11. (Archived by WebCite at https://www.webcitation.org/5utO4iSnR?url=http://www.alglib.net/matrixops/general/svd.php</ref> [[एकवचन मान|एकवचन मानों]] की गणना के लिए क्युआर  एल्गोरिदम के इस संस्करण का वर्णन सबसे पहले किसके द्वारा किया गया था? {{harvtxt|गोलुब |कहन|1965}}. [[LAPACK|लापैक]] सबरूटीन [http://www.netlib.org/lapack/double/dbdsqr.f डीबीडीएसयूआर]  उस मामले को कवर करने के लिए कुछ संशोधनों के साथ इस पुनरावृत्त विधि को कार्यान्वित करता है जहां एकवचन मान बहुत छोटे होते हैं {{harv|डेमेल|कहन|1990}}. हाउसहोल्डर रिफ्लेक्शन और, यदि उपयुक्त हो, क्यूआर अपघटन का उपयोग करते हुए पहले चरण के साथ, यह एकवचन मूल्य अपघटन की गणना के लिए [http://www.netlib.org/lapack/double/dgesvd.f डीजीईएसवीडी] रूटीन बनाता है। क्यूआर एल्गोरिदम को संबंधित अभिसरण परिणामों के साथ अनंत आयामों में भी लागू किया जा सकता है।<ref>{{cite journal |first1=Percy |last1=Deift |first2=Luenchau C. |last2=Li |first3=Carlos |last3=Tomei|title=टोडा अनंत अनेक चरों के साथ बहती है|journal=Journal of Functional Analysis |volume=64 | number=3| pages=358–402 |year=1985 |doi=10.1016/0022-1236(85)90065-5|doi-access=free }}</ref><ref>{{cite journal |first1=Matthew J. |last1=Colbrook |first2=Anders C.|last2=Hansen|title=अनंत-आयामी क्यूआर एल्गोरिदम पर|journal=Numerische Mathematik |volume=143 | number=1| pages=17–83 |year=2019 |doi=10.1007/s00211-019-01047-5|arxiv=2011.08172|doi-access=free }}</ref>




Line 57: Line 61:
== बाहरी संबंध ==
== बाहरी संबंध ==
* {{PlanetMath |urlname=eigenvalueproblem |title=Eigenvalue problem}}
* {{PlanetMath |urlname=eigenvalueproblem |title=Eigenvalue problem}}
* [http://www-users.math.umn.edu/~olver/aims_/qr.pdf Notes on orthogonal bases and the workings of the क्युआर algorithm] by [[Peter J. Olver]]
* [http://www-users.math.umn.edu/~olver/aims_/qr.pdf Notes on orthogonal bases and the workings of the क्युआर algorithm] by [[Peter J. Olver]]
* [https://web.archive.org/web/20081209042103/http://math.fullerton.edu/mathews/n2003/QRMethodMod.html Module for the क्युआर Method]
* [https://web.archive.org/web/20081209042103/http://math.fullerton.edu/mathews/n2003/QRMethodMod.html Module for the क्युआर Method]
* [https://github.com/nom-de-guerre/Matrices C++ Library]
* [https://github.com/nom-de-guerre/Matrices C++ Library]


{{Numerical linear algebra}}
{{Numerical linear algebra}}


{{DEFAULTSORT:Qr Algorithm}}[[Category: संख्यात्मक रैखिक बीजगणित]]
{{DEFAULTSORT:Qr Algorithm}}
 
 


[[Category: Machine Translated Page]]
[[Category:Collapse templates|Qr Algorithm]]
[[Category:Created On 25/07/2023]]
[[Category:Created On 25/07/2023|Qr Algorithm]]
[[Category:Lua-based templates|Qr Algorithm]]
[[Category:Machine Translated Page|Qr Algorithm]]
[[Category:Navigational boxes| ]]
[[Category:Navigational boxes without horizontal lists|Qr Algorithm]]
[[Category:Pages with script errors|Qr Algorithm]]
[[Category:Short description with empty Wikidata description|Qr Algorithm]]
[[Category:Sidebars with styles needing conversion|Qr Algorithm]]
[[Category:Template documentation pages|Documentation/doc]]
[[Category:Templates Vigyan Ready|Qr Algorithm]]
[[Category:Templates generating microformats|Qr Algorithm]]
[[Category:Templates that add a tracking category|Qr Algorithm]]
[[Category:Templates that are not mobile friendly|Qr Algorithm]]
[[Category:Templates that generate short descriptions|Qr Algorithm]]
[[Category:Templates using TemplateData|Qr Algorithm]]
[[Category:Wikipedia metatemplates|Qr Algorithm]]
[[Category:संख्यात्मक रैखिक बीजगणित|Qr Algorithm]]

Latest revision as of 14:14, 11 August 2023

संख्यात्मक रैखिक बीजगणित में, क्युआर एल्गोरिथ्म या क्युआर पुनरावृत्ति आइजेनवैल्यू एल्गोरिथ्म है: अर्थात, आव्युह (गणित) के आइजेनवैल्यू ​​और आइजेनसदिश की गणना करने की प्रक्रिया होती है । तथा क्यूआर एल्गोरिदम को 1950 के दशक के अंत में जॉन जी.एफ. फ्रांसिस और वेरा एन. कुब्लानोव्स्काया द्वारा स्वतंत्र रूप से काम करते हुए विकसित किया गया था।[1][2][3] मूल विचार क्यूआर अपघटन करना है, आव्युह को ऑर्थोगोनल आव्युह और ऊपरी त्रिकोणीय आव्युह के उत्पाद के रूप में लिखना, कारकों को रिवर्स ऑर्डर में गुणा करना और पुनरावृत्त करना है।

व्यावहारिक क्युआर एल्गोरिथ्म

सामान्य रूप से, मान लीजिए कि A वास्तविक आव्युह है जिसके आइजेनवैल्यू ​की गणना हम करना चाहते हैं, और मान लीजिए A0:=A. K-वें चरण पर (k = 0 से प्रारंभ करके), हम क्युआर अपघटन Ak=QkRk की गणना करते हैं जहाँ Qk ऑर्थोगोनल आव्युह है (अर्थात, QT = Q−1) और Rk ऊपरी त्रिकोणीय आव्युह है. फिर हम Ak+1 = RkQk. बनाते हैं ध्यान दें कि


तो सभी Ak समान आव्युह हैं और इसलिए उनके आइजेनवैल्यू ​​​​समान हैं। एल्गोरिथ्म संख्यात्मक स्थिरता है क्योंकि यह ऑर्थोगोनल समानता परिवर्तनों द्वारा आगे बढ़ता है।

कुछ नियमों के अन्तर्गत,[4] आव्युह Ak त्रिकोणीय आव्युह में परिवर्तित हो जाते हैं जहाँ A का स्चुर रूप दर्शाता है । त्रिकोणीय आव्युह के आइजेनवैल्यू ​​​​विकर्ण पर सूचीबद्ध किया जाता हैं, और आइजेनवैल्यू समस्या हल हो जाती है। अभिसरण के परीक्षण में स्पष्ट शून्य की आवश्यकता अव्यावहारिक है, किन्तु गेर्शगोरिन सर्कल प्रमेय त्रुटि पर सीमा प्रदान करता है।

इस क्रूड रूप में पुनरावृत्तियाँ अपेक्षाकृत मूल्यवान हैं। इसे पहले आव्युह A को ऊपरी हेसेनबर्ग रूप में लाकर कम किया जा सकता है (जिसमें हाउसहोल्डर रिडक्शन पर आधारित विधि का उपयोग करके अंकगणितीय संचालन की निवेश होती है), तथा ऑर्थोगोनल समानता परिवर्तनों के सीमित अनुक्रम के साथ, कुछ सीमा तक दो-तरफा क्यूआर अपघटन की तरह।[5][6] (क्यूआर अपघटन के लिए, हाउसहोल्डर रिफ्लेक्टर को केवल बाईं ओर गुणा किया जाता है, किन्तु हेसेनबर्ग स्तिथियाँ के लिए उन्हें बाएं और दाएं दोनों पर गुणा किया जाता है।) ऊपरी हेसेनबर्ग आव्युह निवेश के क्यूआर अपघटन का निर्धारण करने में अंकगणितीय संचालन की निवेश होती है।। इसके अतिरिक्त , क्योंकि हेसेनबर्ग रूप पहले से ही लगभग ऊपरी-त्रिकोणीय है (इसमें प्रत्येक विकर्ण के नीचे केवल गैर-शून्य प्रविष्टि है), इसे प्रारम्भिक बिंदु के रूप में उपयोग करने से क्यूआर एल्गोरिदम के अभिसरण के लिए आवश्यक चरणों की संख्या कम हो जाती है।

यदि मूल आव्युह सममित आव्युह है, तो ऊपरी हेसेनबर्ग आव्युह भी सममित है और इस प्रकार त्रिविकर्ण आव्युह है, और इसी तरह सभी Ak भी हैं. इस प्रक्रिया में हाउसहोल्डर रिडक्शन पर आधारित तकनीक का उपयोग करके अंकगणितीय परिचालन की निवेश आती है।।[5][6] सममित त्रिविकर्ण आव्युह निवेश के क्यूआर अपघटन का निर्धारण संचालन की निवेश आती है।.[7]

अभिसरण की दर आइजेनवैल्यू ​​​​के बीच सेपरेशन पर निर्भर करती है, इसलिए व्यावहारिक एल्गोरिदम सेपरेशन को बढ़ाने और अभिसरण में तेजी लाने के लिए स्पष्ट या अंतर्निहित शिफ्टो का उपयोग करेगा। विशिष्ट सममित क्यूआर एल्गोरिदम केवल या दो पुनरावृत्तियों के साथ प्रत्येक आइजेनवैल्यू को भिन्न करता है (फिर आव्युह के आकार को कम करता है), जिससे यह कुशल और प्रबल हो जाता है।

विज़ुअलाइज़ेशन

चित्र 1: क्यूआर या एलआर एल्गोरिथ्म के एकल पुनरावृत्ति का आउटपुट उसके इनपुट के साथ कैसे भिन्न होता है

मूल क्यूआर एल्गोरिदम की कल्पना उस स्थिति में की जा सकती है जहां A धनात्मक -निश्चित सममित आव्युह है। उस स्थिति में, A को 2 आयामों में दीर्घवृत्त या उच्च आयामों में दीर्घवृत्त के रूप में दर्शाया जा सकता है। एल्गोरिथम के इनपुट और एकल पुनरावृत्ति के बीच संबंध को चित्र 1 (एनीमेशन देखने के लिए क्लिक करें) के रूप में दर्शाया जा सकता है। ध्यान दें कि एलआर एल्गोरिदम को क्यूआर एल्गोरिदम के साथ दर्शाया गया है।

इस प्रकार की एकल पुनरावृत्ति के कारण दीर्घवृत्त x-अक्ष की ओर झुक जाता है या गिर जाता है। तब ऐसी स्थिति में जहां दीर्घवृत्त का बड़ा अर्ध-प्रमुख और अर्ध-लघु अक्ष | अर्ध-अक्ष x-अक्ष के समानांतर होता है, क्युआर का पुनरावृत्ति कुछ नहीं करता है। और स्थिति जहां एल्गोरिदम कुछ नहीं करता है और इसका कारण यह है कि जब बड़ा अर्ध-अक्ष x-अक्ष के अतिरिक्त y-अक्ष के समानांतर होता है। उस घटना में, दीर्घवृत्त को किसी भी दिशा में गिरने में सक्षम हुए बिना अनिश्चित रूप से संतुलन बनाने के रूप में सोचा जा सकता है। तथा दोनों स्थितियों में, आव्युह विकर्ण है। ऐसी स्थिति जहां एल्गोरिथम की पुनरावृत्ति कुछ नहीं करती, उसे निश्चित बिंदु (गणित) कहा जाता है। और एल्गोरिथम द्वारा नियोजित रणनीति निश्चित-बिंदु पुनरावृत्ति या एक निश्चित-बिंदु की ओर पुनरावृत्ति है। ध्यान दें कि निश्चित बिंदु स्थिर है जबकि दूसरा अस्थिर है। यही कारण है कि दीर्घवृत्त को अस्थिर निश्चित बिंदु से बहुत कम मात्रा में झुकाया जाता है, तब क्यूआर के पुनरावृत्ति के कारण दीर्घवृत्त निश्चित बिंदु की ओर झुकने के अतिरिक्त दूर झुक जाता है । चूँकि अंततः एल्गोरिदम भिन्न दुसरे निश्चित बिंदु पर परिवर्तित हो तो जाएगा, किन्तु इसमें लंबा समय लग सकता है ।

आइजनवैल्यू खोजना बनाम आइजेनसदिश खोजना

चित्र 2: जब दो आइजेनवैल्यू ​​​​एक दूसरे के पास आते हैं तो क्युआर या LR के एकल पुनरावृत्ति का आउटपुट कैसे प्रभावित होता है

यह संकेतिक करने योग्य है कि सममित आव्युह का भी आइजनसदिश खोजना गणना योग्य नहीं है अर्थात (गणना योग्य विश्लेषण में परिभाषाओं के अनुसार स्पष्ट वास्तविक अंकगणित में)।[8] यह कठिनाई तब उपस्तिथ होती है जब किसी आव्युह के आइजेनवैल्यू की बहुलताएं जानने योग्य नहीं होती हैं। दूसरी ओर आइजेनवैल्यू ​​​​खोजने के लिए वही समस्या उपस्तिथ नहीं है। आव्युह के आइजेनवैल्यू ​​​​सदैव गणना योग्य होते हैं।

अब हम चर्चा करेंगे कि मूलभूत क्यूआर एल्गोरिदम में ये कठिनाइयाँ कैसे प्रकट होती हैं। इसे चित्र 2 में दर्शाया गया है। (थंबनेल पर क्लिक करना याद रखें)। याद रखें कि दीर्घवृत्त धनात्मक -निश्चित सममित आव्युह का प्रतिनिधित्व करते हैं। जैसे ही इनपुटआव्युह के दो आइगेनवैल्यू एक-दूसरे के समीप आते हैं, तब इनपुट दीर्घवृत्त सर्कल में बदल जाता है। वृत्त पहचान आव्युह के गुणज से मेल खाता है। निकट-वृत्त पहचान आव्युह के निकट-गुणक से मेल खाता है तथा जिसका आइजेनवैल्यू ​​​​आव्युह की विकर्ण प्रविष्टियों के लगभग समान्तर है। इसलिए उस स्तिथियाँ में लगभग आइजेनवैल्यू ​​​​खोजने की समस्या आसान दिखाई देती है। किन्तु ध्यान दें कि दीर्घवृत्त के अर्ध-अक्षों का क्या होता है। क्यूआर (या एलआर) की पुनरावृत्ति अर्ध-अक्षों को कम से कम झुकाती है क्योंकि इनपुट दीर्घवृत्त वृत्त होने के समीप पहुंचता है। आइजनसदिश केवल तभी ज्ञात हो सकते हैं जब अर्ध-अक्ष x-अक्ष और y-अक्ष के समानांतर हों। निकट-समानांतरता प्राप्त करने के लिए आवश्यक पुनरावृत्तियों की संख्या बिना किसी सीमा के बढ़ जाती है क्योंकि इनपुट दीर्घवृत्त अधिक गोलाकार हो जाता है।

चूँकि इच्छा से सममित आव्युह के ईजेंडेकंपोजिशन की गणना करना असंभव हो सकता है, तथा आव्युह को इच्छानुसार छोटी राशि से परेशान करना और परिणामी आव्युह के ईजेंडेकंपोजीशन की गणना करना सदैव संभव होता है। ऐसे स्तिथियाँ में जब आव्युह को निकट-वृत्त के रूप में दर्शाया गया है, आव्युह को उस आव्युह से बदला जा सकता है जिसका चित्रण पूर्ण वृत्त है। उस स्थिति में, आव्युह पहचान आव्युह का गुणक है, और इसका ईजेंडेकंपोजिशन तत्काल है। चूँकि तब भी सावधान रहें कि परिणामी अपना आधार मूल आइजेनबासिस से अधिक दूर हो सकता है।

अंतर्निहित क्यूआर एल्गोरिदम

आधुनिक कम्प्यूटेशनल अभ्यास में, क्यूआर एल्गोरिदम को अंतर्निहित संस्करण में निष्पादित किया जाता है जो कई शिफ्टो के उपयोग को प्रारंभ करना आसान बनाता है।[4] आव्युह को पहले ऊपरी हेसेनबर्ग रूप में लाया जाता है जैसा कि स्पष्ट संस्करण में है; फिर, प्रत्येक चरण पर, का पहला स्तम्भ के पहले स्तम्भ में छोटे आकार के घरेलू समानता परिवर्तन के माध्यम से (या ) रूपांतरित किया गया है, जहाँ , डिग्री का , वह बहुपद है जो स्थानांतरण रणनीति को परिभाषित करता है (प्रायः, जहाँ और ) अनुगामी 2 के आइजेनवैल्यू ​​हैं इसका प्रमुख उपआव्यूह है , तथाकथित अंतर्निहित डबल-शिफ्ट) होती है । फिर आकार का क्रमिक गृहस्वामी परिवर्तन कार्यशील आव्युह को ऊपरी हेसेनबर्ग रूप में वापस करने के लिए किया जाता है । एल्गोरिदम के चरणों के साथ आव्युह की गैर-शून्य प्रविष्टियों के विचित्र आकार के कारण, इस ऑपरेशन को उभार पीछा के रूप में जाना जाता है। पहले संस्करण की तरह, उप-विकर्ण प्रविष्टियों में से के रूप में ही अपस्फीति का प्रदर्शन किया जाता है जहाँ पर्याप्त रूप से छोटा है.

नाम बदलने का प्रस्ताव

चूंकि प्रक्रिया के आधुनिक अंतर्निहित संस्करण में कोई क्यूआर अपघटन स्पष्ट रूप से नहीं किया जाता है, कुछ लेखक, उदाहरण के लिए वॉटकिंस,[9] इसका नाम बदलकर फ्रांसिस एल्गोरिथम रखने का सुझाव दिया जाता है । जीन एच. गोलूब और चार्ल्स एफ. वैन लोन फ्रांसिस क्यूआर स्टेप शब्द का उपयोग करते हैं।

व्याख्या और अभिसरण

क्यूआर एल्गोरिदम को मूलभूत शक्ति पुनरावृत्ति के अधिक परिष्कृत परिवर्तन के रूप में देखा जा सकता है तथा शक्ति आइजेनवैल्यू एल्गोरिदम होती है। अर्थात यह याद रखें कि शक्ति एल्गोरिदम निरंतर सदिश को A से गुणा करता है, जब प्रत्येक पुनरावृत्ति के बाद सामान्य हो जाता है। तब सदिश सबसे बड़े आइजेनवैल्यू के आइजेनसदिश में परिवर्तित हो जाता है। इसके अतिरिक्त ,क्यूआर एल्गोरिदम सदिश के पूर्ण आधार के साथ कार्य करता है, क्यूआर अपघटन का उपयोग करके पुनर्सामान्यीकरण (और ऑर्थोगोनलाइज़) करता है। और सममित आव्युह A के लिए, अभिसरण पर, AQ = QΛ, जहां Λ आइजेनवैल्यू ​​​​का विकर्ण आव्युह है जिसमें A अभिसरण करता है, और जहां Q वहां पहुंचने के लिए आवश्यक सभी ऑर्थोगोनल समानता परिवर्तनों का संयोजन है। इस प्रकार Q के स्तम्भ आइजेनसदिश हैं।

इतिहास

क्यूआर एल्गोरिदम एलआर एल्गोरिदम से पहले था, जो क्यूआर अपघटन के अतिरिक्त एलयू अपघटन का उपयोग करता है। क्यूआर एल्गोरिदम अधिक स्थिर है, इसलिए आजकल एलआर एल्गोरिदम का उपयोग संभवता ही कभी किया जाता है। चूँकि यह क्युआर एल्गोरिदम के विकास में महत्वपूर्ण कदम का प्रतिनिधित्व करता है।

एलआर एल्गोरिथ्म को 1950 के दशक की प्रारंभ में हेंज रूटीशौसर द्वारा विकसित किया गया था, जो उस समय ईटीएच ज्यूरिख में एडवर्ड बूट्स के अनुसंधान सहायक के रूप में कार्य करते थे। तथा स्टिफ़ेल ने सुझाव दिया कि रूटीशौसर A के आइजेनवैल्यू ​​​​को खोजने के लिए क्षणों के अनुक्रम y0TAkx0, k = 0, 1 , … (जहाँ x0 और y0 इच्छानुसार सदिश हैं) का उपयोग करते है। रुतिशौसर ने इस कार्य के लिए अलेक्जेंडर ऐटकेन का एल्गोरिदम लिया और इसे भागफल-अंतर एल्गोरिदम या क्यूडी एल्गोरिदम में विकसित किया। गणना को उपयुक्त आकार में व्यवस्थित करने के बाद, उन्होंने पाया कि क्यूडी एल्गोरिदम वास्तव में पुनरावृत्ति Ak = LkUk (एलयू अपघटन), Ak+1 = UkLk, द्वारा होती है त्रिविकर्ण आव्युह पर प्रयुक्त किया जाता है, जिसमें से एलआर एल्गोरिदम अनुसरण करता है।[10]

अन्य प्रकार

क्यूआर एल्गोरिथ्म का प्रकार, गोलूब-कहान-रीन्स्च एल्गोरिथ्म सामान्य आव्युह को द्विभुजीय आव्युह में कम करने के साथ प्रारंभ होता है।[11] एकवचन मानों की गणना के लिए क्युआर एल्गोरिदम के इस संस्करण का वर्णन सबसे पहले किसके द्वारा किया गया था? अर्थात गोलुब & कहन (1965). लापैक सबरूटीन डीबीडीएसयूआर के द्वारा किया गया था | तब उस स्तिथियाँ को कवर करने के लिए कुछ संशोधनों के साथ इस पुनरावृत्त विधि को कार्यान्वित करता है जहां एकवचन मान बहुत छोटे होते हैं (डेमेल & कहन 1990). हाउसहोल्डर रिफ्लेक्शन और, यदि उपयुक्त हो, क्यूआर अपघटन का उपयोग करते हुए पहले चरण के साथ, यह एकवचन मूल्य अपघटन की गणना के लिए डीजीईएसवीडी रूटीन बनाता है। क्यूआर एल्गोरिदम को संबंधित अभिसरण परिणामों के साथ अनंत आयामों में भी प्रयुक्त किया जा सकता है।[12][13]


संदर्भ

  1. J.G.F. Francis, "The QR Transformation, I", The Computer Journal, 4(3), pages 265–271 (1961, received October 1959). doi:10.1093/comjnl/4.3.265
  2. Francis, J. G. F. (1962). "क्यूआर परिवर्तन, II". The Computer Journal. 4 (4): 332–345. doi:10.1093/comjnl/4.4.332.
  3. Vera N. Kublanovskaya, "On some algorithms for the solution of the complete eigenvalue problem," USSR Computational Mathematics and Mathematical Physics, vol. 1, no. 3, pages 637–657 (1963, received Feb 1961). Also published in: Zhurnal Vychislitel'noi Matematiki i Matematicheskoi Fiziki, vol.1, no. 4, pages 555–570 (1961). doi:10.1016/0041-5553(63)90168-X
  4. 4.0 4.1 Golub, G. H.; Van Loan, C. F. (1996). मैट्रिक्स संगणना (3rd ed.). Baltimore: Johns Hopkins University Press. ISBN 0-8018-5414-8.
  5. 5.0 5.1 Demmel, James W. (1997). अनुप्रयुक्त संख्यात्मक रैखिक बीजगणित. SIAM.
  6. 6.0 6.1 Trefethen, Lloyd N.; Bau, David (1997). संख्यात्मक रैखिक बीजगणित. SIAM.
  7. Ortega, James M.; Kaiser, Henry F. (1963). "The LLT and QR methods for symmetric tridiagonal matrices". The Computer Journal. 6 (1): 99–101. doi:10.1093/comjnl/6.1.99.
  8. "linear algebra - Why is uncomputability of the spectral decomposition not a problem?". MathOverflow. Retrieved 2021-08-09.
  9. Watkins, David S. (2007). The Matrix Eigenvalue Problem: GR and Krylov Subspace Methods. Philadelphia, PA: SIAM. ISBN 978-0-89871-641-2.
  10. Parlett, Beresford N.; Gutknecht, Martin H. (2011), "From qd to LR, or, how were the qd and LR algorithms discovered?" (PDF), IMA Journal of Numerical Analysis, 31 (3): 741–754, doi:10.1093/imanum/drq003, hdl:20.500.11850/159536, ISSN 0272-4979
  11. Bochkanov Sergey Anatolyevich. ALGLIB User Guide - General Matrix operations - Singular value decomposition . ALGLIB Project. 2010-12-11. URL:[1] Accessed: 2010-12-11. (Archived by WebCite at https://www.webcitation.org/5utO4iSnR?url=http://www.alglib.net/matrixops/general/svd.php
  12. Deift, Percy; Li, Luenchau C.; Tomei, Carlos (1985). "टोडा अनंत अनेक चरों के साथ बहती है". Journal of Functional Analysis. 64 (3): 358–402. doi:10.1016/0022-1236(85)90065-5.
  13. Colbrook, Matthew J.; Hansen, Anders C. (2019). "अनंत-आयामी क्यूआर एल्गोरिदम पर". Numerische Mathematik. 143 (1): 17–83. arXiv:2011.08172. doi:10.1007/s00211-019-01047-5.


स्रोत

बाहरी संबंध