24小时热门版块排行榜    

查看: 3647  |  回复: 18

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
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的回帖
相关版块跳转 我要订阅楼主 小小莓邮 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 10/500 2026-08-21 17:32 by chengyan1220
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 18/900 2026-08-21 13:46 by 医学老男孩
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 今天基金会出结果吗?20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[基金申请] 科研孤儿太难了 +15 我4大白菜 2026-08-20 16/800 2026-08-21 11:07 by yming1319
[基金申请] 时间戳今天,20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
信息提示
请填处理意见