24小时热门版块排行榜    

查看: 2071  |  回复: 12

rewind

至尊木虫 (正式写手)

★ ★ ★ ★ ★
robert2020(金币+5):辛苦了~~ 2010-05-14 22:30:18
robert2020:要是能够提供具体程序就更好了~~ 2010-05-14 23:22:06
要是:c=a1*c0*exp(a2*t)+a3*exp(a2*t)+a4*exp(a5*t)

算法1:Levenberg-Marquardt

数                 最佳
a1                 1.04553864358777
a2                 -0.00460021436398113
a3                 2.39604974017078
a4                 -0.2967011653786
a5                 0.0486510014973765

No        c        计算值
1        25        28.2378147
2        24.6        24.5981387
3        21.9        19.8154145
4        14.3        16.1560277
5        9        10.6355689
6        6        1.7225972
7        4        -13.3255453
8        50        54.3762808
9        46        48.0043840
10        40        40.7750346
11        26        35.9899787
12        22        29.4043065
13        17        19.4833304
14        14        3.4813200
15        75        80.5147468
16        73        71.4106294
17        69.7        61.7346548
18        63.6        55.8239296
19        50.9        48.1730441
20        31        37.2440635
21        26.3        20.2881852
22        100        106.6532129
23        97.9        94.8168747
24        92.4        82.6942749
25        83.3        75.6578806
26        58.3        66.9417816
27        41.6        55.0047967
28        30        37.0950504
29        150        158.9301451
30        148        141.6293654
31        142        124.6135152
32        120.7        115.3257825
33        109.5        104.4792568
34        86.8        90.5262631
35        54        70.7087808



算法2:遗传算法

参数                 最佳

a1                 1.26455171001315
a2                 -0.0094630370285552
a3                 -2.49630302207744
a4                 -18.6678330225288
a5                 -5.18673934297126


No        实测值        计算值
1        25        10.4496567
2        24.6        23.2017133
3        21.9        18.4878404
4        14.3        16.5032411
5        9        14.7316811
6        6        13.1502912
7        4        11.7386576
8        50        42.0634495
9        46        48.3925579
10        40        38.5606818
11        26        34.4213395
12        22        30.7263398
13        17        27.4279843
14        14        24.4836947
15        75        73.6772422
16        73        73.5834024
17        69.7        58.6335232
18        63.6        52.3394378
19        50.9        46.7209986
20        31        41.7056774
21        26.3        37.2287318
22        100        105.2910350
23        97.9        98.7742469
24        92.4        78.7063646
25        83.3        70.2575362
26        58.3        62.7156573
27        41.6        55.9833704
28        30        49.9737689
29        150        168.5186205
30        148        149.1559360
31        142        118.8520474
32        120.7        106.0937329
33        109.5        94.7049748
34        86.8        84.5387566
35        54        75.4638431



[ Last edited by rewind on 2010-5-14 at 21:53 ]
在外种地...
11楼2010-05-14 21:52:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rewind

至尊木虫 (正式写手)

★ ★
robert2020(金币+2):鼓励交流经验,如果能够给出源程序就更好了。 2010-05-15 12:16:13
xiexuewang(金币+20):太感谢了 2010-05-19 10:44:15
引用回帖:
Originally posted by xiexuewang at 2010-05-12 21:08:18:
用matlab求出A1-A6的值,给出具体的步骤,先谢了

最近正学最优解 就给您做下吧

您的的拟合式:c=(a1*c0+a2)*(a3*exp(a4*t))+a5*exp(a6*t)

参数                             最优值
a1                 -1.54144909154717
a2                 -3.53252242575036
a3                 -0.678282954176262
a4                 -0.00460021436333272
a5                 -0.296701165791653
a6                 0.0486510014842707



其他拟合式见上贴...祝一切顺利
在外种地...
12楼2010-05-15 09:21:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dayong11

铜虫 (初入文坛)

robert2020:建议使用‘引用回复该贴’回复帖子,方便对方查看。不懂之处也可以PM对方直接请教。欢迎常来仿真模拟版交流。 2010-05-15 12:23:38
我有点看不懂
13楼2010-05-15 12:18:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiexuewang 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见