写一个lstm预测共享单车流量的代码
时间: 2023-12-01 17:31:37 浏览: 74
好的,以下是一个简单的LSTM模型,用于预测共享单车的流量。
首先,我们需要导入所需的库和数据集:
```python
import pandas as pd
import numpy as np
from sklearn.preprocessing import MinMaxScaler
from tensorflow.keras.models import Sequential
from tensorflow.keras.layers import Dense, LSTM
# 读取数据集
data = pd.read_csv('bike_sharing.csv')
```
接下来,我们需要对数据进行一些预处理。我们将使用`MinMaxScaler`来缩放数据,并将其分成训练和测试集。
```python
# 将数据缩放到0和1之间
scaler = MinMaxScaler(feature_range=(0, 1))
scaled_data = scaler.fit_transform(data['cnt'].values.reshape(-1, 1))
# 分割数据集为训练和测试集
training_size = int(len(scaled_data) * 0.7)
test_size = len(scaled_data) - training_size
train_data = scaled_data[0:training_size, :]
test_data = scaled_data[training_size:len(scaled_data), :]
```
接下来,我们将创建一个函数,用于将数据集转换为用于LSTM模型的数据格式。该函数将接受两个参数:数据集和时间步长,它将返回X和Y数据集,其中X是前N天的共享单车流量,Y是第N + 1天的共享单车流量。
```python
def create_dataset(dataset, time_step=1):
X, Y = [], []
for i in range(len(dataset) - time_step - 1):
a = dataset[i:(i + time_step), 0]
X.append(a)
Y.append(dataset[i + time_step, 0])
return np.array(X), np.array(Y)
```
接下来,我们将使用上面的函数来准备我们的训练和测试数据集。
```python
# 创建训练和测试数据集
time_step = 100
X_train, Y_train = create_dataset(train_data, time_step)
X_test, Y_test = create_dataset(test_data, time_step)
# 重塑输入为LSTM的3D格式 [样本数,时间步长,特征数]
X_train = X_train.reshape(X_train.shape[0], X_train.shape[1], 1)
X_test = X_test.reshape(X_test.shape[0], X_test.shape[1], 1)
```
现在,我们可以开始构建我们的LSTM模型。我们将使用一个包含两个LSTM层和一个Dense层的简单模型。
```python
# 创建LSTM模型
model = Sequential()
model.add(LSTM(50, return_sequences=True, input_shape=(time_step, 1)))
model.add(LSTM(50, return_sequences=True))
model.add(LSTM(50))
model.add(Dense(1))
model.compile(loss='mean_squared_error', optimizer='adam')
```
现在,我们可以训练我们的模型。
```python
# 训练模型
model.fit(X_train, Y_train, validation_data=(X_test, Y_test), epochs=100, batch_size=64, verbose=1)
```
最后,我们可以使用测试数据集来评估我们的模型,并进行预测。
```python
# 测试集上的预测
train_predict = model.predict(X_train)
test_predict = model.predict(X_test)
# 将预测结果反缩放为原始值
train_predict = scaler.inverse_transform(train_predict)
Y_train = scaler.inverse_transform([Y_train])
test_predict = scaler.inverse_transform(test_predict)
Y_test = scaler.inverse_transform([Y_test])
# 计算均方根误差
train_score = np.sqrt(mean_squared_error(Y_train[0], train_predict[:, 0]))
print('Train RMSE: %.2f' % (train_score))
test_score = np.sqrt(mean_squared_error(Y_test[0], test_predict[:, 0]))
print('Test RMSE: %.2f' % (test_score))
```
这就是使用LSTM模型预测共享单车流量的代码。
阅读全文