<< Chapter < Page Chapter >> Page >

The following two continuous-time systems are commonlyused in electrical engineering:

differentiator: y ( t ) = d d t x ( t ) integrator: y ( t ) = - t x ( τ ) d τ

For each of these two systems, do the following:

  • Formulate a discrete-time system that approximates the continuous-time function.
  • Write down the difference equation that describes your discrete-time system. Your difference equation should be inclosed form, i.e. no summations.
  • Draw a block diagram of your discrete-time system as in [link] .

Stock market example

One reason that digital signal processing (DSP) techniques are so powerful is thatthey can be used for very different kinds of signals. While most continuous-time systems only process voltage and current signals,a computer can process discrete-time signals which are essentially just sequences of numbers.Therefore DSP may be used in a very wide range of applications. Let's look at an example.

A stockbroker wants to see whether the average value of a certain stock is increasing or decreasing.To do this, the daily fluctuations of the stock values must be eliminated.A popular business magazine recommends three possible methods for computing this average.

avgvalue(today) = 1 3 value(today)+value(yesterday)+value(2 days ago)
avgvalue(today) = 0 . 8 * avgvalue(yesterday) + 0 . 2 * (value(today))
avgvalue(today) = avgvalue(yesterday) + 1 3 (value(today)-(value(3 days ago))

Do the following:

  • For each of the these three methods: 1) write a difference equation,2) draw a system diagram, and 3) calculate the impulse response.
  • Explain why methods [link] and [link] are known as moving averages.

Example discrete-time systems

Write two Matlab functions that will apply the differentiator and integrator systems, designed in the "Example Discrete-time Systems" section, to arbitrary input signals. Then apply the differentiator and integrator to the followingtwo signals for - 10 n 20 .

  • δ ( n ) - δ ( n - 5 )
  • u ( n ) - u ( n - ( N + 1 ) ) with N = 10

Hint: To compute the function u ( n ) for - 10 n 20 , first set n = -10:20 , and then use the Boolean expression u = (n>=0) .

For each of the four cases, use the subplot and stem commands to plot each input and output signal on a single figure.

Submit printouts of your Matlab code and hardcopies containing the input and output signals.Discuss the stability of these systems.

Difference equations

In this section, we will study the effect of two discrete-time filters.The first filter, y = S 1 [ x ] , obeys the difference equation

y ( n ) = x ( n ) - x ( n - 1 )

and the second filter, y = S 2 [ x ] , obeys the difference equation

y ( n ) = 1 2 y ( n - 1 ) + x ( n )

Write Matlab functions to implement each of these filters.

In Matlab, when implementing a difference equation using a loop structure, it is very good practice to pre-define your output vector beforeentering into the loop. Otherwise, Matlab has to resize the output vector at each iteration.For example, say you are using a FOR loop to filter the signal x ( n ) , yielding an output y ( n ) . You can pre-define the output vector by issuing the command y=zeros(1,N) before entering the loop, where N is the final length of y . For long signals, this speeds up the computation dramatically.

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Purdue digital signal processing labs (ece 438). OpenStax CNX. Sep 14, 2009 Download for free at http://cnx.org/content/col10593/1.4
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Purdue digital signal processing labs (ece 438)' conversation and receive update notifications?

Ask