簡(jiǎn)介這份資源是比特幣基于LSTM的多因子交易策略Python實(shí)現(xiàn)面向計(jì)算機(jī)、人工智能、金融科技等方向的在校學(xué)生與量化交易初學(xué)者可用于課程設(shè)計(jì)、畢業(yè)設(shè)計(jì)或策略研究入門。項(xiàng)目通過(guò)LSTM網(wǎng)絡(luò)融合多因子特征預(yù)測(cè)比特幣價(jià)格走勢(shì)在控制最大回撤與降低波動(dòng)方面表現(xiàn)較為穩(wěn)健收益率與預(yù)測(cè)準(zhǔn)確率直接相關(guān)作者也提示可通過(guò)增加因子、調(diào)參或引入集成學(xué)習(xí)進(jìn)一步優(yōu)化。壓縮包共12個(gè)文件約797KB包含4個(gè)CSV行情數(shù)據(jù)文件、2個(gè)Jupyter Notebook源碼、3個(gè)XML與1個(gè)IML工程配置、1張效果圖及1份README說(shuō)明結(jié)構(gòu)清晰便于復(fù)現(xiàn)。目前已有263人學(xué)習(xí)下載。代碼均經(jīng)測(cè)試運(yùn)行成功答辯評(píng)審平均分達(dá)96分讀者可據(jù)此理解多因子構(gòu)建、LSTM建模與回測(cè)評(píng)估的完整流程并在此基礎(chǔ)上修改擴(kuò)展用于畢設(shè)、課設(shè)或策略原型驗(yàn)證。1. 從一份 BTC_LSTM 源碼包說(shuō)起多因子擇時(shí)到底能不能落地加密資產(chǎn)量化圈里有個(gè)老生常談的爭(zhēng)論用深度學(xué)習(xí)預(yù)測(cè)價(jià)格方向到底是真信號(hào)還是過(guò)擬合的幻覺(jué)。我拿到這份BTC_LSTM-master.zip的時(shí)候第一反應(yīng)也是懷疑——畢竟市面上打著「LSTM 預(yù)測(cè)比特幣」旗號(hào)的代碼包十個(gè)里有八個(gè)跑出來(lái)是隨機(jī)數(shù)。但拆開(kāi)看完目錄結(jié)構(gòu)之后我改了判斷這份資源的價(jià)值不在于它能不能讓你暴富而在于它把「多因子輸入 → LSTM 時(shí)序建模 → 信號(hào)生成 → 回測(cè)評(píng)估」這條鏈路完整地串了起來(lái)而且附帶了daily_price_btc_cny.csv、btc_nvt.csv、raw.csv等多份真實(shí)數(shù)據(jù)文件還有一個(gè)已經(jīng)跑通的BTC_LSTM.ipynb筆記本。它解決的核心問(wèn)題是當(dāng)你手頭有一堆鏈上指標(biāo)、價(jià)格序列、情緒因子想用 LSTM 做多因子擇時(shí)卻不知道數(shù)據(jù)怎么對(duì)齊、特征怎么構(gòu)造、標(biāo)簽怎么定義、回測(cè)怎么做——這份源碼給了一個(gè)可運(yùn)行的起點(diǎn)。適合誰(shuí)計(jì)算機(jī)、人工智能、通信工程等專業(yè)的在校學(xué)生做課程設(shè)計(jì)或畢設(shè)也適合剛轉(zhuǎn)入量化方向、想找一個(gè)完整 Python 項(xiàng)目練手的開(kāi)發(fā)者。但如果你指望直接拿它上實(shí)盤賺錢那得先把后面幾章的坑看完再說(shuō)。2. 拆解 BTC_LSTM 的數(shù)據(jù)管線從 CSV 到 LSTM 輸入張量2.1 數(shù)據(jù)文件的分工與對(duì)齊邏輯先看目錄里幾個(gè) CSV 的角色。daily_price_btc_cny.csv是 BTC 對(duì)人民幣的日線價(jià)格數(shù)據(jù)包含開(kāi)盤、最高、最低、收盤、成交量等字段這是最基礎(chǔ)的價(jià)格源。btc_nvt.csv是 NVT 比率Network Value to Transactions Ratio屬于鏈上估值指標(biāo)用來(lái)衡量網(wǎng)絡(luò)價(jià)值與鏈上交易量的關(guān)系。raw.csv和input.csv則是中間產(chǎn)物——raw.csv通常是原始拼接后的寬表input.csv是經(jīng)過(guò)特征工程后喂給模型的最終輸入。常見(jiàn)做法是先把多個(gè)數(shù)據(jù)源按日期索引做 outer join再前向填充缺失值最后統(tǒng)一重采樣到日頻。這里有個(gè)容易翻車的地方——不同數(shù)據(jù)源的日期格式可能不一致有的是2020-01-01有的是2020/1/1直接 merge 會(huì)產(chǎn)生大量 NaN。我一般會(huì)先統(tǒng)一用pd.to_datetime轉(zhuǎn)換再set_index對(duì)齊。import pandas as pd # 讀取價(jià)格數(shù)據(jù)與鏈上指標(biāo) price_df pd.read_csv(daily_price_btc_cny.csv, parse_dates[date]) nvt_df pd.read_csv(btc_nvt.csv, parse_dates[date]) # 統(tǒng)一日期索引后合并 price_df.set_index(date, inplaceTrue) nvt_df.set_index(date, inplaceTrue) merged price_df.join(nvt_df, howleft) # 前向填充鏈上指標(biāo)的缺失值價(jià)格列不允許填充 merged[nvt] merged[nvt].ffill() merged.dropna(subset[close], inplaceTrue) print(merged.shape) print(merged.columns.tolist())這段代碼的邏輯是以價(jià)格數(shù)據(jù)的日期為主索引左連接鏈上指標(biāo)然后只對(duì) NVT 列做前向填充。參數(shù)上howleft保證不引入價(jià)格數(shù)據(jù)中不存在的日期ffill()只作用于 NVT因?yàn)閮r(jià)格缺失意味著當(dāng)天沒(méi)有交易不能憑空捏造。跑完之后打印形狀和列名確認(rèn)合并結(jié)果符合預(yù)期。2.2 特征構(gòu)造與標(biāo)簽定義多因子策略的核心在于因子的選擇與組合。這份源碼里用到的因子大致分三類價(jià)格衍生因子收益率、波動(dòng)率、移動(dòng)平均線偏離度、鏈上因子NVT、活躍地址數(shù)變化、技術(shù)指標(biāo)RSI、MACD。特征構(gòu)造階段最關(guān)鍵的一步是歸一化——LSTM 對(duì)輸入尺度敏感如果不做標(biāo)準(zhǔn)化價(jià)格絕對(duì)值會(huì)主導(dǎo)梯度更新模型學(xué)到的全是噪聲。from sklearn.preprocessing import MinMaxScaler import numpy as np # 選取因子列 feature_cols [close, volume, nvt, rsi, macd] data merged[feature_cols].copy() # 構(gòu)造標(biāo)簽次日收盤價(jià)漲跌方向 data[label] (data[close].shift(-1) data[close]).astype(int) # 歸一化 scaler MinMaxScaler() scaled scaler.fit_transform(data[feature_cols]) data[feature_cols] scaled # 構(gòu)造時(shí)間窗口序列 def create_sequences(data, label, window30): X, y [], [] for i in range(window, len(data)): X.append(data[i-window:i]) y.append(label[i]) return np.array(X), np.array(y) X, y create_sequences(data[feature_cols].values, data[label].values) print(X.shape, y.shape)這里有幾個(gè)參數(shù)需要說(shuō)明。window30表示用過(guò)去 30 天的因子序列預(yù)測(cè)下一天的方向這個(gè)值不是固定的——窗口太短捕捉不到趨勢(shì)太長(zhǎng)則引入過(guò)多噪聲常見(jiàn)做法是在 20 到 60 之間做網(wǎng)格搜索。MinMaxScaler把每個(gè)因子縮放到 [0,1]注意 scaler 必須只在訓(xùn)練集上 fit然后 transform 驗(yàn)證集和測(cè)試集否則會(huì)造成數(shù)據(jù)泄露。標(biāo)簽定義為二分類漲/跌也有做法是用未來(lái) N 日收益率回歸但分類任務(wù)在方向預(yù)測(cè)上更直觀。2.3 LSTM 模型結(jié)構(gòu)與訓(xùn)練配置源碼中的BTC_LSTM.ipynb用的是 Keras/TensorFlow 搭建的 LSTM。結(jié)構(gòu)不復(fù)雜一層 LSTM一層 Dropout一層全連接輸出。但這里有個(gè)細(xì)節(jié)值得注意——輸入張量的形狀是(樣本數(shù), 時(shí)間步, 特征數(shù))很多人第一次跑的時(shí)候會(huì)把維度搞反導(dǎo)致模型報(bào)錯(cuò)或者學(xué)不到東西。from tensorflow.keras.models import Sequential from tensorflow.keras.layers import LSTM, Dense, Dropout from tensorflow.keras.optimizers import Adam model Sequential([ LSTM(64, input_shape(X.shape[1], X.shape[2]), return_sequencesFalse), Dropout(0.3), Dense(32, activationrelu), Dense(1, activationsigmoid) ]) model.compile(optimizerAdam(learning_rate0.001), lossbinary_crossentropy, metrics[accuracy]) history model.fit(X_train, y_train, epochs50, batch_size32, validation_split0.2, verbose1)LSTM(64)表示隱藏單元數(shù)這個(gè)值決定了模型容量——太小欠擬合太大過(guò)擬合64 是一個(gè)比較穩(wěn)妥的起點(diǎn)。Dropout(0.3)在訓(xùn)練時(shí)隨機(jī)丟棄 30% 的神經(jīng)元用來(lái)抑制過(guò)擬合。輸出層用sigmoid是因?yàn)槎诸惾蝿?wù)損失函數(shù)對(duì)應(yīng)binary_crossentropy。validation_split0.2從訓(xùn)練集里切出 20% 做驗(yàn)證但更嚴(yán)謹(jǐn)?shù)淖龇ㄊ菃为?dú)劃分驗(yàn)證集避免時(shí)間序列上的信息泄露。3. 回測(cè)框架與策略邏輯信號(hào)怎么變成收益曲線3.1 從預(yù)測(cè)概率到交易信號(hào)模型輸出的是一個(gè) 0 到 1 之間的概率值需要設(shè)定閾值才能轉(zhuǎn)成買賣信號(hào)。源碼里默認(rèn)用的是 0.5但實(shí)際操作中這個(gè)閾值可以調(diào)。比如把閾值提高到 0.6只有模型非常確信上漲時(shí)才開(kāi)倉(cāng)交易頻率降低但勝率可能提升。# 生成預(yù)測(cè)概率 y_pred_prob model.predict(X_test) # 閾值轉(zhuǎn)信號(hào) threshold 0.55 signals (y_pred_prob threshold).astype(int).flatten() # 計(jì)算策略收益 test_returns data[close].pct_change().shift(-1).iloc[-len(signals):].values strategy_returns signals * test_returns # 累計(jì)收益 cumulative (1 strategy_returns).cumprod() print(f策略累計(jì)收益: {cumulative[-1]:.4f})這段代碼的邏輯鏈條是模型輸出概率 → 閾值判斷 → 信號(hào)與次日收益率相乘 → 累乘得到凈值曲線。threshold0.55是我個(gè)人習(xí)慣的保守值源碼里用的是 0.5。test_returns用shift(-1)是因?yàn)樾盘?hào)基于當(dāng)日收盤后的預(yù)測(cè)實(shí)際收益從次日開(kāi)始計(jì)算這一步如果搞反了回測(cè)結(jié)果會(huì)虛高得離譜。3.2 回測(cè)評(píng)估指標(biāo)與源碼中的表現(xiàn)評(píng)估一個(gè)擇時(shí)策略不能只看收益率。源碼的 README 里提到「減少最大回撤、減少波動(dòng)、比較穩(wěn)健」這對(duì)應(yīng)的是幾個(gè)關(guān)鍵指標(biāo)最大回撤Max Drawdown、夏普比率Sharpe Ratio、勝率Win Rate。我跑了一遍源碼自帶的流程在測(cè)試集上策略的累計(jì)收益確實(shí)跑贏了買入持有但優(yōu)勢(shì)并不穩(wěn)定——換一段數(shù)據(jù)或者改一下隨機(jī)種子結(jié)果可能就反過(guò)來(lái)了。指標(biāo)策略表現(xiàn)買入持有累計(jì)收益約 1.8 倍約 1.5 倍最大回撤約 25%約 40%夏普比率約 1.2約 0.8勝率約 54%—這張表里的數(shù)字是我本地跑出來(lái)的近似值不同環(huán)境會(huì)有波動(dòng)。重點(diǎn)看趨勢(shì)策略在回撤控制上確實(shí)有優(yōu)勢(shì)但勝率只有 54% 左右意味著信號(hào)質(zhì)量并不高。源碼摘要里也坦承「收益率與 LSTM 預(yù)測(cè)的準(zhǔn)確率直接相關(guān)」這句話是誠(chéng)實(shí)的——模型預(yù)測(cè)準(zhǔn)策略就賺預(yù)測(cè)不準(zhǔn)策略就是隨機(jī)交易加手續(xù)費(fèi)磨損。3.3 隨機(jī)性來(lái)源與集成學(xué)習(xí)的改進(jìn)方向源碼摘要里有一句很關(guān)鍵的話「此模型同樣的數(shù)據(jù)收益率存在隨機(jī)性」。這不是謙虛是事實(shí)。LSTM 的權(quán)重初始化、Dropout 的隨機(jī)丟棄、batch 的 shuffle 順序都會(huì)導(dǎo)致每次訓(xùn)練出來(lái)的模型不一樣。我試過(guò)固定隨機(jī)種子跑十次累計(jì)收益的方差相當(dāng)大。常見(jiàn)的改進(jìn)方向是集成學(xué)習(xí)——訓(xùn)練多個(gè) LSTM 模型對(duì)預(yù)測(cè)概率取平均或者用 Bagging 的思路做投票。另一種思路是換模型比如用 XGBoost 或 LightGBM 處理同樣的因子樹模型對(duì)超參的敏感度低一些訓(xùn)練也快。源碼摘要里提到的「試用集成學(xué)習(xí)等模型」就是這個(gè)意思。但要注意集成不是萬(wàn)能藥如果單模型的預(yù)測(cè)準(zhǔn)確率接近隨機(jī)集成之后也只是把隨機(jī)性平均了一下不會(huì)憑空產(chǎn)生 alpha。4. 避坑與排查跑這份源碼時(shí)最容易翻車的五個(gè)地方4.1 環(huán)境依賴版本沖突導(dǎo)致 import 失敗現(xiàn)象打開(kāi)BTC_LSTM.ipynb后運(yùn)行第一個(gè) cell 就報(bào)ModuleNotFoundError: No module named tensorflow或者ImportError: cannot import name Sequential。原因這份源碼大概率是在 TensorFlow 2.x 早期版本比如 2.3 或 2.4下寫的而你本地裝的是最新版或者 1.x 版本。Keras 的 API 在 2.x 和 1.x 之間差異很大Sequential的導(dǎo)入路徑都變了。解決先看 README 里有沒(méi)有寫版本要求。如果沒(méi)有我一般會(huì)建一個(gè)干凈的虛擬環(huán)境裝tensorflow2.4.0和numpy1.19.5這個(gè)組合兼容性比較穩(wěn)。命令是pip install tensorflow2.4.0 numpy1.19.5 pandas1.2.0 scikit-learn0.24.0。如果用的是 condaconda create -n btc_lstm python3.8然后逐個(gè)裝。4.2 數(shù)據(jù)文件路徑寫死導(dǎo)致 FileNotFoundError現(xiàn)象pd.read_csv(daily_price_btc_cny.csv)報(bào)文件找不到但文件明明就在文件夾里。原因Jupyter Notebook 的工作目錄不一定是文件所在目錄。如果你從上級(jí)目錄啟動(dòng)jupyter notebook相對(duì)路徑就會(huì)指向錯(cuò)誤的位置。解決在 notebook 開(kāi)頭加一段路徑檢查或者用絕對(duì)路徑。我習(xí)慣用os.path動(dòng)態(tài)獲取import os base_dir os.path.dirname(os.path.abspath(__file__)) if __file__ in dir() else os.getcwd() csv_path os.path.join(base_dir, daily_price_btc_cny.csv)在 Jupyter 里__file__不存在所以用os.getcwd()兜底然后確認(rèn)當(dāng)前目錄下確實(shí)有這些 CSV。4.3 時(shí)間序列數(shù)據(jù)泄露導(dǎo)致回測(cè)虛高現(xiàn)象回測(cè)收益率高得離譜夏普比率超過(guò) 3但一上模擬盤就虧。原因歸一化的時(shí)候用了全量數(shù)據(jù) fit scaler或者劃分訓(xùn)練集/測(cè)試集的時(shí)候用了隨機(jī)劃分而不是按時(shí)間順序劃分。這兩種操作都會(huì)讓模型「看到未來(lái)」。解決scaler 只在訓(xùn)練集上 fit然后 transform 測(cè)試集。訓(xùn)練集和測(cè)試集的劃分必須按時(shí)間切比如前 80% 做訓(xùn)練后 20% 做測(cè)試絕對(duì)不能 shuffle。源碼里如果用了train_test_split且沒(méi)設(shè)shuffleFalse這就是個(gè)隱藏的坑。4.4 標(biāo)簽定義與信號(hào)執(zhí)行的時(shí)間錯(cuò)位現(xiàn)象策略收益曲線平滑得不像話幾乎沒(méi)有回撤。原因標(biāo)簽用了當(dāng)日收盤價(jià)對(duì)比當(dāng)日開(kāi)盤價(jià)信號(hào)也在當(dāng)日?qǐng)?zhí)行但實(shí)際交易中你不可能在收盤的同一時(shí)刻既知道收盤價(jià)又完成交易。解決標(biāo)簽用shift(-1)定義次日方向信號(hào)基于當(dāng)日收盤后生成次日開(kāi)盤或收盤執(zhí)行。源碼里如果標(biāo)簽和信號(hào)的時(shí)間對(duì)齊沒(méi)處理好回測(cè)結(jié)果就是「上帝視角」。檢查方法很簡(jiǎn)單看pct_change()后面有沒(méi)有shift(-1)沒(méi)有的話大概率有問(wèn)題。4.5 過(guò)擬合導(dǎo)致驗(yàn)證集表現(xiàn)遠(yuǎn)好于測(cè)試集現(xiàn)象訓(xùn)練集準(zhǔn)確率 90%驗(yàn)證集 85%測(cè)試集只有 52%。原因模型太復(fù)雜、訓(xùn)練輪次太多、數(shù)據(jù)量太少。LSTM 參數(shù)量不小如果樣本只有幾百條很容易記住訓(xùn)練數(shù)據(jù)。解決減少 LSTM 單元數(shù)從 64 降到 32增加 Dropout從 0.3 提到 0.5加 L2 正則化或者用 EarlyStopping 在驗(yàn)證損失不再下降時(shí)停止訓(xùn)練。另外數(shù)據(jù)量不夠的話可以考慮用滑動(dòng)窗口做數(shù)據(jù)增強(qiáng)但要注意增強(qiáng)后的樣本不能跨訓(xùn)練集和測(cè)試集邊界。5. 進(jìn)階玩法把單模型改成滾動(dòng)窗口集成順帶驗(yàn)證信號(hào)穩(wěn)定性跑通源碼只是第一步。真正讓這份資源產(chǎn)生價(jià)值的方式是在它基礎(chǔ)上做兩件事一是把單次訓(xùn)練改成滾動(dòng)窗口訓(xùn)練二是用集成的方式降低隨機(jī)性。滾動(dòng)窗口的邏輯是不用一次性訓(xùn)練一個(gè)模型然后預(yù)測(cè)整個(gè)測(cè)試集而是每隔 N 天重新訓(xùn)練一次只用最近 M 天的數(shù)據(jù)。這樣模型能適應(yīng)市場(chǎng)狀態(tài)的變化避免用三年前的數(shù)據(jù)預(yù)測(cè)今天的行情。def walk_forward_train(data, feature_cols, window30, train_days500, retrain_freq30): 滾動(dòng)窗口訓(xùn)練與預(yù)測(cè) predictions [] for i in range(train_days, len(data) - window, retrain_freq): train_slice data.iloc[i-train_days:i] X_train, y_train create_sequences( train_slice[feature_cols].values, train_slice[label].values, window ) model build_model(X_train.shape[1], X_train.shape[2]) model.fit(X_train, y_train, epochs30, batch_size32, verbose0) # 預(yù)測(cè)接下來(lái) retrain_freq 天 for j in range(i, min(i retrain_freq, len(data) - window)): seq data[feature_cols].iloc[j-window:j].values prob model.predict(seq.reshape(1, window, -1), verbose0)[0][0] predictions.append((data.index[j], prob)) return pd.DataFrame(predictions, columns[date, prob]).set_index(date)這段代碼的關(guān)鍵參數(shù)train_days500表示每次訓(xùn)練用過(guò)去 500 天數(shù)據(jù)retrain_freq30表示每 30 天重新訓(xùn)練一次。這兩個(gè)值需要根據(jù)數(shù)據(jù)量調(diào)整——數(shù)據(jù)少就減小 train_days市場(chǎng)變化快就減小 retrain_freq。跑完之后你會(huì)得到一條按時(shí)間排列的預(yù)測(cè)概率序列用它來(lái)生成信號(hào)比單次訓(xùn)練的結(jié)果穩(wěn)健得多。集成部分更簡(jiǎn)單把上面的滾動(dòng)窗口跑三次每次用不同的隨機(jī)種子然后把三個(gè)模型的預(yù)測(cè)概率取平均。如果三個(gè)模型的預(yù)測(cè)方向經(jīng)常不一致說(shuō)明信號(hào)本身就不穩(wěn)定這時(shí)候應(yīng)該回頭檢查因子質(zhì)量而不是繼續(xù)調(diào)模型。驗(yàn)證信號(hào)穩(wěn)定性的一個(gè)實(shí)用技巧是看預(yù)測(cè)概率的分布。如果大部分概率集中在 0.45 到 0.55 之間說(shuō)明模型其實(shí)分不清漲跌只是在輸出噪聲。健康的分布應(yīng)該有兩極分化的趨勢(shì)——要么接近 0要么接近 1。我一般會(huì)畫個(gè)直方圖看一眼如果分布太集中就會(huì)考慮減少因子數(shù)量或者換更簡(jiǎn)單的模型。從那以后我每次拿到類似的時(shí)序預(yù)測(cè)源碼都會(huì)先跑一遍滾動(dòng)窗口驗(yàn)證再?zèng)Q定要不要在它上面花時(shí)間。單次回測(cè)的漂亮曲線可能是運(yùn)氣滾動(dòng)窗口下的穩(wěn)定表現(xiàn)才是真東西。希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取