Python基于numpy灵活定义神经网络结构的方法

发布时间 - 2026-01-11 02:52:15    点击率:

本文实例讲述了Python基于numpy灵活定义神经网络结构的方法。分享给大家供大家参考,具体如下:

用numpy可以灵活定义神经网络结构,还可以应用numpy强大的矩阵运算功能!

一、用法

1). 定义一个三层神经网络:

'''示例一'''
nn = NeuralNetworks([3,4,2]) # 定义神经网络
nn.fit(X,y) # 拟合
print(nn.predict(X)) #预测

说明:
  输入层节点数目:3
  隐藏层节点数目:4
  输出层节点数目:2

2).定义一个五层神经网络:

'''示例二'''
nn = NeuralNetworks([3,5,7,4,2]) # 定义神经网络
nn.fit(X,y) # 拟合
print(nn.predict(X)) #预测

说明:
  输入层节点数目:3
  隐藏层1节点数目:5
  隐藏层2节点数目:7
  隐藏层3节点数目:4
  输出层节点数目:2

二、实现

如下实现方式为本人(@hhh5460)原创。 要点: dtype=object

import numpy as np
class NeuralNetworks(object):
  ''''''
  def __init__(self, n_layers=None, active_type=None, n_iter=10000, error=0.05, alpha=0.5, lamda=0.4):
    '''搭建神经网络框架'''
    # 各层节点数目 (向量)
    self.n = np.array(n_layers) # 'n_layers必须为list类型,如:[3,4,2] 或 n_layers=[3,4,2]'
    self.size = self.n.size # 层的总数
    # 层 (向量)
    self.z = np.empty(self.size, dtype=object) # 先占位(置空),dtype=object !如下皆然
    self.a = np.empty(self.size, dtype=object)
    self.data_a = np.empty(self.size, dtype=object)
    # 偏置 (向量)
    self.b = np.empty(self.size, dtype=object)
    self.delta_b = np.empty(self.size, dtype=object)
    # 权 (矩阵)
    self.w = np.empty(self.size, dtype=object)
    self.delta_w = np.empty(self.size, dtype=object)
    # 填充
    for i in range(self.size):
      self.a[i] = np.zeros(self.n[i]) # 全零
      self.z[i] = np.zeros(self.n[i]) # 全零
      self.data_a[i] = np.zeros(self.n[i]) # 全零
      if i < self.size - 1:
        self.b[i] = np.ones(self.n[i+1])  # 全一
        self.delta_b[i] = np.zeros(self.n[i+1]) # 全零
        mu, sigma = 0, 0.1 # 均值、方差
        self.w[i] = np.random.normal(mu, sigma, (self.n[i], self.n[i+1])) # # 正态分布随机化
        self.delta_w[i] = np.zeros((self.n[i], self.n[i+1])) # 全零

下面完整代码是我学习斯坦福机器学习教程,完全自己敲出来的:

