【Python】最速上升法求函数最大值

最新推荐文章于 2022-02-10 23:49:19 发布

想要变强的瓜瓜

最新推荐文章于 2022-02-10 23:49:19 发布

阅读量1.8k

点赞数 2

分类专栏： Python #微积分文章标签：回归几何学 python

本文链接：https://blog.csdn.net/qq_50143262/article/details/121641110

版权

简介

运用回归设计和响应曲面分析的方法来解决实验工作者通过实验寻找最佳的工艺条件，在实验者所选定的一个小区域上运用回归设计拟合线性回归方程，用最速上升法向最佳区域逼近因子区域，通常，将此法重复使用若干次，就可达到最佳区域。

方法步骤

利用最速上升法寻找最佳区域的步骤为：
(1)在变量的某区域内，通过正交设计，拟合一个线性回归方程；
(2)由拟合的线性回归方程，找出最速上升路线；
(3)沿着这条上升路线进行一系列实验，直到y的值不再明显增大为止；
(4)在(3)中y值不再明显增大的那一点的邻域内重复(1)、(2)、(3)；
(5)当拟合的线性回归方程不再显著，y的曲面具有明显的弯曲时，拟合二次回归，进
而找出最佳条件。

题目

在这里插入图片描述

代码实现

# Method of Steepest Ascent 範例7.5, page 690


import math as m

def f(x,y):                                #定義函數f(x,y)
    return 4*x*y-m.pow(x,4)-m.pow(y,4)+4


def Gf(x,y):                               #定義f(x,y)的一次偏導數
    x1=4*y-4*m.pow(x,3)
    x2=4