next up previous
Next: Convolution (again) Up: No Title Previous: The transfer function and

Solution of linear IC DE's by z-transform

Need for the unilateral z-transform
We have discussed the mathematics of the bilateral z-transform :

displaymath486

In our problems n represents time. We want to solve initial condition problems in which the output is specified prior to a particular time tex2html_wrap_inline836 say and calculated from tex2html_wrap_inline836 onward. For convenience we take tex2html_wrap_inline840 . There is nothing in the bilateral z-transform that makes n=0 special. We make n=0 special by defining a unilateral z-transform :

displaymath763

(subscript u to emphasize unilateral)
To use tex2html_wrap_inline848 rather than F(z) do we need to recalculate our various theorems and tables of z-transforms ? In principle yes we do . However the outputs of the causal systems we study generally have a unit step function factor ie f(n)u(n) and for such functions the values of n<0 make no contribution to the z-transform :

displaymath764

An important case where tex2html_wrap_inline848 differs from F(z) is the shifting theorem.
Shifting theorem for the unilateral z-transform

displaymath765

Substitute r=n-m so n=r+m

displaymath766

displaymath767

displaymath768

We recall that in the corresponding bilateral transform shifting theorem the first term on the right is absent. We use the unilateral transform version of the shifting theorem in the initial condition solutions of DE's. If the output is y(n) , tex2html_wrap_inline848 refers only to the solution y(n) for n>0. The first term on RHS of the theorem is where we specify the initial conditions y(-1) , y(-2) , tex2html_wrap_inline876 , y(-N).

Solution of the DE

displaymath769

Take the z-transform, applying the unilateral shifting theorem :

displaymath770

displaymath771

displaymath772

displaymath773

Collect terms in Y(z) on LHS and arrange in powers of z on RHS :

displaymath774

displaymath775

displaymath776

displaymath777

Write :

displaymath778

The term tex2html_wrap_inline884 contains the initial conditions for x and y .

displaymath779

We see that Y(z) is in the form of sums of polynomial fractions. We then do the inverse transform to find y(n) for n>0 .

Time domain and transform methods

Example solution of a DE by z-transform

displaymath780

with initial conditions y(-1)=y(-2)=1 and an input sequence x(n)=u(n) .

Taking the unilateral z-transform gives :

displaymath781

displaymath782

displaymath783

x(-1)=0 because x(n)=u(n) and y(-1)=y(-2)=1 the given initial conditions.

displaymath784

displaymath785

displaymath786

Multiply by tex2html_wrap_inline906 to generate positive powers of z :

displaymath787

displaymath788

Making a common denominator and multiplying out :

displaymath789

Make a proper fraction by dividing through by z and factorize the denominator :

displaymath790

The PFE expansion then is of the form :

displaymath791

To find tex2html_wrap_inline694 multiply by z+0.80 :

displaymath792

Evaluate this at z=-0.80 to get tex2html_wrap_inline918
Similarly tex2html_wrap_inline920 and tex2html_wrap_inline922

displaymath793

displaymath794

Take the causal function inverse transform :

displaymath795

Check the PFE with MATLAB function residue on tex2html_wrap_inline676 :

b=[0.6 0.8 -0.4]
a=[1 -0.7 -0.7 0.40]
[r,p,k]=residue(b,a)

displaymath796

displaymath797


next up previous
Next: Convolution (again) Up: No Title Previous: The transfer function and

Keith Jones
Tue Oct 27 09:51:00 EST 1998