| 查看: 3269 | 回复: 18 | |||
[求助]
matlab拟合参数显示初始点是一个局部最小值 已有5人参与
|
|
小弟最近刚学习用matlab拟合参数，但编写完运行显示Initial point is a local minimum.得到的k值为k0，求给位大神给指点一下怎么解决这种问题！急求，新人啥也不懂啊！先在此谢过啦！ function zixie clear all; clc; data=[0 8.32 0; 1 7.37 0.211; 2 6.74 0.738; 3 5.92 1.106; 4 5.68 1.296; 5 5.01 1.583; 6 6.09 1.364; 8 5.35 1.84; 10 4.99 2]; x0=[8.32 0]; k0=[10 10 10 10 10 ]; lb=[0 0 0 0 0]; ub=[+inf +inf +inf +inf +inf]; yexp=data(:,2:3); [k,resnorm,residual,exitflag,output,lambda,jacobin]=lsqnonlin(@objFunc1,k0,lb,ub,[],x0,yexp) fprintf('\tk1=%.4\n',k(1)), fprintf('\tk2=%.4\n',k(2)), fprintf('\tk3=%.4\n',k(3)), fprintf('\tk4=%.4\n',k(4)), fprintf('\tk5=%.4\n',k(5)) function f=objFunc1(k,x0,yexp) tspan=[0 1 2 3 4 5 6 8 10]; [t x]=ode45(@funceqs,tspan,x0,[],k); y(:,1)=x(:,1); y(:,2)=x(:,2); f1=y(:,1)-yexp(:,1); f2=y(:,2)-yexp(:,2); f=[f1;f2]; function dxdt=funceqs(t,x,k) dx1dt=-k(1)*k(2)*x(1)*sqrt(k(3)*x(2))/(k(1)+k(2)*x(1))-k(4)*k(5)*x(2); dx2dt=k(1)*k(2)*x(1)*sqrt(k(3)*x(2))/(k(1)+k(2)*x(1))+k(5)*x(2); dxdt=[dx1dt;dx2dt]; Initial point is a local minimum. Optimization completed because the size of the gradient at the initial point is less than the default value of the function tolerance. <stopping criteria details> k = 10 10 10 10 10 Optimization completed: The final point is the initial point. The first-order optimality measure, 0.000000e+00, is less than options.TolFun = 1.000000e-06. Optimization Metric Options relative first-order optimality = 0.00e+00 TolFun = 1e-06 (default) >> |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有86人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 拜求！matlab拟合动力学参数
已经有14人回复
- 关于matlab两个自变量、一个因变量曲线拟合问题
已经有3人回复
- 求助一个最小二乘法matlab拟合程序。
已经有7人回复
- 求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 用Matlab做多变量非线性拟合求参数值，求各位大侠鼎力相助！
已经有3人回复
- 怎样用一组参数同时拟合两个曲线--matlab
已经有5人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- Matlab 含有两个参数的公式拟合
已经有4人回复
- 用MATLAB怎么打开点云并把它拟合到一个平面上？
已经有8人回复
- 模式搜索MATLAB程序
已经有9人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab如何编写共享参数拟合程序
已经有8人回复
- matlab参数拟合
已经有6人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- matlab拟合模型参数
已经有5人回复
- matlab非线性参数拟合问题
已经有7人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】如何使得matlab拟合的曲线强制经过一个点
已经有8人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68131.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-02-26 22:02:02
戴钢盔的猪头
木虫 (知名作家)
- 应助: 200 (高中生)
- 金币: 5083.9
- 散金: 1310
- 红花: 18
- 帖子: 6001
- 在线: 319.9小时
- 虫号: 2854040
- 注册: 2013-12-07
- 性别: GG
- 专业: 波谱分析与成像分析
3楼2014-02-26 22:09:59
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★
感谢参与，应助指数 +1
fegg7502: 金币+1, 鼓励交流 2014-02-27 08:43:32
感谢参与，应助指数 +1
fegg7502: 金币+1, 鼓励交流 2014-02-27 08:43:32
|
x2的起始值为0（0 8.32 0）有问题吧？会导致计算值全部为0； 还有公式还可以化简，否则有过拟合现象。 dx1dt=-k1*k2*x1*sqrt(k3*x2)/(k1+k2*x1)-k4*k5*x2; dx2dt=k1*k2*x1*sqrt(k3*x2)/(k1+k2*x1)+k5*x2; 可简化为： dx1dt=-x1*sqrt(x2)/(k1+k2*x1)-k4*k5*x2; dx2dt=x1*sqrt(x2)/(k1+k2*x1)+k5*x2; |
4楼2014-02-26 23:10:52
pdl9527
专家顾问 (小有名气)
-
专家经验: +8 - 应助: 100 (初中生)
- 金币: 2110.3
- 红花: 23
- 帖子: 282
- 在线: 112小时
- 虫号: 1227333
- 注册: 2011-03-09
- 性别: GG
- 专业: 分离过程
- 管辖: 计算模拟
5楼2014-02-27 00:00:32
6楼2014-02-27 09:45:06
7楼2014-02-27 09:49:39
戴钢盔的猪头
木虫 (知名作家)
- 应助: 200 (高中生)
- 金币: 5083.9
- 散金: 1310
- 红花: 18
- 帖子: 6001
- 在线: 319.9小时
- 虫号: 2854040
- 注册: 2013-12-07
- 性别: GG
- 专业: 波谱分析与成像分析
8楼2014-02-27 10:05:07
9楼2014-02-27 10:05:29
10楼2014-02-27 10:06:54
20