本文章包含以下内容

1、画出最速下降法的算法流程图；

2、MATLAB 编写用数值微分法的梯度计算函数（函数式 M 文件）；

3、MATLAB 编写最速下降法求解无约束优化问题的函数，要求采用黄金分割法精确一维搜索，用数值微分法计算梯度（函数式 M 文件，精度设为 epson 可调）；

4、MATLAB 编写最速下降法求解无约束优化问题的函数，要求采用 Wolfe-Powell 非精确一维搜索，用数值微分法计算梯度（函数式 M 文件，精度设为 epson 可调）；

5、MATLAB 编写程序（命令式 M 文件），分别利用精确搜索和不精确搜索的最速下降法，求解如下问题：

$\min f(x)=100(x_2-x_1^2)^2+(1-x_1)^2$

精度为 0.001，初始点为（-1，1）;

改变初始点为（-1.2，1）重新运行，观察运行结果。

一维搜索，进退法，黄金分割法代码链接

本实验中函数用单独function计算

function y=f(x)
if(length(x)==1)global xk;global pk;x=xk+x*pk;
end
y=100*(x(2)-x(1)^2)^2+(1-x(1))^2;

1.最速下降法的算法流程图

2、MATLAB 编写用数值微分法的梯度计算函数（函数式 M 文件）；

function g=shuzhiweifenfa(x)
for i = 1:length(x)m=zeros(1,length(x));m(i)=(10^-6)/2;g(i)=f(x+m)-f(x-m);
end
g=g/10^-6;

3、最速下降法求解无约束优化问题的函数，采用黄金分割法精确一维搜索，用数值微分法计算梯度（函数式 M 文件，精度设为 epson 可调）；

function x=zuisuxiajiangfa_hjfg(e,x)
%step 1
%没用到k，只存储当前迭代的值。
global xk;
global pk;
while 1%step 2g=shuzhiweifenfa(x);%step 3%范数用的是平方和开根号if sqrt(sum(g.^2))<=ereturn;endpk=-g;xk=x;%这两个函数见之前代码（matlab无约束最优化的一般算法）[a,b,c]=jintuifa(0,0.1);a=huangjinfenge(a,c,10^-4);%step 4x=x+a*pk;
end

4、最速下降法求解无约束优化问题的函数，要求采用 Wolfe-Powell 非精确一维搜索，用数值微分法计算梯度（函数式 M 文件，精度设为 epson 可调）；

