写一个基于lstm的股价预测代码
时间: 2024-01-25 19:03:11 浏览: 76
好的,这是一个基于LSTM的股价预测代码的示例。需要注意的是,股价预测是一个复杂的问题,数据的质量和模型的选择都非常重要,本示例仅供参考。
首先,需要导入必要的库:
```python
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from sklearn.preprocessing import MinMaxScaler
from tensorflow.keras.models import Sequential
from tensorflow.keras.layers import Dense, LSTM
```
接下来,读取股票数据并进行预处理:
```python
# 读取股票数据
df = pd.read_csv('stock_price.csv')
# 只选择收盘价一列,并将数据转换为numpy数组
data = df.filter(['Close']).values
# 将数据分为训练集和测试集
train_size = int(len(data) * 0.8)
train_data = data[:train_size]
test_data = data[train_size:]
# 对数据进行标准化处理
scaler = MinMaxScaler(feature_range=(0, 1))
train_data = scaler.fit_transform(train_data)
test_data = scaler.transform(test_data)
```
接下来,准备训练数据,这里采用时间序列的方法,用前60天的收盘价来预测第61天的收盘价:
```python
# 准备训练数据
look_back = 60
def create_dataset(data):
X, y = [], []
for i in range(look_back, len(data)):
X.append(data[i-look_back:i, 0])
y.append(data[i, 0])
return np.array(X), np.array(y)
train_X, train_y = create_dataset(train_data)
test_X, test_y = create_dataset(test_data)
# 将数据转换为3D格式
train_X = np.reshape(train_X, (train_X.shape[0], train_X.shape[1], 1))
test_X = np.reshape(test_X, (test_X.shape[0], test_X.shape[1], 1))
```
然后,构建LSTM模型,并进行训练:
```python
# 构建LSTM模型
model = Sequential()
model.add(LSTM(50, return_sequences=True, input_shape=(train_X.shape[1], 1)))
model.add(LSTM(50, return_sequences=False))
model.add(Dense(25))
model.add(Dense(1))
# 编译模型
model.compile(optimizer='adam', loss='mean_squared_error')
# 训练模型
model.fit(train_X, train_y, batch_size=1, epochs=1)
```
最后,使用测试数据进行预测,并将预测结果可视化:
```python
# 使用测试数据进行预测
predictions = model.predict(test_X)
predictions = scaler.inverse_transform(predictions)
# 计算均方根误差
rmse = np.sqrt(np.mean(((predictions - test_y) ** 2)))
print(f'RMSE: {rmse}')
# 绘制预测结果
train = df[:train_size]
valid = df[train_size:]
valid['Predictions'] = predictions
plt.figure(figsize=(16, 8))
plt.title('Model')
plt.xlabel('Date', fontsize=18)
plt.ylabel('Close Price USD ($)', fontsize=18)
plt.plot(train['Close'])
plt.plot(valid[['Close', 'Predictions']])
plt.legend(['Train', 'True', 'Predictions'], loc='lower right')
plt.show()
```
这样,一个基于LSTM的股价预测代码就完成了。需要注意的是,这只是一个简单的示例代码,真正的应用需要更加复杂的数据预处理和模型调优。
阅读全文