Изменение\ упрощение кода

@Marat22 · Регистрация: 16.11.2020

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

Постановка задания: нужно матрицу А разложить на 3 матрицы, L - нижнетругольная, D - диагональная, R - верхнетреугольная. А = LDR. Код написан, работает исправно. Однако я создал функцию trans, которая на выходе выдает матрицы LDR, однако при создании функции, внутри нее созданы еще 3, series_d, series_l, series_r, которые считают сумму ряда. Мне нужно вытащить создании данных функций на создание функции trans. Я приложил скрин суммы, которую считают функции и сам код. Я пробовал вытаскивать def series_d и т.д, но пишет, что не найден i и j, т.к я их задаю в цикле, который в trans. Приложенный код работает. Заранее огромное спасибо

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
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
import numpy as np
 
n = int(input('Введите размер матрицы: '))
 
A = np.random.randint(1, 10, size = (n, n))
X = np.random.randint(1, 10, size = (n, 1))
 
 
 
 
 
def trans(A):
    L = np.eye(n)
    D = np.eye(n)
    R = np.eye(n)
 
    def series_d(m):
        summa = 0
        for k in range(1,m):
            res = L[m - 1][k - 1] * D[k - 1][k - 1] * R[k - 1][m - 1]
            summa = summa + res
        return summa
 
    def series_r(m):
        summa = 0
        for k in range(1,m):
            res = L[m - 1][k - 1] * D[k - 1][k - 1] * R[k - 1][j - 1]
            summa = summa + res
        return summa
 
    def series_l(m):
        summa = 0
        for k in range(1,m):
            res = L[i - 1][k - 1] * D[k - 1][k - 1] * R[k - 1][m - 1]
            summa = summa + res
        return summa
    
    for m in range(1, n + 1):
        for i in range(1, n + 1):
            for j in range(1, n + 1):
                D[m - 1][m - 1] = A[m - 1][m - 1] - series_d(m)
                R[m - 1][j - 1] = (A[m - 1][j - 1] - series_r(m))/D[m - 1][m - 1]
                L[i - 1][m - 1] = (A[i - 1][m - 1] - series_l(m))/D[m - 1][m - 1]
    return L, D, R
 
 
 
 
 
 
L , D, R = trans(A)
print(L)
print(D)
print(R)

@Gdez · 25.10.2021, 20:47

Marat22,

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
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
import numpy as np
 
n = int(input('Введите размер матрицы: '))
 
A = np.random.randint(1, 10, size = (n, n))
X = np.random.randint(1, 10, size = (n, 1))
 
 
 
 
 
def trans(A):
    L = np.eye(n)
    D = np.eye(n)
    R = np.eye(n)
    sd, sr, sl = [], [], []
    def series_d(m):
        summa = 0
        for k in range(1,m):
            res = L[m - 1][k - 1] * D[k - 1][k - 1] * R[k - 1][m - 1]
            summa = summa + res
        return summa
 
    def series_r(m):
        summa = 0
        for k in range(1,m):
            res = L[m - 1][k - 1] * D[k - 1][k - 1] * R[k - 1][j - 1]
            summa = summa + res
        return summa
 
    def series_l(m):
        summa = 0
        for k in range(1,m):
            res = L[i - 1][k - 1] * D[k - 1][k - 1] * R[k - 1][m - 1]
            summa = summa + res
        return summa
    
    for m in range(1, n + 1):
        td = series_d(m)
        sd.append(td)
        D[m - 1][m - 1] = A[m - 1][m - 1] - td
        for i in range(1, n + 1):
            tl = series_l(m)
            sl.append(tl)
            L[i - 1][m - 1] = (A[i - 1][m - 1] - tl)/D[m - 1][m - 1]
            for j in range(1, n + 1):
                tr = series_r(m)
                sr.append(tr)
                R[m - 1][j - 1] = (A[m - 1][j - 1] - tr)/D[m - 1][m - 1]
    return L, D, R, sd, sl, sr
 
 
 
 
 
 
L, D, R, sd, sl, sr = trans(A)
print(L, '\n')
print(D, '\n')
print(R, '\n')
 
print(sd, '\n')
print(sl, '\n')
print(sr, '\n')

Новые блоги и статьи Все статьи Все блоги /
Оттенки серого Argus19 18.03.2026 Оттенки серого Нашёл в интернете 3 прекрасных модуля: Модуль класса открытия диалога открытия/ сохранения файла на Win32 API; Модуль класса быстрого перекодирования цветного изображения в оттенки. . .	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 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .