import tensorflow as tf
import numpy as np
import matplotlib.pyplot as plt
import os
os.environ['TF_CPP_MIN_LOG_LEVEL'] = '2'
#图形显示
def moving_average(a,w=10):
if len(a)<w:
return a[:]
return [val if idx <w else sum(a[(idx-w):idx])/w for idx,val in enumerate(a)]
X = tf.placeholder("float")
Y = tf.placeholder("float")
# 占位符
inputsict = {
'x': tf.placeholder("float"),
'y': tf.placeholder("float")
}
# 生成模拟数据
train_X = np.float32(np.linspace(-1, 1, 100))
train_Y=2*train_X+np.random.randn(*train_X.shape)*0.1 #y=2x
plt.plot(train_X,train_Y,'ro',label='data')
plt.show()
#直接定义模型参数
W=tf.Variable(tf.random_normal([1]),name="weight")
b=tf.Variable(tf.zeros([1]),name="bias")
#字典定义学习参数
paradict={
'W':tf.Variable(tf.random_normal([1])),
'b':tf.Variable(tf.zeros([1]))
}
#前向结构
Z=tf.multiply(X,paradict['W'])+paradict['b']
#反向传播优化:
cost=tf.reduce_mean(tf.square(Y-Z))
learning_rate=0.01
optimizer=tf.train.GradientDescentOptimizer(learning_rate).minimize(cost)
#初始化所有变量
init =tf.global_variables_initializer()
#定义参数
training_epochs=20
display_step=2
#启动Session
with tf.Session() as sess:
sess.run(init)
plotdata={"batchsize":[],"loss":[]}
#向模型输入数据
for epoch in range(training_epochs):
for (x,y) in zip(train_X,train_Y):
sess.run(optimizer,feed_dict={X:x,Y:y})
#显示训练中的详细信息
if epoch % display_step==0:
loss=sess.run(cost,feed_dict={X:train_X,Y:train_Y})
print("Epoch:",epoch+1,"cost=",loss,"W=",sess.run(paradict['W']),
"b=",sess.run(paradict['b']))
if not (loss=="NA"):
plotdata["batchsize"].append(epoch)
plotdata["loss"].append(loss)
print("Finished!")
print("cost=",sess.run(cost,feed_dict={X:train_X,Y:train_Y}),
"W=",sess.run(paradict['W']),"b=",sess.run(paradict['b']))
plt.plot(train_X,train_Y,'ro',label='Original data')
plt.plot(train_X,sess.run(paradict['W'])*train_X+sess.run(paradict['b']),label='Fittedline')
plt.legend()
plt.show()
plotdata["avgloss"]=moving_average(plotdata["loss"])
plt.figure(1)
plt.subplot(211)
plt.plot(plotdata["batchsize"],plotdata["avgloss"],'b--')
plt.xlabel('Minibatch number')
plt.ylabel('loss')
plt.title('Minibatch run VS. Training Loss')
plt.show()
""""
#图形显示
def moving_average(a,w=10):
if len(a)<w:
return a[:]
return [val if idx <w else sum(a[(idx-w):idx])/w for idx,val in enumerate(a)]
"""