$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Cette étude n’a pas impliqué de participants humains ni de sujets vertébrés animaux. Toutes les données utilisées sont des séries de prix de matières premières publiquement accessibles par la SMM, qui ne nécessitent pas d’approbation éthique. Par conséquent, aucune approbation éthique n’a été sollicitée ni requise pour cette recherche.
Cette section décrit la conception complète et rigoureuse de la recherche mise en œuvre pour tester empiriquement l’hypothèse centrale. Il fournit une exposition détaillée de la formulation mathématique et des spécificités architecturales des treize modèles d’apprentissage profond évalués, du protocole d’entraînement précis et des métriques formelles d’évaluation. Le flux de travail méthodologique global est résumé visuellement à la Figure 1.

Figure 1 : Aperçu schématique de la méthodologie de recherche. Le diagramme illustre l’ensemble du pipeline expérimental, incluant le partitionnement des données, l’entraînement du modèle exclusivement sur la série de prix Cu, l’évaluation sur l’ensemble de test Cu, et la validation hors échantillon sur les séries indépendantes Al et Zn. La boucle de rétroaction pointillée indique les expériences d’ablation structurée menées pour analyser la contribution des composants architecturaux individuels. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
Le diagramme illustre l’ensemble du pipeline expérimental. Le processus commence par l’utilisation exclusive de la série de prix Cu pour le développement de modèles. Cette série est répartie chronologiquement en ensembles d’entraînement (80 %), de validation (10 %) et de test (10 %). Treize architectures d’apprentissage profond distinctes sont ensuite entraînées et optimisées en hyperparamètre uniquement sur les données d’entraînement Cu, avec un arrêt précoce surveillé via le jeu de validation. Le benchmark principal est l’évaluation de ces modèles sur l’ensemble de test Cu maintenu. De manière cruciale, pour évaluer la généralisabilité, les mêmes modèles entraînés sont appliqués sans modification pour prévoir la série de prix totalement indépendante Al et Zn, représentant un test strict hors échantillon. Enfin, des expériences d’ablation structurée (la boucle de rétroaction pointillée) sont menées pour décomposer et analyser la contribution à la performance des composants architecturaux individuels (par exemple, l’attention, le traitement bidirectionnel et les couches convolutionnelles).
Architectures de modèles et formulation mathématique
Nous avons conçu et mis en œuvre un spectre de 13 modèles DL, augmentant systématiquement en complexité architecturale, allant de simples réseaux récurrents à des hybrides multicomposantes sophistiqués. Tous les modèles partagent le même objectif fondamental : apprendre une application
à partir d’une fenêtre de prix historique Xt = [Pt-L, Pt-L+1,...,Pt-1] de longueur L = 30 vers le prix suivant yt = P t.
Les familles de modèles sont définies comme suit :
Modèles récurrents de référence
GRU : Un réseau récurrent rationalisé utilisant des portes de mise à jour (zt) et de réinitialisation (rt) pour moduler le flux d’information. L’état caché ht est calculé comme suit :
(1)
(2)
(3)
(4)
où X est l’activation sigmoïde,
désigne le produit de Hadamard, et xt est l’entrée au temps t . L’état caché final hL est passé à travers une couche de sortie linéaire. Les équations 1 à 4 sont adaptées de Cho et al.14.
LSTM : Utilise les portes d’entrée (it), d’oubli (ft) et de sortie (ot) pour maintenir un état cellulaire (Ct), offrant un contrôle plus explicite sur la mémoire à long terme.
Modèles bidirectionnels (BiGRU et BiLSTM)
Ces modèles intègrent deux couches récurrentes distinctes qui traitent la séquence en direction directe et inverse. La dernière représentation cachée à chaque étape temporelle est la concaténation
, capturant théoriquement des informations contextuelles du passé et du futur dans la fenêtre d’entrée fixe.
Modèles à Attention Augmentée (GRU–Attention et LSTM–Attention)
Un mécanisme d’attention additive est appliqué à la suite d’états cachés H = [h1,h 2,...,h L] produite par la couche récurrente finale. Le vecteur de contexte est défini comme une somme pondérée :
(5)
(6)
(7)
Ici, α i représente le poids d’attention attribué au i-ième pas de temps historique. Le vecteur de contexte c, qui encapsule un résumé adaptatif de l’histoire pertinente, est transmis à la couche finale de prédiction. Les équations 5 à 7 sont adaptées de Bello et al.47.
CNN–Modèles hybrides (CNN–GRU et CNN–LSTM)
Une couche CNN unidimensionnelle avec activation linéaire redressée (ReLU) est précédée à la couche récurrente.
Modèles hybrides complexes
Ces architectures combinent plusieurs composants (par exemple, CNN–BiGRU–Attention, CNN–BiLSTM–Attention). Ils représentent l’état de l’art en termes de complexité, visant à intégrer l’extraction locale de motifs (CNN), la modélisation contextuelle bidirectionnelle et la pondération temporelle adaptative (attention) dans un seul cadre.
Tous les modèles étaient configurés avec des dimensions d’état caché cohérentes (128 unités pour les couches récurrentes et 64 filtres pour les couches CNN) et étaient finalisés avec une seule couche de sortie linéaire. Cette conception contrôlée garantit que les différences de performance sont attribuables à des choix architecturaux plutôt qu’à des écarts dans l’ajustement de la capacité du modèle. Le nombre de paramètres entraînables a donc considérablement augmenté le long de ce spectre.
Protocole d’entraînement, hyperparamètres et conception de l’étude d’ablation
Le tableau 1 résume la configuration expérimentale unifiée et rigoureuse appliquée à tous les treize modèles afin d’assurer une comparaison équitable et d’atténuer le sur-ajustement. Tous les modèles étaient entraînés de zéro en utilisant uniquement l’ensemble d’entraînement cu Cu. L’optimiseur Adam a été utilisé pour minimiser la perte d’erreur quadratique moyenne (MSE). La technique critique d’arrêt précoce, survue sur le kit de validation Cu, a été appliquée de manière uniforme. Cela garantissait que l’entraînement se terminait au point de généralisation optimale sur des données cu non visibles, empêchant ainsi les modèles de surajuster au bruit d’entraînement.
| Catégorie de paramètres | Spécification / Valeur | Description |
| Tâches principales & Données |
| Objectif de prévision | Prix du jour suivant | Prédiction standard d’un pas en avant. |
| Longueur de la fenêtre d’entrée (L) | 60 jours de bourse | Équilibre un contexte historique suffisant avec la complexité du modèle et la stabilité de l’entraînement. |
| Développement de modèles |
| Ensemble d’entraînement (Cu-only) | Premiers 80 % (~2081 observations) | Utilisé pour apprendre les paramètres du modèle via rétropropagation. |
| Ensemble de validation (Cu-only) | Les 10 % suivants (~260 observations) | Utilisé pour l’ajustement des hyperparamètres et l’arrêt précoce ; Crucial pour éviter le sur-ajustement. |
| Ensemble de tests (Cu-only Set) | 10 % final (~260 observations) | Évaluation finale, tenue, de la performance dans l’échantillon (Cu). |
| Architecture du modèle |
| Unités cachées de la RNN | 128 | Assure une capacité représentative adéquate ; maintenu constant sur tous les modèles basés sur RNN. |
| Filtres CNN | 64 | Nombre de cartes de caractéristiques pour les couches CNN dans les modèles hybrides. |
| Procédure d’entraînement |
| Optimiseur | Adam | Optimiseur adaptatif du taux d’apprentissage pour une convergence stable et efficace. |
| Taux d’apprentissage initial | 1 × 10⁻³ | Taux de départ standard pour Adam. |
| Fonction de perte | Erreur quadratique moyenne (MSE) | Standard pour la régression |
| Taille du lot | 32 | Une formation en mini-lots efficace. |
| Époques maximales | 80 | Limite supérieure pour les itérations d’entraînement. |
| Patience d’arrêt précoce | 10 époques | L’entraînement s’arrête si la perte de validation ne s’améliore pas pendant 20 époques consécutives ; Les poids des modèles de la meilleure époque sont restaurés. |
| Évaluation et validation |
| Indicateurs principaux | MAE, RMSE, R² | Fournir des vues complémentaires de la magnitude et de la variance de l’erreur expliquées. |
| Test de généralisabilité | Prévisions sur les séries complètes Al et Zn (2602 observations chacune) | Les mannequins sont gelés après la formation Cu. Il s’agit d’un test pur et strict hors échantillon sur des produits totalement différents. |
| Conception d’ablation | GRU → BiGRU → BiGRU–Attention → CNN–BiGRU–Attention | Isole systématiquement l’impact de l’ajout de composants bidirectionnels, d’attention et CNN. |
Tableau 1 : Paramètres expérimentaux clés et configuration. Résumé de la configuration expérimentale appliquée à tous les modèles, y compris le partitionnement des données, les paramètres d’architecture du modèle, les paramètres d’entraînement et les métriques d’évaluation.
Pour déconstruire la contribution de chaque composante architecturale, une étude d’ablation structurée a été conçue. En partant de la base la plus performante (GRU), une « chaîne de complexité » progressive a été construite. La figure 2 décrit visuellement cette chaîne de complexité, illustrant l’ajout par étapes des composants. Cette approche progressive permet d’attribuer directement toute modification de la performance de prévision à l’ajout incrémental de la bidirectionnalité, du mécanisme d’attention, et enfin de la couche convolutionnelle du réseau de neurones. Les indicateurs de performance à chaque nœud de cette chaîne fournissent des preuves empiriques claires sur la valeur ou l’inconvénient de chaque composante de complexité pour la tâche spécifique de prévision des prix des métaux.

