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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 常微分方程参数拟合问题，matlab和1stopt均可
41/1返回列表
查看: 622  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

葡萄酱油茶

捐助贵宾 (初入文坛)

爆破鬼才

[求助] 常微分方程参数拟合问题，matlab和1stopt均可

5个参数4个因变量的微分方程组，数据都有，有回复的话数据私聊，matlab和1stopt都可，谢谢！
回复此楼

» 本帖@通知

@dingd @月只蓝

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
一位不愿意透露姓名的大学生
1楼 2015-05-02 23:15:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

葡萄酱油茶

捐助贵宾 (初入文坛)

爆破鬼才
谢谢你哦，不过我之前看过了，自己动手编老是误差太大。另外，你是搞机器人的么？
赞一下 回复此楼
一位不愿意透露姓名的大学生
2楼2015-05-02 23:37:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

葡萄酱油茶

捐助贵宾 (初入文坛)

爆破鬼才
自己写的MATLAB程序算出来不符合实际情况，约束加的不够多
赞一下 回复此楼
一位不愿意透露姓名的大学生
3楼2015-05-03 19:54:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

葡萄酱油茶

捐助贵宾 (初入文坛)

爆破鬼才
S'=-aSI/N-bS;
E'=-E(t-2)+aSI/N;
I'=E(t-2）-cI-dI;
R'=bS+cI;
D'=dI
其中t是自变量，S,E,I,R,D是因变量,a,b,c,d均是大于0的参数，而且b,c,d小于1，E(t-2）代表在t-2时刻e的函数值。N(t)=S(t)+E(t)+I(t)+R(t).
用MATLAB总是得到负的参数值，用1STOPT总是不能大于三个参数。t从1到40一共40组数据。
求大神指点，感激不尽
t        S        E        I        R        D        N
1        2901        99        0        0        0        3000
2        2900        0        85        5        10        2990
3        2875        26        63        12        24        2976
4        2858        42        47        18        35        2965
5        2844        31        56        23        46        2954
6        2828        30        57        28        57        2943
7        2812        32        54        34        68        2932
8        2797        31        54        39        79        2921
9        2782        30        54        45        89        2911
10        2767        30        53        50        100        2900
11        2752        30        52        55        111        2889
12        2737        30        52        60        121        2879
13        2722        30        51        66        131        2869
14        2708        29        51        71        141        2859
15        2694        28        50        76        152        2848
16        2680        28        50        81        161        2839
17        2666        28        49        86        171        2829
18        2653        27        48        91        181        2819
19        2638        28        48        95        191        2809
20        2626        27        47        100        200        2800
21        2613        25        47        105        210        2790
22        2600        26        46        109        219        2781
23        2587        26        45        114        228        2772
24        2574        26        45        118        237        2763
25        2562        25        44        123        246        2754
26        2550        24        44        127        255        2745
27        2538        24        43        132        263        2737
28        2527        23        42        136        272        2728
29        2516        22        42        140        280        2720
30        2504        23        41        144        288        2712
31        2493        23        40        148        296        2704
32        2481        23        40        152        304        2696
33        2470        23        39        156        312        2688
34        2460        21        39        160        320        2680
35        2449        21        38        164        328        2672
36        2439        21        37        168        335        2665
37        2430        19        37        171        343        2657
38        2419        20        36        175        350        2650
赞一下 回复此楼
一位不愿意透露姓名的大学生
4楼2015-05-03 22:43:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 葡萄酱油茶 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