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

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

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

Добрый день, нужно написать код обращения матриц методом окаймления на языке Python. Изначально задается квадратная матрица. Я примерно представляю алгоритм решения на бумаге, и есть алгоритм для программы MatLab'а, но как реализовать код в Python'е не имею ни малейшего представления.

@Datvi · 04.01.2022, 11:34 **[ТС]**

Оставляю здесь, если кому то понадобится решение

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
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
A = np.matrix([
   [2, 4, 3, 8],
   [5, 6, 8, 2],
   [1, 3, 7, 5],
   [7, 7, 2, 1]
])
print("Обратная матрица методом окаймления: ")
print(get_inv(A))
print("Обратная матрица встроенными функциями: ")
print(np.linalg.inv(A))

Результаты вывода:
Обратная матрица методом окаймления:
[[ 2.53 4.33 -5.2 -2.93]
[-2.72 -4.67 5.55 3.32]
[ 0.27 0.67 -0.6 -0.47]
[ 0.75 1. -1.25 -0.75]]
Обратная матрица встроенными функциями:
[[ 2.53 4.33 -5.2 -2.93]
[-2.72 -4.67 5.55 3.32]
[ 0.27 0.67 -0.6 -0.47]
[ 0.75 1. -1.25 -0.75]]

Новые блоги и статьи Все статьи Все блоги /
Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Символьное дифференцирование igorrr37 13.02.2026 / * Программа принимает математическое выражение в виде строки и выдаёт его производную в виде строки и вычисляет значение производной при заданном х Логарифм записывается как: (x-2)log(x^2+2) -. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .
И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»	SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .

@Datvi 0 / 0 / 0 Регистрация: 08.12.2020 Сообщений: 44

	Обращение матриц методом окаймления 20.12.2021, 22:11. Показов 1569. Ответов 1 Метки нет (Все метки) Добрый день, нужно написать код обращения матриц методом окаймления на языке Python. Изначально задается квадратная матрица. Я примерно представляю алгоритм решения на бумаге, и есть алгоритм для программы MatLab'а, но как реализовать код в Python'е не имею ни малейшего представления. 0