pytorch实现梯度下降算法例子

news2024/11/20 14:20:25

如题,利用pytorch,通过代码实现机器学习中的梯度下降算法,求解如下方程:
f ′ ( x , y ) = x 2 + 20 y 2 {f}'(x,y) = x^2 + 20 y^2 f(x,y)=x2+20y2
的最小值。

Latex语法参考:https://blog.csdn.net/ViatorSun/article/details/82826664


import numpy as np
import matplotlib.pyplot as plt
import matplotlib as mpl
import math
from mpl_toolkits.mplot3d import Axes3D
import warnings
import decimal


# 原函数
def Z(x,y):
    return x**2 + 20*(y**2)

# x方向上的梯度
def dx(x):
    return 2*x

# y方向上的梯度
def dy(y):
    return 40*y

# 初始值
X = x_0 = np.float64(3000.0)
Y = y_0 = np.float64(300.0)

# 学习率
alpha = 0.01

# 保存梯度下降所经过的点
globalX = [x_0]
globalY = [y_0]
globalZ = [Z(x_0,y_0)]

cnt = 0

# 迭代30while True:
    lastZ = Z(X,Y)
    temX = X - alpha * dx(X)
    temY = Y - alpha * dy(Y)
    temZ = Z(temX, temY)
    if(lastZ - temZ <= 0.00001):
        break;
    if(cnt >= 1000000):
        break;
    # X,Y 重新赋值
    X = temX
    Y = temY
    
    # 将新值存储起来
    globalX.append(temX)
    globalY.append(temY)
    globalZ.append(temZ)
    cnt = cnt +1
    
    
    
plt.figure(figsize=(8,3))
plt.subplot(111)

plt.ylabel('Y')
plt.xlabel('X')

plt.plot(globalZ)
plt.show()
# 打印结果
print(u"最终结果为:(x,y,z)=(%.5f, %.5f, %.5f),迭代次数:%d" % (X, Y, Z(X,Y),cnt))
print(u"迭代过程中取值")
num = len(globalX)
for i in range(num):
    print(u"第%d次迭代:%.5f, y%d=%.5f, z%d=%.5f" % (i,globalX[i],i,globalY[i],i,globalZ[i]))







输出结果:

