Запитання з тегом «scipy»

1
Коли Ньютон-Крилов не є відповідним вирішувачем?
Нещодавно я порівнював різні нелінійні розв'язувачі з scipy і був особливо вражений прикладом Ньютона-Крилова в кулінарній книзі Scipy, в якому вони вирішують диференціальне рівняння рівняння другого порядку з нелінійним терміном реакції приблизно в 20 рядках коду. Я змінив приклад коду для вирішення нелінійного рівняння Пуассона (його також називають рівнянням Пуассона-Больцмана …

1
Чи може наближений якобійський з скінченними відмінностями викликати нестабільність методу Ньютона?
Я реалізував зворотний рішень-ейлер у python 3 (використовуючи numpy). Для власної зручності і як вправи я також написав невелику функцію, яка обчислює наближення градієнта кінцевою різницею, так що мені не завжди доведеться аналітично визначати якобіанські (якщо це навіть можливо!). Використовуючи описи, наведені в Ascher і Petzold 1998 , я написав …

2
Який найшвидший спосіб обчислити всі власні значення дуже великої і розрідженої матриці суміжності в пітоні?
Я намагаюся розібратися, чи існує швидший спосіб обчислити всі власні значення та власні вектори дуже великої і розрідженої матриці суміжності, ніж використання scipy.sparse.linalg.eigsh Наскільки я знаю, цей метод використовує лише розрідженість і атрибути симетрії матриці. Матриця суміжності також є бінарною, що змушує мене думати, що існує швидший спосіб зробити це. …

1
Алгоритм врівноваження матриці
Я писав набір інструментів системи управління з нуля і суто в Python3 (безсоромний штекер:) harold. У моєму минулому дослідженні я завжди скаржився на вирішувач Ріккаті care.mз причин, які є технічними / неактуальними. Отже, я писав власний набір процедур. Одна річ, яку я не можу знайти, - це отримати високоефективний алгоритм …
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.