<< Chapter < Page Chapter >> Page >

Solving a lccde

In order for a linear constant-coefficient difference equation to be useful in analyzing a LTI system, we must be able tofind the systems output based upon a known input, x t , and a set of initial conditions. Two common methods exist for solving a LCCDE: the direct method and the indirect method , the latter being based on the Laplace-transform. Below we will briefly discussthe formulas for solving a LCCDE using each of these methods.

Direct method

The final solution to the output based on the direct method is the sum of two parts, expressed in the followingequation:

y t y h t y p t
The first part, y h t , is referred to as the homogeneous solution and the second part, y h t , is referred to as particular solution . The following method is very similar to that used to solve many differential equations, so if youhave taken a differential calculus course or used differential equations before then this should seem veryfamiliar.

Homogeneous solution

We begin by assuming that the input is zero, x t 0 .Now we simply need to solve the homogeneous differential equation:

k 0 N a k y t k 0
In order to solve this, we will make the assumption that the solution is in the form of an exponential. We willuse lambda, λ , to represent our exponential terms. We now have to solve thefollowing equation:
k 0 N a k λ t k 0
We can expand this equation out and factor out all of thelambda terms. This will give us a large polynomial in parenthesis, which is referred to as the characteristic polynomial . The roots of this polynomial will be the key to solving the homogeneousequation. If there are all distinct roots, then the general solution to the equation will be as follows:
y h t C 1 λ 1 t C 2 λ 2 t C N λ N t
However, if the characteristic equation contains multiple roots then the above general solution will be slightlydifferent. Below we have the modified version for an equation where λ 1 has K multiple roots:
y h t C 1 λ 1 t C 1 t λ 1 t C 1 t 2 λ 1 t C 1 t K 1 λ 1 t C 2 λ 2 t C N λ N t

Particular solution

The particular solution, y p t , will be any solution that will solve the general differential equation:

k 0 N a k y p t k k 0 M b k x t k
In order to solve, our guess for the solution to y p t will take on the form of the input, x t . After guessing at a solution to the above equation involving the particular solution, one onlyneeds to plug the solution into the differential equation and solve it out.

Indirect method

The indirect method utilizes the relationship between the differential equation and the Laplace-transform, discussed earlier , to find a solution. The basic idea is to convert the differentialequation into a Laplace-transform, as described above , to get the resulting output, Y s . Then by inverse transforming this and using partial-fractionexpansion, we can arrive at the solution.

L d d t y ( t ) = s Y ( s ) - y ( 0 )

This can be interatively extended to an arbitrary order derivative as in Equation [link] .

L d n d t n y ( t ) = s n Y ( s ) - m = 0 n - 1 s n - m - 1 y ( m ) ( 0 )

Now, the Laplace transform of each side of the differential equation can be taken

L k = 0 n a k d k d t k y ( t ) = L x ( t )

which by linearity results in

k = 0 n a k L d k d t k y ( t ) = L x ( t )

and by differentiation properties in

k = 0 n a k s k L y ( t ) - m = 0 k - 1 s k - m - 1 y ( m ) ( 0 ) = L x ( t ) .

Rearranging terms to isolate the Laplace transform of the output,

L y ( t ) = L x ( t ) + k = 0 n m = 0 k - 1 a k s k - m - 1 y ( m ) ( 0 ) k = 0 n a k s k .

Thus, it is found that

Y ( s ) = X ( s ) + k = 0 n m = 0 k - 1 a k s k - m - 1 y ( m ) ( 0 ) k = 0 n a k s k .

In order to find the output, it only remains to find the Laplace transform X ( s ) of the input, substitute the initial conditions, and compute the inverse Laplace transform of the result. Partial fraction expansions are often required for this last step. This may sound daunting while looking at Equation [link] , but it is often easy in practice, especially for low order differential equations. Equation [link] can also be used to determine the transfer function and frequency response.

As an example, consider the differential equation

d 2 d t 2 y ( t ) + 4 d d t y ( t ) + 3 y ( t ) = cos ( t )

with the initial conditions y ' ( 0 ) = 1 and y ( 0 ) = 0 Using the method described above, the Laplace transform of the solution y ( t ) is given by

Y ( s ) = s ( s 2 + 1 ) ( s + 1 ) ( s + 3 ) + 1 ( s + 1 ) ( s + 3 ) .

Performing a partial fraction decomposition, this also equals

Y ( s ) = . 25 1 s + 1 - . 35 1 s + 3 + . 1 s s 2 + 1 + . 2 1 s 2 + 1 .

Computing the inverse Laplace transform,

y ( t ) = ( . 25 e - t - . 35 e - 3 t + . 1 cos ( t ) + . 2 sin ( t ) ) u ( t ) .

