Вопросы с тегом «scipy»

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

1
Может ли приближенный якобиан с конечными разностями вызвать нестабильность в методе Ньютона?
Я реализовал решатель обратной Эйлера в Python 3 (используя Numpy). Для собственного удобства и в качестве упражнения я также написал небольшую функцию, которая вычисляет аппроксимацию градиента с конечной разностью, чтобы мне не всегда приходилось определять якобиан аналитически (если это вообще возможно!). Используя описания, приведенные в Ascher и Petzold 1998 , …

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

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