| 查看: 657 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
champion3568木虫 (职业作家)
|
[求助]
请各位大神帮我看看关于拟牛顿法Maltab代码到底哪里有错？？
|
|
|
代码我是依照一本教程做的，可运行时就出现如下错误： syms x1 x2 >> x0=[1,2]'; >> [x,val,k]=sr1(fun,gfun,x0) ??? Undefined function or variable 'fun'. 这个是代码： function [x,val,k] =sr1(fun,gfun,x0) %功能：用对称秩1算法求解无约束问题：min f(x) %输入：X0是初始点，fun,gfun分别是目标函数及其梯度 %输出：x，val分别是近似最优点和最优值，k是迭代次数. maxk=500; %给出最大迭代次数 rho=0.55;sigma=0.4;epsilon=le-4; k=0; n=length(x0);Hk=eye(n); while(k<maxk) gk=feval(gfun,x0);%计算梯度 dk=-Hk*gk;%计算搜索方向 if(norm(gk)<epsilon),break;end %检验终止准则 m=0;mk=0; while(m<20) %用Armijo搜索求步长 if(feval(fun,x0+rho^m*dk)<feval(fun,x0)+sigma*rho^m*gk'*dk) mk=m;break; end m=m+1; end x=x0+rho^mk*dk; sk=x-x0; yk=feval(gfun,x)-gk; Hk=Hk+(sk-Hk*yk)*(sk-Hk*yk)'/((sk-Hk*yk)'*yk); %秩1校正 k=k+1; x0=x; end val=feval(fun,x0); |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
heaven0419
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2890
- 帖子: 188
- 在线: 214.7小时
- 虫号: 3195508
- 注册: 2014-05-10
- 性别: GG
- 专业: 信号理论与信号处理
3楼2015-11-20 22:31:20
yanshuiting
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1484.8
- 散金: 50
- 红花: 7
- 帖子: 448
- 在线: 68.7小时
- 虫号: 2633300
- 注册: 2013-09-03
- 专业: 计算机科学
2楼2015-11-20 22:28:12
champion3568
木虫 (职业作家)
- 应助: 0 (幼儿园)
- 金币: 3525.6
- 散金: 6880
- 红花: 14
- 沙发: 3
- 帖子: 4187
- 在线: 320.3小时
- 虫号: 4221360
- 注册: 2015-11-15
- 性别: GG
- 专业: 机构学与机器人
4楼2015-11-20 23:45:39