应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab拟合参数显示初始点是一个局部最小值
192/2返回列表上一页12
查看: 3240  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by 小小莓邮 at 2014-02-27 09:49:39
也就是说这造成这种现象是初值的原因是吧！有没有什么方法可以先初步估计一下初值，然后再以这个值为初值进一步估计参数？...

先调用遗传算法预估一个结果，以此结果作为初值，再调用lsqnonlin函数计算。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
11楼2014-02-27 10:10:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小莓邮

铁虫 (小有名气)
引用回帖:
11楼: Originally posted by 月只蓝 at 2014-02-27 10:10:24
先调用遗传算法预估一个结果，以此结果作为初值，再调用lsqnonlin函数计算。...

哦！还有一个问题就是，用ode返回的x的计算值每个点都等于x0，做出来的线是条水平线，这也是k的初值不对造成的吗？谢谢！
赞一下 回复此楼
12楼2014-02-27 10:20:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小莓邮

铁虫 (小有名气)
引用回帖:
11楼: Originally posted by 月只蓝 at 2014-02-27 10:10:24
先调用遗传算法预估一个结果，以此结果作为初值，再调用lsqnonlin函数计算。...

还有就是我随机换了好几个初值，但还是得到k0!
赞一下 回复此楼
13楼2014-02-27 10:45:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
引用回帖:
9楼: Originally posted by 小小莓邮 at 2014-02-27 10:05:29
x2是产物。t=0时是为0啊，不过计算值出来每个点都是x0这个初值。这是什么原因？是x2=0造成的吗？...

x2=0时，微分方程右边每一项都有x2的乘项，会导致全部为0
赞一下 回复此楼
14楼2014-02-27 11:41:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小莓邮

铁虫 (小有名气)
引用回帖:
14楼: Originally posted by dingd at 2014-02-27 11:41:26
x2=0时，微分方程右边每一项都有x2的乘项，会导致全部为0...

额！那岂不是不解？
赞一下 回复此楼
15楼2014-02-27 11:55:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

你把x2=0改为x2=0.1试试，不改的话计算结果会是一条直线（全为0）
赞一下 回复此楼
16楼2014-02-27 12:07:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小莓邮

铁虫 (小有名气)
引用回帖:
16楼: Originally posted by dingd at 2014-02-27 12:07:10
你把x2=0改为x2=0.1试试，不改的话计算结果会是一条直线（全为0）

将x2改成非零后可以拟合出k，但是拟合曲线的趋势跟试验点的趋势差别过大，这是不是就是k0的原因？有没有什么函数可以估算一下k0 ？
赞一下 回复此楼
17楼2014-02-27 13:59:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
公式按前面所说简化，图方便用1stOpt求解：
CODE:
Variable t,x1,x2;
ODEFunction x1'=-x1*sqrt(x2)/(k1+k2*x1)-k4*k5*x2;
            x2'=x1*sqrt(x2)/(k1+k2*x1)+k5*x2;
Data;
0        8.32        0.1
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

均方差(RMSE):0.420626480368569
残差平方和(SSE):2.83082617579601
相关系数(R): 0.852105775635903
相关系数之平方(R^2): 0.726084252872063
确定系数(DC): 0.575698025223075
F统计(F-Statistic): 4.14392332314321

参数                  最佳估算
--------------------        -------------
k1        1.2026577837956
k2        0.369501550975634
k4        0.889293844320815
k5        -1.33346561077959
matlab拟合参数显示初始点是一个局部最小值
c1.jpg


matlab拟合参数显示初始点是一个局部最小值-1
c2.jpg
赞一下 回复此楼
18楼2014-02-27 15:03:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1314168apple

金虫 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
1、可以把初值x0设定为【8.32 1e-18】1e-15等；可能是你的速率方程有误。
2、设定tspan的时候可以将区间中间隔缩小(如（0:0.01:10）)，调用的时候提取需要的时间对应的浓度。否则有可能产生负数的浓度。一般的速率方程比较少见。
很有可能速率方程有问题。
赞一下 回复此楼
为了最终理解你所不理解的，你必须经历一条愚昧无知的道路。为了占有你从未占有的东西，你必须经历被剥夺的道路。为了达到你现在所不在的名位，你必须经历那...
19楼2014-02-27 18:13:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小小莓邮 的主题更新
192/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085404 293求调剂 +3 勇远库爱314 2026-04-08 3/150 2026-04-09 04:24 by 孙小小12457
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 软件工程求调剂22软工296分求调剂，接受跨调 +4 yangchen2017 2026-04-08 5/250 2026-04-08 21:56 by 土木硕士招生
[考研] 一志愿山东大学化学与化工学院材料与化工专硕，360分求调剂 +8 不愿透露姓名的� 2026-04-02 8/400 2026-04-08 20:12 by luoyongfeng
[考研] 材料调剂 +4 hzhahg 2026-04-06 4/200 2026-04-08 16:20 by luoyongfeng
[考研] 求调剂 +14 张zic 2026-04-05 15/750 2026-04-08 16:14 by 一只好果子?
[考研] 一志愿电子科技大学085600材料与化工 329分求调剂 +11 Naiko 2026-04-04 11/550 2026-04-08 14:00 by wutongshun
[考研] 机械工程264学硕求调剂 +3 qiushangxian 2026-04-06 3/150 2026-04-08 01:53 by Linzejun
[考研] 材料工程302分求调剂 +13 zyx上岸！ 2026-04-04 13/650 2026-04-07 11:14 by 诗与自由
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 生物与医药求调剂 +7 heguanhua 2026-04-05 8/400 2026-04-06 18:41 by macy2011
[考研] 315求调剂 +5 ＆123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 326求调剂 +3 顾若浮生 2026-04-05 3/150 2026-04-05 18:32 by 蓝云思雨
[考研] 358求调剂 +7 秋gk 2026-04-04 7/350 2026-04-05 13:29 by huangmoli
[论文投稿] 求文献 5+3 ys879651$ 2026-04-02 3/150 2026-04-04 17:22 by bobvan
[考研] 生物学308分求调剂（一志愿华东师大） +7 相信必会光芒万� 2026-04-02 7/350 2026-04-03 16:48 by rzh123456
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 材料340分调剂 +7 夏夜晚风_long 2026-04-02 9/450 2026-04-02 21:20 by dongzh2009
[考研] 求调剂 +7 Aniyaio 2026-04-02 7/350 2026-04-02 16:42 by zzsw+
信息提示
关闭
请填处理意见
关闭