Optimieren oder korrigieren Sie einen Code für die Interpolation und Parameterschätzung ⇐ Python

[Anonymous]

Dies ist eine einfache Frage, die jedoch etwas lang sein wird, damit ich den gesamten Kontext zum Geschehen angeben kann. Im Wesentlichen handelt es sich hier um eine Frage/ein Problem mit dem Code, und der physische Kontext dahinter dient nur zur Veranschaulichung.
Ich arbeite an der Reproduktion der Abbildungen 1 und 2 aus einem Papier.

Die Abbildungen sind:


und


Für Abb. 2. nur das Hauptdiagramm, nicht die Einfügungen.
Der Artikel lautet wie folgt:

„Der Großteil unserer Daten wird für die Mitte des Energiespektrums angezeigt, wo der Übergang beginnt. Das Modell hat zwei Symmetrieklassen, die durch eine ungerade oder gerade Anzahl von Qubits gekennzeichnet sind, und die Daten werden für eine Symmetrieklasse angegeben. Um statistische Schwankungen zu reduzieren, verwenden wir $5 ≤ ND ≤ 4 × 10^4$ zufällige Realisierungen von Γi und J_{ij} , wie es normalerweise in der Zufallsmatrixtheorie geschieht. Eigenwerte und Eigenvektoren werden durch exakte Diagonalisierung der Hamilton-Matrix für jede Realisierung berechnet. Auf diese Weise beträgt die Gesamtzahl der Abstände 10^4 < NS ≤ 1,6 × 10^5 (NS ∝ ND NH ). Abb. 1. Um die Entwicklung von P(s) mit der Kopplung J zu analysieren, ist es praktisch, den Parameter η = [∫₀^{s₀} P(s)ds - I_W] / [I_P - I_W] zu verwenden, wobei I_P = ∫₀^{s₀} P_P(s)ds, I_W = ∫₀^{s₀} P_W(s)ds, und s₀=0,4729... ist der Schnittpunkt von P_{P}(s) und P_{W}(s). Auf diese Weise entspricht P_{P}(s) η=1 und P_{W}(s) η=0. Wie im Bereich des Quantenchaos üblich, charakterisiert die Variation von η die Entwicklung von P(s). J/∆0 ist in Abb.2 für δ=∆0 dargestellt und zeigt, dass η mit zunehmender Kopplungsstärke tatsächlich steiler zu werden scheint. Der typische J_{c}-Wert, in dessen Nähe der Übergang stattfindet, entspricht den Zwischenwerten von η. Die Abhängigkeit von J_{c} von n ist in Abb.2 angegeben />
Der erwähnte Hamilton-Operator ist:
$H = \sum_i Gamma_i\sigma_i + \sum_{i