Sympy: получить столбец мощности матрицы

Кажется, это простой вопрос. Мой код читает

import sympy as sp
D1 = sp.symbols('D1')
D2 = sp.symbols('D2')
m1 = sp.diag(*[2 * D1, 2 * D1, 2 * D2, 2 * D2])
b = sp.sqrt(m1)

m1.col(0)
b.col(0)

т.е. Я хочу получить столбец квадратного корня матрицы. Однако похоже, что объект MatPow не имеет атрибута b, который есть у обычных матриц, таких как .col.

Выход:

AttributeError                            Traceback (most recent call last)
<ipython-input-55-855c6e53a412> in <module>()
----> 1 b.col(0)

AttributeError: 'MatPow' object has no attribute 'col'

Что я делаю неправильно?

Анализ настроения постов в Twitter с помощью Python, Tweepy и Flair
Анализ настроения постов в Twitter с помощью Python, Tweepy и Flair
Анализ настроения текстовых сообщений может быть настолько сложным или простым, насколько вы его сделаете. Как и в любом ML-проекте, вы можете выбрать...
7 лайфхаков для начинающих Python-программистов
7 лайфхаков для начинающих Python-программистов
В этой статье мы расскажем о хитростях и советах по Python, которые должны быть известны разработчику Python.
Установка Apache Cassandra на Mac OS
Установка Apache Cassandra на Mac OS
Это краткое руководство по установке Apache Cassandra.
Сертификатная программа "Кванты Python": Бэктестер ансамблевых методов на основе ООП
Сертификатная программа "Кванты Python": Бэктестер ансамблевых методов на основе ООП
В одном из недавних постов я рассказал о том, как я использую навыки количественных исследований, которые я совершенствую в рамках программы TPQ...
Создание персонального файлового хранилища
Создание персонального файлового хранилища
Вы когда-нибудь хотели поделиться с кем-то файлом, но он содержал конфиденциальную информацию? Многие думают, что электронная почта безопасна, но это...
Создание приборной панели для анализа данных на GCP - часть I
Создание приборной панели для анализа данных на GCP - часть I
Недавно я столкнулся с интересной бизнес-задачей - визуализацией сбоев в цепочке поставок лекарств, которую могут просматривать врачи и...
1
0
249
1
Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

Вам нужно преобразовать объект MatPow в MutableDenseMatrix с помощью метода as_mutable:

b_as_m = b.as_mutable()
print(b_as_m.col(0))

Выход:

Matrix([[sqrt(2)*sqrt(D1)], [0], [0], [0]])

Спасибо! Я отмечаю этот ответ как принятый, но не могли бы вы объяснить, почему это не делается автоматически?

Dr_Zaszuś 23.05.2019 09:44

@Dr_Zaszuś Думаю, это из-за возможных оптимизаций. Например, если вы попытаетесь вычислить (m1 ** 0.5) ** 2, sympy может вернуть m1 без фактического вычисления temp = m1 ** 0.5 и result = temp ** 2.

sanyassh 23.05.2019 10:13

Верно, но это не объясняет, почему нельзя взять столбец матричного объекта...

Dr_Zaszuś 23.05.2019 13:17

@Dr_Zaszuś есть какая-то ленивая оценка. Матричноподобный объект b не является полностью вычисляемой матрицей - значения его столбцов теперь не определены. Известно лишь, что b = m1 ** 0.5. Чтобы узнать значения, вам нужно сначала оценить их с помощью as_mutable.

sanyassh 23.05.2019 15:12

Другие вопросы по теме