Вычислить выражение матриц

@Getter17 · Регистрация: 16.03.2019

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

Доброго времени суток! Имеется две матрицы:

Python
1
2
3
4
5
6
7
8
MatrixA = np.matrix([[2,2,1,1],[3,4,1,1],
                    [3,2,1,1],[3,4,2,1],
                    [2,4,3,1],[2,4,4,1],
                    [3,3,2,1],[4,3,4,1]])
MatrixB = np.matrix([3, 4,
                     1, 2,
                     4, 3,
                     4, 2])

Нужно найти X = [X, Y, Z, T]
Для его нахождения нужно посчитать выражение: -((A^T * A)^-1) * (A^T * B) Приложил в приложении.
Все бы хорошо, на первое выражение получил ответ, но когда начинаю перемножать две матрицы, выдает ошибку размеров

Python
1
2
3
4
5
print (pow((np.transpose(MatrixA) * MatrixA), -1) * -1)
Вывод: [[-0.29565217 -0.04347826  0.02608696  0.89565217]
       [-0.04347826 -0.23188406  0.07246377  0.71014493]
       [ 0.02608696  0.07246377 -0.11014493 -0.05942029]
       [ 0.89565217  0.71014493 -0.05942029 -4.76231884]]

Python
1
2
print (np.transpose(MatrixA) * MatrixB)
Вывод: ValueError: shapes (4,8) and (1,8) not aligned: 8 (dim 1) != 1 (dim 0)

Попробовал решить на онлайн сайте, выражение решается и выдает следующие результаты: [61 78 52 23]
Не пойму как перемножить.

@Getter17 · 23.12.2019, 16:30 **[ТС]**

Вот само выражение.

@dondublon · 23.12.2019, 19:01

Для начала - "-1" при работе с матрицами обозначает обратную матрицу, а pow(M, -1) - это совсем не то, это поэлементное возведение в степень.
Ну а потом сверяйте размеры, что ж тут ещё посоветовать.

@Getter17 · 24.12.2019, 09:54 **[ТС]**

dondublon, ясно, пробую, пока безрезультатно.

Добавлено через 29 минут
dondublon, сделал так:

Python
1
2
3
4
5
6
print(np.linalg.matrix_power((np.transpose(MatrixA) * MatrixA) * -1, -1))
Вывод:
[[-0.29565217 -0.04347826  0.02608696  0.89565217]
 [-0.04347826 -0.23188406  0.07246377  0.71014493]
 [ 0.02608696  0.07246377 -0.11014493 -0.05942029]
 [ 0.89565217  0.71014493 -0.05942029 -4.76231884]]

А это уже конкретная функция для возведения в степень от np, и где правда, когда оно одинаковое с pow?

@dondublon · 24.12.2019, 11:07

Getter17, то, что вам нужно, называется inverse matrix.

@Getter17 · 24.12.2019, 11:53 **[ТС]**

dondublon, ясно, вообщем потеряв кучу времени, решил.

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
from numpy import matrix
from numpy import linalg
A = matrix([[2,2,1,1],
            [3,4,1,1],
            [3,2,1,1],
            [3,4,2,1],
            [2,4,3,1],
            [2,4,4,1],
            [3,3,2,1],
            [4,3,4,1]])
L = matrix([[3],[4],[1],[2],[4],[3],[4],[2]])
print (-(A.T*A).I * (A.T * L))
вывод:
[[ 0.53043478]
 [-0.63768116]
 [ 0.14927536]
 [-2.59710145]]

Новые блоги и статьи Все статьи Все блоги /
Автозаполнение реквизита при выборе элемента справочника Maks 27.03.2026 Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. При выборе "Спецтехники" (Тип Справочник. Спецтехника), заполняется. . .	Сумматор с применением элементов трёх состояний. Hrethgir 26.03.2026 Тут. https:/ / fips. ru/ EGD/ ab3c85c8-836d-4866-871b-c2f0c5d77fbc Первый документ красиво выглядит, но без схемы. Это конечно не даёт никаких плюсов автору, но тем не менее. . . всё может быть. . .	Автозаполнение реквизитов при создании документа Maks 26.03.2026 Программный код из решения ниже размещается в модуле объекта документа, в процедуре "ПриСозданииНаСервере". Алгоритм проверки заполнения реализован для исключения перезаписи значения реквизита,. . .	Команды формы и диалоговое окно Maks 26.03.2026 1. Команда формы "ЗаполнитьЗапчасти". Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. В качестве источника данных. . .
Кому нужен AOT? DevAlt 26.03.2026 Решил сделать простой ланчер Написал заготовку: dotnet new console --aot -o UrlHandler var items = args. Split(":"); var tag = items; var id = items; var executable = args;. . .	Отправка уведомления на почту при изменении наименования справочника Maks 24.03.2026 Программная отправка письма электронной почты на примере изменения наименования типового справочника "Склады" в конфигурации БП3. Перед реализацией необходимо выполнить настройку системной учетной. . .	модель ЗдравоСохранения 5. Меньше увольнений- больше дохода! anaschu 24.03.2026 Теперь система здравосохранения уменьшает количество увольнений. 9TO2GP2bpX4 a42b81fb172ffc12ca589c7898261ccb/ https:/ / rutube. ru/ video/ a42b81fb172ffc12ca589c7898261ccb/ Слева синяя линия -. . .	Midnight Chicago Blues kumehtar 24.03.2026 Такой Midnight Chicago Blues, знаешь?. . Когда вечерние улицы становятся ночными, а ты не можешь уснуть. Ты идёшь в любимый старый бар, и бармен наливает тебе виски. Ты смотришь на пролетающие. . .

@Getter17 0 / 0 / 0 Регистрация: 16.03.2019 Сообщений: 81
	23.12.2019, 16:30 [ТС]
	Вот само выражение. Изображения 0

@dondublon 4653 / 2073 / 366 Регистрация: 17.03.2012 Сообщений: 10,183 Записей в блоге: 6
	23.12.2019, 19:01
	Для начала - "-1" при работе с матрицами обозначает обратную матрицу, а pow(M, -1) - это совсем не то, это поэлементное возведение в степень. Ну а потом сверяйте размеры, что ж тут ещё посоветовать. 0

@dondublon 4653 / 2073 / 366 Регистрация: 17.03.2012 Сообщений: 10,183 Записей в блоге: 6
	24.12.2019, 11:07
	Getter17, то, что вам нужно, называется inverse matrix. 0