function a=Wolfe_Powell(x,pk)
%step 1
u=0.1;
b=0.5;
a=1;
n=0;
m=10^100;
%step 2
fx=f(x);
g=shuzhiweifenfa(x);
while 1xk=x+a*pk;fxk=f(xk);gk=shuzhiweifenfa(xk);if (fx-fxk)>=(-u*a*g*pk.')%(3-1)if (gk*pk.')>=(b*g*pk.')%(3-2)return;else%step 4n=a;a=min(2*a,(a+m)/2);endelse%step 3m=a;a=(a+n)/2;end
end

function x=zuisuxiajiangfa_Wolfe(e,x)
%step 1
%没用到k，只存储当前迭代的值。
while 1%step 2g=shuzhiweifenfa(x);%step 3%范数用的是平方和开根号if sqrt(sum(g.^2))<=ereturn;endpk=-g;a=Wolfe_Powell(x,pk);%step 4x=x+a*pk;
end

5、分别利用精确搜索和不精确搜索的最速下降法，问题：

$\min f(x)=100(x_2-x_1^2)^2+(1-x_1)^2$

clear
clc
for i=1:2if(i==1)x=[-1,1];fprintf('=========================');fprintf('\nx=%f\t\t%f\n',x(1),x(2));fprintf('=========================\n');elsex=[-1.2,1];fprintf('=========================');fprintf('\nx=%f\t\t%f\n',x(1),x(2));fprintf('=========================\n');endfprintf('精确搜索的最速下降法:\n');x_=zuisuxiajiangfa_hjfg(10^-3,x);fprintf('x*=%f\t%f\n',x_(1),x_(2));fprintf('f(x)=%f\n',f(x_));fprintf('不精确搜索的最速下降法\n');x_=zuisuxiajiangfa_Wolfe(10^-3,x);fprintf('x*=%f\t%f\n',x_(1),x_(2));fprintf('f(x)=%f\n',f(x_));
end

结果：

利用 MATLAB 编程实现最速下降法求解无约束最优化问题相关推荐

利用 MATLAB 编程实现共轭梯度法求解无约束最优化问题
本文章包含以下内容 1.FR 共轭梯度法的算法流程图: 2.MATLAB 编写 FR 共轭梯度法求解无约束优化问题的函数,要求采用黄金分割法精确一维搜索,用数值微分法计算梯度(函数式 M 文件,精度 ...
机器学习之求解无约束最优化问题方法(手推公式版)
文章目录前言 1. 基础知识 1.1 方向导数 1.2 梯度 1.3 方向导数与梯度的关系 1.4 泰勒展开公式 1.5 Jacobian矩阵与Hessian矩阵 1.6 正定矩阵 2. 梯度下降法 ...
matlab粒子群算法求解无约束最小值,pso matlab粒子群算法和遗传是解决约束优化问题，无和多目标的优 259万源代码下载- www.pudn.com...
文件名称: pso下载收藏√ [ 5 4 3 2 1 ] 开发工具: matlab 文件大小: 51 KB 上传时间: 2016-06-01 下载次数: 0 提供者: 孙志勇详细说 ...
牛顿法求解无约束最优化问题
记泰勒公式舍去高阶项在导数为0时候取得极值,令,以上式子可以变形为继续整理可以得到通过观察上式,我们可以发现,当给出一个初始点a的时候,我们可以通过泰勒公式用a点的一阶导数与二阶导数还有a的值 ...
计算智能课程设计(遗传算法求解无约束单目标优化问题)
写在前面前天写完了基于传递闭包的模糊聚类,今天准备写"遗传算法求解无约束单目标优化问题".昨天和npy玩了一下午,去齐白石艺术学院看了画展,一起在最高处看了夕阳,并在落日前接吻. ...
利用MATLAB编程实现系统传递函数的构建以及它们之间进行串联、并联、反馈时的构建方法
本文主要介绍如何利用MATLAB编程实现系统传递函数的构建以及它们之间进行串联.并联.反馈时的构建方法一.传递函数的构建方法首先 ,我们把想要构建的传递函数分子和分母的系数按照阶次从 ...
Matlab：利用Matlab编程实现模拟分子布朗运动的动画展示
Matlab:利用Matlab编程实现模拟分子布朗运动的动画展示目录输出结果实现代码输出结果实现代码 %Brownian motion clf; n=20; s=0.02; x = rand ...
idft重建图像 matlab_利用 MATLAB 编程,打开一幅图像,对其进行 DFT 变换,并置其不同区域内的系数为零,进行 IDFT ,观察其输出效果。_学小易找答案...
[连线题]请对正确的快键键连线 [判断题]板书是指教师在课堂黑板或白板上书写,将教学内容形象.直观.简洁地传授给学生.清晰.流畅.快速的粉笔书写是课堂板书的基本功. [其它]利用 MATLAB 编程, ...
matlab 进行非线性回归,5.利用Matlab编程进行非线性回归分析.doc
5.利用Matlab编程进行非线性回归分析.doc §5. 利用Matlab编程计算非线性回归模型 --以Logistic曲线为例 1.原始数据下表给出了某地区1971-2000年的人口数据(表1) ...

利用 MATLAB 编程实现最速下降法求解无约束最优化问题

本文章包含以下内容

一维搜索，进退法，黄金分割法代码链接

本实验中函数用单独function计算

1.最速下降法的算法流程图

2、MATLAB 编写用数值微分法的梯度计算函数（函数式 M 文件）；

3、最速下降法求解无约束优化问题的函数，采用黄金分割法精确一维搜索，用数值微分法计算梯度（函数式 M 文件，精度设为 epson 可调）；

4、最速下降法求解无约束优化问题的函数，要求采用 Wolfe-Powell 非精确一维搜索，用数值微分法计算梯度（函数式 M 文件，精度设为 epson 可调）；

5、分别利用精确搜索和不精确搜索的最速下降法，问题：

结果：

利用 MATLAB 编程实现最速下降法求解无约束最优化问题相关推荐

最新文章

热门文章

利用 MATLAB 编程实现最速下降法求解无约束最优化问题

本文章包含以下内容

一维搜索，进退法，黄金分割法代码链接

本实验中函数用单独function计算

1.最速下降法的算法流程图

2、MATLAB 编写用数值微分法的梯度计算函数（函数式 M 文件）；

3、最速下降法求解无约束优化问题的函数，采用黄金分割法精确一维搜索，用数值微分法计算梯度（函数式 M 文件，精度设为 epson 可调）；

4、最速下降法求解无约束优化问题的函数，要求采用 Wolfe-Powell 非精确一维搜索，用数值微分法计算梯度（函数式 M 文件， 精度设为 epson 可调）；

5、分别利用精确搜索和不精确搜索 的最速下降法，问题：

结果：

利用 MATLAB 编程实现最速下降法求解无约束最优化问题相关推荐

最新文章

热门文章

4、最速下降法求解无约束优化问题的函数，要求采用 Wolfe-Powell 非精确一维搜索，用数值微分法计算梯度（函数式 M 文件，精度设为 epson 可调）；

5、分别利用精确搜索和不精确搜索的最速下降法，问题：