import numpy as np
'''
参考:http://ufldl.stanford.edu/wiki/index.php/%E7%A5%9E%E7%BB%8F%E7%BD%91%E7%BB%9C
'''
class NeuralNetworks(object):
  ''''''
  def __init__(self, n_layers=None, active_type=None, n_iter=10000, error=0.05, alpha=0.5, lamda=0.4):
    '''搭建神经网络框架'''
    self.n_iter = n_iter # 迭代次数
    self.error = error # 允许最大误差
    self.alpha = alpha # 学习速率
    self.lamda = lamda # 衰减因子 # 此处故意拼写错误!
    if n_layers is None:
      raise '各层的节点数目必须设置!'
    elif not isinstance(n_layers, list):
      raise 'n_layers必须为list类型,如:[3,4,2] 或 n_layers=[3,4,2]'
    # 节点数目 (向量)
    self.n = np.array(n_layers)
    self.size = self.n.size # 层的总数
    # 层 (向量)
    self.a = np.empty(self.size, dtype=object) # 先占位(置空),dtype=object !如下皆然
    self.z = np.empty(self.size, dtype=object)
    # 偏置 (向量)
    self.b = np.empty(self.size, dtype=object)
    self.delta_b = np.empty(self.size, dtype=object)
    # 权 (矩阵)
    self.w = np.empty(self.size, dtype=object)
    self.delta_w = np.empty(self.size, dtype=object)
    # 残差 (向量)
    self.data_a = np.empty(self.size, dtype=object)
    # 填充
    for i in range(self.size):
      self.a[i] = np.zeros(self.n[i]) # 全零
      self.z[i] = np.zeros(self.n[i]) # 全零
      self.data_a[i] = np.zeros(self.n[i]) # 全零
      if i < self.size - 1:
        self.b[i] = np.ones(self.n[i+1])  # 全一
        self.delta_b[i] = np.zeros(self.n[i+1]) # 全零
        mu, sigma = 0, 0.1 # 均值、方差
        self.w[i] = np.random.normal(mu, sigma, (self.n[i], self.n[i+1])) # # 正态分布随机化
        self.delta_w[i] = np.zeros((self.n[i], self.n[i+1])) # 全零
    # 激活函数
    self.active_functions = {
      'sigmoid': self.sigmoid,
      'tanh': self.tanh,
      'radb': self.radb,
      'line': self.line,
    }
    # 激活函数的导函数
    self.derivative_functions = {
      'sigmoid': self.sigmoid_d,
      'tanh': self.tanh_d,
      'radb': self.radb_d,
      'line': self.line_d,
    }
    if active_type is None:
      self.active_type = ['sigmoid'] * (self.size - 1) # 默认激活函数类型
    else:
      self.active_type = active_type
  def sigmoid(self, z):
    if np.max(z) > 600:
      z[z.argmax()] = 600
    return 1.0 / (1.0 + np.exp(-z))
  def tanh(self, z):
    return (np.exp(z) - np.exp(-z)) / (np.exp(z) + np.exp(-z))
  def radb(self, z):
    return np.exp(-z * z)
  def line(self, z):
    return z
  def sigmoid_d(self, z):
    return z * (1.0 - z)
  def tanh_d(self, z):
    return 1.0 - z * z
  def radb_d(self, z):
    return -2.0 * z * np.exp(-z * z)
  def line_d(self, z):
    return np.ones(z.size) # 全一
  def forward(self, x):
    '''正向传播(在线)''' 
    # 用样本 x 走一遍,刷新所有 z, a
    self.a[0] = x
    for i in range(self.size - 1):
      self.z[i+1] = np.dot(self.a[i], self.w[i]) + self.b[i] 
      self.a[i+1] = self.active_functions[self.active_type[i]](self.z[i+1]) # 加了激活函数
  def err(self, X, Y):
    '''误差'''
    last = self.size-1
    err = 0.0
    for x, y in zip(X, Y):
      self.forward(x)
      err += 0.5 * np.sum((self.a[last] - y)**2)
    err /= X.shape[0]
    err += sum([np.sum(w) for w in self.w[:last]**2])
    return err
  def backward(self, y):
    '''反向传播(在线)'''
    last = self.size - 1
    # 用样本 y 走一遍,刷新所有delta_w, delta_b
    self.data_a[last] = -(y - self.a[last]) * self.derivative_functions[self.active_type[last-1]](self.z[last]) # 加了激活函数的导函数
    for i in range(last-1, 1, -1):
      self.data_a[i] = np.dot(self.w[i], self.data_a[i+1]) * self.derivative_functions[self.active_type[i-1]](self.z[i]) # 加了激活函数的导函数
      # 计算偏导
      p_w = np.outer(self.a[i], self.data_a[i+1]) # 外积!感谢 numpy 的强大!
      p_b = self.data_a[i+1]
      # 更新 delta_w, delta_w
      self.delta_w[i] = self.delta_w[i] + p_w
      self.delta_b[i] = self.delta_b[i] + p_b
  def update(self, n_samples):
    '''更新权重参数'''
    last = self.size - 1
    for i in range(last):
      self.w[i] -= self.alpha * ((1/n_samples) * self.delta_w[i] + self.lamda * self.w[i])
      self.b[i] -= self.alpha * ((1/n_samples) * self.delta_b[i])
  def fit(self, X, Y):
    '''拟合'''
    for i in range(self.n_iter):
      # 用所有样本,依次
      for x, y in zip(X, Y):
        self.forward(x) # 前向,更新 a, z;
        self.backward(y) # 后向,更新 delta_w, delta_b
      # 然后,更新 w, b
      self.update(len(X))
      # 计算误差
      err = self.err(X, Y)
      if err < self.error:
        break
      # 整千次显示误差(否则太无聊!)
      if i % 1000 == 0:
        print('iter: {}, error: {}'.format(i, err))
  def predict(self, X):
    '''预测'''
    last = self.size - 1
    res = []
    for x in X:
      self.forward(x)
      res.append(self.a[last])
    return np.array(res)
if __name__ == '__main__':
  nn = NeuralNetworks([2,3,4,3,1], n_iter=5000, alpha=0.4, lamda=0.3, error=0.06) # 定义神经网络
  X = np.array([[0.,0.], # 准备数据
         [0.,1.],
         [1.,0.],
         [1.,1.]])
  y = np.array([0,1,1,0])
  nn.fit(X,y)     # 拟合
  print(nn.predict(X)) # 预测

