1. 功能说明
本代码实现了基于长短期记忆网络(LSTM)的量化交易策略,通过处理时间序列金融数据预测未来价格走势。系统包含数据预处理、特征工程、模型构建、训练验证和实盘接口五个核心模块,支持多维度特征输入和自定义超参数配置。主要风险包括过拟合问题、非平稳时间序列导致的梯度消失、以及市场黑天鹅事件引发的异常波动。
2. 特殊入参设计原理
2.1 时间步长的动态调整机制
LSTM网络的时间步长(timesteps)需要根据交易品种特性进行差异化设置。对于高频交易场景,采用5-30分钟级别的K线数据时,建议设置timesteps=64-128;而日线级别数据可设置为30-60。关键实现在于:
python
def create_sequences(data, timesteps):
"""生成动态时间窗口序列"""
X, y = [], []
for i in range(len(data) - timesteps):
X.append(data[i:(i + timesteps)])
y.append(data[i + timesteps, 3]) # 收盘价作为目标变量
return np.array(X), np.array(y)
2.2 多维特征融合配置
除基础OHLC价格数据外,需集成技术指标作为补充特征。典型配置方案:
| 特征类型 | 示例指标 | 归一化方法 |
|---|---|---|
| 量价特征 | K线形态组合 | MinMaxScaler |
| 动量指标 | RSI(14), MACD(12,26,9) | StandardScaler |
| 波动率指标 | ATR(14), Bollinger Bands | RobustScaler |
| 市场情绪 | VIX指数, 资金流向指标 | Custom Scaling |
python
feature_columns = [
'open', 'high', 'low', 'close',
'rsi_14', 'macd_signal', 'atr_14',
'volume_ma5', 'bb_upper', 'bb_lower'
]
3. 模型架构定制化
3.1 门控机制优化方案
针对金融时间序列的非平稳特性,改进标准LSTM单元结构:
python
from tensorflow.keras.layers import LSTM, Dense, Dropout, Multiply
def attention_lstm(input_shape):
inputs = tf.keras.Input(shape=input_shape)
# 主LSTM层
lstm_out, state_h, state_c = LSTM(
units=64,
return_sequences=True,
return_state=True,
dropout=0.3,
recurrent_dropout=0.2
)(inputs)
# 注意力机制
attention = Dense(1, activation='tanh')(lstm_out)
attention = Flatten()(attention)
attention = Activation('softmax')(attention)
attention = RepeatVector(64)(attention)
attended = Multiply()([lstm_out, attention])
# 状态聚合
last_hidden = attended[:, -1, :]
outputs = Dense(1, activation='linear')(last_hidden)
return tf.keras.Model(inputs, outputs)
3.2 正则化策略实施
为防止过拟合,采用混合正则化方案:
- 递归层使用Zoneout(0.1-0.3)
- 全连接层应用L1/L2正则化(λ=0.001)
- 自适应学习率衰减(ReduceLROnPlateau)
- 早停法(EarlyStopping)监控验证损失
4. 数据预处理规范
4.1 非平稳序列处理流程
- 差分变换:对价格序列进行一阶差分消除趋势
- 小波去噪:使用db4小波分解去除高频噪声
- 标准化:Z-score标准化保持分布一致性
python
def preprocess_series(price_data, lookback=120):
# 计算收益率
returns = np.diff(price_data) / price_data[:-1]
# 小波去噪
coeffs = pywt.wavedec(returns, 'db4', level=3)
sigma = mad(coeffs[-1])
threshold = sigma * np.sqrt(2*np.log(len(returns)))
coeffs[1:] = (pywt.threshold(c, value=threshold, mode='soft') for c in coeffs[1:])
cleaned = pywt.waverec(coeffs, 'db4')
# 标准化
scaler = StandardScaler()
normalized = scaler.fit_transform(cleaned.reshape(-1,1)).flatten()
# 创建序列
X, y = create_sequences(normalized, lookback)
return X, y, scaler
4.2 类别不平衡解决方案
采用SMOTE-Tomek联合采样策略:
- 对少数类样本(下跌行情)进行合成过采样
- 移除多数类中的模糊边界样本
- 权重交叉熵损失函数(class_weight={0:1, 1:3})
5. 回测框架集成要点
5.1 滑点与冲击成本建模
python
def apply_transaction_cost(position, current_price, slippage=0.05):
"""
计算实际交易成本
:param position: 持仓比例(-1到1)
:param current_price: 当前价格
:param slippage: 滑点百分比
:return: 实际成交价格
"""
direction = np.sign(position) if position != 0 else 0
return current_price * (1 + direction * slippage/100)
5.2 风险管理参数绑定
将LSTM输出映射到动态仓位管理系统:
- Volatility Targeting: σ_target = 0.15 × √(252)
- Position Sizing: f(pred) = sign(pred) × min(|pred|, 0.02/σ_t)
- Stop-Loss: ATR(14) × 3倍波动率