$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Questo studio non ha coinvolto partecipanti umani né soggetti animali vertebrati. Tutti i dati utilizzati sono serie di prezzi delle materie prime disponibili pubblicamente dalla SMM, che non richiedono approvazione etica. Pertanto, non è stata richiesta né richiesta alcuna approvazione etica per questa ricerca.
Questa sezione delinea il disegno di ricerca completo e rigoroso implementato per testare empiricamente l'ipotesi centrale. Fornisce un'esposizione dettagliata della formulazione matematica e delle specifiche architettoniche dei tredici modelli di deep learning valutati, del protocollo di addestramento preciso e delle metriche formali di valutazione. Il flusso di lavoro metodologico generale è riassunto visivamente nella Figura 1.

Figura 1: Panoramica schematica della metodologia di ricerca. Il diagramma illustra l'intera pipeline sperimentale, inclusa la partizione dei dati, l'addestramento del modello esclusivamente sulla serie di prezzi Cu, la valutazione sul set di test Cu e la validazione fuori campione su serie indipendenti Al e Zn. Il ciclo di retroazione tratteggiata indica gli esperimenti di ablazione strutturata condotti per analizzare il contributo dei singoli componenti architettonici. Clicca qui per visualizzare una versione più grande di questa figura.
Il diagramma illustra l'intera pipeline sperimentale. Il processo inizia con l'uso esclusivo della serie prezzi Cu per lo sviluppo dei modelli. Questa serie è suddivisa cronologicamente in set di addestramento (80%), validazione (10%) e set di prova (10%). Tredici architetture di deep learning distinte vengono quindi addestrate e ottimizzate per iperparametri esclusivamente sui dati di addestramento Cu, con il early stop monitorato tramite il set di validazione. Il punto di riferimento principale è la valutazione di questi modelli sul set di test Cu tenuto in test. Fondamentalemente, per valutare la generalizzabilità, gli stessi modelli addestrati vengono applicati senza modifiche per prevedere la serie di prezzi completamente indipendente di Al e Zn, rappresentando un rigoroso test fuori campione. Infine, vengono condotti esperimenti di ablazione strutturata (il ciclo di retroazione tratteggiata) per decomporre e analizzare il contributo alle prestazioni dei singoli componenti architettonici (ad esempio, attenzione, elaborazione bidirezionale e strati convoluzionali).
Architetture di modelli e formulazione matematica
Abbiamo progettato e implementato uno spettro di 13 modelli DL, aumentando sistematicamente la complessità architettonica da semplici reti ricorrenti a sofisticati ibridi multicomponenti. Tutti i modelli condividono lo stesso obiettivo fondamentale: apprendere una mappa
da una finestra storica di prezzo Xt = [Pt-L, Pt-L+1,...,Pt-1] di lunghezza L = 30 fino al prezzo successivo yt = Pt.
Le famiglie di modelli sono definite come segue:
Modelli ricorrenti di base
GRU: Una rete ricorrente snellita che impiega porte di aggiornamento (zt) e reset (rt) per modulare il flusso di informazioni. Lo stato nascosto ht si calcola come:
(1)
(2)
(3)
(4)
dove X è l'attivazione sigmoide,
indica il prodotto di Hadamard, e xt è l'ingresso al tempo t . Lo stato nascosto finale hL viene passato attraverso uno strato di uscita lineare. Le equazioni 1–4 sono adattate da Cho et al.14.
LSTM: Utilizza porte di input (it), forget (ft) e output (ot) per mantenere uno stato della cella (Ct), fornendo un controllo più esplicito sulla memoria a lungo termine.
Modelli bidirezionali (BiGRU e BiLSTM)
Questi modelli incorporano due strati ricorrenti separati che elaborano la sequenza in direzioni avanti e indietro. La rappresentazione nascosta finale a ogni passo temporale è la concatenazione
, che teoricamente cattura informazioni contestuali sia dal passato che dal futuro all'interno della finestra di input fissa.
Modelli Potenziati dall'Attenzione (GRU–Attenzione e LSTM–Attenzione)
Un meccanismo di attenzione additiva viene applicato alla sequenza di stati nascosti H = [h 1,h 2,...,h L] prodotta dallo strato finale ricorrente. Il vettore contesto è definito come una somma pesata:
(5)
(6)
(7)
Qui, α i rappresenta il peso di attenzione assegnato al i-esimo passo storico del tempo. Il vettore contesto c, che racchiude un riassunto adattivo della storia rilevante, viene inviato al livello finale di previsione. Le equazioni 5–7 sono adattate da Bello et al.47.
CNN–Modelli ibridi (CNN–GRU e CNN–LSTM)
Un livello CNN unidimensionale con attivazione unitaria lineare rettificata (ReLU) viene preceduto allo strato ricorrente.
Modelli ibridi complessi
Queste architetture combinano più componenti (ad esempio, CNN–BiGRU–Attention, CNN–BiLSTM–Attention). Rappresentano lo stato dell'arte in termini di complessità, mirando a integrare l'estrazione locale di pattern (CNN), la modellazione contestuale bidirezionale e la ponderazione temporale adattativa (attenzione) in un unico quadro.
Tutti i modelli erano configurati con dimensioni di stato nascosti costanti (128 unità per i livelli ricorrenti e 64 filtri per i livelli CNN) e sono stati completati con un unico livello di uscita lineare. Questo design controllato garantisce che le differenze di prestazioni siano attribuibili a scelte architettoniche piuttosto che a discrepanze nella messa a punto della capacità del modello. Il numero di parametri addestrabili è quindi aumentato notevolmente lungo questo spettro.
Protocollo di addestramento, iperparametri e progettazione dello studio di ablazione
La Tabella 1 riassume la configurazione sperimentale unificata e rigorosa applicata a tutti e tredici i modelli per garantire un confronto equo e mitigare l'overfitting. Tutti i modelli venivano addestrati da zero utilizzando solo il set di addestramento designato Cu. L'ottimizzatore Adam veniva utilizzato per minimizzare la perdita di errore quadratico medio (MSE). La tecnica critica dell'arresto precoce, monitorata sul set di validazione del Cu, è stata applicata uniformemente. Ciò garantiva che l'addestramento si concludesse nel punto di generalizzazione ottimale su dati cu non visti, impedendo così ai modelli di sovraadattarsi al rumore di addestramento.
| Categoria dei parametri | Specifiche / Valore | Descrizione |
| Compiti e dati principali |
| Obiettivo di previsione | Prezzo del giorno successivo | La classica previsione a un passo avanti. |
| Lunghezza della finestra di ingresso (L) | 60 giorni di contrattazione | Bilancia un contesto storico sufficiente con la complessità del modello e la stabilità dell'addestramento. |
| Sviluppo del modello |
| Set di addestramento (Solo Cu) | Primo 80% (~2081 osservazioni) | Utilizzato per apprendere i parametri del modello tramite retropropagazione. |
| Set di validazione (Solo Cu) | Il prossimo 10% (~260 osservazioni) | Utilizzato per la regolazione degli iperparametri e il stop anticipato; Fondamentale per prevenire il sovrafitting. |
| Set di test (Solo Cu) | 10% finale (~260 osservazioni) | Valutazione finale, attenta, della performance in campione (Cu). |
| Architettura dei modelli |
| Unità nascoste della RNN | 128 | Fornisce una capacità rappresentativa adeguata; mantenuto costante su tutti i modelli basati su RNN. |
| Filtri CNN | 64 | Numero di feature map per i livelli CNN nei modelli ibridi. |
| Procedura di addestramento |
| Ottimizzatore | Adam | Ottimizzatore adattivo per la velocità di apprendimento per una convergenza stabile ed efficiente. |
| Tasso di apprendimento iniziale | 1 × 10⁻³ | Tasso iniziale standard per Adam. |
| Funzione di perdita | Errore quadratico medio (MSE) | Standard per la regressione |
| Dimensione del lotto | 32 | Addestramento mini-batch efficiente. |
| Epoche Massime | 80 | Limite massimo per le iterazioni di addestramento. |
| Pazienza di fermamento precoce | 10 epoche | L'addestramento si interrompe se la perdita di validazione non migliora per 20 epoche consecutive; I pesi dei modelli dell'epoca migliore sono stati restaurati. |
| Valutazione e Validazione |
| Metriche primarie | MAE, RMSE, R² | Fornire visioni complementari della magnitudo e della varianza dell'errore spiegate. |
| Test di generalizzazione | Previsione sulle serie complete Al e Zn (2602 osservazioni ciascuna) | I modelli vengono congelati dopo l'addestramento Cu. Si tratta di un test puro e rigoroso fuori campione su merci completamente diverse. |
| Progettazione dell'ablazione | GRU → BiGRU → BiGRU–Attenzione → CNN–BiGRU–Attenzione | Isola sistematicamente l'impatto dell'aggiunta di componenti bidirezionali, di attenzione e CNN. |
Tabella 1: Parametri sperimentali chiave e configurazione. Riepilogo della configurazione sperimentale applicata su tutti i modelli, inclusi partizionamento dei dati, parametri dell'architettura del modello, impostazioni di addestramento e metriche di valutazione.
Per decostruire il contributo di ciascun componente architettonico, è stato progettato uno studio strutturato di ablazione. Partendo dalla base di miglior performance (GRU), è stata costruita una progressiva "catena di complessità". La Figura 2 delinea visivamente questa catena di complessità, illustrando l'aggiunta a passo a passo dei componenti. Questo approccio a passo a passo permette l'attribuzione diretta di qualsiasi cambiamento nelle prestazioni di previsione all'aggiunta incrementale della bidirezionalità, al meccanismo di attenzione e infine allo strato della rete neurale convoluzionale. Le metriche di performance in ciascun nodo di questa catena forniscono prove empiriche chiare sul valore o svantaggio di ogni componente di complessità per il compito specifico di previsione dei prezzi dei metalli.

