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

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

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

[复制链接]

该用户从未签到

尚未签到

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

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

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

×
微分方程组如下:
- j8 H& E0 ~7 W) @4 c4 o# Sdy(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;( |* K( F9 ~7 D" q
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;
, X. s5 t6 C, E3 n" n! h" a7 g, o6 k" s% {" r0 U
( U3 r1 M& j- H
两个都是y(1)、y(2)的一阶微分方程,其中第一个方程dy(1)中的I,是一个通过simulink仿真出来的一列数值解,它关于t变化(也能通过曲线拟合得到I的解析式)。y(1)、y(2)的在t=0时的初值都是0;
1 S& A. j6 S6 i% e% u. H% Q% P- c( G哪位能求出上述微分方程组的解析解或者数值解
# I' q. d1 J' v' S0 Z我用ode45函数试过,但求出的数值解全是0,不知事怎么回事。
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

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

    连续签到: 28 天

    [LV.Master]伴坛终老

    累计签到:2770 天
    连续签到:9 天
    发表于 2013-2-25 01:00:06 | 显示全部楼层
    我试着编了一个程序,看看:
    2 P0 }" b" Y* b  q) }3 lfunction dy=myfun(t,y): _  l- D! J6 b+ d$ I, x/ z1 f( ^
    " n8 e9 O. V# Q
    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;
    & W2 {) O8 Y$ Q, y# B! K+ B/ N/ Wdy(2)=(((y(1)-94.25*((y(1)/122243.631033141)^(1/15)-y(2)))/147921.11404)^(1/19)-y(2))*94.25/21.75e-6;$ }7 r9 F/ z$ _6 |$ @: h) h* @
    dy=[dy(1);dy(2)];% C* t7 f1 N5 z' B8 z
    楼主说的I这里模拟成正弦函数
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

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

    连续签到: 28 天

    [LV.Master]伴坛终老

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

    该用户从未签到

    尚未签到

     楼主| 发表于 2013-2-25 16:39:52 | 显示全部楼层
    回复 3# redplum
    $ Q, \4 F# \4 n4 M* N
    5 M" J6 n, |3 z2 {" `+ C2 a& s- J7 |: d3 R" U
        这个I是通过simulink仿真出来的一列关于t的数值,我用曲线拟合拟合为了关于变量t的9九阶多项式,还算比较理想,但不管使用原数值还是这个9阶多项式都没有结果,我待会把这个9阶多项式贴上来,到时你在看下。
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

     楼主| 发表于 2013-2-25 16:51:05 | 显示全部楼层
    回复 5# redplum
    8 k8 e6 Y% Z3 I! p5 u# O# x, Z" S2 g3 Q3 {
    这是我拟合出来的输入电流I关于变量t的9阶多项式函数,* j" y/ c9 Q/ z4 I
    p1 =  2.508e+046  - s) c" n3 m- D2 ?' B
           p2 = -3.116e+042
    1 g# E* p6 z# A# Z& n- a' F       p3 =  1.648e+038 1 d( [' d1 P9 r, I/ r9 W4 ^
           p4 = -4.824e+033 & ^) q0 c5 x, v* v' m, N) P
           p5 =  8.332e+028
    : Z9 B" X% Z9 m, ~) }       p6 = -7.701e+023
    8 K( u. Q. g6 \( v# y7 `       p7 =   3.33e+018 $ r1 w# V8 W5 S- t3 q( X* ]# |- k
           p8 = -9.423e+013
    1 _0 o; n# k: K( O6 \  o       p9 =  1.885e+009
    7 ?" ^8 i, J; S$ i* O+ K  h       p10 =     36.94 ) }$ I$ B2 [9 g$ }' G! {

    " L! o4 i; Z5 S2 f, y% x% o& o3 Gf(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# {7 b1 v; u2 O) y. N! M: F: X
    : y0 p6 e" a% {( v% `5 U/ U
    Goodness of fit:
    3 v1 \9 I% _8 H' ^  SSE: 5195
    ) w# I( P- J1 P2 T) W, R3 l1 I  R-square: 1) _6 o2 Y5 f$ f  P, r! M
      Adjusted R-square: 1
      {! |7 @( q7 e5 h  ?  RMSE: 1.399
    $ |! s4 M. Z% B8 C, V4 S/ b; |你看能不能用上解得出来。
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    愤怒
    2021-6-12 00:00
  • 签到天数: 1657 天

    连续签到: 28 天

    [LV.Master]伴坛终老

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

    连续签到: 28 天

    [LV.Master]伴坛终老

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

    该用户从未签到

    尚未签到

     楼主| 发表于 2013-2-26 21:14:28 | 显示全部楼层
    回复 8# redplum & Y- t" Z% @/ [
    1 h9 y* D7 l) }" z" y5 _
    6 W, a# l) C" m0 ?% N2 `1 r
        是的,那个x就是t,这个微分方程就是从模型中按公式推出来的, 你觉得有什么问题啊?
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

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

    GMT+8, 2025-4-4 15:56

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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