ODE- How to solve this equation??

Question

inbal kashany 2011-6-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/8582-ode-how-to-solve-this-equation

Hi Im trying to solve this differential equation:

dydt=((g0.*L-y)./tau)-((Pin./Esat)*(exp(y)-1));

by ODE45: [T Y]=ode45(@(t,y) fun(t,Y,It,I),[0 3],0);

but it doesnt work... I guess becouse it based on example3 from ODE45 help which looks like this: dydt = -f.*y + g,(y'(t) + f(t)y(t) = g(t))

but here I have the "exp(y)" which I dont know how to hadle with.

g0=R*a*n0*((I/I0)-1); I is a vector dependent on t. L,tau,Pin,Esat,R,a,n0,I0 are constant.

any suggestion may help. thanks

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Arnaud Miege 2011-6-1

Please reformat your code. See http://www.mathworks.co.uk/matlabcentral/answers/6200-tutorial-how-to-ask-a-question-on-answers-and-get-a-fast-answer

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Igor 2011-6-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/8582-ode-how-to-solve-this-equation#answer_11813

be exactly in y(t) definition and remember that y(t) is a COLUMN vector.

also in -- @(t,y) fun(t,Y,It,I) -- you use y and Y!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Laura Proctor 2011-6-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/8582-ode-how-to-solve-this-equation#answer_11822

在 MATLAB Online 中打开

The following code worked for me:

dydt=@(t,y)(((g0*L-y)/tau)-((Pin/Esat)*(exp(y)-1)));
[T Y]=ode45(dydt,[0 3],0);

But, I did have to make some guesses at the constant values, such as g0, L, tau, Pin, and Esat.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

inbal kashany 2011-6-2

g0 is not constant-

g0=R*a*n0*((I/I0)-1);

I=[zeros(10,1);ones(10,1);zeros(10,1);]; %step signal%%

It=linspace(0,30,30); %generate t for I

I=interp1(It,I,t); % Interpolate the data set (It,I) at time t

now I got the error: odearguments at 115 @(T,H)(((G0*L-H)/TAU)-((PIN/ESAT)*(EXP(H)-1))) returns a vector of length 30, but the length of initial conditions vector is 1. The vector returned by @(T,H)(((G0*L-H)/TAU)-((PIN/ESAT)*(EXP(H)-1))) and the initial conditions vector must have the same number of elements.

*h replace y

I'll be happy get an answer.

请先登录，再进行评论。

Answer 3

inbal kashany 2011-6-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/8582-ode-how-to-solve-this-equation#answer_11890

Thanks to both of you for your replay.

g0 is not constant-

g0=R*a*n0*((I/I0)-1);

I=[zeros(10,1);ones(10,1);zeros(10,1);]; %step signal%%

It=linspace(0,30,30); %generate t for I

I=interp1(It,I,t); % Interpolate the data set (It,I) at time t

now I got the error: odearguments at 115 @(T,H)(((G0*L-H)/TAU)-((PIN/ESAT)*(EXP(H)-1))) returns a vector of length 30, but the length of initial conditions vector is 1. The vector returned by @(T,H)(((G0*L-H)/TAU)-((PIN/ESAT)*(EXP(H)-1))) and the initial conditions vector must have the same number of elements.

*h replace y

I'll be happy get an answer.