| 查看: 877 | 回复: 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); |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
-
什么时候出结果，有咨询渠道？？？
已经有3人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有26人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有7人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有8人回复
-
分享一下我之前已中青C的计划书的filecode
已经有4人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有3人回复
-
听说今天filecode变了
已经有49人回复
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
100