$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
מחקר זה לא כלל משתתפים אנושיים או בעלי חיים בעלי חוליות. כל הנתונים המשמשים הם סדרות מחירי סחורות זמינות לציבור מה-SMM, שאינן דורשות אישור אתי. לכן, לא נדרש או נדרש אישור אתי למחקר זה.
סעיף זה מפרט את עיצוב המחקר המקיף והקפדני שיושם כדי לבחון אמפירית את ההשערה המרכזית. הוא מספק הסבר מפורט של הניסוח המתמטי והפרטים האדריכליים של שלושה עשר מודלי הלמידה העמוקה שהוערכו, פרוטוקול ההכשרה המדויק, ומדדי ההערכה הפורמליים. תהליך העבודה המתודולוגי הכולל מסוכם ויזואלית באיור 1.

איור 1: סקירה סכמטית של מתודולוגיית המחקר. הדיאגרמה ממחישה את כל צינור הניסוי, כולל חלוקת נתונים, אימון מודלים בלעדית על סדרת מחירי Cu, הערכה על קבוצת בדיקות Cu, ואימות מחוץ למדגם על סדרות Al ו-Zn עצמאיות. לולאת המשוב המקווקו מצביעה על ניסויי אבלציה מובנית שנערכו לניתוח תרומתם של רכיבים אדריכליים בודדים. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.
התרשים ממחיש את כל צינור הניסוי. התהליך מתחיל בשימוש בלעדי בסדרת המחיר Cu לפיתוח דגמים. סדרה זו מחולקת כרונולוגית לקבוצות אימון (80%), אימות (10%) ובדיקות (10%). שלוש עשרה ארכיטקטורות למידה עמוקה מובחנת מאומנות ומותאמות היפרמטריות אך ורק על נתוני האימון Cu, כאשר ההפסקה המוקדמת מנוטרת באמצעות מערכת הוולידציה. המדד העיקרי הוא הערכת המודלים הללו על סט הבדיקות Cu המוחזק. חשוב לציין, כדי להעריך הכללה, אותם מודלים מאומנים מיושמים ללא שינוי כדי לחזות את סדרת המחירים העצמאית לחלוטין של אלומיניום ואבץ, ומייצגת מבחן מחמיר מחוץ למדגם. לבסוף, מתבצעים ניסויי אבלציה מובנית (לולאת המשוב המקווקוקת) כדי לפרק ולנתח את תרומת הביצועים של רכיבים אדריכליים בודדים (למשל, תשומת לב, עיבוד דו-כיווני ושכבות קונבולוציוניות).
ארכיטקטורות מודלים וניסוח מתמטי
תכננו ויישמנו ספקטרום של 13 מודלים של DL, שגדלו באופן שיטתי במורכבות הארכיטקטונית מרשתות פשוטות וחוזרות להיברידי רב-רכיבי מתוחכמים. כל המודלים חולקים את אותה מטרה מרכזית: ללמוד מיפוי
מחלון מחירים היסטורי Xt = [P t-L,P t-L+1,...,Pt-1] באורך L = 30 למחיר הבא yt = Pt.
משפחות המודל מוגדרות כך:
מודלים חזריים בסיסיים
GRU: רשת חוזרת יעילה המשתמשת בשערי עדכון (zt) ואיפוס (rt) למודולציה של זרימת המידע. המצב הנסתר ht מחושב כך:
(1)
(2)
(3)
(4)
כאשר X הוא ההפעלה הסיגמואידית,
מסמן את מכפלת האדמר, ו-xt הוא הקלט בזמן t . המצב הנסתר הסופי hL עובר דרך שכבת פלט ליניארית. משוואות 1–4 מותאמות מצ'ו ואח' 14.
LSTM: משתמש בשערי קלט (it), שכחה (ft) ושערי פלט (ot) לשמירה על מצב תא (Ct), ומספק שליטה מפורשת יותר על זיכרון לטווח ארוך.
מודלים דו-כיווניים (BiGRU ו-BiLSTM)
מודלים אלו משלבים שתי שכבות חזרתיות נפרדות שמעבדות את הרצף בכיוונים קדימה ואחורה. הייצוג הנסתר הסופי בכל שלב זמן הוא השרשור
, שבאופן תיאורטי לוכד מידע הקשרי מהעבר ומהעתיד בתוך חלון הקלט הקבוע.
מודלים מוגברים על ידי קשב (GRU–קשב ו-LSTM–קשב)
מנגנון קשב מצבור מוחל על רצף המצבים הנסתרים H = [h1,h 2,...,hL] שנוצר על ידי השכבה החוזרת הסופית. וקטור ההקשר מוגדר כסכום משוקלל:
(5)
(6)
(7)
כאן, αi מייצג את משקל תשומת הלב שניתן לצעד הזמן ההיסטורי ה-i. וקטור ההקשר c, המכיל סיכום אדפטיבי של היסטוריה רלוונטית, מוזן לשכבת החיזוי הסופית. משוואות 5–7 מותאמות מ-Bello ואחרים.
CNN–מודלים היברידיים (CNN–GRU ו-CNN–LSTM)
שכבת CNN חד-ממדית עם הפעלה של יחידה ליניארית מתיישרת (ReLU) מצורפת לשכבה החוזרת.
מודלים היברידיים מורכבים
ארכיטקטורות אלו משלבות מספר רכיבים (למשל, CNN–BiGRU–Attention, CNN–BiLSTM–Attention). הם מייצגים את מצב הטכנולוגיה מבחינת מורכבות, ומטרתם לשלב חילוץ דפוסים מקומיים (CNN), מידול דו-כיווני הקשרי ומשקל זמני אדפטיבי (תשומת לב) למסגרת אחת.
כל המודלים הוגדרו עם מידות מצב נסתר עקביות (128 יחידות לשכבות חוזרות ו-64 מסננים לשכבות CNN) והוגדרו עם שכבת פלט ליניארית אחת. עיצוב מבוקר זה מבטיח שהבדלי ביצועים נובעים מבחירות אדריכליות ולא מהבדלים בכיוונון קיבולת הדגם. מספר הפרמטרים הניתנים לאימון גדל משמעותית לאורך הספקטרום הזה.
פרוטוקול אימון, היפרפרמטרים ועיצוב מחקר אבלציה
טבלה 1 מסכמת את התצורה הניסיונית המאוחדת והקפדנית שיושמה על כל שלושה-עשר המודלים כדי להבטיח השוואה הוגנת ולהפחתת התאמת יתר. כל הדגמים אונו מאפס תוך שימוש רק בסט ההדרכה Cu המיועד. אופטימייזר אדם שימש למזעור אובדן השגיאה הממוצעת בריבוע (MSE). הטכניקה הקריטית של עצירה מוקדמת, שנודחה על מערכת האימות Cu, יושמה באופן אחיד. דבר זה הבטיח שהאימון יסתיים בנקודת ההכללה האופטימלית על נתוני Cu בלתי נראים, ובכך מנע מודלים מהתאמת יתר לרעש האימון.
| קטגוריית פרמטרים | מפרט / ערך | תיאור |
| משימה מרכזית ונתונים |
| יעד חיזוי | מחיר למחרת | תחזית סטנדרטית של צעד אחד קדימה. |
| אורך חלון קלט (L) | 60 ימי מסחר | מאזן בין הקשר היסטורי מספק לבין מורכבות המודל ויציבות האימון. |
| פיתוח מודלים |
| ערכת אימון (רק Cu) | 80% הראשונים (~2081 תצפיות) | משמש ללמידת פרמטרי מודל באמצעות הפצה אחורית. |
| סט אימות (רק Cu) | 10% הבאים (~260 תצפיות) | שימש לכיוון היפרפרמטרים ולעצירה מוקדמת; חיוני למניעת התאמת יתר. |
| סט בדיקה (רק Cu) | 10% אחרונים (~260 תצפיות) | הערכה סופית וממושכת של ביצועי המדגם (Cu). |
| ארכיטקטורת מודל |
| יחידות נסתרות של RNN | 128 | מספק יכולת ייצוג מספקת; נשמר קבוע בכל המודלים מבוססי RNN. |
| מסנני CNN | 64 | מספר מפות תכונות לשכבות CNN במודלים היברידיים. |
| הליך ההכשרה |
| אופטימייזר | אדם | מייעל קצב למידה אדפטיבי להתכנסות יציבה ויעילה. |
| קצב הלמידה ההתחלתי | 1 × 10⁻³ | שיעור התחלה סטנדרטי לאדם. |
| פונקציית אובדן | שגיאת ממוצעת בריבוע (MSE) | תקן לרגרסיה |
| גודל אצווה | 32 | אימון מיני-אצווה יעיל. |
| מקסימום אפוקס | 80 | הגבול העליון לאיטרציות אימון. |
| סבלנות עצירה מוקדמת | 10 תקופות | האימון נעצר אם אובדן האימות לא משתפר במשך 20 תקופות רצופות; משקלי הדגמים מהתקופה הטובה ביותר משוחזרים. |
| הערכה ואימות |
| מדדים ראשיים | MAE, RMSE, R² | ספק הסברים משלימים של גודל השגיאה והשונות. |
| מבחן הכללה | תחזית על סדרות מלאות של אלומיניום ואבץ (2602 תצפיות כל אחת) | הדגמים מוקפאים אחרי אימון Cu. זהו מבחן טהור וקפדני מחוץ למדגם על סחורות שונות לחלוטין. |
| עיצוב אבלציה | GRU → BiGRU → BiGRU–Attention → CNN–BiGRU–Attention | מבודדת באופן שיטתי את ההשפעה של הוספת רכיבי דו-כיווניות, תשומת לב ו-CNN. |
טבלה 1: פרמטרים ניסיוניים מרכזיים וקונפיגורציה. סיכום ההגדרה הניסויית שהוחלה על כל המודלים, כולל חלוקת נתונים, פרמטרים של ארכיטקטורת המודלים, הגדרות אימון ומדדי הערכה.
כדי לפרק את התרומה של כל רכיב אדריכלי, תוכנן מחקר אבלציה מובנה. החל מקו הבסיס הטוב ביותר (GRU), נבנתה "שרשרת מורכבות" מתקדמת. איור 2 מתאר ויזואלית את שרשרת המורכבות הזו, וממחיש את החיבור השלבי של הרכיבים. גישה מדורגת זו מאפשרת לייחס ישיר לכל שינוי בביצועי החיזוי לתוספת הדרגתית של דו-כיווניות, מנגנון הקשב, ולבסוף שכבת הרשת העצבית הקונבולוציונית. מדדי הביצועים בכל צומת בשרשרת זו מספקים ראיות אמפיריות ברורות לגבי הערך או החיסרון של כל רכיב מורכבות במשימה הספציפית של חיזוי מחירי מתכות.

