24小时热门版块排行榜    

查看: 2078  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xiexuewang

金虫 (小有名气)

[交流] 【求助】二元函数拟合 已有6人参与

t,c0,c每组35个数据
t=[0 24 48 60 72 84 96 0 24 48 60 72 84 96 0 24 48 60 72 84 96 0 24 48 60 72 84 96 0 24 48 60 72 84 96];
c0=[25 25 25 25 25 25 25 50 50 50 50 50 50 50 75 75 75 75 75 75 75 100 100 100 100 100 100 100 150 150 150 150 150 150 150];
c=[25 24.6 21.9 14.3 9 6 4 50 46 40 26 22 17 14 75 73 69.7 63.6 50.9 31 26.3 100 97.9 92.4 83.3 58.3 41.6 30 150 148 142 120.7 109.5 86.8 54];

拟合函数形式:c=(A(1)*c0+A(2))*(A(3)*exp(A(4)*t)+A(5)*exp(A(6)*t)
求解A1-A6的值

[ Last edited by xiexuewang on 2010-5-13 at 14:05 ]
回复此楼
没有不劳而获,任何成就都是汗水浇灌的结果
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖
查看全部 13 个回答

xiaojun376

金虫 (小有名气)


zzuwangshilei(金币+1):帮助他人解决问题 2010-05-12 19:25:27
你这好像不是拟合,是参数估计,你是想求A的值是吧,那就用lsqnonlin寻找最优值
在交流中共同进步!
2楼2010-05-12 18:05:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiexuewang

金虫 (小有名气)

是啊,我是想求A的值,就是不会用lsqnonlin函数啊,
没有不劳而获,任何成就都是汗水浇灌的结果
3楼2010-05-12 21:04:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiexuewang

金虫 (小有名气)

zzuwangshilei:具体的步骤,一般要自己做的,别人给出的是建议或思路 2010-05-13 15:30:45
用matlab求出A1-A6的值,给出具体的步骤,先谢了
没有不劳而获,任何成就都是汗水浇灌的结果
4楼2010-05-12 21:08:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见