Re: how to solve for periodic solution?

*To*: mathgroup at smc.vnet.net*Subject*: [mg35738] Re: how to solve for periodic solution?*From*: Selwyn Hollis <slhollis at earthlink.net>*Date*: Sun, 28 Jul 2002 03:32:13 -0400 (EDT)*References*: <ahttq2$qjo$1@smc.vnet.net>*Sender*: owner-wri-mathgroup at wolfram.com

It's hard to know for sure what you're expected to do here, so I'll have to guess. First off, I'm going to assume that the problem you're attempting to describe is this: A'(t) = S(t) - k*A(t) A(0) = A0 where k= ln(2)/halflife and S is the period-P extension of a "pulse" of the form q*(UnitStep[t-t1] - UnitStep[t-t2]). Here 0 <= t1 < t2 < t1 + P, and q is the height of the pulse. I'm also guessing that you have numerical values for q, P, A0, and k. The question then is: How wide should the pulse be in order for the solution to be periodic. Here's an approach to get you started toward an "experimental" solution: Make a rough guess at t1 and t2. Use those to define S[t_]:= q*Sum[(UnitStep[t-t1-n*P] - UnitStep[t-t2-n*P]), {n,0,4}] which will define S up to t = 5P. Now compute the solution of the ODE like so: A[t_] = Simplify[E^(-k t)(A0 + Integrate[E^(k x)S[x], {x, 0, t}])] Plot the graph of A[t] over {t,0,5*P}. Then adjust your t1 and t2 until the solution looks periodic. Here's another (better) approach: Up to time t=t1 you have A(t)=A0*e^(-k t). Then for t1 < t < t2 the solution is A1[t_]= A[t]/.First[DSolve[{A'[t] == q - k*A[t], A[t1]==A0*E^(-k t1)}, A[t], t]] Then for t2 < t < t1 + P the solution is A2[t_]= A[t]/.First[DSolve[{A'[t] == - k*A[t], A[t2]==A1[t2]}, A[t], t]] The solution will have period P if A2[P]==A0. Solve that equation for t2. You'll find that t2 can be put in the form of t1 + pulse width. --- Selwyn Hollis slhollis at mac.com http://www.math.armstrong.edu/faculty/hollis tom wrote: > Hi, > > I have no clue as to how to proceed to solve the equation below. > Any ideas on how to solve for the periodic solution of this equation > would be most appreciated!! Thankyou!! > > A(t)=S(t+P)-(ln(2)/halflife)*A(t) > > S(t+P) = S(t) if t1<t<t2 > S(t+P) = 0 otherwise > > where P is the period of the system > t is the time at which a sample is taken > t1 is the beginning of a discontinuity > t2 is the end of a discontinuity > A(t) is the abundance at time t > > I want to fit data to the resulting equation and find t1 and t2 > > Sincerely, > > Tom >