本文实现的原理很简单，优化方法是用的梯度下降。后面有测试结果。

先来看看实现的示例代码：

# coding=utf-8
from math import exp

import matplotlib.pyplot as plt
import numpy as np
from sklearn.datasets.samples_generator import make_blobs


def sigmoid(num):
 '''

 :param num: 待计算的x
 :return: sigmoid之后的数值
 '''
 if type(num) == int or type(num) == float:
  return 1.0 / (1 + exp(-1 * num))
 else:
  raise ValueError, 'only int or float data can compute sigmoid'


class logistic():
 def __init__(self, x, y): 
  if type(x) == type(y) == list:
   self.x = np.array(x)
   self.y = np.array(y)
  elif type(x) == type(y) == np.ndarray:
   self.x = x
   self.y = y
  else:
   raise ValueError, 'input data error'

 def sigmoid(self, x):
  '''

  :param x: 输入向量
  :return: 对输入向量整体进行simgoid计算后的向量结果
  '''
  s = np.frompyfunc(lambda x: sigmoid(x), 1, 1)
  return s(x)

 def train_with_punish(self, alpha, errors, punish=0.0001):
  '''

  :param alpha: alpha为学习速率
  :param errors: 误差小于多少时停止迭代的阈值
  :param punish: 惩罚系数
  :param times: 最大迭代次数
  :return:
  '''
  self.punish = punish
  dimension = self.x.shape[1]
  self.theta = np.random.random(dimension)
  compute_error = 100000000
  times = 0
  while compute_error > errors:
   res = np.dot(self.x, self.theta)
   delta = self.sigmoid(res) - self.y
   self.theta = self.theta - alpha * np.dot(self.x.T, delta) - punish * self.theta # 带惩罚的梯度下降方法
   compute_error = np.sum(delta)
   times += 1

 def predict(self, x):
  '''

  :param x: 给入新的未标注的向量
  :return: 按照计算出的参数返回判定的类别
  '''
  x = np.array(x)
  if self.sigmoid(np.dot(x, self.theta)) > 0.5:
   return 1
  else:
   return 0


def test1():
 '''
 用来进行测试和画图，展现效果
 :return:
 '''
 x, y = make_blobs(n_samples=200, centers=2, n_features=2, random_state=0, center_box=(10, 20))
 x1 = []
 y1 = []
 x2 = []
 y2 = []
 for i in range(len(y)):
  if y[i] == 0:
   x1.append(x[i][0])
   y1.append(x[i][1])
  elif y[i] == 1:
   x2.append(x[i][0])
   y2.append(x[i][1])
 # 以上均为处理数据，生成出两类数据
 p = logistic(x, y)
 p.train_with_punish(alpha=0.00001, errors=0.005, punish=0.01) # 步长是0.00001，最大允许误差是0.005，惩罚系数是0.01
 x_test = np.arange(10, 20, 0.01)
 y_test = (-1 * p.theta[0] / p.theta[1]) * x_test
 plt.plot(x_test, y_test, c='g', label='logistic_line')
 plt.scatter(x1, y1, c='r', label='positive')
 plt.scatter(x2, y2, c='b', label='negative')
 plt.legend(loc=2)
 plt.title('punish value = ' + p.punish.__str__())
 plt.show()


if __name__ == '__main__':
 test1()

运行结果如下图

总结

以上就是这篇文章的全部内容了，希望本文的内容对大家的学习或者工作能带来一定的帮助，如果有疑问大家可以留言交流，谢谢大家对的支持。

