| 查看: 3763 | 回复: 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) >> |
回复此楼» 猜你喜欢
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有66人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有120人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 拜求！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
- 金币: 68177.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.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
