深度学习之前馈神经网络

1.导入常用工具包

#在终端中输入以下命令就可以安装工具包
pip install numpy
pip install pandas
Pip install matplotlib

注：
numpy是科学计算基础包
pandas能方便处理结构化数据和函数
matplotlib主要用于绘制图表。

#导包的代码：
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt

2.导入数据集

2.1.数据集的下载与转换

数据集下载地址：http://yann.lecun.com/exdb/mnist/

注：下载这四个文件，但由于不是csv的格式，所以我们要将这四个文件转换成csv格式。

import struct 
import numpy as np 
import pandas as pd

#以二进制读取模式。struct.unpack('>IIII', f.read(16))从文件中读取前16个字节，并按照大端字节序解析出魔数、图像数量、行数和列数。
def read_mnist_image(filename):
    with open(filename, 'rb') as f:
        magic_number, num_images, num_rows, num_cols = struct.unpack('>IIII', f.read(16))
        image_data = np.fromfile(f, dtype=np.uint8).reshape(num_images, num_rows * num_cols)
    return image_data

#读取MNIST数据集中的标签数据。与读取图像数据类似，它打开文件，解析魔数和标签数量，然后读取剩余的数据，将其转换为NumPy数组。
def read_mnist_labels(filename):
    with open(filename, 'rb') as f:
        magic_number, num_labels = struct.unpack('>II', f.read(8))
        label_data = np.fromfile(f, dtype=np.uint8)
    return label_data

# 读取图像和标签文件
image_filename = 'train-images.idx3-ubyte'
label_filename = 'train-labels.idx1-ubyte'
images = read_mnist_image(image_filename)
labels = read_mnist_labels(label_filename)

# 将图像和标签合并为一个DataFrame
train_Data = pd.DataFrame(images)
train_Data['label'] = labels

# 保存为CSV文件
train_Data.to_csv('mnist_train.csv', index=False)

# 对测试数据进行相同的操作
image_filename = 't10k-images.idx3-ubyte'
label_filename = 't10k-labels.idx1-ubyte'
images = read_mnist_image(image_filename)
labels = read_mnist_labels(label_filename)
test_Data = pd.DataFrame(images)
test_Data['label'] = labels
test_Data.to_csv('mnist_test.csv', index=False)

数据集图片：

2.2.数据观察

import pandas as pd
train_Data = pd.read_csv('mnist_train.csv',header = None)
test_Data = pd.read_csv('mnist_test.csv',header = None)
print("Train data:")
train_Data.info()
print("Test data:")
test_Data.info()
#继续观察训练数据前五行
train_Data.head(5)

注：可以发现训练数据中包含60000个数据样本，维度785，包括标签信息与784个特征维度；测试数据中包含10000个样本，维度785，包括标签信息与784个特征维度。
运行结果：

2.3.读取第一行数据

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
#取第一行数据
x = train_Data.iloc[0]
#标签信息
y = x[0]
#将1*784转换成28*28
img = x[1:].values.reshape(28,28)
#画图
plt.imshow(img)
plt.title('label = ' + str(y))
plt.show()

注：这段代码的主要作用是从训练数据集中取出第一行数据，将其中的图像数据转换为28x28的二维数组，并使用matplotlib库显示这个图像。同时，它还展示了图像对应的标签（即手写数字的类别）。
运行结果：

2.4.从sklearn中导入数据并观察

from sklearn.datasets import fetch_openml
mnist = fetch_openml('mnist_784', version=1)
mnist.keys()

运行结果：

2.5.打印维度信息

data, label = mnist["data"], mnist["target"]
print("数据维度:", data.shape)
print("标签为度:", label.shape)

注：这段代码从之前使用fetch_openml函数获取的MNIST数据集字典中提取出图像数据（data）和标签（label），并打印它们的维度信息。

2.6.显示第一幅图像及其对应的标签

x = data.iloc[0]
y = label[0] #从label数组中取出第一个元素，即第一幅图像的标签
img = x.values.reshape(28,28)
plt.imshow(img)
plt.title('label = ' + str(y))
plt.show()

运行结果：

3.数据预处理

X = train_Data.iloc[:,1:].values#样本数据
y = train_Data.iloc[:,0].values#样本标签
print("数据X中最大值:",X.max())
print("数据X中最小值:",X.min())

运行结果：

3.1.对X进行归一化处理

#归一化
X = X/255
#此时将数值大小缩小在[0，1]范围内，重现观察数据中的最大、最小值
print("数据X中最大值:",X.max())
print("数据X中最小值:",X.min())

注：这段代码的目的是对数据集X进行归一化处理，并将处理后的数据范围缩放到[0, 1]之间。
运行结果：

3.2.分类

X_valid, X_train = X[:5000], X[5000:]
y_valid, y_train = y[:5000], y[5000:]

X_test,y_test = test_Data.iloc[:,1:].values/255, test_Data.iloc[:,0].values

注：将数据集分割为训练集、验证集和测试集，并对这些数据集进行归一化处理。

4.前馈全连接神经网络（Sequential模型）

4.1.创建模型

import tensorflow as tf
from tensorflow import keras
model = keras.models.Sequential([
    keras.layers.Flatten(input_shape=[784]),#输入层784个神经元
    keras.layers.Dense(300, activation="relu"),#隐藏层300个神经元
    keras.layers.Dense(100, activation="relu"),#隐藏层100个神经元
    keras.layers.Dense(10, activation="softmax")#输入层10个神经元
])
model.layers[1]
weight_l,bias_l = model.layers[1].get_weights()
print(weight_l.shape)
print(bias_l.shape)

