python实现最速下降法
作者：zcc_TPJH 发布时间：2023-08-10 18:19:51
标签：python,最速下降法
本文实例为大家分享了python实现最速下降法的具体代码，供大家参考，具体内容如下
代码：
from sympy import *
import numpy as np
def backtracking_line_search(f,df,x,x_k,p_k,alpha0):
rho=0.5
c=10**-4
alpha=alpha0
replacements1=zip(x,x_k)
replacements2=zip(x,x_k+alpha*p_k)
f_k=f.subs(replacements1)
df_p=np.dot([df_.subs(replacements1) for df_ in df],p_k)
while f.subs(replacements2)>f_k+c*alpha*df_p:
alpha=rho*alpha
replacements2 = zip(x, x_k +alpha * p_k)
return alpha
def stepest_line_search(f,x,x0,alpha0):
df = [diff(f, x_) for x_ in x]
x_k=x0
alpha=alpha0
replacements=zip(x,x_k)
len_df = sqrt(np.sum([df_.subs(replacements) ** 2 for df_ in df]))
while len_df>1e-6:
p_k=-1*np.array([df_.subs(replacements) for df_ in df])
alpha = backtracking_line_search(f, df, x, x_k, p_k, alpha)
x_k=x_k+alpha*p_k
replacements = zip(x, x_k)
len_df=np.sum([df_.subs(replacements)**2 for df_ in df])
return x_k
if __name__=="__main__":
init_printing(use_unicode=True)
x1 = symbols("x1")
x2 = symbols("x2")
x = np.array([x1, x2])
f = 100 * (x2 - x1 ** 2)**2 + (1 - x1) ** 2
ans=stepest_line_search(f, x, np.array([1.2, 1]), 1)
print "the minimal value in point:",ans
分析：
这个采用的是backtracking line search来寻找alpha。
来源：https://blog.csdn.net/weixin_39881922/article/details/80540799
0
投稿猜你喜欢
- 如何实现刷新当前页面呢？借助js你将无所不能。1，reload 方法，该方法强迫浏览器刷新当前页面。语法：location.reload([
- 1998年，W3C发布HTML 4.0 Specification，里面清清楚楚的写了每个标签的用法和语义。搜索引擎的算法参考了W3C的语义
- 代码如下:set fso=server.createobject("scripting.filesystemobject"
在用户体验这个行业，经常会听到，可用性，可访问性这样专业的名词，但是，事实上在很多产品实现过程里都忽略了这一点！WHY？举个很简单的例子，用- 如何用SysOjects来获知数据库的信息？SysObjects中就保存了数据库中所有对象的信息，如：SELECT * FROM SysOb
- SQL Server导出表到EXCEL文件的存储过程：*--数据导出EXCEL导出表中的数据到Excel,包含字段名,文件为真正的Excel
- PHP mysqli_select_db() 函数更改连接的默认数据库：删除数据库<?php // 假定数据库用户名：root，密码：
- HTTP头中一般断点下载时才用到Range和Content-Range实体头，Range用户请求头中，指定第一个字节的位置和最后一个字节的位
- SQL Server TEXT、NTEXT字段拆分的问题引用的内容：SET NOCOUNT ON CREATE
- 翻译：ShiningRay @ Nirvana Studio作者：Douglas Crockford来源：http://www.crockf
- 上期回顾:亚马逊购物用户体验分析 (一)“查找内部”功能书是在亚马逊最常被购买的产品之一，所以毋庸置疑亚马逊的开发小组已经建立了一个关于“查
- pre标签会原样保留HTML内容的格式，可是如果宽度过大会把页面撑坏，这时候需要自动换行来帮忙：Making preformate
- 字符串索引示意图字符串切片也就是截取字符串，取子串Python中字符串切片方法字符串[开始索引：结束索引：步长]切取字符串为开始索引到结束索
- 如何用组件实现自动发送电子邮件？我想做一个能够自动发送电子邮件的程序，该如何做？ 这就要用到w3 upl
- 客户用的数据库是mysql，而研发好的产品支持oracle，为了让客户掏腰包，我们必须把数据库环境从oracle转向mysql。我们在转换的
- 北京时间2月15日据国外媒体报道，美国知名sns网站Facebook全球活跃用户量已突破1.75亿大关。数据显示，全球20％的网民都使用Fa
- 具体代码和实现方法见下：第一个办法，这个程序可以进行万亿以下的货币金额转换（够用的了吧），其中汉字与数字均按一位计：Function&nbs
- 代码如下:---找出促销活动中销售额最高的职员 ---你刚在一家服装销售公司中找到了一份工作，此时经理要求你根据数据库中的两张表
- 可控制的滚动新闻不同于自动的滚动条，它是通过按钮控制移动的，当你把鼠标放在按钮上时，新闻内容就会向上或
- 使用SQL Server事件探查器工具，你可用一个捕获到的跟踪来收集有关服务器的重要信息。与索引优化向导(Index Tuning Wiza
