【Python】【matLab】模拟退火算法求二元高次函数最小值

一、目标函数

求二元高次函数的最小值。目标函数选择:

用于测试算法的简单的目标函数:

二、Python代码实现

python 复制代码
import numpy as np


# 目标函数(2变量)
def objective_function(x):
    return x[0] ** 2 + 2 * x[0] - 15 + 4 * 4 * 2 * x[1] + 4 * x[1] ** 2
    # 测试:return x[0] ** 2 + x[1] ** 2


# 模拟退火
def simulated_annealing(objective_func,  # 目标函数
                        initial_solution=np.array([0, 0]),  # 初始解
                        temperature=100,  # 初始温度
                        min_temperature=0.1,  # 最小温度值
                        cooling_rate=0.90,  # 冷却率(乘法系数)
                        iter_max=100,  # 最大迭代次数
                        seed=0  # 随机种子
                        ):
    np.random.seed(seed)
    current_solution = initial_solution
    best_solution = current_solution

    while temperature > min_temperature:
        for j in range(iter_max):
            # 生成新的解
            new_solution = current_solution + np.random.uniform(-1, 1, len(current_solution))

            # 计算新解与当前解之间的目标函数值差异
            current_cost = objective_func(current_solution)
            new_cost = objective_func(new_solution)
            cost_diff = new_cost - current_cost

            # 判断是否接受新解
            if cost_diff < 0 or np.exp(-cost_diff / temperature) > np.random.random():
                current_solution = new_solution

            # 更新最优解
            if objective_func(current_solution) < objective_func(best_solution):
                best_solution = current_solution

        # 降低温度
        temperature *= cooling_rate

    return best_solution


# 调用退火算法求解最小值
best_solution = simulated_annealing(objective_function)

print(f"函数最小值: {objective_function(best_solution)} 自变量取值:{best_solution}")

代码参考博客:利用Python 实现 模拟退火算法

三、程序输出

测试函数输出:

目标函数输出:

四、MatLab验证程序

参考博客:MATLAB求解二元(多元)函数极值

先定义目标函数(位于fun2_3.m中):

matlab 复制代码
function f = fun2_3(x)
f = x(1) ^ 2 + 2 * x(1) - 15 + 32 * x(2) + 4 * x(2) ^ 2;

模拟退火算法求极值:

matlab 复制代码
clc, clear
[x, y]=meshgrid(-10:0.3:10,-10:0.3:10);
z = x.^2 + 2 * x -15 + 32 * y + 4 * y.^2;
figure(1)
surf(x,y,z)
xlabel('X');
ylabel('Y');
zlabel('Z');

figure(2)
contour(x,y,z)
xlabel('X');
ylabel('Y');
grid on;

x0=[-3,-3];
% [x,fmin]=fminsearch(@fun2_3,x0)
[x,fmin] = fminunc(@fun2_3,x0)

程序输出:

可见,两种方法的求解结果基本相同。

相关推荐
jieyu111915 分钟前
Python 实战:内网渗透中的信息收集自动化脚本(2)
python·网络安全·脚本开发
码界筑梦坊3 小时前
171-基于Flask的笔记本电脑数据可视化分析系统
python·信息可视化·flask·毕业设计·echarts
hui函数7 小时前
Flask电影投票系统全解析
后端·python·flask
闲人编程9 小时前
Python第三方库IPFS-API使用详解:构建去中心化应用的完整指南
开发语言·python·去中心化·内存·寻址·存储·ipfs
计算机编程小咖10 小时前
《基于大数据的农产品交易数据分析与可视化系统》选题不当,毕业答辩可能直接挂科
java·大数据·hadoop·python·数据挖掘·数据分析·spark
zhangfeng113311 小时前
以下是基于图论的归一化切割(Normalized Cut)图像分割工具的完整实现,结合Tkinter界面设计及Python代码示
开发语言·python·图论
flashlight_hi11 小时前
LeetCode 分类刷题:2529. 正整数和负整数的最大计数
python·算法·leetcode
Ashlee_code12 小时前
香港券商櫃台系統跨境金融研究
java·python·科技·金融·架构·系统架构·区块链
Jia-Hui Su12 小时前
GDSFactory环境配置(PyCharm+Git+KLayout)
git·python·pycharm
凌晨7点13 小时前
控制建模matlab练习14:线性状态反馈控制器-③极点配置
matlab