Преобразование кода метода Якоби в метод Зейделя

@blackjstyle · Регистрация: 20.09.2016

Студворк — интернет-сервис помощи студентам

Что нужно добавить в код чтобы это был метод Зейделя

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
import numpy as np
 
def jacobi(A, B, tol=1e-6, max_iter=500):
    n = A.shape[0]
    B = np.reshape(B, (n, -1))
    C = np.eye(n) - A / np.outer(np.diag(A), np.ones(n))
    D = B / np.outer(np.diag(A), np.ones(B.shape[1]))
    P, X = D.copy(), D.copy()
    for i in range(max_iter):
        X = C @ X + D
        if np.amax(np.absolute(X-P)) <= tol:
            return X, i + 1
        P = X.copy()
    print('Maximum number of iterations exceeded.')
    return X, i + 1
 
if __name__ == '__main__':
    n = 5
    A = np.random.normal(0, 100, (n, n))
    for i in range(n):
        A[i,i] = np.sum(np.absolute(A[i,:])) + 1
    B = np.random.normal(0, 100, (n, 3))
    X, it = jacobi(A, B)
    print(X, it)

@dondublon · 02.05.2018, 13:14

Функцию zeidel (или как-то так), которая будет решать систему методом Зейделя.

(Каков вопрос - таков ответ.)

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Desktop (MinGW): Рисуем цветные прямоугольники с помощью рисовальщика SDL3 на Си и C++ 8Observer8 17.03.2026 Содержание блога Финальные проекты на Си и на C++: finish-rectangles-sdl3-c. zip finish-rectangles-sdl3-cpp. zip	Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая ссылка» (hard link),. . .	[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .	делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .

@dondublon 4653 / 2073 / 366 Регистрация: 17.03.2012 Сообщений: 10,183 Записей в блоге: 6
	02.05.2018, 13:14
	Функцию zeidel (или как-то так), которая будет решать систему методом Зейделя. (Каков вопрос - таков ответ.) 0