在这里插入图片描述
最终结果为:(x,y,z)=(397.85867, 0.00000, 158291.51945),迭代次数:100
迭代过程中取值
第0次迭代:3000.00000, y0=300.00000, z0=10800000.00000
第1次迭代:2940.00000, y1=180.00000, z1=9291600.00000
第2次迭代:2881.20000, y2=108.00000, z2=8534593.44000
第3次迭代:2823.57600, y3=64.80000, z3=8056562.22778
第4次迭代:2767.10448, y4=38.88000, z4=7687100.29124
第5次迭代:2711.76239, y5=23.32800, z5=7364539.17367
第6次迭代:2657.52714, y6=13.99680, z6=7066368.72182
第7次迭代:2604.37660, y7=8.39808, z7=6784188.02823
第8次迭代:2552.28907, y8=5.03885, z8=6514687.28512
第9次迭代:2501.24329, y9=3.02331, z9=6256400.78564
第10次迭代:2451.21842, y10=1.81399, z10=6008537.55665
第11次迭代:2402.19405, y11=1.08839, z11=5770559.95657
第12次迭代:2354.15017, y12=0.65303, z12=5542031.55767
第13次迭代:2307.06717, y13=0.39182, z13=5322561.98712
第14次迭代:2260.92582, y14=0.23509, z14=5111786.68892
第15次迭代:2215.70731, y15=0.14106, z15=4909359.27237
第16次迭代:2171.39316, y16=0.08463, z16=4714948.40627
第17次迭代:2127.96530, y17=0.05078, z17=4528236.36337
第18次迭代:2085.40599, y18=0.03047, z18=4348918.17242
第19次迭代:2043.69787, y19=0.01828, z19=4176701.00164
第20次迭代:2002.82392, y20=0.01097, z20=4011303.63796
第21次迭代:1962.76744, y21=0.00658, z21=3852456.01246
第22次迭代:1923.51209, y22=0.00395, z22=3699898.75384
第23次迭代:1885.04185, y23=0.00237, z23=3553382.76300
第24次迭代:1847.34101, y24=0.00142, z24=3412668.80552
第25次迭代:1810.39419, y25=0.00085, z25=3277527.12080
第26次迭代:1774.18631, y26=0.00051, z26=3147737.04681
第27次迭代:1738.70258, y27=0.00031, z27=3023086.65975
第28次迭代:1703.92853, y28=0.00018, z28=2903372.42802
第29次迭代:1669.84996, y29=0.00011, z29=2788398.87987
第30次迭代:1636.45296, y30=0.00007, z30=2677978.28423
第31次迭代:1603.72390, y31=0.00004, z31=2571930.34417
第32次迭代:1571.64942, y32=0.00002, z32=2470081.90254
第33次迭代:1540.21643, y33=0.00001, z33=2372266.65920
第34次迭代:1509.41210, y34=0.00001, z34=2278324.89950
第35次迭代:1479.22386, y35=0.00001, z35=2188103.23348
第36次迭代:1449.63938, y36=0.00000, z36=2101454.34543
第37次迭代:1420.64660, y37=0.00000, z37=2018236.75335
第38次迭代:1392.23366, y38=0.00000, z38=1938314.57792
第39次迭代:1364.38899, y39=0.00000, z39=1861557.32064
第40次迭代:1337.10121, y40=0.00000, z40=1787839.65074
第41次迭代:1310.35919, y41=0.00000, z41=1717041.20057
第42次迭代:1284.15200, y42=0.00000, z42=1649046.36903
第43次迭代:1258.46896, y43=0.00000, z43=1583744.13281
第44次迭代:1233.29958, y44=0.00000, z44=1521027.86515
第45次迭代:1208.63359, y45=0.00000, z45=1460795.16169
第46次迭代:1184.46092, y46=0.00000, z46=1402947.67329
第47次迭代:1160.77170, y47=0.00000, z47=1347390.94543
第48次迭代:1137.55627, y48=0.00000, z48=1294034.26399
第49次迭代:1114.80514, y49=0.00000, z49=1242790.50714
第50次迭代:1092.50904, y50=0.00000, z50=1193576.00305
第51次迭代:1070.65886, y51=0.00000, z51=1146310.39333
第52次迭代:1049.24568, y52=0.00000, z52=1100916.50176
第53次迭代:1028.26077, y53=0.00000, z53=1057320.20829
第54次迭代:1007.69555, y54=0.00000, z54=1015450.32804
第55次迭代:987.54164, y55=0.00000, z55=975238.49505
第56次迭代:967.79081, y56=0.00000, z56=936619.05064
第57次迭代:948.43499, y57=0.00000, z57=899528.93624
第58次迭代:929.46629, y58=0.00000, z58=863907.59036
第59次迭代:910.87697, y59=0.00000, z59=829696.84979
第60次迭代:892.65943, y60=0.00000, z60=796840.85453
第61次迭代:874.80624, y61=0.00000, z61=765285.95669
第62次迭代:857.31011, y62=0.00000, z62=734980.63281
第63次迭代:840.16391, y63=0.00000, z63=705875.39975
第64次迭代:823.36063, y64=0.00000, z64=677922.73392
第65次迭代:806.89342, y65=0.00000, z65=651076.99366
第66次迭代:790.75555, y66=0.00000, z66=625294.34471
第67次迭代:774.94044, y67=0.00000, z67=600532.68866
第68次迭代:759.44163, y68=0.00000, z68=576751.59419
第69次迭代:744.25280, y69=0.00000, z69=553912.23106
第70次迭代:729.36774, y70=0.00000, z70=531977.30671
第71次迭代:714.78039, y71=0.00000, z71=510911.00536
第72次迭代:700.48478, y72=0.00000, z72=490678.92955
第73次迭代:686.47509, y73=0.00000, z73=471248.04394
第74次迭代:672.74558, y74=0.00000, z74=452586.62140
第75次迭代:659.29067, y75=0.00000, z75=434664.19119
第76次迭代:646.10486, y76=0.00000, z76=417451.48922
第77次迭代:633.18276, y77=0.00000, z77=400920.41025
第78次迭代:620.51911, y78=0.00000, z78=385043.96200
第79次迭代:608.10872, y79=0.00000, z79=369796.22111
第80次迭代:595.94655, y80=0.00000, z80=355152.29075
第81次迭代:584.02762, y81=0.00000, z81=341088.26004
第82次迭代:572.34707, y82=0.00000, z82=327581.16494
第83次迭代:560.90013, y83=0.00000, z83=314608.95081
第84次迭代:549.68212, y84=0.00000, z84=302150.43636
第85次迭代:538.68848, y85=0.00000, z85=290185.27908
第86次迭代:527.91471, y86=0.00000, z86=278693.94202
第87次迭代:517.35642, y87=0.00000, z87=267657.66192
第88次迭代:507.00929, y88=0.00000, z88=257058.41851
第89次迭代:496.86910, y89=0.00000, z89=246878.90514
第90次迭代:486.93172, y90=0.00000, z90=237102.50049
第91次迭代:477.19309, y91=0.00000, z91=227713.24147
第92次迭代:467.64922, y92=0.00000, z92=218695.79711
第93次迭代:458.29624, y93=0.00000, z93=210035.44354
第94次迭代:449.13032, y94=0.00000, z94=201718.03998
第95次迭代:440.14771, y95=0.00000, z95=193730.00560
第96次迭代:431.34475, y96=0.00000, z96=186058.29738
第97次迭代:422.71786, y97=0.00000, z97=178690.38880
第98次迭代:414.26350, y98=0.00000, z98=171614.24940
第99次迭代:405.97823, y99=0.00000, z99=164818.32513
第100次迭代:397.85867, y100=0.00000, z100=158291.51945

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/790665.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

