24小时热门版块排行榜    

查看: 1705  |  回复: 14

wodemiss

铁虫 (小有名气)

[求助] 一阶微分方程求解,并优化参数值

紧急求助,非常感谢!!!
我有一系列y-t值,如何计算dy/dt=a*y*(b-y)/(c*d+d*y+e*(b-y)+y*(b-y)).
或者说我已知a,b,c,d,e的值如何根据时间t的变化求出Y的值,然后依据求得的y-t值来拟合y-t值求得最佳的参数值a,c,d,e。
十分感谢,现在对这个相当困惑。希望大家指导一下。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

一蓑烟雨任平生
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wodemiss: 金币+20, ★★★很有帮助 2012-05-07 15:31:12
微分方程拟合可以试下1stOpt,简单易上手。有数据的话贴上来可以帮你试试。

» 本帖已获得的红花(最新10朵)

2楼2012-05-07 11:57:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodemiss

铁虫 (小有名气)

送鲜花一朵
引用回帖:
2楼: Originally posted by dingd at 2012-05-07 11:57:38:
微分方程拟合可以试下1stOpt,简单易上手。有数据的话贴上来可以帮你试试。

嗯,我现在只是这个思路,还没有数据,而且我以后还会加参数,进行拟合,会很麻烦,唉,我学生物的,做反应动力学,但是不懂计算,我在网上下到了1.5破解版,是不是没有这个功能,我以后一两个月会经常用,也不好一直让你弄吧。。。如果买1stopt4.0或5.0版需要多少钱?
一蓑烟雨任平生
3楼2012-05-07 14:43:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

六本木

新虫 (小有名气)

★ ★
zhangguangping: 金币+2, 谢谢提示! 2012-05-14 19:10:42
初步看了一下,你这个其实不难,用matlab,mathtype,或者是自己写一个小程序(具体你可以参考任何一本数值分析的书,如果不需要太高的精度的话,你就用欧拉方法)就可以解决你的第一个问题。
4楼2012-05-10 10:19:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wjx0162

银虫 (著名写手)


zhangguangping: 金币+1, 谢谢提醒! 2012-05-14 19:11:19
想将y-t值做数值微分得到dy/dt-t的曲线,然后用最小二乘法拟合该曲线,得到abcde。
5楼2012-05-10 11:03:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodemiss

铁虫 (小有名气)

引用回帖:
5楼: Originally posted by wjx0162 at 2012-05-10 11:03:11:
想将y-t值做数值微分得到dy/dt-t的曲线,然后用最小二乘法拟合该曲线,得到abcde。

嗯,谢谢,版主已经帮我了,谢谢。
一蓑烟雨任平生
6楼2012-05-10 16:10:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodemiss

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by 六本木 at 2012-05-10 10:19:09:
初步看了一下,你这个其实不难,用matlab,mathtype,或者是自己写一个小程序(具体你可以参考任何一本数值分析的书,如果不需要太高的精度的话,你就用欧拉方法)就可以解决你的第一个问题。

嗯,谢谢,版主已经帮我解决了,谢谢
一蓑烟雨任平生
7楼2012-05-10 16:11:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodemiss

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2012-05-07 11:57:38:
微分方程拟合可以试下1stOpt,简单易上手。有数据的话贴上来可以帮你试试。

你好,你能帮我跑一下程序吗?非常感谢!!!
方程dy/dt = V*y*(2-y)/(Kia*Kmb+Kmb*y+Kma*(2-y)+y*(2-y))
数据1:
0        0.129421222
0.167        0.127813505
0.333        0.126366559
0.5        0.124758842
0.667        0.123633441
0.833        0.122508039
1        0.121061093
1.167        0.119614148
1.333        0.118649518
1.5        0.117363344
1.667        0.116237942
1.833        0.114951768
2        0.113826367
2.167        0.112861736
2.333        0.111897106
2.5        0.110771704
2.667        0.109485531
2.833        0.1085209
3        0.10755627
3.167        0.106430868
3.333        0.10562701
3.5        0.104501608
3.667        0.103697749
3.833        0.102572347
4        0.101768489
4.167        0.100803859
4.333        0.1
4.5        0.098874598
4.667        0.09807074
4.833        0.097106109
5        0.096463023
5.167        0.095659164
5.333        0.094533762
5.5        0.093729904
5.667        0.092926045
5.833        0.092122186
6        0.091318328
6.167        0.090514469
6.333        0.089710611
6.5        0.089067524
6.667        0.088263666
6.833        0.087459807
7        0.086655949
7.167        0.086012862
7.333        0.085209003
7.5        0.084405145
7.667        0.083762058
7.833        0.082958199
8        0.082154341
8.167        0.081511254
8.333        0.080868167
8.5        0.080064309
8.667        0.079421222
8.833        0.078617363
9        0.078135048
9.167        0.077491961
9.333        0.076848875
9.5        0.076045016
9.667        0.075401929
9.833        0.074758842
10        0.074115756
10.167        0.073472669
10.333        0.072829582
10.5        0.072347267
10.667        0.07170418
10.833        0.071061093
11        0.070418006
11.167        0.06977492
11.333        0.069292605
11.5        0.068649518
11.667        0.068167203
11.833        0.067524116
12        0.066881029
12.167        0.066398714
12.333        0.065916399
12.5        0.065273312
12.667        0.064630225
12.833        0.06414791
13        0.063665595
13.167        0.06318328
13.333        0.062540193
13.5        0.062057878
13.667        0.061575563
13.833        0.061093248
14        0.060610932
14.167        0.059967846
14.333        0.059485531
14.5        0.059003215
14.667        0.0585209
14.833        0.058199357
15        0.05755627
15.167        0.057073955
15.333        0.056752412
15.5        0.056109325
15.667        0.055787781
15.833        0.055144695
16        0.054823151
16.167        0.054501608
16.333        0.053858521
16.5        0.053376206
16.667        0.053054662
16.833        0.052572347
17        0.052090032
17.167        0.051768489
17.333        0.051286174
17.5        0.050803859
17.667        0.050482315
17.833        0.05
18        0.049678457
18.167        0.049196141
18.333        0.048874598
18.5        0.048392283
18.667        0.04807074
18.833        0.047749196
19        0.047427653
19.167        0.046945338
19.333        0.046463023
19.5        0.046141479
19.667        0.045819936
19.833        0.045498392
20        0.045176849
20.167        0.044694534
20.333        0.04437299
20.5        0.044051447
20.667        0.043729904
20.833        0.04340836
21        0.042926045
21.167        0.042765273
21.333        0.04244373
21.5        0.041961415
21.667        0.041800643
21.833        0.041318328
22        0.040996785
22.167        0.040836013
22.333        0.040514469
22.5        0.040192926
22.667        0.039871383
22.833        0.039549839
23        0.039228296
23.167        0.038906752
23.333        0.038745981
23.5        0.038424437
23.667        0.038102894
23.833        0.03778135
24        0.037459807
24.167        0.037299035
24.333        0.036977492
24.5        0.036655949
24.667        0.036495177
24.833        0.036173633
25        0.03585209
25.167        0.035691318
25.333        0.035530547
25.5        0.035209003
25.667        0.035048232
25.833        0.03488746
26        0.034565916
26.167        0.034565916
26.333        0.034083601
26.5        0.034083601
26.667        0.033762058
26.833        0.033440514
27        0.033279743
27.167        0.033118971
27.333        0.032797428
27.5        0.032636656
27.667        0.032315113
27.833        0.032154341
28        0.031993569
28.167        0.031672026
28.333        0.031672026
28.5        0.031350482
28.667        0.031189711
28.833        0.030868167
29        0.030707395
29.167        0.030546624
29.333        0.03022508
29.5        0.03022508
29.667        0.029903537
29.833        0.029742765
30        0.029581994
30.167        0.02926045
30.333        0.029099678
30.5        0.029099678
30.667        0.028778135
30.833        0.028778135
31        0.028456592
31.167        0.02829582
31.333        0.028135048
31.5        0.027974277
31.667        0.027813505
31.833        0.027491961
32        0.027491961
32.167        0.02733119
32.333        0.027170418
32.5        0.027009646
32.667        0.026848875
32.833        0.026688103
33        0.026527331
33.167        0.026366559
33.333        0.026366559
33.5        0.026045016
33.667        0.025884244
33.833        0.025723473
34        0.025723473
34.167        0.025401929
34.333        0.025241158
34.5        0.025241158
34.667        0.025080386
34.833        0.024919614
35        0.024919614
35.167        0.024758842
35.333        0.024598071
35.5        0.024437299
35.667        0.024276527
35.833        0.024115756
36        0.024115756
36.167        0.023954984
36.333        0.023794212
36.5        0.023633441
36.667        0.023633441
36.833        0.023472669
37        0.023311897
37.167        0.023151125
37.333        0.023151125
37.5        0.022990354
37.667        0.022829582
37.833        0.02266881
38        0.02266881
38.167        0.022508039
38.333        0.022508039
38.5        0.022347267
38.667        0.022186495
38.833        0.022186495
39        0.022025723
39.167        0.021864952
39.333        0.021864952
39.5        0.021543408
39.667        0.021382637
39.833        0.021382637
40        0.021382637
数据2
0        0.065755627
0.167        0.064469453
0.333        0.062861736
0.5        0.061414791
0.667        0.059967846
0.833        0.058681672
1        0.057395498
1.167        0.056270096
1.333        0.054983923
1.5        0.054019293
1.667        0.052733119
1.833        0.051607717
2        0.050482315
2.167        0.049678457
2.333        0.048553055
2.5        0.047588424
2.667        0.046623794
2.833        0.045659164
3        0.044694534
3.167        0.043890675
3.333        0.042926045
3.5        0.042122186
3.667        0.041157556
3.833        0.040353698
4        0.039549839
4.167        0.038745981
4.333        0.037942122
4.5        0.037299035
4.667        0.036495177
4.833        0.035691318
5        0.03488746
5.167        0.034244373
5.333        0.033601286
5.5        0.032958199
5.667        0.032315113
5.833        0.031511254
6        0.031028939
6.167        0.030385852
6.333        0.029742765
6.5        0.02926045
6.667        0.028617363
6.833        0.027974277
7        0.027491961
7.167        0.027009646
7.333        0.026527331
7.5        0.025884244
7.667        0.025401929
7.833        0.024919614
8        0.024437299
8.167        0.023794212
8.333        0.023472669
8.5        0.022990354
8.667        0.022508039
8.833        0.022025723
9        0.021543408
9.167        0.021221865
9.333        0.020900322
9.5        0.020418006
9.667        0.020096463
9.833        0.01977492
10        0.019292605
10.167        0.018971061
10.333        0.018649518
10.5        0.018327974
10.667        0.017845659
10.833        0.017684887
11        0.017363344
11.167        0.017041801
11.333        0.016720257
11.5        0.016398714
11.667        0.01607717
11.833        0.015916399
12        0.015594855
12.167        0.015434084
12.333        0.01511254
12.5        0.014790997
12.667        0.014630225
12.833        0.014308682
13        0.01414791
13.167        0.013987138
13.333        0.013665595
13.5        0.013504823
13.667        0.01318328
13.833        0.013022508
14        0.012861736
14.167        0.012861736
14.333        0.012540193
14.5        0.012379421
14.667        0.01221865
14.833        0.011897106
15        0.011736334
15.167        0.011575563
15.333        0.011414791
15.5        0.011414791
15.667        0.011254019
15.833        0.011093248
16        0.010932476
16.167        0.010771704
16.333        0.010771704
16.5        0.010450161
16.667        0.010450161
16.833        0.010289389
17        0.010289389
17.167        0.010128617
17.333        0.009967846
17.5        0.009967846
17.667        0.009807074
17.833        0.009807074
18        0.009646302
18.167        0.009485531
18.333        0.009485531
18.5        0.009485531
18.667        0.009324759
18.833        0.009163987
19        0.009163987
19.167        0.009163987
19.333        0.009003215
19.5        0.008842444
19.667        0.008842444
19.833        0.008842444
20        0.008842444
20.167        0.0085209
20.333        0.008681672
20.5        0.008681672
20.667        0.0085209
20.833        0.0085209
21        0.0085209
21.167        0.008360129
21.333        0.008360129
21.5        0.008360129
21.667        0.008199357
21.833        0.008199357
22        0.008199357
22.167        0.008038585
22.333        0.008199357
22.5        0.008038585
22.667        0.008038585
22.833        0.008038585
23        0.008038585
23.167        0.008038585
23.333        0.007877814
23.5        0.007877814
23.667        0.007717042
23.833        0.007877814
24        0.007877814
24.167        0.007877814
24.333        0.007877814
24.5        0.007717042
24.667        0.007717042
24.833        0.007717042
25        0.007717042
25.167        0.007717042
25.333        0.007717042
25.5        0.007717042
25.667        0.007717042
25.833        0.007717042
26        0.00755627
26.167        0.00755627
26.333        0.00755627
26.5        0.00755627
26.667        0.00755627
26.833        0.00755627
27        0.00755627
27.167        0.00755627
27.333        0.00755627
27.5        0.00755627
27.667        0.00755627
27.833        0.00755627
28        0.00755627
28.167        0.00755627
28.333        0.00755627
28.5        0.00755627
28.667        0.00755627
28.833        0.00755627
29        0.00755627
29.167        0.00755627
29.333        0.00755627
29.5        0.00755627
29.667        0.00755627
29.833        0.00755627
30        0.00755627
一蓑烟雨任平生
8楼2012-05-14 15:29:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodemiss

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2012-05-07 11:57:38:
微分方程拟合可以试下1stOpt,简单易上手。有数据的话贴上来可以帮你试试。

