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

@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 **[ТС]**

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

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): сборка C/C++ проекта из консоли 8Observer8 30.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .	Установка Emscripten SDK (emsdk) и CMake на Windows для сборки C и C++ приложений в WebAssembly (Wasm) 8Observer8 30.01.2026 Чтобы скачать Emscripten SDK (emsdk) необходимо сначало скачать и уставить Git: Install for Windows. Следуйте стандартной процедуре установки Git через установщик. Система контроля версиями Git. . .	Подключение Box2D v3 к SDL3 для Android: физика и отрисовка коллайдеров 8Observer8 29.01.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами. Версия v3 была полностью переписана на Си, в. . .	Инструменты COM: Сохранение данный из VARIANT в файл и загрузка из файла в VARIANT bedvit 28.01.2026 Сохранение базовых типов COM и массивов (одномерных или двухмерных) любой вложенности (деревья) в файл, с возможностью выбора алгоритмов сжатия и шифрования. Часть библиотеки BedvitCOM Использованы. . .
Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 28.01.2026 Содержание блога SDL3 имеет собственные средства для загрузки и отображения PNG-файлов с альфа-каналом и базовой работы с ними. В этой инструкции используется функция SDL_LoadPNG(), которая. . .	Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL3_image 8Observer8 27.01.2026 Содержание блога SDL3_image - это библиотека для загрузки и работы с изображениями. Эта пошаговая инструкция покажет, как загрузить и вывести на экран смартфона картинку с альфа-каналом, то есть с. . .	Влияние грибов на сукцессию anaschu 26.01.2026 Бифуркационные изменения массы гриба происходят тогда, когда мы уменьшаем массу компоста в 10 раз, а скорость прироста биомассы уменьшаем в три раза. Скорость прироста биомассы может уменьшаться за. . .	Воспроизведение звукового файла с помощью SDL3_mixer при касании экрана Android 8Observer8 26.01.2026 Содержание блога SDL3_mixer - это библиотека я для воспроизведения аудио. В отличие от инструкции по добавлению текста код по проигрыванию звука уже содержится в шаблоне примера. Нужно только. . .

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

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

Решение