【最优化】scipy.optimize.fmin.
来源:互联网 发布:图书管理系统c语言链表 编辑:程序博客网 时间:2024/05/16 16:01
文章首发于本人的个人博客
这代码借鉴了《Python科学计算》,进行了改动
import scipy.optimize as optimport numpy as nppoints=[]def obj_func(p): x,y=p z=(1-x)**2+100*(y-x**2)**2 points.append((x,y,z)) return z#偏导数,有些优化方法用得到,有些用不到def fprime(p): x,y=p dx=-2+2*x-400*x*(y-x**2) dy=200*y-200*x**2 return np.array([dx,dy])init_point=(-2,-2)#这两种优化方法没用到偏导#result=opt.fmin(obj_func,init_point)#result=opt.fmin_powell(obj_func,init_point)#用到偏导的:#result=opt.fmin_cg(obj_func,init_point,fprime=fprime)#result=opt.fmin_bfgs(obj_func,init_point,fprime=fprime)result=opt.fmin_tnc(obj_func,init_point,fprime=fprime)#result=opt.fmin_l_bfgs_b(obj_func,init_point,fprime=fprime)#其它#result=opt.fmin_cobyla(obj_func,init_point,[])print(result)### 绘图import pylab as plp=np.array(points)xmin,xmax=np.min(p[:,0])-1,np.max(p[:,0])+1ymin,ymax=np.min(p[:,1]),np.max(p[:,1])+1Y,X=np.ogrid[ymin:ymax:500j,xmin:xmax:500j]Z=np.log10(obj_func((X,Y)))zmin,zmax=np.min(Z),np.max(Z)pl.imshow(Z,extent=(xmin,xmax,ymin,ymax),origin="bottom",aspect="auto")pl.plot(p[:,0],p[:,1])pl.scatter(p[:,0],p[:,1],c=range(len(p)))pl.xlim(xmin,xmax)pl.ylim(ymin,ymax)pl.show()
下面这2个优化方法的展示:
阅读全文
1 0
- 【最优化】scipy.optimize.fmin.
- Scipy教程 - 优化和拟合库scipy.optimize
- scipy.optimize.minimize 的优化算法(2): Powell
- 使用scipy.optimize进行目标函数中的矩阵参数优化
- 【解方程】scipy.optimize.solve.
- scipy.optimize.minimize 的优化算法(1): Nelder–Mead Simplex
- scipy.optimize.minimize 的优化算法(2): Powell中的line search算法
- scipy.optimize.minimize 的优化算法(3): DFP和BFGS的证明
- python scipy.optimize curve_fit 多高斯拟合
- Python multiprocessing & scipy optimize leastsq线程安全问题
- optimize 优化mysql数据表
- optimize table 优化表
- python科学计算之scipy——optimize
- MySQL中optimize优化表
- optimize table 表优化问题
- optimize table优化mysql例子
- Python导入Scipy子模块时出现问题(DLL load failed with scipy.optimize)
- 用 OPTIMIZE TABLE 优化 mysql 表
- 在IGBT的开启过程中,IGBT的电压降低,电流上升,在IGBT的关断过程中IGBT的电压上升,电流下降,在一段时间内,电压和电流均不为0,由于功率等于电压乘以电流,即P=U×I,因此将产生损耗,开
- 对人生的一点思考
- C++拷贝构造函数详解
- 微信小程序 View布局
- flex中的Repeater的用法,类似Java中的for循环
- 【最优化】scipy.optimize.fmin.
- ubuntu 16.04 安装opencv 2.4.13
- dcloud-mui 仿ios 向下滑动页面,回弹效果
- java常见面试题
- 浏览器加载和渲染html的顺序-css渲染效率的探究
- 关于mysql_query()的注意点
- iOS开发—调用第三方地图导航
- Python简单递归
- 全面理解Java内存模型(JMM)及volatile关键字