这篇文章记录我所看的视频课【tensorflow笔记】的代码，这些代码并不是课上讲的代码，而是我自己听了课之后根据自己的想法敲写的，所以不是100%还原老师的代码。
1.前向传播
一个很简单很简单的前向传播例子
import numpy as np
import tensorflow as tf
#生成输入数据
x_data=np.random.rand(100)
y_data=0.2*x_data+0.1
#前向传播
k=0.5
b=0.3
x=tf.placeholder(tf.float32)
y=k*x+b
#初始化
init=tf.global_variables_initializer()
with tf.Session() as sess:
sess.run(init)
print(sess.run(y,feed_dict={x:x_data}))2.反向传播
一个很简单很简单的反向传播的例子
import numpy as np
import tensorflow as tf
#生成输入数据
x_data=np.random.rand(100)
y_data=0.2*x_data+0.1
#前向传播
k=tf.Variable(0.5)
b=tf.Variable(0.3)
y=tf.Variable(0)
y=k*x_data+b
#反向传播
error=y-y_data
loss=tf.reduce_mean(tf.square(error))
optimizer=tf.train.GradientDescentOptimizer(0.2).minimize(loss)
#初始化
init=tf.global_variables_initializer()
with tf.Session() as sess:
sess.run(init)
for i in range(200):
sess.run(optimizer)
if i%20==0:
print(sess.run([k,b,loss]))下面也是反向传播，但这个代码比上面那个有难度提升,因为我们开始使用权重矩阵
import numpy as np
import tensorflow as tf
#输入数据
x_data= np.random.rand(32,2)
y_data=[[x1+x2+(rmd.rand()/10.0-0.05)] for (x1,x2) in x_data]
"""这行代码与下句等同
y=[]
for [x1,x2] in x:
y.append([x1+x2+(np.random.rand()/10-0.05)])
"""
#前向传播
x=tf.placeholder(tf.float32,shape=(None,2))
y=tf.placeholder(tf.float32,shape=(None,1))
w1= tf.Variable(tf.random_normal([2, 1]))#正态分布
z= tf.matmul(x, w1)#点积
#反向传播
loss=tf.reduce_mean(tf.square(z-y))
optimizer=tf.train.GradientDescentOptimizer(0.002).minimize(loss)
#初始化
init=tf.global_variables_initializer()
with tf.Session() as sess:
sess.run(init)
for i in range(20000):
sess.run(optimizer,feed_dict={x:x_data,y:y_data})
if i%5000==0:
print(f"第{i}次训练结果：")
print(sess.run(w1))
print("最终结果：\n",sess.run(w1))3.损失函数
损失函数很多，你可以自己写，也可以用经典的损失函数（线性回归、sigmoid、softmax）
经典的线性回归上面两个代码我们已经用了，这次写一下自己定义的损失函数：
import tensorflow as tf
w=tf.Variable(5.)
loss=tf.square(w+1)
optimizer=tf.train.GradientDescentOptimizer(0.2).minimize(loss)
inti=tf.global_variables_initializer()
with tf.Session() as sess:
sess.run(inti)
for i in range(50):
sess.run(optimizer)
print(sess.run([w,loss]))可以看到，loss是自己定义的函数，而不是经典的三大金刚。loss可以随便定义，你想要的啥就写啥，这里方便理解所以写了个最简单的loss=（w+1）^2，在实际问题中肯定不会这么简单。
4.学习率
这个就不写代码了，就是在使用各种优化器的时候要设定学习率，就是反向传播一次更新的程度大小。
学习率的设置也有讲究，小了更新慢，大了会爆炸，所以我们说一下怎么设置学习率：
在实际编码时，我们通常使用【指数衰减学习率】：
这是指数衰减学习率的公式，你可以不太懂它的含义，只要你会敲代码就可以
代码如下：
globe_step=tf.Variable(0，trainable=False) #计数器，用来记录当前运行了几轮batch_size，该参数不参与训练，所以trainable=False
learning_rate=tf.train.exponential_decay(
LEARNING_RATE_BASE, #学习率初始值
globe_step, #计数器
LEARNING_RATE_STEP, #学习率更新周期，值为总样本数÷batch_size样本数
LEARNING_RATE_DECAY, #学习率衰减率
staircase=True #True表示以不连续的间隔衰减学习速率，最后曲线是锯齿状
False表示以连续的间隔衰减学习速率，最后曲线是光滑状
)
在训练时，要求在tf.train.GradientDescentOptimizer(learning_rate).minimize(loss,global_step)的minimize里写上global_step
实际应用如下：
import tensorflow as tf
#设置learning_rate
global_step=tf.Variable(0，trainable=False)
LEARNING_RATE_BASE=0.1
LEARNING_RATE_STEP=1
LEARNING_RATE_DECAY=0.99
learning_rate=tf.train.exponential_decay(LEARNING_RATE_BASE,global_step,LEARNING_RATE_STEP,LEARNING_RATE_DECAY,staircase=True)
#设置训练数据
w=tf.Variable(5.)
#设置损失函数
loss=tf.square(w+1)
#设置优化器
optimize=tf.train.GradientDescentOptimizer(learning_rate).minimize(loss,global_step)
#初始化
init=tf.global_variables_initializer()
with tf.Session() as sess:
sess.run(init)
for i in range(40):
sess.run(optimize)
learning_rate_val=sess.run(learning_rate)
global_step_val=sess.run(global_step)
w_val=sess.run(w)
loss_val=sess.run(loss)
print(f"经过{i}步，当前global_step为{ global_step_val}，w为{w_val}，learning_rate为{learning_rate_val}，loss为{loss_val}")一定多多留意独特之处！
5.正则化
我们用正则化来缓解泛化效果差的问题，就是给权重一个矩阵，加在loss上。
如图所示，正则化后的loss，是原来的loss再加上正则化项，正则化项就是正则化参数*loss（w），这个
，跟吴恩达讲的一模一样。
正则化项的表示代码为tf.contrib.layers.l2_regularizer（lambda）(w)，其中lambda是正则化参数，w是权重
tf.add_to_collection("losses",tf.contrib.layers.l2_regularizer(lambda)(w)) #把正则化项放入一个名为"losses"的集合器里面，我们打算用这个losses来表示整个正则化项了。
loss=cem+tf.add_n(tf.get_collection("losses")) #cem是原来的loss，现在把正则化项和之前的最原始的loss加在一起，就得到了正则化后的loss
ok，你可能看不懂，问“为什么要这么麻烦？直接把正则化项加进cem不就行了吗？”是的，我也不懂为什么，但就是这样的一环套一环，你还就真得按照这个代码去操作，不然你做不出来。
下面看代码：
import tensorflow as tf
import numpy as np
import matplotlib.pyplot as plt
#设置batch（临时样本）
step=30
#产生输入数据
x_data=np.random.randn(300,2)
#产生输入数据集的标签（正确答案）
y_data=[int(x0*x0+x1*x1<2) for (x0,x1) in x_data]
#遍历y_data中的每个元素，1赋值red，其余赋值blue，这样可视化显示时容易区分
y_s=[["red" if y else "blue"] for y in y_data]
#数据变形，按行堆叠成一个新数组
x_data=np.vstack(x_data).reshape(-1,2)#numpy中，shape里写-1，相当于tensorflow里写None
y_data=np.vstack(y_data).reshape(-1,1)
print(x_data)
print(y_data)
print(y_s)
#用plt.scatter画出数据集x_data各行中第0列元素和第一列元素的点即各行的（x0,x1），用各行y_s对应的值表示颜色
plt.scatter(x_data[:,0],x_data[:,1],c=np.squeeze(y_s))
plt.show()
#定义神经网络的输入、参数和输出，定义前向传播
def get_weight(shape,regularizer):
w=tf.Variable(tf.random_normal(shape),dtype=tf.float32)
tf.add_to_collection("losses",tf.contrib.layers.l2_regularizer(regularizer)(w))
return w
def get_bias(shape):
b=tf.Variable(tf.constant(0.01,shape=shape))
return b
#输入层输入
x=tf.placeholder(tf.float32,shape=(None,2))
y_=tf.placeholder(tf.float32,shape=(None,1))
#输入层到隐藏层
w1=get_weight([2,11],0.01)
b1=get_bias([11])
y1=tf.nn.relu(tf.matmul(x,w1)+b1)
#隐藏层到输出层
w2=get_weight([11,1],0.01)
b2=get_bias([1])
y=tf.matmul(y1,w2)+b2#输出层不激活（仅限这个代码不激活）
#定义损失函数
loss_mse=tf.reduce_mean(tf.square(y-y_))
loss=loss_mse + tf.add_n(tf.get_collection("losses"))
#定义反向传播
optimizer=tf.train.AdamOptimizer(0.0001).minimize(loss)
#初始化
init=tf.global_variables_initializer()
with tf.Session() as sess:
sess.run(init)
for i in range(40000):
start=(i*step)%300
end=start+step
sess.run(optimizer,feed_dict={x:x_data[start:end],y_:y_data[start:end]})
if i%2000==0:
loss_mse_val=sess.run(loss_mse,feed_dict={x:x_data,y_:y_data})
print(f"经过{i}步，loss现为：{loss_mse_val}")
# xx在-3到3之间以步长为0.01，yy在-3到3之间以步长0.01,生成二维网格坐标点
xx, yy = np.mgrid[-3:3:.01, -3:3:.01]
# 将xx , yy拉直，并合并成一个2列的矩阵，得到一个网格坐标点的集合
grid = np.c_[xx.ravel(), yy.ravel()]
# 将网格坐标点喂入神经网络 ，probs为输出
probs = sess.run(y, feed_dict={x: grid})
# probs的shape调整成xx的样子
probs = probs.reshape(xx.shape)
print("w1:\n", sess.run(w1))
print("b1:\n", sess.run(b1))
print("w2:\n", sess.run(w2))
print("b2:\n", sess.run(b2))
plt.scatter(x_data[:, 0], x_data[:, 1], c=np.squeeze(y_s))
plt.contour(xx, yy, probs, levels=[.5])
plt.show()
扫一扫
207
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