איור 2: שרשרת המורכבות ששימשה במחקר האבלציה. הדיאגרמה ממחישה את התוספת הדרגתית של רכיבים אדריכליים, המתקדמים מ-GRU ל-BiGRU, BiGRU–Attention, ו-CNN–BiGRU–Attention. רצף זה מייצג את העלייה השיטתית במורכבות המודל המשמשת להערכת השפעת כל רכיב על ביצועי התחזית. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.
מדדי הערכת ביצועים
ביצועי המודל נמדדו בקפדנות באמצעות שלושה מדדי רגרסיה סטנדרטיים, שהציעו תובנות משלימות לגבי דיוק החיזוי וכוח ההסבר.
שגיאה ממוצעת מוחלטת (MAE)
מודד את גודל השגיאות הממוצע, ומספק סולם סטייה איתן וקל לפרש.
(8)
שגיאת ריבוע ממוצעת שורש (RMSE)
מדגיש שגיאות גדולות יותר בגלל פעולת הריבוע, מה שהופך אותו לרגיש יותר לחריגים ולשגיאות גדולות.
(9)
מקדם קביעה (R2)
מייצג את שיעור השונות במשתנה היעד שניתן לחיזוי מהמודל.
(10)
כאשר
הוא הממוצע של הערכים האמיתיים. ערך R2 הקרוב יותר ל-1 מציין מודל שמסביר את רוב השונות בנתונים. משוואות 8–10 הן מדדי רגרסיה סטנדרטיים48. ההערכה בוצעה בשני שלבים נפרדים, רציפים, כדי להעריך בנפרד ביצועי מדד במדגם והכללות מחוץ למדגם. (1) שלב 1 (מדד בסיסי): כל שלושה-עשר המודלים, לאחר אימון והפסקה מוקדמת בנתוני Cu, הוערכו על סט בדיקות Cu שהוחזק. (2) שלב 2 (מבחן הכללה): אותם מודלים בדיוק, עם פרמטריהם הוקפאים, הופעלו כדי לייצר תחזיות לסדרות המחירים המלאות והעצמאיות של אלומיניום ואבץ. לא בוצעה הכשרה מחדש או התאמה.
שחזוריות: הגדרות ניסיוניות מפורטות
מחירי ספוט יומיים (CNY/ton) עבור דרגת A Cu, Alum, ו-Zn התקבלו מהפלטפורמה הציבורית של SMM (https://www.smm.cn/), מ-5 בינואר 2015 ועד 12 בספטמבר 2025. הנתונים הגולמיים והמעובדים זמינים במאגר ציבורי (DOI: 10.5281/zenodo.19976985). קובץ הנתונים מכיל את העמודות תאריך, Cu, Al, ו-Zn. התאריכים מומרים לפורמט datetime וממוינים בסדר עולה. ערכים חסרים מטופלים על ידי מילוי קדמי ואחריו מילוי אחורה. התכונות מתוקננות באמצעות סולם ניקוד z המותאם רק על סט האימון (ממוצע μ j, סטיית תקן σj ,
); אותו דבר ומוחל על ערכות האימות והבדיקה ללא התאמה. משתנה היעד (Cu, Al, או Zn) מותקן בנפרד באמצעות סטטיסטיקות סט אימון משלו.
רצפי קלט–פלט נבנים באמצעות חלון הזזה עם אורך קלט L = 30 ימי מסחר ואופק תחזית h = 1 (תחזית ליום הבא). עבור אינדקס יעד (0 = Cu, 1 = Al, 2 = Zn), כל דגימה מוגדרת כ-Xi = V[ t - L : t, : ] (צורה 30 × 3) ו-yi = V[ t + h , k] (סקלר). לא נעשה ערבוב לשימור הסדר הזמני. מערך הנתונים מחולק כרונולוגית ללא אקראיות: ההכשרה כוללת אינדקסים 0–2080 (2,081 תצפיות, 80%), מדדי אימות 2081–2340 (260 תצפיות, 10%), ומדדי בדיקה 2341–2601 (261 תצפיות, 10%). גבולות התאריכים המתאימים הם 5 בינואר 2015 עד 31 ביולי 2023 (אימון), 1 באוגוסט 2023 עד 19 באוקטובר 2023 (אימות), ו-20 באוקטובר 2023 עד 12 בספטמבר 2025 (מבחן); הקובץ במאגר מספק פרטים מדויקים.
זרעים אקראיים קבועים כך: seed ניסוי ראשי = 42, ו-seeds Python, NumPy ו-TensorFlow מוגדרים כולם על 42. אתחול משקל משתמש באחידות גלורוט לגרעיני קלט, אורתוגונל לגרעינים חוזרים, ואפסים להטיות. סביבת התוכנה כוללת את Python 3.10.19, TensorFlow 2.20.0/Keras, NumPy 1.26.4, pandas 2.3.3, scikit-learn 1.7.2, ו-Matplotlib 3.10.6. הניסויים בוצעו על מחשב Windows 11 עם Intel Core i7 (2.20 GHz) ו-32 GB RAM; לא השתמשו ב-GPU.
אופטימייזר אדם משמש עם learning_rate = 1×10-3, β1 = 0.9, β2 = 0.999,
ו-weight_decay = 0. פונקציית האובדן היא MSE. מתזמן ReduceLROnPlateau עוקב אחרי אובדן אימות עם גורם 0.5, סבלנות 5, וקצב למידה מינימלי 1 ×-10-5. עצירה מוקדמת מיושמת עם מוניטור = val_loss, סבלנות = 10, restore_best_weights = נכון, ו-min_delta = 0. כל תקופת אימון כוללת מעבר קדימה על קבוצת האימון, חישוב אובדן MSE, הפצה אחורית ועדכון פרמטר אדם. לאחר כל תקופה, מחושב אובדן אימות; הפסקה מוקדמת והפחתת קצב הלמידה מופעלים בהתאם לערך זה. המודל עם אובדן האימות הנמוך ביותר משוחזר לבדיקה. גודל האצווה הוא 32, והדגימות מוזנות בסדר כרונולוגי ללא ערבוב (ערבוב = שגוי).
במודלים היברידיים של CNN, משתמשים בשכבת Conv1D אחת עם 64 מסננים, kernel_size = 3, stride = 1, ריפוד = 'אותו דבר', והפעלה ליניארית מתוקנת (ReLU), ואחריה MaxPooling1D(pool_size = 2) ו-Dropout (0.15). במודלים מוגברים על ידי תשומת לב, הרשת העצבית החוזרת מחזירה את הרצף הנסתר המלא H עם צורה B × T × C. שכבה צפופה עם יחידה אחת יוצרת ניקוד, וסופטמקס לאורך זמן ממיר את הציונים הללו למשקלי תשומת לב , כאשר וקטור ההקשר מוגדר כ-c = ∑t αt ht. לאחר מכן יש שכבה צפופה עם 64 יחידות והפעלה של ReLU, Dropout (0.15), ושכבת הפלט הצפופה. מודלים דו-כיווניים מחברים מצבים נסתרים קדימה ואחורית (כל אחד 64 יחידות), מה שמוביל ל-128 ממדים; כאשר משתמשים בתשומת לב, return_sequences = אמת שומרת על הרצף המלא (B × T × 128).
הערכה משתמשת בחיזוי צעד אחד ישיר (לא רקורסיבי). כל התחזיות עוברות הפוכה לסולם המחירים המקורי לפני חישוב MAE, RMSE ו-R2 בסולם זה. לבדיקות הכללה על אלומיניום ואבץ, סולן הקלט המותקן על Cu משמש מחדש ללא שינוי, בעוד שלכל מתכת מטרה יש סולם מטרות משלה המותקן על מטרות האימון שלה. במחקר האבלציה, כל הפרמטרים הלא-אדריכליים (נתונים, פיצול, קנה מידה, זרע אקראי, מספר תקופות, גודל אצווה, אופטימייזר, קצב למידה, פונקציית אובדן, עצירה מוקדמת, מתזמן, ניתוק) נשמרים זהים לאורך השרשרת; רק הארכיטקטורה משתנה. קוד המקור המלא והוראות השכפול זמינים לציבור ב-Zenodo (10.5281/zenodo.19976985). כל הדמויות נוצרו באמצעות Matplotlib 3.10.6 עם הסקריפט שסופק; הפלטים נשמרים כ-PDF, SVG ו-PNG ברזולוציה גבוהה (600 dpi). כל הדגמים עברו הכשרה עם מקסימום של 80 אפוקים. הפסקה מוקדמת עם סבלנות = 10 (מנוטרת עם אובדן אימות) הופעלה בכל מודל לפני שהגיעו למגבלת התקופה. לדוגמה, מודל ה-GRU נעצר בתקופה 37 (אפוק הטוב ביותר 27, אובדן אימות מיטבי 0.0040), בעוד שהמודל ההיברידי המורכב ביותר CNN–BiLSTM–Attention עצר בתקופה 23 (אפוק הטוב ביותר 13, אובדן אימות מיטבי 0.0072). הרשימה המלאה של התקופות שנעצרו, התקופות הטובות ביותר ואובדן האימות הטוב ביותר לכל 13 המודלים זמינה במאגר זנודו, ומבטיחה שקיפות מלאה ושחזוריות מבלי להעמיס טבלה על הטקסט הראשי. הדיונים התיאורטיים (גבולות ליפשיץ, מורכבות מדגם, מורכבות רדמאכר, פירוק הטיה-שונות, אנטרופיית קשב ומידע הדדי) הם הסברים מושגיים לתוצאות האמפיריות ואינם משנים את מטרת האימון או את יישום המודל. לבסוף, להערכת יציבות, הניסויים המרכזיים חזרו על עצמם עם חמישה זרעים אקראיים (1, 7, 21, 42 ו-2024); הממוצע והסטיית התקן של RMSE בין ריצות אלו מדווחות בסעיף התוצאות , שם ה-GRU שמר על RMSE ממוצע תחרותי עם שונות נמוכה, ותמך בביצועים יציבים באתחולים אקראיים.