【雕爷学编程】Arduino动手做(87)---ULN2003步进电机模组2

37款传感器与执行器的提法&#xff0c;在网络上广泛流传&#xff0c;其实Arduino能够兼容的传感器模块肯定是不止这37种的。鉴于本人手头积累了一些传感器和执行器模块&#xff0c;依照实践出真知&#xff08;一定要动手做&#xff09;的理念&#xff0c;以学习和交流为目的&am…

openAI API使用教程、openai.Completion.create() 详细解释一下

openAI提供了几种不同场景的模型,主要有text completion、code completion、chat completion、image completion,例如chat completion,则调用方式为。而且请求的token和回复的token数会被加一起计费,例如说输入了10个token,openAI回复了20个token,那么最终收费是按照30个…

微信登录账户文件、聊天信息存储转移(存储空间足够、想保留历史聊天记录、文件)

1、打开电脑版微信、点击左下角的三根横线 2、点击左侧的“设置” 3、弹出层左侧点击“文件管理” 4、记录历史存储路径验证使用 历史默认存储位置“./我的文档/WeChat Files” 5、在其他存储空间充足盘符创建存储路径&#xff08;也可在 步骤7过程创建&#xff09; 具体位…

SpringBoot开发小技巧使用(DEBUG、启动图标修改、Lombok、devtools、Spring Initializr)

目录 1. 通过DEBUG查看自动配置的组件2. springboot启动图标修改3. Lombok4. devtools5. 通过IDEA的Spring Initializr快速创建新项目 1. 通过DEBUG查看自动配置的组件 在resources/application.properties中添加如下&#xff0c;开启DEBUG功能&#xff1a; debugtrue然后启动…

.Net Core 6.0依赖注入

.Net Core 6.0依赖注入 往期文章&#xff1a; .ner Core实现接口限流.net Core程序发布到IIS(Window Server 2019) 文章目录 .Net Core 6.0依赖注入前言一、ICO 和DI和DLICO [控制反转]DI [依赖注入]DL [依赖查找] 二、.net Core 中的依赖注入【Autofac】瞬时模式作用域模式单…

中文输入法开发-关键代码

续上篇介绍了嵌入式Linux下中文输入法&#xff0c; 嵌入式Linux下开发中文输入法_嵌入式输入法_小刚学長的博客-CSDN博客 本篇继续介绍核心关键功能 展现效果图如下&#xff1a; 1、如何跟应用关联起来&#xff0c;比如说&#xff0c;希望当LineEdit 输入状态激活后&#xff0…

怎么清空回收站?3个方法轻松搞定!

有没有大佬知道该怎么清空回收站呀&#xff1f;想把回收站清空了不知道该如何操作&#xff01;求一个清空方法&#xff01;感谢大家啦&#xff01; 电脑的回收站会为我们保存一些删除的数据&#xff0c;如果我们删除文件后意识到这些文件仍然是有用的&#xff0c;可以利用回收站…

报考CDGP一定要先有CDGA证书吗?

据DAMA中国官方网站消息&#xff0c;2023年度第三期DAMA中国CDGA和CDGP认证考试定于2023年9月23日举行&#xff0c;报名通道现已开启。 最近好多想报名CDGP认证考试的朋友都来咨询考试报名事宜&#xff0c;其中有一个重合度很高的问题“报考CDGP一定要先有CDGA证书吗&#xff…

【leetcode难题】2569. 更新数组后处理求和查询【线段树实现01翻转和区间求和模版】

题目截图 题目分析 关键就是记录每次操作2时&#xff0c;nums1中的1的个数这就需要实现线段树进行区间反转以及区间求和 ac code class Solution:def handleQuery(self, nums1: List[int], nums2: List[int], queries: List[List[int]]) -> List[int]:n len(nums1)m le…

校企合作谱新篇︱西藏农牧学院水土学院与鹏业软件签署校企合作协议