Figure 2 : Chaîne de complexité utilisée dans l’étude d’ablation. Le diagramme illustre l’ajout progressive de composants architecturaux, allant de GRU à BiGRU, BiGRU–Attention et CNN–BiGRU–Attention. Cette séquence représente l’augmentation systématique de la complexité du modèle utilisée pour évaluer l’impact de chaque composant sur la performance des prévisions. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
Indicateurs d’évaluation de la performance
La performance du modèle a été rigoureusement quantifiée à l’aide de trois métriques de régression standard, offrant des perspectives complémentaires sur la précision des prédictions et le pouvoir explicatif.
Erreur absolue moyenne (MAE)
Mesure l’ampleur moyenne des erreurs, fournissant une échelle robuste et facilement interprétable de déviation.
(8)
Erreur quadratique moyenne de la racine (RMSE)
Met l’accent sur les erreurs plus importantes dues à l’opération de carré, ce qui la rend plus sensible aux valeurs aberrantes et aux grosses erreurs.
(9)
Coefficient de détermination (R2)
Représente la proportion de variance dans la variable cible qui est prévisible à partir du modèle.
(10)
où
est la moyenne des vraies valeurs. Une valeur deR 2 plus proche de 1 indique un modèle qui explique la plupart des variances dans les données. Les équations 8 à 10 sont les mesures de régressionstandard 48. L’évaluation a été réalisée en deux phases distinctes et séquentielles afin d’évaluer séparément la performance des benchmarks dans l’échantillon et la généralisabilité hors échantillon. (1) Phase 1 (Benchmark primaire) : Les treize modèles, après entraînement et arrêt précoce des données Cu, ont été évalués sur l’ensemble de test Cu maintenu. (2) Phase 2 (Test de généralisabilité) : Les mêmes modèles, avec leurs paramètres gelés, ont été déployés pour générer des prévisions pour la série complète et indépendante de prix Al et Zn. Aucune réentraînement ni adaptation n’a été effectué.
Reproductibilité : Paramètres expérimentaux détaillés
Les prix au comptant quotidiens (CNY/tonne) pour les Grade A Cu, Al et Zn ont été obtenus sur la plateforme publique de la SMM (https://www.smm.cn/), du 5 janvier 2015 au 12 septembre 2025. Les données brutes et traitées sont disponibles dans un dépôt public (DOI : 10.5281/zenodo.19976985). Le fichier de données contient les colonnes date, Cu, Al et Zn. Les dates sont converties au format datetime et triées par ordre croissant. Les valeurs manquantes sont gérées par un remplissage en avant suivi d’un remplissage vers l’arrière. Les caractéristiques sont standardisées à l’aide d’un scaleur z-score installé uniquement sur l’ensemble d’entraînement (moyenne μj, écart-type σ j, ) ;
les mêmes et sont appliquées aux ensembles de validation et de test sans réajustement. La variable cible (Cu, Al ou Zn) est mise à l’échelle séparément à l’aide de ses propres statistiques d’ensemble d’entraînement respectives.
Les séquences entrée-sortie sont construites à l’aide d’une fenêtre glissante avec la longueur d’entrée L = 30 jours de négociation et l’horizon de prévision h = 1 (prédiction du lendemain). Pour un indice cible (0 = Cu, 1 = Al, 2 = Zn), chaque échantillon est défini comme Xi = V[ t - L : t, : ] (forme 30 × 3) et yi = V[ t + h, k] (scalaire). Aucun mélange n’est appliqué pour préserver l’ordre temporel. L’ensemble de données est divisé chronologiquement sans aléatoire : l’entraînement comprend les indices 0–2080 (2 081 observations, 80 %), les indices de validation 2081–2340 (260 observations, 10 %) et les indices de test 2341–2601 (261 observations, 10 %). Les limites de dates correspondantes sont du 5 janvier 2015 au 31 juillet 2023 (formation), du 1er août 2023 au 19 octobre 2023 (validation), et du 20 octobre 2023 au 12 septembre 2025 (test) ; Le fichier dans le dépôt fournit des détails exacts.
Les graines aléatoires sont fixées comme suit : graine de l’expérience principale = 42, et graines Python, NumPy et TensorFlow sont toutes réglées à 42. L’initialisation des poids utilise un Glorot uniforme pour les noyaux d’entrée, orthogonal pour les noyaux récurrents, et zéros pour les polarisations. L’environnement logiciel se compose de Python 3.10.19, TensorFlow 2.20.0/Keras, NumPy 1.26.4, pandas 2.3.3, scikit-learn 1.7.2 et Matplotlib 3.10.6. Les expériences ont été menées sur un PC Windows 11 équipé d’un Intel Core i7 (2,20 GHz) et de 32 Go de RAM ; aucun GPU n’a été utilisé.
L’optimiseur d’Adam est utilisé avec learning_rate = 1×10-3, β1 = 0,9, β2 = 0,999,
, et weight_decay = 0. La fonction de perte est MSE. Un planificateur ReduceLROnPlateau surveille la perte de validation avec un facteur de 0,5, une patience de 5 et un taux d’apprentissage minimum de 1 × 10-5. L’arrêt précoce s’applique avec moniteur = val_loss, patience = 10, restore_best_weights = Vrai, et min_delta = 0. Chaque époque d’entraînement comprend un passage en avant sur le lot d’entraînement, le calcul de perte MSE, la rétropropagation et la mise à jour des paramètres Adam. Après chaque époque, une perte de validation est calculée ; L’arrêt précoce et la réduction du taux d’apprentissage sont déclenchés en fonction de cette valeur. Le modèle avec la perte de validation la plus faible est restauré pour les tests. La taille du lot est de 32, et les échantillons sont distribués dans l’ordre chronologique sans mélange (mélange = Faux).
Pour les modèles hybrides CNN, une couche Conv1D avec 64 filtres, kernel_size = 3, foulée = 1, bourrage = « identique » et activation unité linéaire rectificée (ReLU) est utilisée, suivie de MaxPooling1D(pool_size = 2) et Dropout(0,15). Dans les modèles à attention augmentée, le réseau neuronal récurrent renvoie la suite cachée complète H avec la forme B × T × C. Une couche dense d’une unité produit un score, et un softmax convertit ces scores en pondérations d’attention , avec le vecteur de contexte défini comme c = ∑t αt ht. Cela est suivi par une couche dense avec 64 unités et l’activation ReLU, Dropout(0,15), et la couche dense de sortie. Les modèles bidirectionnels concaténent des états cachés avant et arrière (chacun de 64 unités), ce qui donne 128 dimensions ; lorsque l’attention est utilisée, return_sequences = Vrai préserve la suite complète (B × T × 128).
L’évaluation utilise une prédiction directe à un pas d’avance (non récursive). Toutes les prédictions sont transformées inversement à l’échelle de prix originale avant de calculer MAE, RMSE et R 2 sur cette échelle. Pour les tests de généralisabilité sur Al et Zn, le scaleur d’entrée monté sur Cu est réutilisé sans modification, tandis que chaque métal cible possède son propre scaleur de cible installé sur ses propres cibles d’entraînement. Dans l’étude d’ablation, tous les paramètres non architecturaux (données, division, mise à l’échelle, graine aléatoire, nombre d’époques, taille du lot, optimiseur, taux d’apprentissage, fonction de perte, arrêt anticipé, planificateur, abandon) sont maintenus identiques sur toute la chaîne ; seule l’architecture change. Le code source complet et les instructions de réplication sont disponibles publiquement sur Zenodo (10.5281/zenodo.19976985). Toutes les figures ont été générées avec Matplotlib 3.10.6 avec le script fourni ; les sorties sont enregistrées en PDF, SVG et PNG haute résolution (600 dpi). Tous les modèles ont été entraînés avec un maximum de 80 époques. Un arrêt précoce avec patience = 10 (surveillé lors de la perte de validation) était déclenché pour chaque modèle avant d’atteindre la limite d’époque. Par exemple, le modèle GRU s’est arrêté à l’époque 37 (meilleure époque 27, meilleure perte de validation 0,0040), tandis que le modèle hybride le plus complexe CNN–BiLSTM–Attention s’est arrêté à l’époque 23 (meilleure époque 13, meilleure perte de validation 0,0072). La liste complète des époques arrêtées, meilleures époques et meilleures pertes de validation pour les 13 modèles est disponible dans le dépôt Zenodo, garantissant une transparence et une reproductibilité totales sans surcharger le texte principal avec un tableau. Les discussions théoriques (bornes de Lipschitz, complexité de l’échantillon, complexité de Rademacher, décomposition biais–variance, entropie de l’attention et information mutuelle) sont des explications conceptuelles des résultats empiriques et ne modifient pas l’objectif d’entraînement ni la mise en œuvre du modèle. Enfin, pour évaluer la stabilité, les expériences principales ont été répétées avec cinq graines aléatoires (1, 7, 21, 42 et 2024) ; la moyenne et l’écart-type du RMSE sur ces séries sont rapportées dans la section Résultats, où le GRU a maintenu une moyenne RMSE compétitive avec faible variance, soutenant une performance stable sur des initialisations aléatoires.