How Do You Solve a Second Order Linear ODE with Non-Constant Coefficients?

pat666
Messages
703
Reaction score
0

Homework Statement



22e^{2t}=y''+8y'-95

Homework Equations


The Attempt at a Solution



I've been reading a textbook on this and think that I should use "method of undetermined coefficients" I know r(x)=ke^{\gamma*x} so y_p(x)=Ce^{\gamma*x}
The trouble is after reading the entire chapter I don't know what to do with that? I'm using Advanced Engineering Mathematics by Kreyszig and it is the worst textbook I have ever used there is huge gaps in methods and no easy to understand descriptions of the maths! Anyone know of any good textbooks for DE's??

Thanks for any help
 
Physics news on Phys.org
I cannot recommand you a good book, but here is a guide to the solution of your particular problem.

First you should bring all y terms on one side and the rest on the other, so you have
y''(t) + 8 y'(t) = 22 \exp(2t)+95. Now you see that there is no y(t) term, only derivatives, so you can substitute a new function z(t) = y'(t) to simplify. That's completely general up to now.

In your specific case you can guess the form of z(t) (it should have two unknown constants) and then substitute into the differential equation to determine those constants. A last integration then brings you to a specific (particular) solution y(t).

Since your differential equation is linear, you can add to this particular solution any solution of the homogeneous equation (the one where you get rid of the right hand side). Do you see why?

Try, and if you're still stuck, come back, show what you tried and ask for more :)
 
Well, how I've been taught involves finding the complementary function and the particular integral and adding them to get the complete solution. I'd follow these steps -:

1. Write in terms of the D operator. Here it would be -: (D² + 8D)y = 22exp(t) + 95exp(0).
I wrote 95 as 95exp(0) so that I can have a recognizable particular integral for it.

2. Find the auxiliary equation and its roots. The auxiliary equation here would be -: D² + 8D = 0. Surely you can find the roots by quadratic concepts.

3. Now as the roots would be real and different here(as is obvious), if the roots be A and B, the complementary function is of the form -:
c1eAx + c2eBx + ... so on as many number of roots there are. Here there are two. The c's are arbitrary constants.

4. To find the particular integral, note these inverse D operator formulae as 1/f(D) of the RHS is the particular integral -:

gif.latex?\frac{1}{D}T%20=%20\int%20Tdt.gif


gif.latex?\frac{1}{(D%20-%20a)}T%20=%20e^{aT}\int%20Te^{-aT}dt.gif


This isn't multiplication, this is the application of the inverse D operator on the RHS to find the particular integral. Note that T is a function of t(T being the RHS, then). Our particular integral here would be of the form -:

gif.latex?\frac{1}{(D^{2}%20+%208D)}(22e^{2t}%20+%2095e^{0}).gif


5. This is getting long, but we're almost there. Now to get the PI of a function like e^(t), we have to put the value of the coefficient of t from e^(something t) into D and check if the expression evaluates to 0(which it may in this question). For example, if it were e^(t), I'd put 1 into D² + 8D. Operate on e^(2t) and e^(0t) separately and add their PIs.

6. IF it does evaluate to zero (and it might here ;) ), we differentiate the expression and then put it again to check. We keep doing this until it doesn't evaluate to zero. The cases where it evaluates to zero are called cases of failure. We may or may not require those inverse D formulae I wrote, depending on the question. Just thought you should know 'em, though. They do come in handy :)

In the end,
Complete Solution = Complementary Function + Particular Integral. It isn't as hard as it seems, its just hard to explain over the internet!
Btw, this is from my college-supplied book...I can't really pinpoint you to a good book as I don't know your syllabus.
Note : Simply differentiating the expression is done specifically in the case of e^(x) functions. If you want, PM me and I'll send you some of my notes for all 4 cases.
 

Attachments

  • gif.latex?\frac{1}{D}T%20=%20\int%20Tdt.gif
    gif.latex?\frac{1}{D}T%20=%20\int%20Tdt.gif
    543 bytes · Views: 427
Last edited:
Doesn't this have a direct first integral?

<br /> 11e^{2x}=y&#039;+8y-95t<br />
This is now solveable via the integrating factor method.
 
Hey thanks for the responses:
so far I have:
roots are 0 or -8
the complementary equation is
c_1e^{0x} + c_2e^{-8x}=C_1+c_2e^{-8x}

Then 1/(d^2+8D)(22e^{2t}+95)
I am getting a bit confused here, would the coefficient you are talking about for e^(2t) be 1 or 2? none of them make D^2+8D=0 Also what does PI stand for/mean I know we've used it in these questions but I can't remember what exactly it means? also is a=A or where did a=-8D come from?

Thanks
 
For e^(2t), the coefficient would be 2. e^(2t) won't make it zero, so just put 2 into D^2 + 8D and you get the particular integral as e^(2t) * 22/20. However, in 95e^(0t), it will be zero, so differentiate till it will be non zero on putting D = 0, OR until you get a constant after differentiation. PI means particular integral.
Those inverse D formulae are not required for this case!
 
Last edited:
Note that your equation can be written
y&#039;&#039;+ 8y&#039;= 2e^{2t}+ 95
so its characteristic equation is r^2+ 8r= 0 which is easy to solve. Since r= 0 is a root of that equation, so a constant is already a solution to the associated homogeneous equation, you then seek a particular solution of the form y(t)= Ae^{2t}+ Bt

y&#039;= 2Ae^{2t}+ B and y&#039;&#039;= 4Ae^{2t}. Put those into the differential equation and solve for A and B.
 
Last edited by a moderator:
r^2+ 8r[= 0/itex] which is easy to solve. Since r= 0 is a root of that equation, so a constant is already a solution to the associated homogeneous equation, you then seek a particular solution of the form y(t)= Ae^{2t}+ Bt<br /> <br /> only posted because the original didn&#039;t show up in my browser
 
See also https://www.physicsforums.com/showpost.php?p=3206657&postcount=7

There I am talking about particular integrals and assuming known:

SVXX said:
Well, how I've been taught involves finding the complementary function and the particular integral and adding them to get the complete solution.

In the end,
Complete Solution = Complementary Function + Particular Integral. It isn't as hard as it seems, its just hard to explain over the internet!
 
Last edited:
  • #10
pat666 said:
r^2+ 8r[= 0/itex] which is easy to solve. Since r= 0 is a root of that equation, so a constant is already a solution to the associated homogeneous equation, you then seek a particular solution of the form y(t)= Ae^{2t}+ Bt<br /> <br /> only posted because the original didn&#039;t show up in my browser
<br /> Thanks. I&#039;ve fixed my post. Forgot a &quot;[&quot;!
 
  • #11
Hey guys,
I have managed to really confuse myself now!(
HallsofIvy you said to sub y' and y'' to solve for A and B which yields
4A*e^{2t} +8(2Ae^{2t}+B)=2e^{2t}+ 95 But my problem here is that i have 3 unknowns but 1 equation? Is A=0 and B=-8 from svxx not valid for some reason? I also have confused myself so much that I can't see a linear sequence of events that have led me anywhere, I seem to be back at the beginning?
Not complaining, very grateful for the help from everyone!

thanks
 
Back
Top