| 查看: 455 | 回复: 1 | |||
[交流]
matlab牛顿法如何调用问题，调用部分总是出错 已有1人参与
|
|
牛顿法的程序： function [minx,minf]=example3g(f,x0,var,eps) %目标函数：f %初始点：x0 %自变量向量var %精度：eps %目标函数取最小值时的自变量值：minx; %目标函数的最小值：minf format long; if nargin==4 eps=1.0e-6; end tol=1; x0=transpose(x0); while tol>eps gradutef=jacobian(f,var); %梯度方向 %其大概意思是：函数jacobian就是雅可比矩阵, %命令jacobian（f，v）是计算标量函数或矢量函数f关于矢量v的雅可比矩阵。 %其结果就相当于df/dv。需要注意的是，当f是标量时，就相当于是求f的梯度。 jacf=Jacobian(gradutef,var); %雅克比矩阵 v=Funval(gradutef,var,x0); tol=norm(v); pv=Funval(jacf,var,x0); p=-inv(pv)*transpose(v); %搜索方向 x1=x0+p; x0=x1; end minx=x1; minf=Funval(f,var,minx); format short; 想优化的函数： function f=op(x) f=561+7.92*x(1)+0.001562*x(1)^2+310+7.85*x(2)+0.00194*x(2)^2+78+7.97*x(3)+0.00482*x(3)^2+x(4)*(800-x(1)-x(2)-x(3)); 调用部分，不知道var部分如何初始化 syms x x0; x0=[300 200 300 10]; [minx,minf]=example3g(op,x0,[x(1) x(2) x(3) x(4)],1e-6) |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
-
什么时候出结果，有咨询渠道？？？
已经有3人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有26人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有7人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有8人回复
-
分享一下我之前已中青C的计划书的filecode
已经有4人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有3人回复
-
听说今天filecode变了
已经有49人回复
48038074
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 147
- 散金: 2
- 帖子: 37
- 在线: 38.4小时
- 虫号: 4032798
- 注册: 2015-08-21
- 性别: GG
- 专业: 电力系统
2楼2015-11-23 00:40:35
100