| 查看: 723 | 回复: 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); |
回复此楼» 猜你喜欢
346求调剂[0856]
已经有7人回复
-
一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！
已经有3人回复
-
310求调剂
已经有3人回复
-
机械专硕299求调剂至材料
已经有4人回复
-
070300化学319求调剂
已经有7人回复
-
08工科 320总分 求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有12人回复
-
307求调剂
已经有3人回复
-
265求调剂
已经有3人回复
-
301求调剂
已经有10人回复
champion3568
木虫 (职业作家)
- 应助: 0 (幼儿园)
- 金币: 3525.6
- 散金: 6880
- 红花: 14
- 沙发: 3
- 帖子: 4187
- 在线: 320.3小时
- 虫号: 4221360
- 注册: 2015-11-15
- 性别: GG
- 专业: 机构学与机器人
4楼2015-11-20 23:45:39
yanshuiting
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1484.8
- 散金: 50
- 红花: 7
- 帖子: 448
- 在线: 68.7小时
- 虫号: 2633300
- 注册: 2013-09-03
- 专业: 计算机科学
2楼2015-11-20 22:28:12
heaven0419
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2890
- 帖子: 188
- 在线: 214.7小时
- 虫号: 3195508
- 注册: 2014-05-10
- 性别: GG
- 专业: 信号理论与信号处理
3楼2015-11-20 22:31:20
