Оптимизируйте или исправьте код для интерполяции и оценки параметров. ⇐ Python

[Anonymous]

Это простой вопрос, но он будет немного длинным, чтобы я мог объяснить весь контекст происходящего. По сути, это вопрос/проблема кода, и физический контекст служит лишь для иллюстрации.
Я работаю над воспроизведением рисунков 1 и 2 из статьи.

Эти рисунки:


и


Для Рис. 2. только основной сюжет, а не вставки.
В статье говорится следующее:

"Большая часть наших данных отображается для середины энергетического спектра, где начинается переход. Модель имеет два класса симметрии, характеризующиеся нечетным или четным количеством кубитов вверх, и данные приведены для одного класса симметрии. Чтобы уменьшить статистические флуктуации, мы используем $5 ≤ ND ≤ 4 × 10^4$ случайных реализаций Γi и J_{ij} , как это обычно делается в теории случайных матриц. Собственные значения и собственные векторы вычисляются путем точной диагонализации матрицы Гамильтона для каждой реализации. Таким образом, общее количество интервалов составляет 10^4
Упомянутый гамильтониан:
$H = \sum_i Gamma_i\sigma_i + \sum_{i

Подробнее здесь: https://stackoverflow.com/questions/798 ... estimation