Dieser Beitrag ist ein Artikel für den 11. Dezember 2019 in Adventskalender für numerische Berechnungen 2019, aber es ist fast ein Gedicht ...
Ist ein geordneter Vektorraum ($ \ deg V_k = k $). Zu diesem Zeitpunkt sollte ** Poancare-Reihe ** $ P (V; t) $ von $ V $ die folgende formale Potenzreihe sein.
Sei $ V = \ mathbb {C} [x] $ (komplexer polymorpher polymorpher Ring mit komplexem Koeffizienten) und sei $ V_k, \ (k = 0, 1, 2 ...) $ das gesamte Polypoly der Ordnung der Ordnung Versuchen wir es als Vektorraum. Dann kann $ \ mathbb {C} [x] $ als ein Vektorraum mit Ordnung betrachtet werden, aber da es nur eine Variable gibt, ist das synchrone polymorphe Polymorph $ k $ Grad $ a_k x ^ k (a_k \ in \ mathbb {C}) Es gibt nur eine in Form von)
Betrachten Sie nun den polymorphen Ring mit zwei Variablen $ V = \ mathbb {C} [x, y] $. Die $ k $ -Bestellungskomponente $ V_k $ von $ V $ ist der Vektorraum des gesamten $ k $ -Ordnungspolypolys wie oben. In diesem Fall ändert sich $ \ dim V_k $ abhängig vom Wert von $ k $. Blick von der kleineren
\begin{align}
V_0 &= \text{Span}\langle 1 \rangle, \\
V_1 &= \text{Span}\langle x, y\rangle, \\
V_2 &= \text{Span}\langle x^2, xy, y^2\rangle, \\
V_3 &= \text{Span}\langle x^3, x^2y, xy^2, y^3\rangle, \\
&...
\end{align}
Und so weiter können Sie sehen, dass $ \ dim V_k = k + 1
\begin{align}
P(\mathbb{C}[x, y]; t) &= \sum_{k = 0}^\infty \frac{d}{dt}t^{k + 1} \\
&= \frac{d}{dt} \left( t \sum_{k = 0} t^{k} \right) = \frac{d}{dt} \frac{t}{1 - t} = \frac{1}{(1 - t)^2}
\end{align}
Sie können in geschlossener Form schreiben.
Was ist in ähnlicher Weise, wenn $ V = \ mathbb {C} [x, y, z] $, $ V_k $ ein Unterraum ist, der das gesamte Polypoly $ k $ order-order umfasst? $ V_k $ basiert auf einem sequentiellen Ausdruck der Form $ x ^ p y ^ q z ^ r $, $ p + q + r = k $.
\frac{(k + 2)!}{k!2!} = \dbinom{k + 2}{2}
Sie können sehen, dass es verschiedene Typen gibt
Nun möchte ich durch die Berechnung des ersten numerischen Terms der Potenzreihe bestätigen, dass diese (★) Formel wirklich zu gelten scheint. Was soll ich tun?
In diesem Artikel werde ich versuchen, durch Programmieren mit Python zu berechnen.
Zunächst sollten Sie den Binomialkoeffizienten auf der rechten Seite berechnen. Die Python-Bibliothek SciPy hat eine Funktion "Kamm", die Binomialkoeffizienten berechnet. Wenn Sie dies verwenden, um bis zu den ersten 10 Begriffen zu berechnen,
from scipy.special import comb
combinations = [comb(k + 2, k, exact=True) for k in range(10)]
print(combinations)
Alles was Sie tun müssen, ist das Programm auszuführen. Ergebnis ist
[1, 3, 6, 10, 15, 21, 28, 36, 45, 55]
ist geworden.
Was ist nun mit der linken Seite?
Allgemeiner Ausdruck für die Taylor-Expansion um $ t = 0 $
f(t) = f(0) + f'(0)t + \frac{1}{2!}f''(0)t^2 + \frac{1}{3!}f'''(0)t^3 + ... = \sum_{k=0}^\infty \frac{1}{k!} \frac{d^kf(0)}{dt^k}t^k
Unter Hinweis darauf ist der Koeffizient der Differenzkoeffizient der Ordnung $ k $ bei $ t = 0 $ auf der linken Seite geteilt durch $ k! $. Berechnen wir dies durch Programmieren mit Python. Solche mathematischen Berechnungen können mit einer Bibliothek namens SymPy durchgeführt werden. Wenn Sie sich entscheiden, bis zu den ersten 10 Termen zu berechnen, müssen Sie das Differential 10. Ordnung (!) ** berechnen, was durch manuelle Berechnung sehr kompliziert ist, aber Sie können es leicht finden, indem Sie ein Programm schreiben. Ich werde. Das eigentliche Programm ist wie folgt. (Ich habe "scipy" verwendet, um den Multiplikator zu berechnen.)
from sympy import *
from scipy.special import factorial
t = Symbol('t')
P = 1 / (1 - t) ** 3
coeffs = []
for k in range(10):
dk = diff(P, t, k).subs([(t, 0)])
coeff = dk / factorial(k, exact=True)
coeffs.append(coeff)
print(coeffs)
Das Ausführungsergebnis ist noch
[1, 3, 6, 10, 15, 21, 28, 36, 45, 55]
Es wurde bestätigt, dass die Koeffizienten auf der linken und rechten Seite mit den ersten 10 Termen der Potenzreihenerweiterung übereinstimmen.
Ich dachte vage, dass ich eine solche Kombinationstheorie mit einem Computer berechnen wollte, also schrieb ich sie, obwohl sie sehr einfach ist.
Recommended Posts