- 博客(6)
- 收藏
- 关注
原创 0.618法matlab
0.618法建立golds.m文件function [s,phis,k,G,E]=golds(phi,a,b,delta,epsilon)%输入:phi是目标函数,a,b是搜索区间的两个端点% delta,epsilon分别是自变量和函数值的容许误差%输出:s,phis分别是近似极小点和极小值,G是nx4矩阵% 其第k行分贝时a,p,q,b的第k次迭代值[ak,pk,qk,bk]% E=[ds,dphi],分别是s和phis的误差限t=(sqrt(5)-1)/2;h=
2021-05-21 17:14:04
1967
原创 乘子法matlab
乘子法建立multphr.m文件function [x,mu,lambda,output]=multphr(fun,hf,gf,dfun,dhf,dgf,x0)%功能:用乘子法一般解约束问题;min f(x),s.t.h(x)=0,g(x)>=0%输入:x0是初始点,fun,dfun分别是目标函数及其梯度;%hf,dhf分别是等式约束(向量)函数及其Jacobi矩阵的转置;%gf,dgf分别是不等式约束(向量)函数及其Jacobi矩阵的转置;%输出:x是近似最优点,mu,lambda分别
2021-05-21 17:12:48
2223
1
原创 共轭梯度法matlab
共轭梯度法建立frcg.m文件function [x,val,k]=frcg(fun,gfun,x0)%功能:用FR共轭梯度法求解无约束问题:min f(x)%输入:x0是初始点,fun,gfun分别是目标函数和梯度%输出:x,val分别是近似最优点和最优值,k是迭代次数maxk=5000; %最大迭代次数rho=0.6;sigma=0.4;k=0; epsilon=1e-4;n=length(x0);while(k<maxk)g=feval(gfun,x0); %计算梯度it
2021-05-21 17:11:50
1779
1
原创 内点法matlab
内点法建立neidianfa.m文件function [x,tao,output]=neidianfa(fun,gf,dfun,dgf,ddfun,ddgf,x0,tao)maxk=500;%最大迭代次数ganma=0.7;%算法中的实参数k=0;ink=0;%k,ink分别是外迭代和内迭代次数epsilon=1e-5;%终止条件x=x0;gi=feval(gf,x);n=length(x);m=length(gi);btak=10;while(abs(btak)>epsilon
2021-05-21 17:10:03
3217
原创 最优化拟牛顿法matlab
最优化备忘录(拟牛顿法,matlab)建立bfgs.m文件function[x,val,k]=bfgs(fun,gfun,x0,varargin)%功能:用BFGS算法求解无约束问题:min f(x)%输入:x0是初始点,fun,gfun分别是目标函数及其梯度;%varargin是输入的可变参数变量,简单调用bfgs时可以忽略它%但若其他程序循环调用该程序时发挥重要的作用%输出:x,val分别是近似最优点和最优值,是迭代次数maxk=500;%给出最大迭代次数rho=0.55;sigma=
2021-05-21 17:08:10
1101
原创 最优化随笔(外罚函数)
最优化随笔(外罚函数)通过matlab运行建立waifahanshu.m文件function [x,sigma,output]=waifahanshu(fun,hf,gf,dfun,dhf,dgf,x0,sigma)maxk=500;%最大迭代次数ganma=2;%算法中的实参数k=0;ink=0;%k,ink分别是外迭代和内迭代次数epsilon=1e-5;%终止条件x=x0;he=feval(hf,x);gi=feval(gf,x);n=length(x);l=length(he)
2021-05-21 17:04:44
263
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人