c++分治法求最大最小值实现_最优化计算与matlab实现(13)——非线性最小二乘优化问题——修正G-N法...

6fa308af866cf6f05b577b451e0dab2d.png

参考资料

《精通MATLAB最优化计算(第二版)》

编程工具

Matlab 2019a

目录

石中居士:最优化计算与Matlab实现——目录​zhuanlan.zhihu.com
8e2c6c11a2adc1402c25218a3d0c00a5.png

修正G-N法

  • 原理

为了克服G-N的缺点,有以下两种修正方案:

(1)在已经得到一个近似极小点

后,计算
,其中
的大小用一维无约束优化方法求解
确定,由于每一步都要进行一维搜索,计算量很大;

(2)由于

为目标函数的下降方向,选取很小的正数
,使得
此种方案计算量比较小,可以用比较简单的方法来确定
的取值。

因此,根据第二种方案可编制下面的修正G-N算法。

  • 算法步骤

用修正G-N法求解非线性最小二乘优化问题

的算法过程如下:

【1】给定初始点

,及精度
,置

【2】计算

【3】计算

【4】计算

【5】解方程

【6】令

【7】

,令
,转【8】;否则令
,转【7】;

【8】终止条件是否满足,不满足则令

,转【2】。
  • Matlab代码与试算

修正G-N法求解非线性最小二乘优化问题实例。用修正G-N法求下面的优化问题:

,其中初始点取

test.m

syms 

Modified_Gauss_Newton_Method.m

function

命令行窗口

x_optimization =

   -0.6941


f_optimization =

    0.7140

图像

6f6d783ef633f84785cc63ed7d4d70ca.png

由于目标函数是单谷函数,只有一个极值点,因此修正G-N法求出的极小点也是目标函数的最小点。


b5296ea20864dbc8ea0fe4fab63126aa.png

f2f000b01fc22aab68cc2bbb5ea8d647.png

目录

石中居士:最优化计算与Matlab实现——目录​zhuanlan.zhihu.com
60cafd8378734a7f4e038ab1df4a38c5.png

分治法是一种递归的思想,将问题分成若干个子问题,解决子问题后合并得到原问题的解。对于最大最小值问题,也可以使用分治法来解决。 具体思路如下: 1. 将序列分成两个子序列。 2. 对两个子序列分别递归最大最小值。 3. 将两个子序列的最大最小值合并,得到原序列的最大最小值代码实现如下: ``` #include <stdio.h> #include <stdlib.h> struct result { int max; // 最大int min; // 最小值 }; struct result find_max_min(int* arr, int left, int right) { struct result res, left_res, right_res; int mid; if (left == right) { // 数组只有一个元素 res.max = arr[left]; res.min = arr[left]; return res; } if (left + 1 == right) { // 数组有两个元素 if (arr[left] > arr[right]) { res.max = arr[left]; res.min = arr[right]; } else { res.max = arr[right]; res.min = arr[left]; } return res; } // 数组有三个以上的元素 mid = left + (right - left) / 2; left_res = find_max_min(arr, left, mid); right_res = find_max_min(arr, mid + 1, right); // 合并两个子序列的最大最小值 if (left_res.max > right_res.max) { res.max = left_res.max; } else { res.max = right_res.max; } if (left_res.min < right_res.min) { res.min = left_res.min; } else { res.min = right_res.min; } return res; } int main() { int arr[] = { 3, 7, 1, 4, 2, 8, 5, 9, 6 }; int len = sizeof(arr) / sizeof(arr[0]); struct result res = find_max_min(arr, 0, len - 1); printf("Max: %d, Min: %d\n", res.max, res.min); return 0; } ``` 运行结果: ``` Max: 9, Min: 1 ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值