恳请帮助,谢谢
一蓑烟雨任平生
9楼2012-05-14 16:44:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodemiss

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2012-05-07 11:57:38:
微分方程拟合可以试下1stOpt,简单易上手。有数据的话贴上来可以帮你试试。

V=[-10,0],Kia=[0, 1],kmb=[0,20],kma=[0,1]
这是约束条件。
恳请帮助,谢谢
一蓑烟雨任平生
10楼2012-05-14 17:52:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wodemiss 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 论文已接收,但发现修改稿传的是旧版该怎么办? +3 zsq54321 2024-06-15 6/300 2024-06-17 19:45 by wjykycg
[基金申请] 有人中过人文社科类的博后特助吗? +3 outsider1986 2024-06-16 4/200 2024-06-17 19:35 by outsider1986
[基金申请] 最新消息:2024国自然下载文件名变了 +15 dream200606 2024-06-16 28/1400 2024-06-17 19:04 by 6543yes
[考博] 34岁读博士晚吗 +30 emitdne 2024-06-13 30/1500 2024-06-17 17:08 by 白米粒向日葵
[有机交流] 车间生产,真空度很高,温度很高,但减压蒸馏速度很慢。 10+12 召唤鬼泣lL 2024-06-13 38/1900 2024-06-17 16:46 by 科研虫子小助手
[硕博家园] 硕博巨婴,也许才刚刚开始 +32 SNaiL1995 2024-06-12 84/4200 2024-06-17 16:04 by 20171214
[论文投稿] 审稿人含糊拒稿,还需要回复吗?如何回复? 20+3 BruceChum 2024-06-15 19/950 2024-06-17 12:03 by 投必得科研顾问
[找工作] 江西双非一本和四川双一流高校如何选择? 5+8 寒山敲钟 2024-06-12 25/1250 2024-06-16 22:05 by zhaojiang427
[基金申请] 博后基金,博管会会提前知道消息吗? +4 yuyiang 2024-06-13 4/200 2024-06-16 11:40 by yangyuzhong4
[基金申请] 博士后创新人才支持计划公示 +9 aishida144 2024-06-14 15/750 2024-06-16 09:52 by msjy
[基金申请] 博士后基金需要结题吗? +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[基金申请] 工材E口JQ有消息了吗 +4 babyduck 2024-06-11 4/200 2024-06-14 17:23 by firepick
[基金申请] E12面上申请 +4 汉风之遗 2024-06-13 4/200 2024-06-14 15:28 by 天外飞去来
[硕博家园] 关于硕博连读的一些疑问? +4 Lwenter 2024-06-14 4/200 2024-06-14 14:32 by ou0551
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
[硕博家园] 申博 +3 悦悦小小鱼 2024-06-12 3/150 2024-06-12 15:11 by chen5805
信息提示
请填处理意见