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

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

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

[复制链接]

该用户从未签到

尚未签到

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

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

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

×
微分方程组如下:9 n% a; P( l6 m
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 l6 p; R; g+ w0 ?dy(2)=(((y(1)-94.25*((y(1)/122243.631033141)^(1/15)-y(2)))/147921.11404)^(1/19)-y(2))*94.25/21.75e-6;2 p* P; w/ s4 O  V$ Q0 P8 ~% A
9 Z" H$ O7 J% @6 M' f3 V5 H/ J
( p( K& t- B" k3 W' C
两个都是y(1)、y(2)的一阶微分方程,其中第一个方程dy(1)中的I,是一个通过simulink仿真出来的一列数值解,它关于t变化(也能通过曲线拟合得到I的解析式)。y(1)、y(2)的在t=0时的初值都是0;/ Z, Q& n0 {  ]$ s& [5 x1 _" l. \
哪位能求出上述微分方程组的解析解或者数值解
% ~' z3 g4 b1 p4 [, `' o我用ode45函数试过,但求出的数值解全是0,不知事怎么回事。
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

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

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到:2787 天
    连续签到:2 天
    发表于 2013-2-25 01:00:06 | 显示全部楼层
    我试着编了一个程序,看看:
    % d! _" Y4 t0 N  W% f8 v% |function dy=myfun(t,y)
    3 _9 h8 L5 S2 |0 h5 Q$ O  k5 @$ M" g8 E. q* l1 }1 q/ o
    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;7 z* b6 J) p" [" [/ k
    dy(2)=(((y(1)-94.25*((y(1)/122243.631033141)^(1/15)-y(2)))/147921.11404)^(1/19)-y(2))*94.25/21.75e-6;
    & p9 [) z$ n) n# I2 \! qdy=[dy(1);dy(2)];9 @* _2 a5 l6 n, [6 |: D; l
    楼主说的I这里模拟成正弦函数
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

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

    连续签到: 28 天

    [LV.Master]伴坛终老

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

    该用户从未签到

    尚未签到

     楼主| 发表于 2013-2-25 16:39:52 | 显示全部楼层
    回复 3# redplum
    ' g, w7 I: l3 n
    3 a# r+ E& [% ^1 v" Z: N6 A" v
    " ~. ?- w# `! ~9 J  n% \, M6 I    这个I是通过simulink仿真出来的一列关于t的数值,我用曲线拟合拟合为了关于变量t的9九阶多项式,还算比较理想,但不管使用原数值还是这个9阶多项式都没有结果,我待会把这个9阶多项式贴上来,到时你在看下。
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

     楼主| 发表于 2013-2-25 16:51:05 | 显示全部楼层
    回复 5# redplum 6 u: H8 V" ?% E/ P7 g

    ( Y4 y9 Y! [" a+ J' X这是我拟合出来的输入电流I关于变量t的9阶多项式函数,. Q/ }+ U3 G- \! o) N$ Y
    p1 =  2.508e+046  
    8 {3 M7 J8 }+ e$ @( q) x       p2 = -3.116e+042 & T5 d1 J" R! S
           p3 =  1.648e+038
    ( Y* M' a9 W8 p* }9 g       p4 = -4.824e+033 5 `4 j. z) x2 ^' l
           p5 =  8.332e+028   f; J" M) |. M8 d
           p6 = -7.701e+023
    3 D1 `: [  G: _+ J       p7 =   3.33e+018
    , R: J% O) u9 E       p8 = -9.423e+013
    ' W3 A. t8 |6 I0 |/ W       p9 =  1.885e+009 ; v7 U3 v: m; [( D
           p10 =     36.94 7 F# S! d. d3 [, y4 C' T

    9 a- o  \, B* Z9 u1 |f(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  o0 o2 x, i5 J* c+ v6 g
    7 }1 Q1 c. r' r/ P& F
    Goodness of fit:
    ' _9 t; A' X2 w5 W! \- M! Y( }  SSE: 5195
    , E+ [$ h: x: h1 k+ X5 D  R-square: 12 A* y0 g# h# P. r# B" t* ~
      Adjusted R-square: 1
    , B. X  ~# Q  o' M! x) Y  RMSE: 1.3995 M( z( {3 u; m% ^; i9 P
    你看能不能用上解得出来。
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

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

    连续签到: 28 天

    [LV.Master]伴坛终老

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

    该用户从未签到

    尚未签到

     楼主| 发表于 2013-2-26 21:14:28 | 显示全部楼层
    回复 8# redplum
    9 ?, ?1 }4 F) P# d1 C' l
    0 }4 n  y: L& `6 n1 F% A
    ) O9 }* V. W1 O. i2 p9 @    是的,那个x就是t,这个微分方程就是从模型中按公式推出来的, 你觉得有什么问题啊?
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

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

    GMT+8, 2025-4-24 22:07

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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