2023年7月17日&#xff0c;西藏农牧学院水利土木工程学院&#xff08;以下简称“水土学院”与成都鹏业软件股份有限公司&#xff08;以下简称“鹏业软件”&#xff09;在西藏林芝签署校企合作协议并揭牌。水土学院副书记、院长王培清、副院长孙海波、土木教研室主任宗永臣、土木…

【SpringBoot应用篇】SpringBoot+MybatisPlus集成国产DM8(达梦)数据库

【SpringBoot应用篇】SpringBootMybatisPlus集成国产DM8&#xff08;达梦&#xff09;数据库 简介和安装基本概念介绍SpringBootMP整合DM8pomymlAddressAddressMapper启动类测试类增删改查 yml配置与实体类TableName表映射问题 简介和安装 DM8 是达梦数据库有限公司推出的新一…

AI加速游戏开发 亚马逊云科技适配3大场景,打造下一代游戏体验

随着疫情的消散&#xff0c;中国游戏产业正在快速前进。在伴随着游戏产业升级的同时&#xff0c;整个行业都在面临着新的挑战与新的诉求。亚马逊云科技游戏研发解决方案和服务&#xff0c;覆盖端到端3大场景&#xff0c;为游戏公司与游戏开发人员赋能。 场景1&#xff1a;AI辅助…

memcached最大的优势是什么?

在Java中使用Memcached的最大优势之一是它提供了高性能的分布式缓存解决方案。Memcached是一个简单、快速、开源的内存缓存系统&#xff0c;可以帮助加速应用程序的访问速度&#xff0c;减轻数据库的负担&#xff0c;提高应用的扩展性和吞吐量。它常用于缓存频繁读取的数据&…

仿微信开发

介绍 一款仿微信界面的即时通讯APP, 采用uniapp开发(uview1.x版本),后端采用java.(springbootwebsocket) ####功能支持 1:添加好友,单聊,群聊群组的发起 2:支持文字,表情,语音,图片,文件信息的收发,离线消息推送仅限安卓(ios没证书,安卓需要自行申请unipush) 3:支持发送红包…

高压开关柜无线测温系统设计与实现 安科瑞 许敏

摘要&#xff1a;文章提出了一种用于小车式开关柜的无线测温系统&#xff0c;包括设于小车动触头套管上的温度传感器、环绕设于小车动触头上的电流传感器&#xff0c;温度传感器及电流传感器均分别连接有数据接收处理器&#xff0c;数据接收处理器电连接有设于小车式开关柜上的…

04 - 慎重使用正则表达式

在讲 String 对象优化时&#xff0c;提到了 Split() 方法&#xff0c;该方法使用的正则表达式可能引起回溯问题&#xff0c;今天就来深入了解下&#xff0c;这究竟是怎么回事&#xff1f; 开始之前&#xff0c;我们先来看一个案例&#xff0c;可以帮助你更好地理解内容。 在一…

uniapp兼容微信小程序和支付宝小程序遇到的坑

1、支付宝不支持v-show 改为v-if。 2、v-html App端和H5端支持 v-html &#xff0c;微信小程序会被转为 rich-text&#xff0c;其他端不支持 v-html。 解决方法&#xff1a;去插件市场找一个支持跨端的富文本组件。 3、导航栏处有背景色延伸至导航栏外 兼容微信小程序和支…

【装饰器设计模式详解】C/Java/JS/Go/Python/TS不同语言实现

简介 装饰器模式&#xff08;Decorator Pattern&#xff09;是一种结构型设计模式。将对象放入到一个特殊封装的对象中&#xff0c;为这个对象绑定新的行为&#xff0c;具备新的能力&#xff0c;同时又不改变其原有结构。 如果你希望在无需修改代码的情况下即可使用对象&…

《水经注地图服务》发布的卫星影像数据在OpenLayers中调用

OpenLayers是一个用于开发WebGIS客户端的JavaScript包。 OpenLayers 支持的地图来源包括Google Maps、Yahoo、 Map、微软Virtual Earth 等&#xff0c;用户还可以用简单的图片地图作为背景图&#xff0c;与其他的图层在OpenLayers 中进行叠加&#xff0c;在这一方面OpenLayers…

Linux6.15 Docker 私有仓库(harbor)

文章目录 计算机系统5G云计算第四章 LINUX Docker 私有仓库&#xff08;harbor&#xff09;一、搭建本地私有仓库二、Docker--harbor私有仓库部署与管理1.Harbor 简介1&#xff09;什么是Harbor2&#xff09;Harbor的特性3&#xff09;Harbor的构成 2.Harbor 部署1&#xff09;…