$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Este estudo não envolveu participantes humanos nem vertebrados animais. Todos os dados utilizados são séries de preços de commodities publicamente disponíveis da SMM, que não requerem aprovação ética. Portanto, nenhuma aprovação ética foi solicitada ou exigida para essa pesquisa.
Esta seção delineia o desenho abrangente e rigoroso de pesquisa implementado para testar empiricamente a hipótese central. Ele fornece uma exposição detalhada da formulação matemática e das especificidades arquitetônicas dos treze modelos de aprendizado profundo avaliados, do protocolo preciso de treinamento e das métricas formais de avaliação. O fluxo de trabalho metodológico geral é visualmente resumido na Figura 1.

Figura 1: Visão geral esquemática da metodologia de pesquisa. O diagrama ilustra todo o pipeline experimental, incluindo particionamento de dados, treinamento de modelos exclusivamente na série de preços de, avaliação no conjunto de teste e validação fora da amostra em séries independentes Al e Zn. O ciclo de retroalimentação tracejado indica os experimentos estruturados de ablação realizados para analisar a contribuição de componentes arquitetônicos individuais. Por favor, clique aqui para ver uma versão ampliada desta figura.
O diagrama ilustra todo o pipeline experimental. O processo começa com o uso exclusivo da série de preços para o desenvolvimento de modelos. Essa série é dividida cronologicamente em conjuntos de treinamento (80%), validação (10%) e testes (10%). Treze arquiteturas distintas de deep learning são então treinadas e otimizadas por hiperparâmetros exclusivamente com os dados de treinamento de, com a parada precoce monitorada pelo conjunto de validação. O principal parâmetro é a avaliação desses modelos no conjunto de teste que está estendido. Crucialmente, para avaliar a generalizabilidade, os mesmos modelos treinados são aplicados sem modificações para prever a série de preços totalmente independente de Al e Zn, representando um teste rigoroso fora da amostra. Por fim, experimentos de ablação estruturada (o ciclo de retroalimentação tracejada) são realizados para decompor e analisar a contribuição de desempenho de componentes arquitetônicos individuais (por exemplo, atenção, processamento bidirecional e camadas convolucionais).
Arquiteturas de Modelos e Formulação Matemática
Projetamos e implementamos um espectro de 13 modelos DL, aumentando sistematicamente a complexidade arquitetônica, desde redes recorrentes simples até híbridos multicomponentes sofisticados. Todos os modelos compartilham o mesmo objetivo central: aprender um mapeamento
a partir de uma janela histórica de preços Xt = [Pt-L, Pt-L+1,...,Pt-1] de comprimento L = 30 para o preço subsequente yt = Pt.
As famílias de modelos são definidas da seguinte forma:
Modelos Recorrentes de Referência
GRU: Uma rede recorrente simplificada que utiliza portas de atualização (zt) e reset (rt) para modular o fluxo de informação. O estado oculto ht é calculado como:
(1)
(2)
(3)
(4)
onde X é a ativação sigmoide,
denota o produto de Hadamard, e xt é a entrada no tempo t . O estado oculto final hL é passado por uma camada de saída linear. As equações 1–4 são adaptadas de Cho et al.14.
LSTM: Utiliza portas de entrada (it), forget (ft) e saída (ot) para manter um estado de célula (Ct), proporcionando controle mais explícito sobre a memória de longo prazo.
Modelos Bidirecionais (BiGRU e BiLSTM)
Esses modelos incorporam duas camadas recorrentes separadas que processam a sequência em direções para frente e para trás. A representação oculta final em cada passo de tempo é a concatenação
, teoricamente capturando informações contextuais tanto do passado quanto do futuro dentro da janela de entrada fixa.
Modelos Aumentados pela Atenção (GRU–Atenção e LSTM–Atenção)
Um mecanismo aditivo de atenção é aplicado à sequência de estados ocultos H = [h1,h 2,...,h L] produzida pela camada final recorrente. O vetor de contexto é definido como uma soma ponderada:
(5)
(6)
(7)
Aqui, α i representa o peso de atenção atribuído ao i-ésimo passo histórico de tempo. O vetor de contexto c, que encapsula um resumo adaptativo da história relevante, é alimentado para a camada final de previsão. As equações 5–7 são adaptadas de Bello et al.47.
CNN–Modelos Híbridos (CNN–GRU e CNN–LSTM)
Uma camada CNN unidimensional com ativação linear retificada (ReLU) é precedida à camada recorrente.
Modelos Híbridos Complexos
Essas arquiteturas combinam múltiplos componentes (por exemplo, CNN–BiGRU–Atenção, CNN–BiLSTM–Atenção). Eles representam o estado da arte em termos de complexidade, visando integrar extração local de padrões (CNN), modelagem contextual bidirecional e ponderação temporal adaptativa (atenção) em um único framework.
Todos os modelos foram configurados com dimensões de estado oculto consistentes (128 unidades para camadas recorrentes e 64 filtros para camadas CNN) e finalizados com uma única camada linear de saída. Esse projeto controlado garante que as diferenças de desempenho sejam atribuídas a escolhas arquitetônicas, e não a discrepâncias na afinação da capacidade do modelo. O número de parâmetros treináveis, consequentemente, aumentou substancialmente ao longo desse espectro.
Protocolo de Treinamento, Hiperparâmetros e Projeto do Estudo de Ablação
A Tabela 1 resume a configuração experimental unificada e rigorosa aplicada a todos os treze modelos para garantir uma comparação justa e mitigar o sobreajuste. Todos os modelos foram treinados do zero usando apenas o conjunto de treinamento designado. O otimizador Adam foi usado para minimizar a perda por erro quadrático médio (MSE). A técnica crítica de parada precoce, monitorada no conjunto de validação de, foi aplicada de forma uniforme. Isso garantia que o treinamento fosse concluído no ponto de generalização ótima sobre dados de não vistos, prevenindo assim que os modelos superajustassem ao ruído de treinamento.
| Categoria de parâmetros | Especificação / Valor | Descrição |
| Tarefa Principal & Dados |
| Meta de Previsão | Preço do dia seguinte | Previsão padrão de um passo à frente. |
| Comprimento da janela de entrada (L) | 60 dias de negociação | Equilibra contexto histórico suficiente com complexidade do modelo e estabilidade no treinamento. |
| Desenvolvimento de Modelos |
| Conjunto de Treinamento (Somente) | Primeiros 80% (~2081 observações) | Usado para aprender parâmetros do modelo por meio de retropropagação. |
| Conjunto de Validação (Somente) | Próximos 10% (~260 observações) | Usado para ajuste de hiperparâmetros e parada precoce; Crucial para evitar o excesso de ajuste. |
| Conjunto de Teste (Somente) | 10% final (~260 observações) | Avaliação final, mantida, do desempenho em amostra (). |
| Arquitetura de Modelos |
| Unidades Ocultas da RNN | 128 | Oferece capacidade representativa adequada; mantido constante em todos os modelos baseados em RNN. |
| Filtros CNN | 64 | Número de mapas de características para camadas CNN em modelos híbridos. |
| Procedimento de Treinamento |
| Otimizador | Adam | Otimizador de taxa de aprendizado adaptativo para convergência estável e eficiente. |
| Taxa de Aprendizagem Inicial | 1 × 10⁻³ | Taxa inicial padrão para Adam. |
| Função de Perda | Erro quadrático médio (MSE) | Padrão para regressão |
| Tamanho do lote | 32 | Treinamento eficiente em mini-lote. |
| Épocas Máximas | 80 | Limite superior para iterações de treinamento. |
| Parando Antecemente a Paciência | 10 épocas | O treinamento é interrompido se a perda de validação não melhorar por 20 épocas consecutivas; Os pesos dos modelos da melhor época são restaurados. |
| Avaliação e Validação |
| Métricas Primárias | MAE, RMSE, R² | Forneça visões complementares da magnitude e variância do erro explicadas. |
| Teste de Generalização | Previsão para séries completas Al e Zn (2602 observações cada) | Os modelos são congelados após o treinamento de. Este é um teste puro e rigoroso fora da amostra em commodities totalmente diferentes. |
| Design de Ablação | GRU → BiGRU → BiGRU–Atenção → CNN–BiGRU–Atenção | Isola sistematicamente o impacto da adição de componentes bidirecionais, atenção e CNN. |
Tabela 1: Parâmetros experimentais chave e configuração. Resumo da configuração experimental aplicada em todos os modelos, incluindo particionamento de dados, parâmetros da arquitetura do modelo, configurações de treinamento e métricas de avaliação.
Para desconstruir a contribuição de cada componente arquitetônico, foi elaborado um estudo estruturado de ablação. Partindo da linha de base de melhor desempenho (GRU), foi construída uma "cadeia de complexidade" progressiva. A Figura 2 descreve visualmente essa cadeia de complexidade, ilustrando a adição de componentes passo a passo. Essa abordagem passo a passo permite a atribuição direta de qualquer alteração no desempenho da previsão à adição incremental da bidirecionalidade, ao mecanismo de atenção e, finalmente, à camada convolucional de rede neural. As métricas de desempenho em cada nó dessa cadeia fornecem evidências empíricas claras sobre o valor ou desvantagem de cada componente de complexidade para a tarefa específica da previsão de preços dos metais.

