जैकोबी विधि: Difference between revisions

From Vigyanwiki
(Created page with "{{Short description|Iterative method used to solve a linear system of equations}} {{Distinguish|Jacobi eigenvalue algorithm}} संख्यात्मक रैखिक...")
 
No edit summary
Line 2: Line 2:
{{Distinguish|Jacobi eigenvalue algorithm}}
{{Distinguish|Jacobi eigenvalue algorithm}}


[[संख्यात्मक रैखिक बीजगणित]] में, जैकोबी विधि (उर्फ जैकोबी पुनरावृति विधि) रैखिक समीकरणों के विकर्ण रूप से प्रभावी मैट्रिक्स प्रणाली के समाधान का निर्धारण करने के लिए एक पुनरावृत्त एल्गोरिथम है। प्रत्येक विकर्ण तत्व के लिए हल किया जाता है, और एक अनुमानित मान प्लग इन किया जाता है। प्रक्रिया तब तक दोहराई जाती है जब तक कि यह अभिसरण न हो जाए। यह एल्गोरिथम [[जैकोबी ईजेनवेल्यू एल्गोरिथम]] का एक स्ट्रिप्ड-डाउन संस्करण है। विधि का नाम [[कार्ल गुस्ताव जैकब जैकोबी]] के नाम पर रखा गया है।
[[संख्यात्मक रैखिक बीजगणित]] में, '''जैकोबी विधि''' रैखिक समीकरणों के एक सख्ती से  विकर्णतः प्रभावी प्रणाली के समाधान का निर्धारण करने के लिए एक पुनरावृत्ति एल्गोरिथ्म है। प्रत्येक विकर्ण तत्व के लिए हल किया जाता है, और एक अनुमानित मान प्लग इन किया जाता है। प्रक्रिया तब तक दोहराई जाती है जब तक कि यह अभिसरण न हो जाए। यह एल्गोरिथम [[जैकोबी ईजेनवेल्यू एल्गोरिथम]] का एक स्ट्रिप्ड-डाउन संस्करण है। विधि का नाम [[कार्ल गुस्ताव जैकब जैकोबी]] के नाम पर रखा गया है।


== विवरण ==
== विवरण ==
Line 19: Line 19:


== एल्गोरिथम ==
== एल्गोरिथम ==
  इनपुट: {{nowrap|initial guess ''x''<sup>(0)</sup> to the solution}}, (विकर्ण प्रभावी) मैट्रिक्स A, दाएँ हाथ की ओर सदिश b, अभिसरण मानदंड
  '''Input:''' initial guess ''x''<sup>(0)</sup> to the solution, (diagonal dominant) matrix ''A'', right-hand side vector ''b'', convergence criterion
  'आउटपुट:' {{nowrap|solution when convergence is reached}}
  '''Output:''' solution when convergence is reached
  टिप्पणियाँ: उपरोक्त तत्व-आधारित सूत्र के आधार पर स्यूडोकोड
  '''Comments:''' pseudocode based on the element-based formula above
   
   
{{nowrap|1=''k'' = 0}}
''k'' = 0
  जबकि अभिसरण नहीं हुआ है
  '''while''' convergence not reached '''do'''
     ''i'' के लिए := 1 चरण तक n करें       
     '''for''' ''i'' := 1 '''step until''' n '''do'''
{{nowrap|1=''σ'' = 0}}
        ''σ'' = 0
         for ''j'' := 1 कदम तक n करते हैं
         '''for''' ''j'' := 1 '''step until''' n '''do'''
             अगर ''जे'' ≠ ''मैं'' तो               
             '''if''' ''j'' ≠ ''i'' '''then'''
{{nowrap|1=''σ'' = ''σ'' + ''a''<sub>''ij''</sub> ''x''<sub>''j''</sub><sup>(''k'')</sup>}}
                ''σ'' = ''σ'' + ''a<sub>ij</sub>'' ''x<sub>j</sub>''<sup>(''k'')</sup>
             अंत
             '''end'''
         अंत       
         '''end'''
{{nowrap|1=''x''<sub>''i''</sub><sup>(''k''+1)</sup> = (''b''<sub>''i''</sub> ''σ'') / ''a''<sub>''ii''</sub>}}
        ''x<sub>i</sub>''<sup>(''k''+1)</sup> = (''b<sub>i</sub>'' − ''σ'') / ''a<sub>ii</sub>''
     अंत
     '''end'''
     वेतन वृद्धि ''के''
     increment ''k''
  अंत
  '''end'''


