num_epochs=2000

如果您想要训练2000个epoch的神经网络模型，可以按以下方式进行修改：

import numpy as np
import pandas as pd

# 定义sigmoid激活函数
def sigmoid(x):
    return 1 / (1 + np.exp(-x))

# 定义神经网络模型
def neural_network(X, W1, W2, b1, b2):
    z1 = np.dot(X, W1) + b1
    a1 = sigmoid(z1)
    z2 = np.dot(a1, W2) + b2
    return z2

# 加载训练数据
train_data = pd.read_csv('train_data.csv')
X_train = np.array(train_data.iloc[:, :-1])
y_train = np.array(train_data.iloc[:, -1])

# 对训练数据进行预处理
X_train = (X_train - np.mean(X_train, axis=0)) / np.std(X_train, axis=0)

# 初始化模型参数
input_size = X_train.shape[1]
hidden_size = 10
output_size = 1
W1 = np.random.randn(input_size, hidden_size)
W2 = np.random.randn(hidden_size, output_size)
b1 = np.zeros((1, hidden_size))
b2 = np.zeros((1, output_size))

# 定义超参数
learning_rate = 0.01
num_epochs = 2000

# 训练模型
for epoch in range(num_epochs):
    # 前向传播，计算预测值
    z1 = np.dot(X_train, W1) + b1
    a1 = sigmoid(z1)
    z2 = np.dot(a1, W2) + b2
    y_pred = z2.ravel()
    
    # 计算损失函数值
    loss = np.mean((y_pred - y_train) ** 2)
    
    # 反向传播，更新模型参数
    delta2 = (y_pred - y_train).reshape(-1, 1)
    delta1 = delta2.dot(W2.T) * sigmoid(z1) * (1 - sigmoid(z1))
    dW2 = np.dot(a1.T, delta2)
    db2 = np.sum(delta2, axis=0, keepdims=True)
    dW1 = np.dot(X_train.T, delta1)
    db1 = np.sum(delta1, axis=0)
    W2 -= learning_rate * dW2
    b2 -= learning_rate * db2
    W1 -= learning_rate * dW1
    b1 -= learning_rate * db1
    
    # 打印训练过程中的损失函数值
    if epoch % 100 == 0:
        print('Epoch [{}/{}], Loss: {:.4f}'.format(epoch+1, num_epochs, loss))

# 保存模型参数
np.save('W1.npy', W1)
np.save('W2.npy', W2)
np.save('b1.npy', b1)
np.save('b2.npy', b2)

在这段代码中，我们将 `num_epochs` 的值设置为 2000。在训练过程中，每经过100个epoch，就打印一次当前的损失函数值。训练完成后，将模型参数保存到相应的文件中。请注意，在实际应用中，需要根据数据集的大小和模型的复杂度来调整 `num_epochs` 的值，以避免过拟合或欠拟合等问题。