| 查看: 359 | 回复: 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) |
回复此楼» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
48038074
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 147
- 散金: 2
- 帖子: 37
- 在线: 38.4小时
- 虫号: 4032798
- 注册: 2015-08-21
- 性别: GG
- 专业: 电力系统
2楼2015-11-23 00:40:35