找回密码
 注册
关于网站域名变更的通知
查看: 459|回复: 6
打印 上一主题 下一主题

已知含有9个微分的微分方程组,带参数,根据已知数据,拟合曲线,得参数最优解

[复制链接]
  • TA的每日心情

    2019-11-19 15:29
  • 签到天数: 1 天

    [LV.1]初来乍到

    跳转到指定楼层
    1#
    发表于 2020-12-8 11:21 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

    EDA365欢迎您登录!

    您需要 登录 才可以下载或查看,没有帐号?注册

    x
    % syms s(t) e(t) i1(t) p(t) h(t) f(t) iph(t) %定义符号变量  I% X- q  `! b  m3 _# O4 l& R
    % syms Beta IoTa Beta_2 Kappa Rho_1 Rho_2 Gamma_a Gamma_i Gamma_r Delta_i Delta_p Delta_h 希腊字母参数
      ?6 w2 c* s( ?N = 11000000/250;/ v- c7 {. j# w$ P: F; T7 W
    %dS/dt/ @0 U: O/ ?0 H# w$ t+ l* L+ F: E$ ]
    dy(1)=-Beta*(y(3)/N)*y(1)-Iota*Beta*(y(6)/N)*y(1)-Beta_2*(y(4)/N)*y(1);/ Q# _) {' j4 B- d% T2 ~% V
    %dE/dt
    + Z. {. q* U- D) @  O; M8 M4 v/ Qdy(2)=Beta*(y(3)/N)*y(1)+Iota*Beta*(y(6)/N)*y(1)+Beta_2*(y(4)/N)*y(1)-Kappa*y(2);
    0 Y/ G+ u" h5 j- V% _1 C7 x%dI/dt, t4 ^) [/ |2 H2 r' y2 _
    dy(3)=Kappa*Rho_1*y(2)-(Gamma_a+Gamma_i)*y(3)-Delta_i*y(3);: q7 H! f* c8 Y" b3 v# {  n4 t
    %dP/dt6 h/ G$ Z( b- X6 u
    dy(4)=Kappa*Rho_2*y(2)-(Gamma_a+Gamma_i)*y(4)-Delta_p*y(4);: o( n5 n5 J- y& B5 d3 v4 ~
    %dA/dt
    # H8 M+ ]3 k. b+ ndy(5)=Kappa*(1-Rho_1-Rho_2)*y(2);
    9 _5 o0 }. L/ h+ K. u% F%dH/dt7 N4 b) e  Z  }, G7 Z6 W
    dy(6)=Gamma_a*(y(3)+y(4))-Gamma_r*y(6)-Delta_h*y(6);
    $ N3 Q5 J& o8 D5 x1 J) ]* [%dR/dt
    3 ~  W: q0 y5 {dy(7)=Gamma_i*(y(3)+y(4))+Gamma_r*y(6);: c' ?2 W" L6 ^9 V2 j
    %dF/dt3 _$ X( J+ E7 N: j
    dy(8)=Delta_i*y(3)+Delta_p*y(4)+Delta_h*y(6);%
    ' x- T- Z# J) t%d(I+P+H)/dt/ a) {" y. S, B. O# M- v% F
    dy(9)=dy(3)+dy(4)+dy(6);+ I9 f$ P6 a5 N, z# i7 a
    %y(9)数据
    3 ~$ E3 b0 c, q. oy9data1 = [6,  12,  19,  25,  31,  38,  44,  60,  80,  131,  131,  259,  467,  688,  776 ...,
    . |5 t& @  M' A1 ^+ x7 M$ m+ ^7 L0 C         1776,  1460,  1739,  1984,  2101,  2590,  2827,  3233,  3892,  3697,  3151 ...,8 M: S7 z- H7 k1 U# B
             3387,  2653,  2984,  2473,  2022,  1820,  1998,  1506,  1278,  2051,  1772 ...,, e, q/ m, l) p! V4 C
             1891,  399,  894,  397,  650,  415,  518,  412,  439,  441,  435,  579,  206 ...,
    ( w& c8 J6 f9 B& @% P1 K$ J         130,  120,  143,  146,  102,  46,  45,  20,  31,  26,  11,  18,  27,  29,  39,  39];0 Q6 U/ Y3 F% `  F- Y
    %y(8)数据
    * u0 B! v6 G: e7 x+ n/ X$ Ty8data1= [0,  0,  0,  0,  0,  0,  0,  0,  4,  4,  4,  8,  15,  15,  25,  26,  26 ...,
    0 E/ X2 f- B5 `4 G( h) q3 K    38,  43,  46,  45,  57,  64,  66,  73,  73,  86,  89,  97,  108,  97,  254 ...,3 U7 Q# y1 }; V7 y+ [7 ~9 M
        121,  121,  142,  106,  106,  98,  115,  118,  109,  97,  150,  71,  52,  29 ...,4 c8 K7 P' b% p! w
        44,  37,  35,  42,  31,  38,  31,  30,  28,  27,  23,  17,  22,  11,  7,  14 ...,* ]+ [2 x. k$ O0 C$ C  v
        10,  14,  13,  13];
      G' n  y( o% C: [# U7 h( q- x1 @8 q
  • TA的每日心情
    开心
    2022-6-1 15:04
  • 签到天数: 81 天

    [LV.6]常住居民II

    2#
    发表于 2020-12-8 12:37 | 只看该作者
    不怎么了解,有待学习

    该用户从未签到

    4#
    发表于 2020-12-8 13:47 | 只看该作者
    这是让我们看你写的程序呀

    该用户从未签到

    5#
    发表于 2020-12-8 13:48 | 只看该作者
    你这个看着应该属于微分方程拟合问题吧,不知道t数据在哪儿?

    点评

    是的,t=1:1:66.我在网上收了很多相关数据,目前找到很多拟合的方法,但就对这个微分方程组难处理  详情 回复 发表于 2020-12-8 14:08
  • TA的每日心情

    2019-11-19 15:29
  • 签到天数: 1 天

    [LV.1]初来乍到

    6#
     楼主| 发表于 2020-12-8 14:08 | 只看该作者
    pulbieup 发表于 2020-12-8 13:481 ?; z0 v( _1 L
    你这个看着应该属于微分方程拟合问题吧,不知道t数据在哪儿?
    3 }' r) h# ~  O8 d0 n  W9 G8 M6 c% n
    是的,t=1:1:66.我在网上收了很多相关数据,目前找到很多拟合的方法,但就对这个微分方程组难处理" K/ K# v  y+ b! C/ y

    点评

    用1stOpt试试,应该可以解决  详情 回复 发表于 2020-12-8 14:10

    该用户从未签到

    7#
    发表于 2020-12-8 14:10 | 只看该作者
    adogpp 发表于 2020-12-8 14:08
    2 x, H% _" \4 m. b是的,t=1:1:66.我在网上收了很多相关数据,目前找到很多拟合的方法,但就对这个微分方程组难处理
    ( a$ u2 }3 K1 l6 [7 M) Q
    用1stOpt试试,应该可以解决
    4 w. Q- H$ ~* r
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-11-1 02:12 , Processed in 0.140625 second(s), 28 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表