Python手写实现梯度下降算法(核心是求极值,必须是凸函数),使用matplotlib绘制过程图

1、什么是梯度下降:

梯度就是函数在一个点的斜率,梯度下降,就是在求函数极值的过程中,让斜率减小,从而找到极大值或极小值点,因为在几何坐标系里,极值点的斜率是零

# 梯度下降 == 导数值下降
import matplotlib.pyplot as plt
import numpy as np

# f(x) = (x-10)**2, 目标函数
# f'(x) = 2*x - 20  梯度函数:一阶导数函数

# 梯度下降算法是一个方法,是帮助我们找极值点的方法cost # 凸函数
def targetFunc(x):
    return (x-10)**2
    pass

def gradientFunc(x):
    return 2*x - 20
    pass

listx = []

# 猜测的过程
def gradientCal(initX, targetFunc, gradientFunc, rating=0.1, tolent=0.000000001, times = 500000):
    '''
    :param initX: 猜测的点
    :param targetFunc: 目标函数
    :param gradientFunc: 梯度函数
    :param rating: 步进系数
    :param tolent: 收敛条件
    :return: 返回极值点x值
    '''
    result = targetFunc(initX)           # 计算出initX这个点的实际值
    gradientResult = gradientFunc(initX) # 计算出initX这个点的导数,也是斜率,梯度
    listx.append(initX)
    newX = initX - rating*gradientFunc(initX)
    newResult = targetFunc(newX)

    reResult = np.abs(result - newResult)
    t = 0
    while reResult > tolent and t < times:
        t += 1
        initX = newX
        result = newResult
        listx.append(initX)
        newX = newX - rating * gradientFunc(newX)
        newResult = targetFunc(newX)
        reResult = np.abs(result - newResult)
        pass
    return initX

    pass

if __name__ == "__main__":
    print(gradientCal(-10, targetFunc, gradientFunc))
    x = np.arange(-10, 10+1, 1)
    y = (x - 10)**2

    plt.plot(x, y)
    plt.grid(linestyle='--')
    plt.scatter(np.array(listx), (np.array(listx) - 10)**2 , s=20)
    plt.show()

    print(listx)
    pass

2、收敛过程如下图所示:

  • 5
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值