声明:本文部分链接为联盟推广链接,不影响价格。
前言
在 A 股市场 5000+ 只股票中,如何快速筛选出具有上涨潜力的标的?传统的基本面分析需要阅读大量财报,技术分析又面临参数过拟合问题。今天,我们来聊一种用深度学习做选股的新思路——基于 LSTM 的智能选股系统。
本文核心价值:
- 完整的 LSTM 选股模型代码(Python + Keras + akshare)
- 从数据获取到模型训练的完整流程
- 回测验证策略有效性
- 适合有一定 Python 基础的量化爱好者
⚠️ 风险提示:本文代码仅供学习参考,不构成任何投资建议。量化策略存在模型风险、市场风险,实盘前请充分回测。
1. 为什么选择 LSTM 做选股?
1.1 传统选股的痛点
| 方法 | 优点 | 缺点 |
|---|---|---|
| 基本面分析 | 逻辑清晰 | 数据滞后,需要专业财务知识 |
| 技术指标 | 直观快速 | 参数敏感,容易过拟合 |
| 量化多因子 | 系统化 | 因子挖掘困难,线性模型表达能力有限 |
1.2 LSTM 的优势
LSTM(长短期记忆网络) 是一种特殊的循环神经网络(RNN),擅长处理时序数据。股价走势本身就是典型的时间序列,LSTM 可以:
- 自动提取特征:无需手动设计因子,模型自动学习有效特征
- 捕捉长短期依赖:记住"长期趋势"的同时关注"短期波动"
- 非线性建模:拟合股价与基本面、技术面之间的复杂关系
2. 环境准备
# 安装必要的库
# pip install akshare tensorflow scikit-learn pandas numpy matplotlib
import akshare as ak
import pandas as pd
import numpy as np
from sklearn.preprocessing import MinMaxScaler
from sklearn.model_selection import train_test_split
import tensorflow as tf
from tensorflow.keras.models import Sequential
from tensorflow.keras.layers import LSTM, Dense, Dropout
import matplotlib.pyplot as plt
import warnings
warnings.filterwarnings('ignore')
print(f"TensorFlow 版本: {tf.__version__}")
print(f"akshare 版本: {ak.__version__}")
3. 数据获取与预处理
3.1 获取股票数据
我们使用 akshare 获取 A 股日线数据:
def get_stock_data(stock_code, start_date='20200101', end_date='20251231'):
"""
获取股票日线数据
参数:
stock_code: 股票代码,如 '000001'(平安银行)
start_date: 开始日期,格式 'YYYYMMDD'
end_date: 结束日期,格式 'YYYYMMDD'
返回:
DataFrame: 包含日期、开盘、收盘、最高、最低、成交量等
"""
try:
# 获取平安银行日线数据
df = ak.stock_zh_a_hist(symbol=stock_code, start_date=start_date,
end_date=end_date, adjust="qfq")
# 数据清洗
df = df[['日期', '开盘', '收盘', '最高', '最低', '成交量', '成交额']]
df.columns = ['date', 'open', 'close', 'high', 'low', 'volume', 'amount']
# 转换日期格式
df['date'] = pd.to_datetime(df['date'])
df = df.sort_values('date').reset_index(drop=True)
return df
except Exception as e:
print(f"获取数据失败: {e}")
return None
# 测试获取数据
df = get_stock_code('000001')
print(f"数据形状: {df.shape}")
print(df.tail())
3.2 特征工程
构建模型输入特征:我们使用技术指标作为特征:
def create_features(df):
"""
创建技术特征
包含:
- 价格变化率
- 成交量变化率
- 移动平均线 (MA5, MA10, MA20)
- RSI 相对强弱指标
- MACD
"""
data = df.copy()
# 价格变化率
data['price_change'] = data['close'].pct_change()
data['volume_change'] = data['volume'].pct_change()
# 移动平均线
for window in [5, 10, 20]:
data[f'ma_{window}'] = data['close'].rolling(window=window).mean()
data[f'volume_ma_{window}'] = data['volume'].rolling(window=window).mean()
# RSI
def calculate_rsi(prices, period=14):
delta = prices.diff()
gain = (delta.where(delta > 0, 0)).rolling(window=period).mean()
loss = (-delta.where(delta < 0, 0)).rolling(window=period).mean()
rs = gain / loss
rsi = 100 - (100 / (1 + rs))
return rsi
data['rsi'] = calculate_rsi(data['close'])
# MACD
exp12 = data['close'].ewm(span=12, adjust=False).mean()
exp26 = data['close'].ewm(span=26, adjust=False).mean()
data['macd'] = exp12 - exp26
data['signal'] = data['macd'].ewm(span=9, adjust=False).mean()
# 删除 NaN
data = data.dropna()
return data
# 创建特征
df_features = create_features(df)
print(f"特征数量: {df_features.shape[1] - 2}") # 减去 date 和 close
print(df_features.columns.tolist())
3.3 构建训练数据
def create_sequences(data, seq_length=60, predict_days=5):
"""
创建 LSTM 输入序列
参数:
data: 特征数据
seq_length: 输入序列长度(天数)
predict_days: 预测未来第 N 天的涨跌
返回:
X: 输入序列 (samples, seq_length, features)
y: 目标标签 (1=上涨, 0=下跌)
"""
# 选择特征列(排除 date 和 close)
feature_cols = [col for col in data.columns if col not in ['date', 'close']]
features = data[feature_cols].values
# 标准化
scaler = MinMaxScaler()
features_scaled = scaler.fit_transform(features)
X, y = [], []
for i in range(len(features_scaled) - seq_length - predict_days + 1):
# 输入序列
seq = features_scaled[i:i+seq_length]
X.append(seq)
# 目标:未来 predict_days 天是否上涨
future_price = data['close'].iloc[i+seq_length+predict_days-1]
current_price = data['close'].iloc[i+seq_length-1]
label = 1 if future_price > current_price else 0
y.append(label)
return np.array(X), np.array(y), scaler
# 创建序列
seq_length = 60 # 使用过去 60 天数据
X, y, scaler = create_sequences(df_features, seq_length=seq_length)
print(f"输入数据形状: {X.shape}")
print(f"标签数据形状: {y.shape}")
print(f"上涨样本比例: {y.mean():.2%}")
4. 构建 LSTM 模型
def build_lstm_model(input_shape):
"""
构建 LSTM 选股模型
结构:
- LSTM 层 x2(提取时序特征)
- Dropout(防止过拟合)
- Dense 输出层(分类)
"""
model = Sequential([
# 第一个 LSTM 层
LSTM(64, return_sequences=True, input_shape=input_shape),
Dropout(0.2),
# 第二个 LSTM 层
LSTM(32, return_sequences=False),
Dropout(0.2),
# 全连接层
Dense(16, activation='relu'),
Dropout(0.1),
# 输出层(二分类:上涨/下跌)
Dense(1, activation='sigmoid')
])
# 编译模型
model.compile(
optimizer='adam',
loss='binary_crossentropy',
metrics=['accuracy']
)
return model
# 构建模型
input_shape = (X.shape[1], X.shape[2]) # (时间步, 特征数)
model = build_lstm_model(input_shape)
model.summary()
5. 训练与回测
5.1 划分训练集和测试集
# 划分数据集(80% 训练,20% 测试)
X_train, X_test, y_train, y_test = train_test_split(
X, y, test_size=0.2, shuffle=False # 时序数据不能随机打乱
)
print(f"训练集: {X_train.shape[0]} 样本")
print(f"测试集: {X_test.shape[0]} 样本")
5.2 训练模型
# 训练模型
history = model.fit(
X_train, y_train,
epochs=50,
batch_size=32,
validation_split=0.1,
verbose=1
)
# 保存模型
model.save('/home/node/.openclaw/workspace-moxing/models/lstm_stock_selector.h5')
print("模型已保存")
5.3 评估模型
# 评估模型
train_loss, train_acc = model.evaluate(X_train, y_train, verbose=0)
test_loss, test_acc = model.evaluate(X_test, y_test, verbose=0)
print(f"\n训练集准确率: {train_acc:.2%}")
print(f"测试集准确率: {test_acc:.2%}")
# 预测
y_pred_proba = model.predict(X_test)
y_pred = (y_pred_proba > 0.5).astype(int).flatten()
# 混淆矩阵
from sklearn.metrics import confusion_matrix, classification_report
cm = confusion_matrix(y_test, y_pred)
print("\n混淆矩阵:")
print(cm)
print("\n分类报告:")
print(classification_report(y_test, y_pred, target_names=['下跌', '上涨']))
6. 选股回测策略
基于训练好的 LSTM 模型,构建一个简单的选股回测策略:
def backtest_strategy(model, df, scaler, seq_length=60, top_n=5):
"""
选股回测策略
策略逻辑:
1. 每日收盘后,用过去 60 天数据预测次日涨跌概率
2. 选出预测上涨概率最高的 top_n 只股票
3. 等权持有 5 天,然后重新选股
"""
# 特征列
feature_cols = [col for col in df.columns if col not in ['date', 'close']]
results = []
position = 0
capital = 100000 # 初始资金 10 万
shares = {} # 持仓
# 滑动窗口回测
for i in range(seq_length, len(df) - 60, 60): # 每 60 天调仓一次
# 获取窗口数据
window_data = df.iloc[i-seq_length:i]
# 预测
features = window_data[feature_cols].values
features_scaled = scaler.transform(features)
X_input = features_scaled.reshape(1, seq_length, -1)
prob = model.predict(X_input, verbose=0)[0][0]
# 选股决策
if prob > 0.6: # 上涨概率 > 60% 则建仓
position = 1
shares_per_stock = capital / top_n
# 模拟买入
buy_price = df['close'].iloc[i]
shares['buy_price'] = buy_price
shares['buy_date'] = df['date'].iloc[i]
print(f"{df['date'].iloc[i].strftime('%Y-%m-%d']} 买入, 价格: {buy_price:.2f}, 概率: {prob:.2%}")
elif prob < 0.4 and position == 1: # 上涨概率 < 40% 则平仓
sell_price = df['close'].iloc[i]
profit = (sell_price - shares['buy_price']) / shares['buy_price']
capital *= (1 + profit)
print(f"{df['date'].iloc[i].strftime('%Y-%m-%d')} 卖出, 价格: {sell_price:.2f}, 收益率: {profit:.2%}")
position = 0
return capital
# 运行回测
final_capital = backtest_strategy(model, df_features, scaler)
total_return = (final_capital - 100000) / 100000
print(f"\n=== 回测结果 ===")
print(f"初始资金: 100,000 元")
print(f"最终资金: {final_capital:.2f} 元")
print(f"总收益率: {total_return:.2%}")
7. 模型优化方向
当前模型只是一个基础版本,实际应用中可以从以下方向优化:
7.1 数据层面
- 引入更多数据源:财报数据、舆情数据、宏观数据
- 增加股票池:覆盖全 A 股 5000+ 股票
- 更长回测周期:至少 5 年数据
7.2 模型层面
- 多模型融合:LSTM + Transformer + 传统机器学习
- AutoML:使用 KerasTuner 自动调参
- 特征选择:使用 SHAP 值解释特征重要性
7.3 策略层面
- 仓位管理:根据预测概率动态调整仓位
- 风控体系:设置止损线、最大回撤控制
- 实盘对接:对接券商 API 实现程序化交易
8. 总结
本文详细介绍了如何用 Python + LSTM 构建一个智能选股系统:
- 数据获取:使用 akshare 获取 A 股日线数据
- 特征工程:构建技术指标特征(MA、RSI、MACD)
- 模型训练:构建双层 LSTM 模型,进行二分类训练
- 回测验证:验证策略有效性
核心代码行数:约 150 行 测试集准确率:约 55-65%(取决于数据周期) 相对基准:跑赢随机选股(50%)
⚠️ 再次提醒:本文代码仅供学习参考,不构成投资建议。量化策略存在模型过拟合风险、市场风格变化风险,实盘前请充分回测并控制仓位。
参考资料
- akshare 文档 - A 股数据获取
- TensorFlow Keras LSTM - LSTM 层文档
- 《Python量化交易实战》 - 量化交易入门书籍
- 《动手学深度学习》 - 深度学习教材
相关文章推荐:
- 《量化交易"智能刹车系统":用 Python 构建动态仓位管理框架》
- 《Python 量化回测框架全景对比:2026 年主流框架实战评测》
- 《量化因子"雷达站":用 3 层信号过滤网,选股胜率提升 42%》