One can check that this satisfies that this satisfies both the differential equation and the initial conditions.

Summary

One of the most important concepts of DSP is to be able to properly represent the input/output relationship to a given LTIsystem. A linear constant-coefficient difference equation (LCCDE) serves as a way to express just this relationship in a discrete-time system. Writing the sequenceof inputs and outputs, which represent the characteristics of the LTI system, as a difference equation helps in understandingand manipulating a system.

Questions & Answers

differentiate between demand and supply giving examples
Lambiv Reply
differentiated between demand and supply using examples
Lambiv
what is labour ?
Lambiv
how will I do?
Venny Reply
how is the graph works?I don't fully understand
Rezat Reply
information
Eliyee
devaluation
Eliyee
t
WARKISA
hi guys good evening to all
Lambiv
multiple choice question
Aster Reply
appreciation
Eliyee
explain perfect market
Lindiwe Reply
In economics, a perfect market refers to a theoretical construct where all participants have perfect information, goods are homogenous, there are no barriers to entry or exit, and prices are determined solely by supply and demand. It's an idealized model used for analysis,
Ezea
What is ceteris paribus?
Shukri Reply
other things being equal
AI-Robot
When MP₁ becomes negative, TP start to decline. Extuples Suppose that the short-run production function of certain cut-flower firm is given by: Q=4KL-0.6K2 - 0.112 • Where is quantity of cut flower produced, I is labour input and K is fixed capital input (K-5). Determine the average product of lab
Kelo
Extuples Suppose that the short-run production function of certain cut-flower firm is given by: Q=4KL-0.6K2 - 0.112 • Where is quantity of cut flower produced, I is labour input and K is fixed capital input (K-5). Determine the average product of labour (APL) and marginal product of labour (MPL)
Kelo
yes,thank you
Shukri
Can I ask you other question?
Shukri
what is monopoly mean?
Habtamu Reply
What is different between quantity demand and demand?
Shukri Reply
Quantity demanded refers to the specific amount of a good or service that consumers are willing and able to purchase at a give price and within a specific time period. Demand, on the other hand, is a broader concept that encompasses the entire relationship between price and quantity demanded
Ezea
ok
Shukri
how do you save a country economic situation when it's falling apart
Lilia Reply
what is the difference between economic growth and development
Fiker Reply
Economic growth as an increase in the production and consumption of goods and services within an economy.but Economic development as a broader concept that encompasses not only economic growth but also social & human well being.
Shukri
production function means
Jabir
What do you think is more important to focus on when considering inequality ?
Abdisa Reply
any question about economics?
Awais Reply
sir...I just want to ask one question... Define the term contract curve? if you are free please help me to find this answer 🙏
Asui
it is a curve that we get after connecting the pareto optimal combinations of two consumers after their mutually beneficial trade offs
Awais
thank you so much 👍 sir
Asui
In economics, the contract curve refers to the set of points in an Edgeworth box diagram where both parties involved in a trade cannot be made better off without making one of them worse off. It represents the Pareto efficient allocations of goods between two individuals or entities, where neither p
Cornelius
In economics, the contract curve refers to the set of points in an Edgeworth box diagram where both parties involved in a trade cannot be made better off without making one of them worse off. It represents the Pareto efficient allocations of goods between two individuals or entities,
Cornelius
Suppose a consumer consuming two commodities X and Y has The following utility function u=X0.4 Y0.6. If the price of the X and Y are 2 and 3 respectively and income Constraint is birr 50. A,Calculate quantities of x and y which maximize utility. B,Calculate value of Lagrange multiplier. C,Calculate quantities of X and Y consumed with a given price. D,alculate optimum level of output .
Feyisa Reply
Answer
Feyisa
c
Jabir
the market for lemon has 10 potential consumers, each having an individual demand curve p=101-10Qi, where p is price in dollar's per cup and Qi is the number of cups demanded per week by the i th consumer.Find the market demand curve using algebra. Draw an individual demand curve and the market dema
Gsbwnw Reply
suppose the production function is given by ( L, K)=L¼K¾.assuming capital is fixed find APL and MPL. consider the following short run production function:Q=6L²-0.4L³ a) find the value of L that maximizes output b)find the value of L that maximizes marginal product
Abdureman
types of unemployment
Yomi Reply
What is the difference between perfect competition and monopolistic competition?
Mohammed
Got questions? Join the online conversation and get instant answers!
Jobilize.com Reply

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Signals and systems. OpenStax CNX. Aug 14, 2014 Download for free at http://legacy.cnx.org/content/col10064/1.15
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Signals and systems' conversation and receive update notifications?

Ask