注：通过打印权重和偏置的形状，可以确认模型的第一个隐藏层是否正确地连接到输入层，并且可以了解该层的参数数量。这对于调试和理解模型的结构非常重要。在实际应用中，这些权重和偏置会在模型训练过程中通过反向传播算法自动调整，以最小化预测误差。
运行结果：

4.2.通过.summary()观察神经网络的整体情况

model.summary()

注：model.summary()是Keras模型的一个方法，它用于打印出模型的概述信息。
运行结果：

4.3.训练网格

#编译网络
model.compile(loss="sparse_categorical_crossentropy",optimizer="sgd",metrics=["accuracy"])

注：是Keras中的一个重要步骤，它用于编译刚刚创建的神经网络模型。编译过程定义了模型训练时需要使用的损失函数、优化器和评估指标。

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
h = model.fit(X_train,y_train,batch_size=32,epochs=30,validation_data=(X_valid,y_valid))

注：epochs=30：训练过程中遍历整个训练数据集的次数。每个时期包含一次完整的训练数据遍历。
运行结果：

4.4.将Keras的History对象转换为Pandas的DataFrame

pd.DataFrame(h.history)

注：h.history属性是一个字典，其中包含了训练过程中的损失值和评估指标值。
运行结果：

4.5.绘图

#绘图
pd.DataFrame(h.history).plot(figsize=(8, 5))
plt.grid(True)
plt.gca().set_ylim(0, 1)#set the vertical range to [0-1]
plt.show()

注：
Loss（损失）：损失函数的值越低，表示模型的预测越接近实际值。
Accuracy（准确率）：准确率是指模型正确预测的样本数与总样本数之间的比例。
Val Loss（验证损失）：如果训练损失持续下降，但验证损失开始上升，这可能表明模型出现了过拟合。
Val Accuracy（验证准确率）：用于评估模型的泛化能力，并且是模型性能的一个重要指标。
运行结果：

4.6.识别准确率

model.evaluate(X_test, y_test, batch_size = 1)

注：使用Keras模型的evaluate方法来评估模型在测试集上的性能。evaluate方法会计算并返回模型在给定测试数据上的损失和评估指标。
运行结果：

4.7.对样本进行预测

x_sample, y_sample = X_test[11:12], y_test[11]
y_prob = model.predict(x_sample).round(2)
y_prob

img = x_sample.reshape(28,28)
plt.imshow(img)
plt.title('label = ' + str(np.argmax(y_prob)))
plt.show()

注：从测试数据集中选取索引为11的单个样本，并将其特征和标签分别存储在x_sample和y_sample中。这里使用切片[11:12]来确保x_sample是一个二维数组，符合模型的输入要求。
运行结果：

5.使用Sequential()方法，对鸢尾花数据集进行分类

5.1划分

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from sklearn.datasets import load_iris
from sklearn.model_selection import train_test_split

iris = load_iris()  #鸢尾花数据集

x_train,x_test,y_train,y_test = train_test_split(iris.data,iris.target,test_size=0.2,random_state=23)

X_train,X_valid,y_train,y_valid = train_test_split(x_train,y_train,test_size=0.2,random_state=12)

print(X_valid.shape)
print(X_train.shape)

注：从Scikit-learn的model_selection模块中导入train_test_split函数，用于分割数据集。
运行结果：

5.2. 构建模型

import tensorflow as tf
from tensorflow import keras

model = keras.models.Sequential([
    keras.layers.Flatten(input_shape=[4]),
    keras.layers.Dense(16,activation='relu'),
    keras.layers.Dense(16,activation='relu'),
    keras.layers.Dense(16,activation='relu'),
    keras.layers.Dense(16,activation='relu'),
    keras.layers.Dense(16,activation='relu'),
    keras.layers.Dense(16,activation='relu'),
    keras.layers.Dense(16,activation='relu'),
    keras.layers.Dense(16,activation='relu'),
    keras.layers.Dense(16,activation='relu'),
    keras.layers.Dense(16,activation='relu'),
    keras.layers.Dropout(rate=0.2),
    keras.layers.Dense(3,activation='softmax'),
])

model.summary()

注：这段代码使用TensorFlow和Keras库创建了一个神经网络模型，用于分类鸢尾花数据集。模型的结构是顺序的，包含了多个全连接层（Dense layers）和一个Dropout层。
运行结果：

5.3.提高准确率添加方式：keras.layers.Dropout(rate=0.2)

model.layers[1]

注：从之前定义的Keras模型中获取第二个层的对象。在Keras模型中，层是按照它们添加到模型中的顺序存储在一个列表中的，索引从0开始。因此，model.layers[1]将返回模型中第一个隐藏层的对象。

weight_1,bias_1 = model.layers[1].get_weights()

print(weight_1.shape)
print(bias_1.shape)

运行结果：

注：从之前定义的Keras模型中获取第一个隐藏层的权重和偏置，并打印它们的形状。

model.compile(loss='sparse_categorical_crossentropy',
             optimizer='sgd',metrics=["accuracy"])

h = model.fit(X_train,y_train,batch_size=10,epochs=50,validation_data=(X_valid,y_valid))

注：optimizer=‘sgd’：这是模型训练时使用的优化器。sgd代表随机梯度下降（Stochastic Gradient Descent），它是一种简单的优化算法，用于在训练过程中更新模型的权重。
运行结果：

pd.DataFrame(h.history)

运行结果：

pd.DataFrame(h.history).plot(figsize=(8,5))
plt.grid(True)
plt.gca().set_ylim(0,1)
plt.show()

运行结果：

model.evaluate(x_test,y_test,batch_size = 1)

注：
x_test：测试数据集的特征，通常是NumPy数组或TensorFlow张量。
y_test：测试数据集的标签，与x_test中的每个样本相对应。
batch_size = 1：评估过程中每次前向传播所使用的数据样本数量。
运行结果：