算法的开始
import numpy as np
import math
import random
x = np.array([[2104], [1600], [2400], [1416], [3000]])
x1 = np.array([[3],[3],[3],[2],[4]])
t = np.array([[400], [330], [369], [232], [540]])
a = 0.00100
error = 0.00001
dert = random.random()
dert1 = random.random()
dert2 = random.random()
sit = random.random()
sit2 = random.random()
sit3 = random.random()
while(dert>error or dert1>error or dert2>error):
for i in range(0, 5):
dert = (sit+sit2*x[i]+sit3*x1[i]-t[i])*1
dert1 = (sit + sit2 * x[i] + sit3 * x1[i] - t[i]) * x[i]
dert2 = (sit + sit2 * x[i] + sit3 * x1[i] - t[i]) * x1[i]
sit = sit - dert * a
sit2 = sit2 - dert1*a
sit3 = sit3 - dert2*a
print(dert)
print(dert1)
print(dert2)
print(sit)
print(sit2)
print(sit3)
x3 = 3600
x4 = 5
h = sit + sit2*x3 + sit3*x4
print(h)