TensorFlow系列（4）——基于MNIST数据集的CNN实现

本文主要是尝试搭建一个简单的卷积神经网络（CNN）模型，并用它来训练MNIST数据集。

1. CNN简介

卷积神经网络（Convolutional neural network）属于人工神经网络的一种，它的权值共享的网络结构显著降低了模型的复杂度，减少了权值的数量，是目前语音分析和图像识别领域研究的热点。

和传统神经网络相比，卷积神经网络的特点在于隐藏层分为卷积层和池化层（pooling layer,又叫下采样层）。
（1）卷积层通过一块块卷积核在原始图像上平移来提取特征，每一个特征就是一个特征映射。
（2）池化层通过汇聚特征后稀疏参数来减少学习的参数，降低网络的复杂度，池化层最常见的包括最大值池化（max pooling）和平均值池化（average pooling）。

2. 基于MNIST数据集的CNN实现

实现代码如下所示：

# -*- coding: utf-8 -*-
"""
Created on Mon Jul 24 21:48:37 2017

@author: Administrator
"""

import tensorflow as tf
import numpy as np
from tensorflow.examples.tutorials.mnist import input_data


'''
1.构建模型
'''

#定义输入数据并预处理数据
mnist = input_data.read_data_sets("MNIST_data",one_hot=True)
trX,trY,teX,teY = mnist.train.images,mnist.train.labels,mnist.test.images,mnist.test.labels

trX = trX.reshape(-1,28,28,1)
teX = teX.reshape(-1,28,28,1)

X = tf.placeholder("float",[None,28,28,1])
Y = tf.placeholder("float",[None,10])

#初始化权重与定义网络结构
def init_weights(shape):
    return tf.Variable(tf.random_normal(shape,stddev=0.01))

w = init_weights([3,3,1,32])
w2 = init_weights([3,3,32,64])
w3 = init_weights([3,3,64,128])
w4 = init_weights([128*4*4,625])

w_o = init_weights([625,10])

'''
神经网络模型的构造函数，传入以下参数：
X:输入数据
Y:每一层的权重
p_keep_conv,p_keep_hidden:dropout 要保留的神经元比例
'''

def model(X,w,w2,w3,w4,w_o,p_keep_conv,p_keep_hidden):
    #第一层卷积层及池化层，
    l1a = tf.nn.relu(tf.nn.conv2d(X,w,strides=[1,1,1,1],padding='SAME'))

    l1 = tf.nn.max_pool(l1a,ksize=[1,2,2,1],strides=[1,2,2,1],padding='SAME')

    l1 = tf.nn.dropout(l1,p_keep_conv)

    #第二组卷积层及池化层，最后dropput一些神经元
    l2a = tf.nn.relu(tf.nn.conv2d(l1,w2,strides=[1,1,1,1],padding='SAME'))

    l2 = tf.nn.max_pool(l2a,ksize=[1,2,2,1],strides=[1,2,2,1],padding='SAME')

    l2 = tf.nn.dropout(l2,p_keep_conv)

    #第三组卷积层及池化层，最后dropput一些神经元
    l3a = tf.nn.relu(tf.nn.conv2d(l2,w3,strides=[1,1,1,1],padding='SAME'))

    l3 = tf.nn.max_pool(l3a,ksize=[1,2,2,1],strides=[1,2,2,1],padding='SAME')

    l3 = tf.reshape(l3,[-1,w4.get_shape().as_list()[0]])

    l3 = tf.nn.dropout(l3,p_keep_conv)

    #全连接层，最后dropout一些神经元
    l4 = tf.nn.relu(tf.matmul(l3,w4))
    l4 = tf.nn.dropout(l4,p_keep_hidden)

    #输出层
    pyx = tf.matmul(l4,w_o)
    return pyx

p_keep_conv = tf.placeholder("float")
p_keep_hidden = tf.placeholder("float")
#得到预测值
py_x = model(X,w,w2,w3,w4,w_o,p_keep_conv,p_keep_hidden)

#定义损失函数
cost = tf.reduce_mean(tf.nn.softmax_cross_entropy_with_logits(logits=py_x,labels=Y))
train_op = tf.train.RMSPropOptimizer(0.001,0.9).minimize(cost)
predict_op = tf.argmax(py_x,1)

'''
2.训练模型和评估模型
'''

#定义训练和评估批次大小
batch_size = 128
test_size = 256

#在会话中启动图
with tf.Session() as sess:
    tf.global_variables_initializer().run()

    for i in range(100):
        training_batch = zip(range(0,len(trX),batch_size),range(batch_size,len(trX)+1,batch_size))
        for start,end in training_batch:
            sess.run(train_op,feed_dict={X:trX[start:end],Y:trY[start:end],
                                         p_keep_conv:0.8,p_keep_hidden:0.5})

        test_indices = np.arange(len(teX))
        np.random.shuffle(test_indices)

        test_indices = test_indices[0:test_size]

        print(i,np.mean(np.argmax(teY[test_indices],axis=1) == 
                        sess.run(predict_op,feed_dict={X:teX[test_indices],
                                                       p_keep_conv:1.0,
                                                       p_keep_hidden:1.0})))

运行部分结果如下图所示：