# python实现逻辑回归  # 逻辑回归  # python  # 实例  # pytorch使用nn.Moudle实现逻辑回归  # 基于Pytorch实现逻辑回归  # PyTorch实现多维度特征输入逻辑回归  # pytorch实现逻辑回归  # PyTorch线性回归和逻辑回归实战示例  # python sklearn库实现简单逻辑回归的实例代码  # python编写Logistic逻辑回归  # python代码实现逻辑回归logistic原理  # Pytorch实现逻辑回归分类  # 迭代  # 均为  # 很简单  # 这篇文章  # 谢谢大家  # 两类  # 先来  # 如下图  # 计算出  # 进行测试  # 有疑问  # ValueError  # return  # data  # compute  # param  # float  # int  # type  # raise 

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

相关推荐： 长沙做网站要多少钱,长沙国安网络怎么样？  悟空识字怎么关闭自动续费_悟空识字取消会员自动扣费步骤  Laravel如何使用Gate和Policy进行权限控制_Laravel权限判定与策略规则配置  Laravel如何集成Inertia.js与Vue/React？（安装配置）  Laravel如何使用Eloquent进行子查询  头像制作网站在线观看,除了站酷，还有哪些比较好的设计网站？  edge浏览器无法安装扩展 edge浏览器插件安装失败【解决方法】  作用域操作符会触发自动加载吗_php类自动加载机制与::调用【教程】  Laravel如何清理系统缓存命令_Laravel清除路由配置及视图缓存的方法【总结】  Laravel Sail是什么_基于Docker的Laravel本地开发环境Sail入门  企业在线网站设计制作流程,想建设一个属于自己的企业网站，该如何去做？  儿童网站界面设计图片,中国少年儿童教育网站-怎么去注册？  jimdo怎样用html5做选项卡_jimdo选项卡html5实现与切换效果【指南】  如何快速启动建站代理加盟业务？  大同网页,大同瑞慈医院官网？  Laravel Artisan命令怎么自定义_创建自己的Laravel命令行工具完全指南  Laravel如何处理CORS跨域请求？（配置示例）  Laravel怎么配置S3云存储驱动_Laravel集成阿里云OSS或AWS S3存储桶【教程】  Laravel Octane如何提升性能_使用Laravel Octane加速你的应用  如何生成腾讯云建站专用兑换码？  javascript事件捕获机制【深入分析IE和DOM中的事件模型】  网站制作免费,什么网站能看正片电影？  如何用PHP快速搭建高效网站？分步指南  Android利用动画实现背景逐渐变暗  LinuxCD持续部署教程_自动发布与回滚机制  历史网站制作软件,华为如何找回被删除的网站？  深圳网站制作培训,深圳哪些招聘网站比较好？  在线制作视频网站免费,都有哪些好的动漫网站？  图册素材网站设计制作软件,图册的导出方式有几种？  如何在Windows虚拟主机上快速搭建网站？  制作企业网站建设方案,怎样建设一个公司网站？  Laravel如何使用Scope本地作用域_Laravel模型常用查询逻辑封装技巧【手册】  Laravel如何与Docker（Sail）协同开发？（环境搭建教程）  Laravel定时任务怎么设置_Laravel Crontab调度器配置  Java Adapter 适配器模式（类适配器，对象适配器）优缺点对比  JS碰撞运动实现方法详解  手机怎么制作网站教程步骤,手机怎么做自己的网页链接？  Laravel如何设置定时任务（Cron Job）_Laravel调度器与任务计划配置  如何快速生成橙子建站落地页链接？  DeepSeek是免费使用的吗 DeepSeek收费模式与Pro版本功能详解  如何在万网开始建站？分步指南解析  canvas 画布在主流浏览器中的尺寸限制详细介绍  米侠浏览器网页背景异常怎么办 米侠显示修复  如何快速生成ASP一键建站模板并优化安全性？  Claude怎样写约束型提示词_Claude约束提示词写法【教程】  详解Nginx + Tomcat 反向代理 如何在高效的在一台服务器部署多个站点  如何在IIS7上新建站点并设置安全权限？  如何在香港服务器上快速搭建免备案网站？  用yum安装MySQLdb模块的步骤方法  LinuxShell函数封装方法_脚本复用设计思路【教程】