| 查看: 719 | 回复: 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); |
回复此楼» 猜你喜欢
22408 344分 求调剂 一志愿 华电计算机技术
已经有4人回复
-
材料专硕英一数二306
已经有7人回复
-
南京大学化学376求调剂
已经有6人回复
-
南昌大学材料专硕311分求调剂
已经有5人回复
-
295求调剂
已经有6人回复
-
一志愿南昌大学，327分，材料与化工085600
已经有9人回复
-
294求调剂材料与化工专硕
已经有15人回复
-
308求调剂
已经有3人回复
-
321求调剂
已经有11人回复
-
一志愿中海洋材料工程专硕330分求调剂
已经有8人回复
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
champion3568
木虫 (职业作家)
- 应助: 0 (幼儿园)
- 金币: 3525.6
- 散金: 6880
- 红花: 14
- 沙发: 3
- 帖子: 4187
- 在线: 320.3小时
- 虫号: 4221360
- 注册: 2015-11-15
- 性别: GG
- 专业: 机构学与机器人
4楼2015-11-20 23:45:39
10