Prev: scaling in negative axis
Next: plotyy and while loop= :(
From: Dave on 1 Dec 2009 01:13
I want to solve an ODE in matlab or simulink, its a 3-element Windkessel model and I have already derived the ode:



(1+R1/R2)*(Q(t)+CR1*Q'(t)=(P(t)/R2)+C(P'(t))

I have all the constants (R1, R2, and C) and Also I have an input of P(t) (its a square wave input)

My question is how can I manipulate equation to make it solvable by ode23? I can't isolate what I'm solving for (q(t)) because I have its derivative present also. Can someone please help? The same question goes for simulink.. can I do it in there also?

Thanks in advance.

Dave
From: Torsten Hennig on 30 Nov 2009 17:08
> I want to solve an ODE in matlab or simulink, its a
> 3-element Windkessel model and I have already derived
> the ode:
>
>
>
> (1+R1/R2)*(Q(t)+CR1*Q'(t)=(P(t)/R2)+C(P'(t))
>
> I have all the constants (R1, R2, and C) and Also I
> have an input of P(t) (its a square wave input)
>
> My question is how can I manipulate equation to make
> it solvable by ode23? I can't isolate what I'm
> solving for (q(t)) because I have its derivative
> present also. Can someone please help? The same
> question goes for simulink.. can I do it in there
> also?
>
> Thanks in advance.
>
> Dave

So P(t) (and thus P'(t) (obtainable by a numerical
approximation)) are known ?
Then you can isolate Q'(t) and you will get a
standard ordinary differential equation to solve for
Q(t).

Best wishes
Torsten.
 | 
Pages: 1
Prev: scaling in negative axis
Next: plotyy and while loop= :(