== अभिसरण ==
== अभिसरण ==
<!-- [[Matrix splitting]] links here.  Please do not change. -->
मानक अभिसरण स्थिति (किसी पुनरावृत्त विधि के लिए) तब होती है जब पुनरावृत्ति मैट्रिक्स का [[वर्णक्रमीय त्रिज्या]] 1 से कम होता है:
मानक अभिसरण स्थिति (किसी पुनरावृत्त विधि के लिए) तब होती है जब पुनरावृत्ति मैट्रिक्स का [[वर्णक्रमीय त्रिज्या]] 1 से कम होता है:


Line 256: Line 255:


=== पायथन उदाहरण ===
=== पायथन उदाहरण ===
<वाक्यविन्यास प्रकाश लैंग = संख्यात्मक रेखा = 1>
import numpy as np
Numpy को np के रूप में आयात करें
<blockquote>
 
  ITERATION_LIMIT = 1000
ITERATION_LIMIT = 1000
</blockquote>
 
initialize the matrix
# मैट्रिक्स को इनिशियलाइज़ करें
<blockquote>
= एनपी। सरणी (10।, -1।, 2।, 0।],
A = np.array([[10., -1., 2., 0.],
</blockquote>
               [-1., 11., -1., 3.],
               [-1., 11., -1., 3.],
               [2., -1., 10., -1.],
               [2., -1., 10., -1.],
               [0.0, 3., -1., 8.)
               [0.0, 3., -1., 8.]])
# आरएचएस वेक्टर को इनिशियलाइज़ करें
बी = एनपी। सरणी ([6।, 25।, -11।, 15।])
 
# सिस्टम प्रिंट करता है
प्रिंट (सिस्टम:)
for i in range(A.shape[0]):
    पंक्ति = [एफ {ए [आई, जे]} * एक्स {जे + 1} फॉर जे इन रेंज (ए आकार [1])]
    प्रिंट (एफ '{ + .जॉइन (पंक्ति)} = {बी [i]}')
प्रिंट ()


x = np.zeros_like ()
initialize the RHS vector
इसके लिए_गिनती श्रेणी में (ITERATION_LIMIT):
<blockquote>
     अगर it_count! = 0:
b = np.array([6., 25., -11., 15.])
         प्रिंट (एफ पुनरावृत्ति {it_count}: {x})
</blockquote>
prints the system
<blockquote>
print("System:") for i in range(A.shape[0]):
</blockquote>
    row = [f"{A[i, j]}*x{j + 1}" for j in range(A.shape[1])]
    print(f'{" + ".join(row)} = {b[i]}')
<blockquote>
print() x = np.zeros_like(b) for it_count in range(ITERATION_LIMIT):
</blockquote>
     if it_count != 0:
         print(f"Iteration {it_count}: {x}")
     x_new = np.zeros_like(x)
     x_new = np.zeros_like(x)
     for i in range(A.shape[0]):
     for i in range(A.shape[0]):
         s1 = np.dot(A[i, :i], x[:i])
         s1 = np.dot(A[i, :i], x[:i])
         s2 = np.dot(A[i, i + 1:], x[i + 1:])
         s2 = np.dot(A[i, i + 1:], x[i + 1:])
         x_new[i] = (b[i] - s1 - s2) / A[i, i]
         x_new[i] = (b[i] - s1 - s2) / A[i, i]
         अगर x_new[i] == x_new[i-1]:
         if x_new[i] == x_new[i-1]:
           तोड़ना
           break
 
     if np.allclose(x, x_new, atol=1e-10, rtol=0.):
     अगर np.allclose(x, x_new, atol=1e-10, rtol=0.):
         break
         तोड़ना
     x = x_new
 
<blockquote>
     एक्स = x_new
print("Solution: ") print(x) error = np.dot(A, x) - b print("Error:") print(error)
 
</blockquote>
प्रिंट (समाधान:)
प्रिंट (एक्स)
त्रुटि = एनपी डॉट (, एक्स) - बी
प्रिंट (त्रुटि:)
प्रिंट (त्रुटि)
</वाक्यविन्यास हाइलाइट>


== भारित जैकोबी विधि ==
== भारित जैकोबी विधि ==

Revision as of 23:18, 25 May 2023

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

विवरण

होने देना n रैखिक समीकरणों की एक वर्ग प्रणाली हो, जहाँ:

कब और जाने जाते हैं, और अज्ञात है, हम अनुमान लगाने के लिए जैकोबी पद्धति का उपयोग कर सकते हैं . सदिश के लिए हमारे प्रारंभिक अनुमान को दर्शाता है (अक्सर के लिए ). हम निरूपित करते हैं के-वें सन्निकटन या पुनरावृत्ति के रूप में , और का अगला (या k+1) पुनरावृत्ति है .

मैट्रिक्स आधारित सूत्र

तब A को एक विकर्ण मैट्रिक्स घटक D, एक निचला त्रिकोणीय भाग L और एक ऊपरी त्रिकोणीय भाग U में विघटित किया जा सकता है: