$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Diese Studie umfasste weder menschliche Teilnehmer noch Wirbeltiere als Testpersonen. Alle verwendeten Daten sind öffentlich zugängliche Warenkursreihen der SMM, die keine ethische Genehmigung erfordern. Daher wurde für diese Forschung keine ethische Genehmigung eingeholt oder erforderlich.
Dieser Abschnitt beschreibt das umfassende und rigorose Forschungsdesign, das zur empirischen Überprüfung der Kernhypothese durchgeführt wurde. Sie bietet eine detaillierte Darstellung der mathematischen Formulierung und architektonischen Details der dreizehn bewerteten Deep-Learning-Modelle, des präzisen Trainingsprotokolls und der formalen Bewertungsmetriken. Der übergeordnete methodische Workflow wird in Abbildung 1 visuell zusammengefasst.

Abbildung 1: Schematische Übersicht der Forschungsmethodik. Das Diagramm illustriert die vollständige experimentelle Pipeline, einschließlich Datenpartitionierung, Modelltraining ausschließlich auf der Cu-Preisreihe, Bewertung auf dem Cu-Testsatz und Validierung außerhalb der Stichprobe auf unabhängigen Al- und Zn-Reihen. Die gestrichelte Rückkopplungsschleife zeigt die strukturierten Ablationsexperimente an, die durchgeführt wurden, um den Beitrag einzelner architektonischer Komponenten zu analysieren. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Das Diagramm zeigt die vollständige experimentelle Pipeline. Der Prozess beginnt mit der exklusiven Nutzung der Cu-Preisreihe für die Modellentwicklung. Diese Serie ist chronologisch in Trainingssets (80 %), Validierung (10 %) und Testsets (10 %) unterteilt. Dreizehn verschiedene Deep-Learning-Architekturen werden dann ausschließlich auf den Cu-Trainingsdaten trainiert und hyperparametrisiert, wobei das frühzeitige Stoppen über das Validierungsset überwacht wird. Der primäre Benchmark ist die Bewertung dieser Modelle auf dem zurückgehaltenen Cu-Testset. Entscheidend ist, dass zur Bewertung der Generalisierbarkeit exakt dieselben trainierten Modelle ohne Modifikation angewendet werden, um die vollständig unabhängige Al- und Zn-Preisreihen vorherzusagen, was einen strikten Out-of-Sample-Test darstellt. Schließlich werden strukturierte Ablationsexperimente (die gestrichelte Rückkopplungsschleife) durchgeführt, um den Leistungsbeitrag einzelner architektonischer Komponenten (z. B. Aufmerksamkeit, bidirektionale Verarbeitung und Faltungsschichten) zu zerlegen und zu analysieren.
Modellarchitekturen und mathematische Formulierung
Wir entwarfen und implementierten ein Spektrum von 13 DL-Modellen und erhöhten systematisch die architektonische Komplexität von einfachen rekurrenten Netzwerken bis hin zu ausgefeilten Multikomponenten-Hybriden. Alle Modelle teilen dasselbe Kernziel: eine Abbildung
aus einem historischen Preisfenster Xt = [P t-L,P t-L+1,...,Pt-1] der Länge L = 30 zum nachfolgenden Preis yt = Pt zu lernen.
Die Modellfamilien sind wie folgt definiert:
Basismodelle für rekurrente Modelle
GRU: Ein schlankes rekurrierendes Netzwerk, das Update- (z-t) und Reset-(rt)-Gatter verwendet, um den Informationsfluss zu modulieren. Der verborgene Zustand ht wird berechnet als:
(1)
(2)
(3)
(4)
wobei X die sigmoidale Aktivierung ist,
das Hadamard-Produkt bezeichnet und xt der Eingang zum Zeitpunkt t ist. Der endgültige verborgene Zustand hL wird durch eine lineare Ausgangsschicht geleitet. Die Gleichungen 1–4 sind von Cho et al.14 übernommen.
LSTM: Nutzt Eingangsgatter (it), Forget (ft) und Ausgabegatter (ot), um einen Zellzustand (Ct) aufrechtzuerhalten, was eine explizite Kontrolle über das Langzeitgedächtnis ermöglicht.
Bidirektionale Modelle (BiGRU und BiLSTM)
Diese Modelle integrieren zwei separate rekurrente Schichten, die die Sequenz in Vorwärts- und Rückwärtsrichtung verarbeiten. Die letzte versteckte Darstellung in jedem Zeitschritt ist die Verkettung
, die theoretisch kontextuelle Informationen aus Vergangenheit und Zukunft innerhalb des festen Eingabefensters erfasst.
Aufmerksamkeitsaugmentierte Modelle (GRU–Attention und LSTM–Attention)
Ein additiver Aufmerksamkeitsmechanismus wird auf die Folge der verborgenen Zustände H = [h1,h 2,...,hL] angewendet, die von der finalen rekurrenten Schicht erzeugt werden. Der Kontextvektor ist als gewichtete Summe definiert:
(5)
(6)
(7)
Hier steht αi für das Aufmerksamkeitsgewicht, das dem i-ten historischen Zeitschritt zugewiesen wird. Der Kontextvektor c, der eine adaptive Zusammenfassung der relevanten Geschichte kapselt, wird der finalen Vorhersageschicht zugeführt. Die Gleichungen 5–7 sind von Bello et al.47 adaptiert.
CNN–Hybridmodelle (CNN–GRU und CNN–LSTM)
Eine eindimensionale CNN-Schicht mit rektifizierter linearer Einheitsaktivierung (ReLU) wird der rekurrenten Schicht vorangebracht.
Komplexe Hybridmodelle
Diese Architekturen kombinieren mehrere Komponenten (z. B. CNN–BiGRU–Attention, CNN–BiLSTM–Attention). Sie repräsentieren den Stand der Technik in Bezug auf Komplexität und zielen darauf ab, lokale Musterextraktion (CNN), bidirektionale kontextuelle Modellierung und adaptive zeitliche Gewichtung (Aufmerksamkeit) in einem einzigen Rahmen zu integrieren.
Alle Modelle wurden mit konsistenten versteckten Zustandsdimensionen konfiguriert (128 Einheiten für rekurrente Schichten und 64 Filter für CNN-Schichten) und wurden mit einer einzigen linearen Ausgangsschicht finalisiert. Dieses kontrollierte Design stellt sicher, dass Leistungsunterschiede auf architektonische Entscheidungen und nicht auf Abweichungen bei der Modellkapazitätsoptimierung zurückzuführen sind. Die Anzahl der trainierbaren Parameter stieg folglich entlang dieses Spektrums erheblich an.
Trainingsprotokoll, Hyperparameter und Ablationsstudiendesign
Tabelle 1 fasst die einheitliche und rigorose experimentelle Konfiguration zusammen, die auf alle dreizehn Modelle angewendet wurde, um einen fairen Vergleich zu gewährleisten und Überanpassungen zu mindern. Alle Modelle wurden von Grund auf ausschließlich mit dem vorgesehenen Cu-Trainingsset trainiert. Der Adam-Optimierer wurde verwendet, um den mittleren Quadratfehler (MSE) zu minimieren. Die kritische Technik des frühen Stopps, überwacht auf dem Cu-Validierungsset, wurde einheitlich angewendet. Dies stellte sicher, dass das Training am Punkt der optimalen Verallgemeinerung auf unsichtbaren Cu-Daten abgeschlossen wurde, wodurch verhindert wurde, dass Modelle überangepasst auf das Training von Rauschen übergingen.
| Parameterkategorie | Spezifikation / Wert | Beschreibung |
| Kernaufgabe & Daten |
| Vorhersageziel | Preis für den nächsten Tag | Standard-Ein-Schritt-voraus-Vorhersage. |
| Eingabefensterlänge (L) | 60 Handelstage | Balanciert ausreichend historischen Kontext mit Modellkomplexität und Trainingsstabilität. |
| Modellentwicklung |
| Trainingsset (nur Cu) | Erste 80 % (~2081 Beobachtungen) | Verwendet zum Lernen von Modellparametern über Backpropagation. |
| Validierungsmenge (nur Cu) | Nächste 10 % (~260 Beobachtungen) | Verwendet für Hyperparameter-Abstimmung und frühes Stoppen; Entscheidend zur Verhinderung von Überanpassungen. |
| Testset (nur Cu) | Endgültige 10 % (~260 Beobachtungen) | Abschließende, vorgehaltene Bewertung der In-Sample-(Cu)-Leistung. |
| Modellarchitektur |
| RNN Versteckte Einheiten | 128 | Bietet eine angemessene Vertretungskapazität; in allen RNN-basierten Modellen konstant gehalten. |
| CNN-Filter | 64 | Anzahl der Merkmalskarten für CNN-Schichten in hybriden Modellen. |
| Ausbildungsverfahren |
| Optimierer | Adam | Adaptiver Lernrate-Optimierer für stabile und effiziente Konvergenz. |
| Anfängliche Lernrate | 1 × 10⁻³ | Standard-Startrate für Adam. |
| Verlustfunktion | Mittlerer quadratischer Fehler (MSE) | Standard für Regression |
| Chargengröße | 32 | Effizientes Mini-Batch-Training. |
| Maximale Epochen | 80 | Obere Grenze für Trainings-Iterationen. |
| Frühes Aufhören von Geduld | 10 Epochen | Das Training stoppt, wenn sich der Validierungsverlust für 20 aufeinanderfolgende Epochen nicht verbessert; Modellgewichte aus der besten Epoche werden restauriert. |
| Bewertung und Validierung |
| Primäre Kennzahlen | MAE, RMSE, R² | Bieten Sie ergänzende Sichtweisen auf erklärte Fehlergrößen und Varianz an. |
| Generalisierbarkeitstest | Vorhersage auf vollständigen Al- und Zn-Serien (jeweils 2602 Beobachtungen) | Die Modelle werden nach dem Cu-Training eingefroren. Dies ist ein reiner, strenger Test außerhalb der Stichprobe mit völlig anderen Produkten. |
| Ablationsdesign | GRU → BiGRU → BiGRU–Aufmerksamkeit → CNN–BiGRU–Aufmerksamkeit | Isoliert systematisch die Auswirkungen des Hinzufügens von bidirektionaler Komponenten, Aufmerksamkeit und CNN. |
Tabelle 1: Wichtige experimentelle Parameter und Konfiguration. Zusammenfassung des experimentellen Setups, das auf alle Modelle angewendet wurde, einschließlich Datenpartitionierung, Modellarchitekturparametern, Trainingseinstellungen und Evaluierungsmetriken.
Um den Beitrag jeder architektonischen Komponente zu dekonstruieren, wurde eine strukturierte Ablationsstudie entwickelt. Ausgehend vom best-performenden Baseline (GRU) wurde eine progressive "Komplexitätskette" konstruiert. Abbildung 2 zeigt diese Komplexitätskette visuell und veranschaulicht die schrittweise Addition von Komponenten. Dieser schrittweise Ansatz erlaubt es, jede Änderung der Prognoseleistung direkt auf die inkrementelle Hinzufügung von bidirektionaler Fähigkeit, dem Aufmerksamkeitsmechanismus und schließlich der Schicht des konvolutionalen neuronalen Netzwerks zurückzuführen. Die Leistungskennzahlen an jedem Knoten dieser Kette liefern klare empirische Belege für den Wert oder Nachteil jeder Komplexitätskomponente für die spezifische Aufgabe der Metallpreisprognose.

Abbildung 2: Komplexitätskette verwendet in der Ablationsstudie. Das Diagramm zeigt die schrittweise Hinzufügung architektonischer Komponenten, die von GRU zu BiGRU, BiGRU–Attention und CNN–BiGRU–Attention übergehen. Diese Sequenz stellt die systematische Erhöhung der Modellkomplexität dar, die verwendet wird, um die Auswirkungen jeder Komponente auf die Prognoseleistung zu bewerten. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Leistungsbewertungsmetriken
Die Modellleistung wurde streng mit drei Standardregressionsmetriken quantifiziert, die ergänzende Einblicke in die Vorhersagegenauigkeit und Erklärungskraft boten.
Mittlerer absoluter Fehler (MAE)
Misst die durchschnittliche Fehlergröße und bietet eine robuste und leicht interpretierbare Abweichungsskala.
(8)
Wurzelmittlerer Quadratfehler (RMSE)
Betont größere Fehler aufgrund des Quadrierens, was sie empfindlicher gegenüber Ausreißern und großen Fehlern macht.
(9)
Bestimmungskoeffizient (R2)
Repräsentiert den Anteil der Varianz in der Zielvariable, der aus dem Modell vorhersehbar ist.
(10)
wobei
der Mittelwert der wahren Werte ist. Ein R2-Wert näher an 1 zeigt ein Modell an, das den Großteil der Varianz in den Daten erklärt. Die Gleichungen 8–10 sind Standardregressionsmetriken48. Die Bewertung wurde in zwei getrennten, sequentiellen Phasen durchgeführt, um die Leistung der Benchmarks innerhalb der Stichprobe und die Generalisierbarkeit außerhalb der Stichprobe separat zu bewerten. (1) Phase 1 (Primärer Benchmark): Alle dreizehn Modelle wurden nach Training und frühzeitigem Stoppen der Cu-Daten auf dem zurückgehaltenen Cu-Testsatz bewertet. (2) Phase 2 (Generalisierbarkeitstest): Exakt dieselben Modelle, deren Parameter eingefroren waren, wurden eingesetzt, um Prognosen für die vollständige, unabhängige Preisreihe von Al und Zn zu erstellen. Es wurde kein Nachtraining oder Anpassung durchgeführt.
Reproduzierbarkeit: Detaillierte experimentelle Settings
Die täglichen Spotpreise (CNY/Tonne) für Grade A Cu, Al und Zn wurden von der SMM-Plattform (https://www.smm.cn/ übernommen und erstrecken sich vom 5. Januar 2015 bis zum 12. September 2025. Die Roh- und Verarbeitungsdaten sind in einem öffentlichen Repository verfügbar (DOI: 10.5281/zenodo.19976985). Die Datendatei enthält die Spalten Datum, Cu, Al und Zn. Daten werden in das Datum-Zeit-Format umgewandelt und in aufsteigender Reihenfolge sortiert. Fehlende Werte werden durch Vorwärtsfüllung gefolgt von Rückwärtsfüllung behandelt. Die Funktionen werden standardisiert, indem ein Z-Score-Skalierer verwendet wird, der nur auf dem Trainingssatz angepasst wird (Mittelwert μj, Standardabweichung σj ,
); dieselben und werden auf die Validierungs- und Testsätze angewendet, ohne dass sie neu angepasst werden. Die Zielvariable (Cu, Al oder Zn) wird separat mit ihren jeweiligen Trainingsmengenstatistiken skaliert.
Eingabe-Ausgabe-Sequenzen werden mithilfe eines gleitenden Fensters mit Eingabelänge L = 30 Handelstagen und Prognosehorizont h = 1 (Vorhersage für den nächsten Tag) konstruiert. Für einen Zielindex (0 = Cu, 1 = Al, 2 = Zn) ist jede Probe definiert als Xi = V[ t - L : t, : ] (Form 30 × 3) und yi = V[ t + h , k] (Skalar). Es wird kein Mischen angewendet, um die zeitliche Ordnung zu bewahren. Der Datensatz ist chronologisch ohne Zufälligkeit aufgeteilt: Das Training umfasst die Indizes 0–2080 (2.081 Beobachtungen, 80 %), Validierungsindizes 2081–2340 (260 Beobachtungen, 10 %) und Testindizes 2341–2601 (261 Beobachtungen, 10 %). Die entsprechenden Datumsgrenzen sind vom 5. Januar 2015 bis zum 31. Juli 2023 (Training), vom 1. August 2023 bis zum 19. Oktober 2023 (Validierung) und vom 20. Oktober 2023 bis zum 12. September 2025 (Test); Die Datei im Repository liefert genaue Details.
Zufällige Seeds sind wie folgt festgelegt: Hauptexperiment-Seed = 42, und Python-, NumPy- und TensorFlow-Seeds sind alle auf 42 gesetzt. Die Gewichtsinitialisierung verwendet Glorot-Uniform für Eingabekerne, orthogonal für rekurrente Kernel und Nullstellen für Verzerrungen. Die Softwareumgebung besteht aus Python 3.10.19, TensorFlow 2.20.0/Keras, NumPy 1.26.4, pandas 2.3.3, scikit-learn 1.7.2 und Matplotlib 3.10.6. Die Experimente wurden auf einem Windows-11-PC mit einem Intel Core i7 (2,20 GHz) und 32 GB RAM durchgeführt; es wurde keine GPU verwendet.
Der Adam-Optimierer wird verwendet mit learning_rate = 1×10-3, β1 = 0,9, β2 = 0,999
und weight_decay = 0. Die Verlustfunktion ist MSE. Ein ReduceLROnPlateau-Scheduler überwacht den Validierungsverlust mit dem Faktor 0,5, Geduld 5 und einer Mindestlernrate von 1 × 10-5. Frühes Stoppen erfolgt mit Monitor = val_loss, Geduld = 10, restore_best_weights = Wahr und min_delta = 0. Jede Trainingsepoche besteht aus einem Vorwärtspass auf dem Trainingsbatch, MSE-Verlustberechnung, Backpropagation und Adam-Parameteraktualisierung. Nach jeder Epoche wird der Validierungsverlust berechnet; Frühzeitiges Stoppen und die Reduktion der Lernrate werden auf Basis dieses Wertes ausgelöst. Das Modell mit dem geringsten Validierungsverlust wird zum Testen wiederhergestellt. Die Chargengröße beträgt 32, und die Proben werden in chronologischer Reihenfolge ohne Mischen (Mischen = Falsch) eingegeben.
Für CNN-Hybridmodelle wird eine Conv1D-Schicht mit 64 Filtern verwendet, kernel_size = 3, Stride = 1, Padding = 'gleich' und Aktivierung der rektifizierten linearen Einheit (ReLU), gefolgt von MaxPooling1D(pool_size = 2) und Dropout(0,15). In aufmerksamkeitsaugmentierten Modellen gibt das rekurrente neuronale Netzwerk die vollständige verborgene Sequenz H mit Form B × T × C zurück. Eine dichte Schicht mit einer Einheit erzeugt einen Score, und ein Softmax wandelt diese Werte über die Zeit in Aufmerksamkeitsgewichte um, wobei der Kontextvektor definiert ist als c = ∑ t αtht. Darauf folgen eine dichte Schicht mit 64 Einheiten und ReLU-Aktivierung, Dropout (0,15) und die Output-Dichteschicht. Bidirektionale Modelle verketten vorwärts und rückwärts versteckte Zustände (jeweils 64 Einheiten), was zu 128 Dimensionen führt; wenn Aufmerksamkeit verwendet wird, bewahrt return_sequences = Wahr die vollständige Folge (B × T × 128).
Die Bewertung verwendet direkte Ein-Schritt-voraus-Vorhersage (nicht-rekursiv). Alle Vorhersagen werden invers auf die ursprüngliche Preisskala umgewandelt, bevor MAE, RMSE und R2 auf dieser Skala berechnet werden. Für Generalisierbarkeitstests auf Al und Zn wird der auf Cu montierte Eingangsskaler ohne Modifikation wiederverwendet, während jedes Zielmetall seinen eigenen Zielskaler auf seinen eigenen Trainingszielen besitzt. In der Ablationsstudie bleiben alle nicht-architektonischen Parameter (Daten, Split, Skalierung, Zufallsseed, Anzahl der Epochen, Chargengröße, Optimierer, Lernrate, Verlustfunktion, frühes Stoppen, Scheduler, Dropout) über die gesamte Kette hinweg identisch; Nur die Architektur ändert sich. Der vollständige Quellcode und die Replikationsanleitung sind öffentlich bei Zenodo verfügbar (10.5281/zenodo.19976985). Alle Figuren wurden mit Matplotlib 3.10.6 und dem bereitgestellten Skript generiert; Ausgaben werden als PDF, SVG und hochauflösende PNG (600 dpi) gespeichert. Alle Modelle wurden mit maximal 80 Epochen trainiert. Frühes Stoppen mit Geduld = 10 (überwacht bei Validierungsverlust) wurde für jedes Modell vor Erreichen der Epochengrenze ausgelöst. Zum Beispiel stoppte das GRU-Modell bei Epoche 37 (beste Epoche 27, bester Validierungsverlust 0,0040), während der komplexeste Hybrid CNN–BiLSTM–Attention bei Epoche 23 endete (beste Epoche 13, bester Validierungsverlust 0,0072). Die vollständige Liste der gestoppten Epochen, besten Epochen und besten Validierungsverluste für alle 13 Modelle ist im Zenodo-Repository verfügbar und gewährleistet volle Transparenz und Reproduzierbarkeit, ohne den Haupttext mit einer Tabelle zu überladen. Die theoretischen Diskussionen (Lipschitz-Grenzen, Stichprobenkomplexität, Rademacher-Komplexität, Bias-Varianz-Dekomposition, Aufmerksamkeitsentropie und gegenseitige Informationen) sind konzeptionelle Erklärungen für die empirischen Ergebnisse und verändern weder das Trainingsziel noch die Modellimplementierung. Schließlich wurden zur Bewertung der Stabilität die Hauptexperimente mit fünf zufälligen Samen (1, 7, 21, 42 und 2024) wiederholt; der Mittelwert und die Standardabweichung des RMSE über diese Durchläufe hinweg werden im Results-Abschnitt berichtet, wo der GRU einen wettbewerbsfähigen Mittelwert RMSE mit geringer Varianz beibehielt und stabile Leistung über zufällige Initialisierungen unterstützte.