जैकोबी विधि: 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: | ||
== एल्गोरिथम == | == एल्गोरिथम == | ||
'''Input:''' initial guess ''x''<sup>(0)</sup> to the solution, (diagonal dominant) matrix ''A'', right-hand side vector ''b'', convergence criterion | |||
' | '''Output:''' solution when convergence is reached | ||
'''Comments:''' pseudocode based on the element-based formula above | |||
''k'' = 0 | |||
'''while''' convergence not reached '''do''' | |||
''i'' | '''for''' ''i'' := 1 '''step until''' n '''do''' | ||
''σ'' = 0 | |||
for ''j'' := 1 | '''for''' ''j'' := 1 '''step until''' n '''do''' | ||
'''if''' ''j'' ≠ ''i'' '''then''' | |||
''σ'' = ''σ'' + ''a<sub>ij</sub>'' ''x<sub>j</sub>''<sup>(''k'')</sup> | |||
'''end''' | |||
'''end''' | |||
''x<sub>i</sub>''<sup>(''k''+1)</sup> = (''b<sub>i</sub>'' − ''σ'') / ''a<sub>ii</sub>'' | |||
'''end''' | |||
increment ''k'' | |||
'''end''' | |||
== अभिसरण == | == अभिसरण == | ||
मानक अभिसरण स्थिति (किसी पुनरावृत्त विधि के लिए) तब होती है जब पुनरावृत्ति मैट्रिक्स का [[वर्णक्रमीय त्रिज्या]] 1 से कम होता है: | मानक अभिसरण स्थिति (किसी पुनरावृत्त विधि के लिए) तब होती है जब पुनरावृत्ति मैट्रिक्स का [[वर्णक्रमीय त्रिज्या]] 1 से कम होता है: | ||
| Line 256: | Line 255: | ||
=== पायथन उदाहरण === | === पायथन उदाहरण === | ||
< | import numpy as np | ||
<blockquote> | |||
ITERATION_LIMIT = 1000 | |||
ITERATION_LIMIT = 1000 | </blockquote> | ||
initialize the matrix | |||
<blockquote> | |||
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.]]) | ||
x = np.zeros_like ( | initialize the RHS vector | ||
<blockquote> | |||
b = np.array([6., 25., -11., 15.]) | |||
</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] | ||
if x_new[i] == x_new[i-1]: | |||
break | |||
if np.allclose(x, x_new, atol=1e-10, rtol=0.): | |||
break | |||
x = x_new | |||
<blockquote> | |||
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 में विघटित किया जा सकता है: