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

zhuce

微分方程组如下：
: c( v0 O- F# ~' Xdy(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;
/ u" G* R+ a( S$ @4 @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;

( p" {2 k- n* }. ]
- d& D% z8 v) X6 M! ~两个都是y（1）、y（2）的一阶微分方程，其中第一个方程dy(1)中的I，是一个通过simulink仿真出来的一列数值解，它关于t变化（也能通过曲线拟合得到I的解析式）。y（1）、y（2）的在t=0时的初值都是0；
哪位能求出上述微分方程组的解析解或者数值解
我用ode45函数试过，但求出的数值解全是0，不知事怎么回事。

redplum

I与t的关系是确定的函数关系吗

redplum

我试着编了一个程序，看看：
function dy=myfun(t,y)
/ Q1 e. r) ]& D0 x0 b. w
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;
" o: l; M- k0 \! kdy(2)=(((y(1)-94.25*((y(1)/122243.631033141)^(1/15)-y(2)))/147921.11404)^(1/19)-y(2))*94.25/21.75e-6;
8 W" ?* W. p. B: r( o' @  Udy=[dy(1);dy(2)];
楼主说的I这里模拟成正弦函数

redplum

不过好像执行不了，楼主说用ode45,试过用ode23了吗，另外你说的I究竟为定值还是常值

redplum

我觉得是楼主的微分方程的问题，解不出来

zhuce

回复 3# redplum

' E; y. \/ W- H( t3 P. W1 w" A
! r5 N6 X! L5 Q/ F* i) t1 P7 \9 n+ g    这个I是通过simulink仿真出来的一列关于t的数值，我用曲线拟合拟合为了关于变量t的9九阶多项式，还算比较理想，但不管使用原数值还是这个9阶多项式都没有结果，我待会把这个9阶多项式贴上来，到时你在看下。

zhuce

回复 5# redplum

9 g& F$ x. w& V$ n4 E9 d$ b! W& J这是我拟合出来的输入电流I关于变量t的9阶多项式函数，
p1 =  2.508e+046  
       p2 = -3.116e+042
       p3 =  1.648e+038
       p4 = -4.824e+033
       p5 =  8.332e+028
       p6 = -7.701e+023
- g# Z' I9 g; n$ Y- L       p7 =   3.33e+018
. H6 J/ j: H0 M3 k       p8 = -9.423e+013
3 w( V: S6 k+ |3 i' V2 H6 w       p9 =  1.885e+009
# M' g. `" ~( [. s0 L7 Q       p10 =     36.94

6 H/ K0 w" y7 ]% bf(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

) O& E9 m; @" J- ^/ |4 [8 k" qGoodness of fit:
5 T$ K9 b' H5 i. D* u4 L" n  SSE: 5195
8 @8 b8 s0 @: f6 x# [7 o- j( ]6 y  R-square: 1
' @) T# I3 X& k. N) d. c$ |0 G  Adjusted R-square: 1
' {, C3 N' W/ S+ s$ Q  RMSE: 1.399
/ R: d: `5 c4 E1 L. t' }你看能不能用上解得出来。

redplum

你那里的x就是t吧

redplum

还是解不出来，我觉得还是微分方程有问题

zhuce

回复 8# redplum
- k8 V6 u) W2 h( a. t
2 d0 W% W, Q* \+ m& ~. F
: ~: L* ]. J: Z# Z4 g& }( z5 t    是的，那个x就是t，这个微分方程就是从模型中按公式推出来的， 你觉得有什么问题啊？