更多关于Python相关内容可查看本站专题:《Python数学运算技巧总结》、《Python数据结构与算法教程》、《Python函数使用技巧总结》、《Python字符串操作技巧汇总》及《Python入门与进阶经典教程》

希望本文所述对大家Python程序设计有所帮助。


# Python  # numpy  # 定义  # 神经网络结构  # Python利用numpy实现三层神经网络的示例代码  # Numpy实现卷积神经网络(CNN)的示例  # numpy实现神经网络反向传播算法的步骤  # 纯用NumPy实现神经网络的示例代码  # Python使用numpy实现BP神经网络  # numpy创建神经网络框架  # 正态分布  # 一遍  # 斯坦福  # 进阶  # 均值  # 相关内容  # 还可以  # 数据结构  # 给大家  # 随机化  # 更多关于  # 所述  # 程序设计  # 使用技巧  # 五层  # 前向  # 后向  # 操作技巧  # 敲出  # 迭代 


相关栏目: 【 网站优化151355 】 【 网络推广146373 】 【 网络技术251813 】 【 AI营销90571


相关推荐: Laravel Eloquent:优雅地将关联模型字段扁平化到主模型中  网站制作企业,网站的banner和导航栏是指什么?  Laravel如何使用Gate和Policy进行权限控制_Laravel权限判定与策略规则配置  成都品牌网站制作公司,成都营业执照年报网上怎么办理?  如何在 Telegram Web View(iOS)中防止键盘遮挡底部输入框  HTML透明颜色代码怎么让图片透明_给img元素加透明色的技巧【方法】  如何用AI帮你把自己的生活经历写成一个有趣的故事?  进行网站优化必须要坚持的四大原则  Laravel Vite是做什么的_Laravel前端资源打包工具Vite配置与使用  JavaScript如何实现继承_有哪些常用方法  微信小程序 HTTPS报错整理常见问题及解决方案  郑州企业网站制作公司,郑州招聘网站有哪些?  Laravel如何使用Service Provider注册服务_Laravel服务提供者配置与加载  如何挑选最适合建站的高性能VPS主机?  Laravel怎么返回JSON格式数据_Laravel API资源Response响应格式化【技巧】  详解jQuery中基本的动画方法  在Oracle关闭情况下如何修改spfile的参数  如何用已有域名快速搭建网站?  Windows家庭版如何开启组策略(gpedit.msc)?(安装方法)  再谈Python中的字符串与字符编码(推荐)  javascript事件捕获机制【深入分析IE和DOM中的事件模型】  Python文件异常处理策略_健壮性说明【指导】  Laravel怎么创建自己的包(Package)_Laravel扩展包开发入门到发布  深圳网站制作培训,深圳哪些招聘网站比较好?  如何续费美橙建站之星域名及服务?  香港网站服务器数量如何影响SEO优化效果?  Laravel如何创建和注册中间件_Laravel中间件编写与应用流程  香港代理服务器配置指南:高匿IP选择、跨境加速与SEO优化技巧  Laravel的HTTP客户端怎么用_Laravel HTTP Client发起API请求教程  百度输入法ai面板怎么关 百度输入法ai面板隐藏技巧  Laravel怎么防止CSRF攻击_Laravel CSRF保护中间件原理与实践  Laravel事件监听器怎么写_Laravel Event和Listener使用教程  Laravel Eloquent性能优化技巧_Laravel N+1查询问题解决  如何在橙子建站上传落地页?操作指南详解  详解jQuery停止动画——stop()方法的使用  Win11关机界面怎么改_Win11自定义关机画面设置【工具】  Laravel如何配置中间件Middleware_Laravel自定义中间件拦截请求与权限校验【步骤】  音乐网站服务器如何优化API响应速度?  制作电商网页,电商供应链怎么做?  企业网站制作这些问题要关注  Laravel如何创建自定义Artisan命令?(代码示例)  英语简历制作免费网站推荐,如何将简历翻译成英文?  Laravel如何操作JSON类型的数据库字段?(Eloquent示例)  ,南京靠谱的征婚网站?  家族网站制作贴纸教程视频,用豆子做粘帖画怎么制作?  如何快速搭建高效WAP手机网站?  用v-html解决Vue.js渲染中html标签不被解析的问题  php后缀怎么变mp4格式错误_修改扩展名提示格式不对怎么办【技巧】  如何获取上海专业网站定制建站电话?  Laravel中的Facade(门面)到底是什么原理