设为首页收藏本站|繁體中文 快速切换版块

 找回密码
 立即加入
搜索
查看: 2706|回复: 10

哪位知道怎么用Matlab求解如下simulink模型中导出的微分方程组

[复制链接]

该用户从未签到

尚未签到

发表于 2013-2-24 11:15:07 | 显示全部楼层 |阅读模式

马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!

您需要 登录 才可以下载或查看,没有账号?立即加入

×
微分方程组如下:: ]+ D, G, r# t
dy(1)=(I-((y(1)-94.25*(y(1)/122243.631033141)^(1/15)-y(2))/147921.11404)^(1/19)-(y(1)/122243.631033141)^(1/15))/6.897e-11;
. @2 N4 o/ x3 I& m# d% [2 sdy(2)=(((y(1)-94.25*((y(1)/122243.631033141)^(1/15)-y(2)))/147921.11404)^(1/19)-y(2))*94.25/21.75e-6;' l* h8 t: S: n2 @
: V# r! F  T0 x* h8 y' x

9 z2 J3 ?+ V9 V  N/ h两个都是y(1)、y(2)的一阶微分方程,其中第一个方程dy(1)中的I,是一个通过simulink仿真出来的一列数值解,它关于t变化(也能通过曲线拟合得到I的解析式)。y(1)、y(2)的在t=0时的初值都是0;
8 y& a, u% J6 ?  \: A哪位能求出上述微分方程组的解析解或者数值解/ J& y+ l, r  H4 p. K3 w4 C* ^% ^1 `
我用ode45函数试过,但求出的数值解全是0,不知事怎么回事。
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到:3057 天
    连续签到:1 天
    发表于 2013-2-24 23:38:54 | 显示全部楼层
    I与t的关系是确定的函数关系吗
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到:3057 天
    连续签到:1 天
    发表于 2013-2-25 01:00:06 | 显示全部楼层
    我试着编了一个程序,看看:/ }" s1 _  J# ?
    function dy=myfun(t,y)) N' C! \8 Z2 H- y. P5 C
    + P( A6 f4 _- J0 U& q7 a
    dy(1)=(sin(2*pi*4*t)-((y(1)-94.25*(y(1)/122243.631033141)^(1/15)-y(2))/147921.11404)^(1/19)-(y(1)/122243.631033141)^(1/15))/6.897e-11;
    + M! g: _# O! M7 r* U) {- Idy(2)=(((y(1)-94.25*((y(1)/122243.631033141)^(1/15)-y(2)))/147921.11404)^(1/19)-y(2))*94.25/21.75e-6;4 Y: \. f  @( G
    dy=[dy(1);dy(2)];: ^$ f) X8 S9 q5 S) j, _
    楼主说的I这里模拟成正弦函数
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到:3057 天
    连续签到:1 天
    发表于 2013-2-25 01:04:03 | 显示全部楼层
    不过好像执行不了,楼主说用ode45,试过用ode23了吗,另外你说的I究竟为定值还是常值
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到:3057 天
    连续签到:1 天
    发表于 2013-2-25 16:15:11 | 显示全部楼层
    我觉得是楼主的微分方程的问题,解不出来
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

     楼主| 发表于 2013-2-25 16:39:52 | 显示全部楼层
    回复 3# redplum 0 {5 @; q# T* ]! K( ]
      ]4 G% x1 l5 y8 E" |( |8 ^
    & c7 L. @3 |# A5 p* g% z* W3 d7 f) r
        这个I是通过simulink仿真出来的一列关于t的数值,我用曲线拟合拟合为了关于变量t的9九阶多项式,还算比较理想,但不管使用原数值还是这个9阶多项式都没有结果,我待会把这个9阶多项式贴上来,到时你在看下。
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

     楼主| 发表于 2013-2-25 16:51:05 | 显示全部楼层
    回复 5# redplum ) z+ a0 M# ?# i" u  s

    . C" y2 I  S" V/ q/ O这是我拟合出来的输入电流I关于变量t的9阶多项式函数,' `# L+ U, A& g" k0 d
    p1 =  2.508e+046  
    , x# a0 z6 ~% N       p2 = -3.116e+042
    " g( g7 `) E4 F       p3 =  1.648e+038 5 b% ~7 x2 W* I* D8 B
           p4 = -4.824e+033
    5 W3 O* `6 C! ]- k/ G. c+ ^$ Y       p5 =  8.332e+028 / \. J, l/ S. h
           p6 = -7.701e+023
    2 V7 D7 `7 A. z1 d0 q7 f       p7 =   3.33e+018 3 e, `! k9 L" |1 T; s+ g/ F
           p8 = -9.423e+013
    ! C& e1 x& i% D9 @; c, ^       p9 =  1.885e+009 - ^7 E& n, }# ?/ n9 [
           p10 =     36.94
    4 G+ l. c) a! H1 {1 T
    ( a2 N4 b3 z8 |$ u4 u) B6 n# Nf(x)=p1*x^9+p2*x^8+p3*x^7+p4*x^6+p5*x^5+p6*x^4+p7*x^3+p8*x^2+p9*x+p10  @! @0 u5 y: Y' ^

    2 n) L/ ?; f: P6 O$ kGoodness of fit:
    0 C6 _; ]  O  t. o  SSE: 5195
    2 b8 D# p4 Q7 p4 q( o+ w5 m& r, m  R-square: 1
    % c8 R$ q5 R; l8 f0 d2 `+ F  Adjusted R-square: 1% e% D! D; |* j
      RMSE: 1.399# L# ?9 _" x# L: C
    你看能不能用上解得出来。
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到:3057 天
    连续签到:1 天
    发表于 2013-2-26 00:19:00 | 显示全部楼层
    你那里的x就是t吧
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到:3057 天
    连续签到:1 天
    发表于 2013-2-26 01:11:54 | 显示全部楼层
    还是解不出来,我觉得还是微分方程有问题
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

     楼主| 发表于 2013-2-26 21:14:28 | 显示全部楼层
    回复 8# redplum , H0 R" T2 N& Z/ f( w

    9 H3 T5 w# @/ u7 Z5 Z  ^! m* x
    6 [! W2 X7 U; g* N6 y    是的,那个x就是t,这个微分方程就是从模型中按公式推出来的, 你觉得有什么问题啊?
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

    小黑屋|手机版|APP下载(beta)|Archiver|电力研学网 ( 赣ICP备12000811号-1|赣公网安备36040302000210号 )|网站地图

    GMT+8, 2026-4-30 15:02

    Powered by Discuz! X3.5 Licensed

    © 2001-2026 Discuz! Team.

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