Figura 2: Cadeia de complexidade usada no estudo de ablação. O diagrama ilustra a adição gradual de componentes arquitetônicos, progredindo de GRU para BiGRU, BiGRU–Atenção e CNN–BiGRU–Atenção. Essa sequência representa o aumento sistemático da complexidade do modelo usado para avaliar o impacto de cada componente no desempenho da previsão. Por favor, clique aqui para ver uma versão ampliada desta figura.
Métricas de Avaliação de Desempenho
O desempenho do modelo foi rigorosamente quantificado usando três métricas padrão de regressão, oferecendo insights complementares sobre a precisão da previsão e o poder explicativo.
Erro Absoluto Médio (EMA)
Mede a magnitude média dos erros, fornecendo uma escala robusta e facilmente interpretável de desvio.
(8)
Erro Quadrático Médio Raiz (RMSE)
Enfatiza erros maiores devido à operação de quadrado, tornando-a mais sensível a outliers e erros grandes.
(9)
Coeficiente de Determinação (R 2)
Representa a proporção de variância na variável-alvo que é previsível a partir do modelo.
(10)
onde
é a média dos valores verdadeiros. Um valor de R2 mais próximo de 1 indica um modelo que explica a maior parte da variância nos dados. As equações 8–10 são métricas padrão de regressão48. A avaliação foi realizada em duas fases distintas e sequenciais para avaliar separadamente o desempenho em benchmarks dentro da amostra e a generalização fora da amostra. (1) Fase 1 (Benchmark Primário): Todos os treze modelos, após treinamento e paragem precoce dos dados de, foram avaliados no conjunto de teste de que foi mantido em aberto. (2) Fase 2 (Teste de Generalizabilidade): Os mesmos modelos, com seus parâmetros congelados, foram implantados para gerar previsões para a série completa e independente de preços Al e Zn. Nenhuma reeducação ou adaptação foi realizada.
Reprodutibilidade: Configurações Experimentais Detalhadas
Os preços à vista diários (CNY/tonelada) para, Al e Zn de Grau A foram obtidos da plataforma pública da SMM (https://www.smm.cn/), abrangendo de 5 de janeiro de 2015 a 12 de setembro de 2025. Os dados brutos e processados estão disponíveis em um repositório público (DOI: 10.5281/zenodo.19976985). O arquivo de dados contém as colunas date,, Al e Zn. As datas são convertidas para o formato de data-hora e ordenadas em ordem crescente. Valores ausentes são tratados por preenchimento direto seguido por preenchimento para trás. As características são padronizadas usando um escalador de pontuação z ajustado apenas no conjunto de treinamento (média μj, desvio padrão σ j, );
as mesmas e são aplicadas aos conjuntos de validação e teste sem necessidade de reajuste. A variável-alvo (, Al ou Zn) é escalada separadamente usando suas próprias estatísticas do conjunto de treinamento.
As sequências de entrada–saída são construídas usando uma janela deslizante com comprimento de entrada L = 30 dias de negociação e horizonte de previsão h = 1 (previsão do dia seguinte). Para um índice-alvo (0 =, 1 = Al, 2 = Zn), cada amostra é definida como Xi = V[ t - L : t, : ] (forma 30 × 3) e yi = V[ t + h , k] (escalar). Não é aplicado embaralhamento para preservar a ordem temporal. O conjunto de dados é dividido cronologicamente sem aleatoriedade: o treinamento compreende índices 0–2080 (2.081 observações, 80%), índices de validação 2081–2340 (260 observações, 10%) e índices de teste 2341–2601 (261 observações, 10%). Os limites correspondentes para datas são de 5 de janeiro de 2015 a 31 de julho de 2023 (treinamento), 1º de agosto de 2023 a 19 de outubro de 2023 (validação) e 20 de outubro de 2023 a 12 de setembro de 2025 (teste); O arquivo no repositório fornece detalhes exatos.
Sementes aleatórias são fixadas da seguinte forma: semente do experimento principal = 42, e sementes Python, NumPy e TensorFlow estão todas definidas para 42. A inicialização de pesos usa uniforme de Glorot para núcleos de entrada, ortogonal para núcleos recorrentes e zeros para vieses. O ambiente de software consiste em 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. Os experimentos foram realizados em um PC com Windows 11 e um Intel Core i7 (2,20 GHz) e 32 GB de RAM; nenhuma GPU foi usada.
O otimizador de Adam é empregado com learning_rate = 1×10-3, β1 = 0,9, β2 = 0,999,
, e weight_decay = 0. A função de perda é MSE. Um escalonador ReduceLROnPlateau monitora a perda de validação com fator 0,5, paciência 5 e taxa mínima de aprendizado 1 × 10-5. Parada antecipada é aplicada com monitor = val_loss, paciência = 10, restore_best_weights = Verdadeiro e min_delta = 0. Cada época de treinamento consiste em uma passagem direta no lote de treinamento, cálculo de perda MSE, retropropagação e atualização do parâmetro Adam. Após cada época, a perda de validação é calculada; Parar precocemente e reduzir a taxa de aprendizado são acionados com base nesse valor. O modelo com menor perda de validação é restaurado para testes. O tamanho do lote é 32, e as amostras são alimentadas em ordem cronológica sem embaralhar (embaralhamento = Falso).
Para modelos híbridos CNN, usa-se uma camada Conv1D com 64 filtros, kernel_size = 3, stride = 1, padding = 'igual' e ativação de unidade linear retificada (ReLU), seguida por MaxPooling1D(pool_size = 2) e Dropout(0,15). Em modelos aumentados pela atenção, a rede neural recorrente retorna a sequência oculta completa H com a forma B × T × C. Uma camada densa com uma unidade produz uma pontuação, e um softmax ao longo do tempo converte essas pontuações em pesos de atenção , com o vetor de contexto definido como c = ∑t αth t. Em seguida, vem uma camada densa com 64 unidades e ativação ReLU, Dropout(0,15) e a camada densa de saída. Modelos bidirecionais concatenam estados ocultos para frente e para trás (cada um com 64 unidades), resultando em 128 dimensões; quando a atenção é usada, return_sequences = Verdadeiro preserva a sequência completa (B × T × 128).
A avaliação utiliza previsão direta de um passo à frente (não recursiva). Todas as previsões são transformadas inversamente para a escala de preços original antes de calcular MAE, RMSE e R 2 nessa escala. Para testes de generalização em Al e Zn, o escalador de entrada instalado em é reutilizado sem modificações, enquanto cada metal alvo possui seu próprio escalador de alvo instalado em seus próprios alvos de treinamento. No estudo de ablação, todos os parâmetros não arquitetônicos (dados, divisão, escalonamento, semente aleatória, número de épocas, tamanho do lote, otimizador, taxa de aprendizado, função de perda, parada antecipada, agendador, desistência) são mantidos idênticos ao longo da cadeia; Apenas a arquitetura muda. O código-fonte completo e as instruções de replicação estão disponíveis publicamente em Zenodo (10.5281/zenodo.19976985). Todas as figuras foram geradas usando Matplotlib 3.10.6 com o script fornecido; as saídas são salvas como PDF, SVG e PNG de alta resolução (600 dpi). Todos os modelos foram treinados com um máximo de 80 épocas. Paradas precoces com paciência = 10 (monitoradas na perda de validação) foram acionadas para cada modelo antes de atingir o limite da época. Por exemplo, o modelo GRU parou na época 37 (melhor época 27, melhor perda de validação 0,0040), enquanto o híbrido mais complexo CNN–BiLSTM–Atenção parou na época 23 (melhor época 13, melhor perda de validação 0,0072). A lista completa de épocas paradas, melhores épocas e melhores perdas de validação para todos os 13 modelos está disponível no repositório Zenodo, garantindo total transparência e reprodutibilidade sem sobrecarregar o texto principal com uma tabela. As discussões teóricas (limites de Lipschitz, complexidade da amostra, complexidade de Rademacher, decomposição por viés-variância, entropia da atenção e informação mútua) são explicações conceituais para os resultados empíricos e não alteram o objetivo de treinamento ou a implementação do modelo. Por fim, para avaliar a estabilidade, os experimentos principais foram repetidos com cinco sementes aleatórias (1, 7, 21, 42 e 2024); a média e o desvio padrão do RMSE nessas execuções são relatados na seção de Resultados, onde o GRU manteve RMSE média competitiva com baixa variância, apoiando desempenho estável em inicializações aleatórias.