Figura 2: Catena di complessità utilizzata nello studio di ablazione. Il diagramma illustra l'aggiunta graduale di componenti architettonici, proseguendo da GRU a BiGRU, BiGRU–Attenzione e CNN–BiGRU–Attenzione. Questa sequenza rappresenta l'aumento sistematico della complessità del modello utilizzato per valutare l'impatto di ciascun componente sulle prestazioni di previsione. Clicca qui per visualizzare una versione più grande di questa figura.
Metriche di Valutazione delle Prestazioni
Le prestazioni del modello sono state rigorosamente quantificate utilizzando tre metriche standard di regressione, offrendo approfondimenti complementari sull'accuratezza della previsione e sul potere esplicativo.
Errore Assoluto Medio (MAE)
Misura la grandezza media degli errori, fornendo una scala robusta e facilmente interpretabile di deviazione.
(8)
Errore quadratico medio radice (RMSE)
Enfatizza errori maggiori dovuti all'operazione di squadratura, rendendola più sensibile a outlier e errori grandi.
(9)
Coefficiente di determinazione (R2)
Rappresenta la proporzione di varianza nella variabile target che è prevedibile dal modello.
(10)
dove
è la media dei veri valori. Un valore diR 2 più vicino a 1 indica un modello che spiega la maggior parte della varianza nei dati. Le equazioni 8–10 sono metriche standard di regressione48. La valutazione è stata condotta in due fasi distinte e sequenziali per valutare separatamente le prestazioni del benchmark in campione e la generalizzazione fuori campione. (1) Fase 1 (Benchmark Primario): Tutti e tredici i modelli, dopo l'addestramento e l'arresto anticipato sui dati di Cu, sono stati valutati sul set di test con il Cu offerto. (2) Fase 2 (Test di Generalizzabilità): Gli stessi esatti modelli, con i parametri congelati, sono stati impiegati per generare previsioni per l'intera serie di prezzi indipendenti di Al e Zn. Non è stato effettuato alcun ritraining o adattamento.
Riproducibilità: Impostazioni sperimentali dettagliate
I prezzi spot giornalieri (CNY/tonnellata) per Cu, Al e Zn di Grado A sono stati ottenuti dalla piattaforma pubblica SMM (https://www.smm.cn/), dal 5 gennaio 2015 al 12 settembre 2025. I dati grezzi e processati sono disponibili in un repository pubblico (DOI: 10.5281/zenodo.19976985). Il file dati contiene le colonne date, Cu, Al e Zn. Le date vengono convertite in formato datetime e ordinate in ordine crescente. I valori mancanti vengono gestiti dal riempimento in avanti seguito da un riempimento all'indietro. Le caratteristiche sono standardizzate utilizzando uno scaler z-score installato solo sul set di addestramento (media μj, deviazione standard σj, );
le stesse e vengono applicate ai set di validazione e test senza riadattamento. La variabile target (Cu, Al o Zn) viene scalata separatamente utilizzando le rispettive statistiche del set di addestramento.
Le sequenze input-output sono costruite utilizzando una finestra scorrevole con lunghezza di input L = 30 giorni di negoziazione e orizzonte di previsione h = 1 (previsione del giorno successivo). Per un indice target (0 = Cu, 1 = Al, 2 = Zn), ogni campione è definito come Xi = V[ t - L : t, : ] (forma 30 × 3) e yi = V[ t + h, k] (scalare). Non viene applicato alcun mescolare per preservare l'ordine temporale. Il dataset è suddiviso cronologicamente senza casualità: l'addestramento comprende indici 0–2080 (2.081 osservazioni, 80%), indici di validazione 2081–2340 (260 osservazioni, 10%) e indici di test 2341–2601 (261 osservazioni, 10%). I limiti corrispondenti delle date sono dal 5 gennaio 2015 al 31 luglio 2023 (addestramento), dal 1° agosto 2023 al 19 ottobre 2023 (validazione) e dal 20 ottobre 2023 al 12 settembre 2025 (test); Il file nel repository fornisce dettagli precisi.
I semi casuali sono fissati come segue: seed dell'esperimento principale = 42, e seed Python, NumPy e TensorFlow sono tutti impostati a 42. L'inizializzazione dei pesi utilizza l'uniforme di Glorot per i nuclei di ingresso, ortogonale per i nuclei ricorrenti e zeri per i bias. L'ambiente software è composto da Python 3.10.19, TensorFlow 2.20.0/Keras, NumPy 1.26.4, pandas 2.3.3, scikit-learn 1.7.2 e Matplotlib 3.10.6. Gli esperimenti sono stati eseguiti su un PC con Windows 11 e un Intel Core i7 (2,20 GHz) e 32 GB di RAM; non è stata utilizzata alcuna GPU.
L'ottimizzatore di Adam viene impiegato con learning_rate = 1×10-3, β1 = 0,9, β2 = 0,999,
e weight_decay = 0. La funzione di perdita è MSE. Uno scheduler ReduceLROnPlateau monitora la perdita di validazione con fattore 0,5, pazienza 5 e tasso minimo di apprendimento 1 × 10-5. L'arresto anticipato viene applicato con monitor = val_loss, pazienza = 10, restore_best_weights = Vero e min_delta = 0. Ogni epoca di addestramento consiste in un passaggio in avanti sul batch di addestramento, il calcolo delle perdite MSE, la retropropagazione e l'aggiornamento dei parametri Adam. Dopo ogni epoca, viene calcolata la perdita di validazione; La fermata precoce e la riduzione della velocità di apprendimento sono attivate in base a questo valore. Il modello con la perdita di validazione più bassa viene ripristinato per i test. La dimensione del lotto è 32, e i campioni vengono distribuiti in ordine cronologico senza mescolare (mescolamento = Falso).
Per i modelli ibridi CNN, viene utilizzato un layer Conv1D con 64 filtri, kernel_size = 3, stride = 1, padding = 'stesso' e attivazione rettificata di unità lineari (ReLU), seguito da MaxPooling1D(pool_size = 2) e Dropout(0.15). Nei modelli ad attenzione aumentata, la rete neurale ricorrente restituisce la sequenza nascosta completa H con forma B × T × C. Uno strato denso con un'unità produce un punteggio, e un softmax nel tempo converte questi punteggi in pesi di attenzione , con il vettore contestuale definito come c = ∑t αt ht. Segue uno strato denso con 64 unità e l'attivazione ReLU, Dropout(0.15) e lo strato denso di output. I modelli bidirezionali concatenano stati nascosti avanti e indietro (ciascuno con 64 unità), risultando in 128 dimensioni; quando si usa attenzione, return_sequences = Vero preserva la sequenza completa (B × T × 128).
La valutazione utilizza una previsione diretta a un passo avanti (non ricorsiva). Tutte le previsioni vengono trasformate inversamente alla scala di prezzo originale prima di calcolare MAE, RMSE e R2 su quella scala. Per i test di generalizzazione su Al e Zn, lo scaler di input montato su Cu viene riutilizzato senza modifiche, mentre ogni metallo bersaglio ha il proprio scalatore di bersaglio montato sui propri bersagli di addestramento. Nello studio di ablazione, tutti i parametri non architettonici (dati, divisione, scala, seed casuale, numero di epoche, dimensione del lotto, ottimizzatore, tasso di apprendimento, funzione di perdita, stop anticipato, scheduler, dropout) sono mantenuti identici lungo tutta la catena; cambia solo l'architettura. Il codice sorgente completo e le istruzioni di replica sono disponibili pubblicamente su Zenodo (10.5281/zenodo.19976985). Tutte le figure sono state generate utilizzando Matplotlib 3.10.6 con lo script fornito; gli output sono salvati come PDF, SVG e PNG ad alta risoluzione (600 dpi). Tutti i modelli furono addestrati con un massimo di 80 epoche. Una fermata anticipata con pazienza = 10 (monitorata in caso di perdita di validazione) veniva attivata per ogni modello prima di raggiungere il limite dell'epoca. Ad esempio, il modello GRU si è fermato all'epoca 37 (migliore epoca 27, migliore perdita di validazione 0,0040), mentre l'ibrido più complesso CNN–BiLSTM–Attenzione si è fermato all'epoca 23 (migliore epoca 13, migliore perdita di validazione 0,0072). L'elenco completo delle epoche ferme, delle migliori epoche e delle migliori perdite di validazione per tutti e 13 i modelli è disponibile nel repository Zenodo, garantendo completa trasparenza e riproducibilità senza sovraccaricare il testo principale con una tabella. Le discussioni teoriche (limiti di Lipschitz, complessità del campione, complessità di Rademacher, decomposizione bias–varianza, entropia dell'attenzione e informazione reciproca) sono spiegazioni concettuali per i risultati empirici e non alterano l'obiettivo di addestramento o l'implementazione del modello. Infine, per valutare la stabilità, gli esperimenti principali sono stati ripetuti con cinque semi casuali (1, 7, 21, 42 e 2024); la media e la deviazione standard dell'RMSE in queste corse sono riportate nella sezione Risultati , dove il GRU ha mantenuto una media RMSE competitiva con bassa varianza, supportando prestazioni stabili tra inizializzazioni casuali.