| 查看: 395 | 回复: 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) |
回复此楼» 猜你喜欢
材料学硕297已过四六级求调剂推荐
已经有8人回复
-
本人考085602 化学工程 专硕
已经有21人回复
-
265求调剂
已经有12人回复
-
求调剂
已经有5人回复
-
一志愿中南化学337求调剂
已经有6人回复
-
350求调剂
已经有3人回复
-
295复试调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有22人回复
-
287求调剂
已经有7人回复
-
329求调剂
已经有7人回复
48038074
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 147
- 散金: 2
- 帖子: 37
- 在线: 38.4小时
- 虫号: 4032798
- 注册: 2015-08-21
- 性别: GG
- 专业: 电力系统
2楼2015-11-23 00:40:35
