[问答]

DFP与BFGS算法的原理是什么？

850 DFP

问答对人有帮助，内容完整，我也想知道答案 0 DFP算法原理是什么？ BFGS算法原理是什么？ 0
2021-11-19 07:26:36　　评论淘帖0 邀请回答您可以邀请以下用户，快速回答问题 × heks 该类别下有 50 个回答。邀请回答 hgimtk 该类别下有 44 个回答。邀请回答新星之火12138 该类别下有 42 个回答。邀请回答 wang21cj 该类别下有 41 个回答。邀请回答 chm5 该类别下有 40 个回答。邀请回答 hjfjsdgfjdsf 该类别下有 37 个回答。邀请回答 fdjslkjd 该类别下有 35 个回答。邀请回答 uvysdfydad 该类别下有 35 个回答。邀请回答 werywer 该类别下有 35 个回答。邀请回答 h1654155957.9520 该类别下有 35 个回答。邀请回答 jenny042 该类别下有 34 个回答。邀请回答 dfzvzs 该类别下有 34 个回答。邀请回答维生素B2 该类别下有 34 个回答。邀请回答凤毛麟角该类别下有 34 个回答。邀请回答 hfgdzc 该类别下有 33 个回答。邀请回答 ggfvxv 该类别下有 33 个回答。邀请回答刘洋该类别下有 33 个回答。邀请回答 meihuacg 该类别下有 33 个回答。邀请回答储蓄叛逆该类别下有 32 个回答。邀请回答尼克wo 该类别下有 32 个回答。邀请回答举报发生的方式相关推荐 • 请问哪位有Keil.STM32F4xx_DFP.2.8.0.pack的离线文件吗 8717 • 导入Nuvoton.NuMicro_DFP.1.1.0.pack时出错是什么原因？ 1515 • 请问Keil.WCH57x_DFP.1.2.0.pack要去哪下载？ 1016 • stm32f0x的DFP文件导入不成功 3885 • HDMI信号的接收端和发送端与PD的DFP和UFP有对应关系吗？ 333 • keil 5.14.0.0导入Nuvoton.NuMicro_DFP.1.1.0.pack时出错的原因？ 224 • 请问Keil.STM32F1xx_DFP.1.0.5代码无法烧录怎么解决？ 835 • M6，RHEL 6.6，无法连接将EDID监视器连接到有效的显示设备 2121 • 关于type c接口充电的疑问，CC引脚检测 26746 • 使用keil5下载了STM32F3xx_DFP的包，，该怎么用keil5编写stm32f301k6tx的程序？ 214 1个回答

答案对人有帮助，有参考价值 0 DFP算法原理由于博主使用WPS编辑的文本，公式无法赋值粘贴，这里以截图的方法给出了推导过程。博主会上传该DOC文档。牛顿法阻尼牛顿法 BFGS算法原理 matlab代码(DFP) syms x1 x2 f=@(x1,x2) x1.^2+x2.^2-x1x2-10x1-4x2+60; X=DFP(f,[0 0],1e-8,100) function x=DFP(f,x0,eps,k) %目标函数f %初始迭代点x0 %迭代精度eps %迭代次数K x0=x0'; TiDu=matlabFunction(gradient(sym(f))); m=length(x0); H=eye(m);%构造初始海塞矩阵 %写出函数值，梯度值的表达式 f_cal='f(x0(1),x0(2))'; %初始点函数值 tidu_cal='TiDu(x0(1),x0(2))'; f1_cal='f(x_1(1),x_1(2))';%下一点函数值 %% 从第一个点计算到第二个点 f_x0=eval(f_cal);%计算初始点函数值 tidu_x0=eval(tidu_cal);%计算初始点梯度值 if norm(tidu_x0) < eps%判断是否满足终止条件 x=x0; return; end d=-Atidu_x0;% 下降方向 syms alfa %定义步长 x_1=x0+alfad;%更新迭代点位置 f_x1=eval(f1_cal);%计算迭代点的函数表达式 df_x1=diff(f_x1);%迭代点的梯度表达式 dalfa=double(solve(df_x1));%求解，得到步长alfa x0=x0+dalfad;%更新初始点 tidu_x1=eval(tidu_cal);%计算该点梯度 n=1; while n < k && norm(tidu_x1) > eps delta_x=dalfad; %计算sK delta_g=tidu_x1-tidu_x0; %计算yk delta_H=delta_xdelta_x'/(delta_x'delta_g)- Hdelta_gdelta_g'H/(delta_g'Hdelta_g); %计算delta_Dk H=H+delta_H; %dfp迭代公式 tidu_x0=tidu_x1; %将该点梯度作为新的初始点继续迭代 d=-Atidu_x0;% 下降方向 syms alfa x_1=x0+alfad; f_x1=eval(f1_cal); df_x1=diff(f_x1); dalfa=double(solve(df_x1)); x0=x0+dalfad; tidu_x1=eval(tidu_cal); end x=x0; end matlab代码(BFGS) %% BFGS算法与DFP算法过程类似，只是迭代函数不同 clc clear syms x1 x2 f=@(x1,x2) x1.^2+x2.^2-x1x2-10x1-4x2+60; X=BFGS(f,[0 0],1e-8,100) function x=BFGS(f,x0,eps,k) x0=x0'; TiDu=matlabFunction(gradient(sym(f))); m=length(x0); A=eye(m);%构造初始矩阵 f_cal='f(x0(1),x0(2))'; tidu_cal='TiDu(x0(1),x0(2))'; f1_cal='f(x_1(1),x_1(2))'; %从第一个点计算到第二个点 f_x0=eval(f_cal); tidu_x0=eval(tidu_cal); if norm(tidu_x0) < eps x=x0; return; end d=-Atidu_x0;% 下降方向 syms alfa x_1=x0+alfad; f_x1=eval(f1_cal); df_x1=diff(f_x1); dalfa=double(solve(df_x1)); x0=x0+dalfad; tidu_x1=eval(tidu_cal); n=1; while n < k && norm(tidu_x1) > eps delta_x=dalfad; delta_g=tidu_x1-tidu_x0; delta_A=-Adelta_xdelta_x'A/(delta_x'Adelta_x)+delta_gdelta_g'/(delta_g'delta_x); A=A+delta_A; tidu_x0=tidu_x1; d=-Atidu_x0;% 下降方向 syms alfa x_1=x0+alfad; f_x1=eval(f1_cal); df_x1=diff(f_x1); dalfa=double(solve(df_x1)); x0=x0+dalfad; tidu_x1=eval(tidu_cal); end x=x0; end

2021-11-19 14:30:23 评论举报李淑兰