Обращение матриц методом окаймления

@Datvi · Регистрация: 08.12.2020

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

Добрый день, у меня есть код для обращения матрицы методом окаймления, если я задаю матрицу вручную то всё прекрасно работает, сейчас пытаюсь заполнять матрицу NxN размерности с помощью рандома, сама матрица создается нормально, встроенная функция спокойно ее обращает, но на мою функцию ругается, кто-нибудь может подсказать как это решить?

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
import numpy as np
np.set_printoptions(precision=2, suppress=True)
 
def get_inv(A, depth = 0):
    n = len(A)
    k = n - 1
 
    if n == 1:
        return np.matrix([[ 1 / A[0, 0] ]])
 
    Ap = A[:k, :k]
    V, U = A[k, :k], A[:k, k]
 
    Ap_inv = get_inv(Ap, depth + 1)
    
    alpha = 1 / (A[k, k] - V * Ap_inv * U).item()
    Q = -V * Ap_inv * alpha
    P = Ap_inv - Ap_inv * U * Q
    R = - Ap_inv * U * alpha
 
    A_inv = np.matrix([[0.0] * n for _ in range(n)])
    A_inv[:k, :k] = P
    A_inv[k, :k] = Q[0]
 
    A_inv[:k, k] = R[:, 0]
    A_inv[k, k] = alpha
 
    return A_inv
    
n = input("Введите размерность матрицы: ")
n = int(n)
a1 = np.squeeze((np.random.randint(-10,10, size=(n,n)).astype(np.float)))
 
print("Исходная матрица: ")
print(a1)
print("Обратная матрица методом окаймления: ")
print(get_inv(a1))
print("Обратная матрица встроенными функциями: ")
print(np.linalg.inv(a1))

@Gdez · 19.01.2022, 22:01

Datvi, 12-я строчка - нужно:

Python
1
V, U = A[k, :k], A[:k, k].reshape(-1,1)

@Datvi · 20.01.2022, 10:27 **[ТС]**

Спасибо большое за помощь.

Новые блоги и статьи Все статьи Все блоги /
Контроль заполнения и очистка дат в зависимости от значения перечислений Maks 12.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .	Архитектура слоя интернета для сервера-слоя. Hrethgir 11.04.2026 В продолжение https:/ / www. cyberforum. ru/ blogs/ 223907/ 10860. html Знаешь что я подумал? Раз мы все источники пишем в голове ветки, то ничего не мешает добавить в голову такой источник, который сам. . .	Подстановка значения реквизита справочника в табличную часть документа Maks 10.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: при выборе сотрудника (справочник Сотрудники) в ТЧ документа. . .	Очистка реквизитов документа при копировании Maks 09.04.2026 Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .
модель ЗдравоСохранения 8. Подготовка к разному выполнению заданий anaschu 08.04.2026 https:/ / github. com/ shumilovas/ med2. git main ветка * содержимое блока дэлэй из старой модели теперь внутри зайца новой модели 8ATzM_2aurI	Блокировка документа от изменений, если он открыт у другого пользователя Maks 08.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа, разработанного в конфигурации КА2. Задача: запретить редактирование документа, если он открыт у другого пользователя. / / . . .	Система безопасности+живучести для сервера-слоя интернета (сети). Двойная привязка. Hrethgir 08.04.2026 Далее были размышления о системе безопасности. Сообщения с наклонным текстом - мои. А как нам будет можно проверить, что ссылка наша, а не подделана хулиганами, которая выбросит на другую ветку и. . .	Модель ЗдрввоСохранения 7: больше работников, больше ресурсов. anaschu 08.04.2026 работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20% kYBz3eJf3jQ

@Datvi 0 / 0 / 0 Регистрация: 08.12.2020 Сообщений: 44
	20.01.2022, 10:27 [ТС]
	Спасибо большое за помощь. 0

Обращение матриц методом окаймления

Решение