Purdue University: ECE438 - Digital Signal Processing with Applications

ECE438 - Laboratory 3: Frequency Analysis

Introduction

In this experiment, we will use Fourier series and Fourier transforms to analyze continuoustime and discrete-time signals and systems. The Fourier representations of signals involve the decomposition of the signal in terms of complex exponential functions. These decompositions are very important in the analysis of linear time-invariant (LTI) systems, due to the property that the response of an LTI system to a complex exponential input is a complex exponential of the same frequency! Only the amplitude and phase of the input signal are changed. Therefore, studying the frequency response of an LTI system gives complete insight into its behavior. In this experiment and others to follow, we will use the Simulink extension to Matlab. Simulink is an icon-driven dynamic simulation package that allows the user to represent a system or a process by a block diagram. Once the representation is completed, Simulink may be used to digitally simulate the behavior of the continuous or discrete-time system. Simulink inputs can be Matlab variables from the workspace, or waveforms or sequences generated by Simulink itself. These Simulink-generated inputs can represent continuoustime or discrete-time sources. The behavior of the simulated system can be monitored using Simulinks version of common lab instruments, such as scopes, spectrum analyzers and network analyzers.

Questions or comments concerning this laboratory should be directed to Prof. Charles A. Bouman, School of Electrical and Computer Engineering, Purdue University, West Lafayette IN 47907; (765) 4940340; [email protected]

Purdue University: ECE438 - Digital Signal Processing with Applications

Background Exercises

INLAB REPORT: Submit these background exercises with the lab report.

2.1

Synthesis of Periodic Signals

Each signal given below represents one period of a periodic signal with period T0 . 1. Period T0 = 2. For t [0, 2]: 2. Period T0 = 1. For t [ 1 , 1 ]: 2 2 s(t) = rect(2t) For each of these two signals, do the following: i. Compute the Fourier series expansion in the form

1 s(t) = rect(t ) 2

1 2

s(t) = a0 +
k=1

Ak sin(2kf0 t + k )

where f0 = 1/T0 . Hint: You may want to use one of the following references: Sec. 4.1 of Digital Signal Processing, by Proakis and Manolakis, 1996; Sec. 4.2 of Signals and Systems, by A. Oppenheim and A. Willsky, 1983; Sec. 3.3 of Signals and Systems, A. Oppenheim and A. Willsky, 1997. Note that in the expression above, the function in the summation is sin(2kf0 t + k ), rather than a complex sinusoid. The formulas in the above references must be modied to accommodate this. You can compute the cos/sin version of the Fourier series, then convert the coecients. ii. Sketch the signal on the interval [0, T0 ].

2.2

Magnitude and Phase of Discrete-Time Systems

For the discrete-time system described by the following dierence equation, y(n) = 0.9y(n 1) + 0.3x(n) + 0.24x(n 1) i. Compute the impulse response.

Purdue University: ECE438 - Digital Signal Processing with Applications ii. Draw a system diagram.

iii. Take the Z-transform of the dierence equation using the linearity and the time shifting properties of the Z-transform. iv. Find the transfer function, dened as H(z)
Y (z) X(z)

v. Use Matlab to compute and plot the magnitude and phase responses, |H(ej )| and H(ej ), for < < . You may use Matlab commands phase and abs.

Purdue University: ECE438 - Digital Signal Processing with Applications

Full Block Library Sine Wave Simulink Blocks Scope Spectrum Analyzer Network Analyzer

Synthesizer Experiment 1

Modulator Experiment 2

CT System Analysis using a Network Analyzer Experiment 3

CT System Analysis using a Unit Step Experiment 4

DT System Analysis Experiment 5

Figure 1: Simulink utilities for lab 3.

Getting Started with Simulink

In this section, we will learn the basics of Simulink and build a simple system.

Help on Simulink Down load Lab3Utilities.zip

To get the library of Simulink functions for this laboratory, download the le Lab3Utilities.zip . Once Matlab is started, type Lab3 to bring up the library of Simulink components shown in Fig. 1. This library contains a full library of Simulink blocks, a spectrum analyzer and network analyzer designed for this laboratory, a sine wave generator, a scope, and pre-design systems for each of the experiments that you will be running.

Sine Wave

Scope

Spectrum Analyzer

Figure 2: Simulink model for the introductory example. In order to familiarize yourself with Simulink, you will rst build the system shown in Fig. 2. This system consists of a sine wave generator that feeds a scope and a spectrum analyzer.

Purdue University: ECE438 - Digital Signal Processing with Applications

1. Open a window for a new system by using the New option from the File pull-down menu, and select Model. 2. Drag the Sine Wave, Scope, and Spectrum Analyzer blocks from the Lab3 window into the new window you created. 3. Now you need to connect these three blocks. With the left mouse button, click on the output of the Sine Wave and drag it to the input of the Scope. Now use the right button to click on the line you just created, and drag to the input of the Spectrum Analyzer block. Your system should now look like Fig. 2. 4. Double click on the Scope block to make the plotting window for the scope appear. 5. Set the simulation parameters by selecting Conguration Parameters from the Simulation pull-down menu. Under the Solver tab, set the Stop time to 50, and the Max step size to 0.02. Then select OK. This will allow the Spectrum Analyzer to make a more accurate calculation. 6. Start the simulation by using the Start option from the Simulation pull-down menu. A standard Matlab gure window will pop up showing the output of the Spectrum Analyzer. 7. Change the frequency of the sine wave to 5*pi rad/sec by double clicking on the Sine Wave icon and changing the number in the Frequency eld. Restart the simulation. Observe the change in the waveform and its spectral density. If you want to change the time scaling in the plot generated by the spectrum analyzer, from the Matlab prompt use the subplot(2,1,1) and axis() commands. 8. When you are done, close the system window you created by using the Close option from the File pull-down menu.

Continuous-Time Frequency Analysis

Help on printing gures in Simulink

In this section, we will study the use and properties of the continuous-time Fourier transform with Simulink. The Simulink package is especially useful for continuous-time systems because it allows the simulation of their behavior on a digital computer.

4.1

Synthesis of Periodic Signals

Double click the icon labeled Synthesizer to bring up a model as shown in Fig. 3. This system may be used to synthesize periodic signals by adding together the harmonic components of a

Purdue University: ECE438 - Digital Signal Processing with Applications

0 Constant

Sine Wave1

Sine Wave2

Sine Wave3

Spectrum Analyzer Sine Wave4 Sum

Sine Wave5

Scope1

Sine Wave6

Sine Wave7

Mux Scope2 Mux

Figure 3: Simulink model for the synthesizer experiment. Fourier series expansion. Each Sin Wave block can be set to a specic frequency, amplitude and phase. The initial settings of the Sin Wave blocks are set to generate the Fourier series expansion 13 4 sin(2kt) . x(t) = 0 + k=1 k
k odd

These are the rst 8 terms in the Fourier series of the periodic square wave shown in Fig. 4. Run the model by selecting Start under the Simulation menu. A graph will pop up that shows the synthesized square wave signal and its spectrum. This is the output of the Spectrum Analyzer. After the simulation runs for a while, the Spectrum Analyzer element will update the plot of the spectral energy and the incoming waveform. Notice that the energy is concentrated in peaks corresponding to the individual sine waves. Print the output of the Spectrum Analyzer. You may have a closer look at the synthesized signal by double clicking on the Scope1 icon. You can also see a plot of all the individual sine waves by double clicking on the Scope2 icon.

Purdue University: ECE438 - Digital Signal Processing with Applications

Figure 4: The desired waveform for the synthesizer experiment. Synthesize the two periodic waveforms listed in section 2.1 of the background exercises. Do this by setting the frequency, amplitude, and phase of each sinewave generator to the proper values. For each case, print the output of the Spectrum Analyzer. INLAB REPORT: Hand in plots of the Spectrum Analyzer output for each of the three synthesized waveforms. For each case, comment on how the synthesized waveform diers from the desired signal, and on the structure of the spectral density.

4.2

Modulation Property

Amplitude

Repeating Sequence Product

Scope

Sine Wave Spectrum Analyzer

Figure 5: Simulink model for the modulation experiment. Double click the icon labeled Modulator to bring up a system as shown in Fig. 5. This system modulates a triangular pulse signal with a sine wave. You can control the duration and duty cycle of the triangular envelope and the frequency of the modulating sine wave. The system also contains a spectrum analyzer which plots the modulated signal and its spectrum. Generate the following signals by adjusting the Time values and Output values of the Repeating Sequence block and the Frequency of the Sine Wave. The Time values vector contains entries spanning one period of the repeating signal. The Output values vector contains the values of the repeating signal at the times specied in the Time values vector.

Purdue University: ECE438 - Digital Signal Processing with Applications

Note that the Repeating Sequence block does NOT create a discrete time signal. It creates a continuous time signal by connecting the output values with line segments. Print the output of the Spectrum Analyzer for each signal. 1. Triangular pulse duration of 1 sec; period of 2 sec; modulating frequency of 10 Hz (initial settings of the experiment). 2. Triangular pulse duration of 1 sec; period of 2 sec; modulating frequency of 15 Hz. 3. Triangular pulse duration of 1 sec; period of 3 sec; modulating frequency of 10 Hz. 4. Triangular pulse duration of 1 sec; period of 6 sec; modulating frequency of 10 Hz. Notice that the spectrum of the modulated signal consists of a comb of impulses in the frequency domain, arranged around a center frequency. INLAB REPORT: Hand in plots of the output of the Spectrum Analyzer for each signal. Answer following questions: What eect does changing the modulating frequency have on the spectral density? Why does the spectrum have a comb structure and what is the spectral distance between impulses? Why? What would happen to the spectral density if the period of the triangle pulse were to increase toward innity? (in the limit)

4.3

System Analysis

Scope

Network Analyzer

Analog Butterworth LP Filter

Figure 6: Simulink model for the continuous-time system analysis experiment using a network analyzer. Double click the icon labeled CT System Analysis using a Network Analyzer to bring up a system as shown in Fig. 6. This system includes a Network Analyzer model for measuring the frequency response of a system. The Network Analyzer works by generating a weighted chirp signal (shown on the Scope) as an input to the system-under-test. The analyzer measures

Purdue University: ECE438 - Digital Signal Processing with Applications

the frequency response of the input and output of the system and computes the transfer function. By computing the inverse Fourier transform, it then computes the impulse response of the system. Use this setup to compute the frequency and impulse response of the given fourth-order Butterworth lter with a cut-o frequency of 1Hz. Print the gure showing the magnitude response, the phase response and the impulse response of the system. To use the tall mode to obtain a larger printout, type orient(tall); directly before you print.

du/dt Step Analog Butterworth LP Filter Derivative

Scope

Figure 7: Simulink model for the continuous-time system analysis experiment using a unit step. An alternative method for computing the impulse response is to input a step function into the system and then to compute the derivative of the output. The model for doing this is given in the CT System Analysis using a Unit Step block. Double click on this icon and compute the impulse response of the lter using this setup (Fig. 7). Make sure that the characteristics of the lter are the same as in the previous setup. After running the simulation, print the graph of the impulse response from the scope. INLAB REPORT: Hand in the printout of the output of the Network Analyzer (magnitude and phase of the frequency response, and the impulse response) and the plot of the impulse response obtained using a unit step. What are the advantages and disadvantages of each method?

Discrete-Time Frequency Analysis

In this section of the laboratory, we will study the use of the discrete-time Fourier transform.

5.1

Discrete-Time Fourier Transform

The DTFT (Discrete-Time Fourier Transform) is the Fourier representation used for nite energy discrete-time signals. For a discrete-time signal, x(n), we denote the DTFT as the function X(ej ) given by the expression

X(ej ) =
n=

x(n)ejn .

Purdue University: ECE438 - Digital Signal Processing with Applications

10

Since X(ej ) is a periodic function of with a period of 2, we need only to compute X(ej ) for < < . Write a Matlab function X = DTFT(x,n0,dw) that computes the DTFT of the discretetime signal x. Here n0 is the time index corresponding to the 1st element of the x vector, and dw is the spacing between the samples of the Matlab vector X. For example, if x is a vector of length N , then its DTFT is computed by
N

X(w) =
n=1

x(n)ejw(n+n01)

where w is a vector of values formed by w=(-pi:dw:pi). Hint: In Matlab, j or i is dened as 1. However, you may also compute this value using the Matlab expression i = sqrt(-1). For the following signals use your DTFT function to i. Compute X(ej ) ii. Plot the magnitude and the phase of X(ej ) in a single plot using the subplot command. Hint: Use the abs() and angle() commands. 1. x(n) = (n) 2. x(n) = (n 5) 3. x(n) = (0.5)n u(n)

INLAB REPORT: Hand in a printout of your Matlab function. Also hand in plots of the DTFTs magnitude and phase for each of the three signals.

5.2

System Analysis

Help on printing Simulink system windows

Double click the icon labeled DT System Analysis to bring up an incomplete block diagram as shown in Fig. 8. It is for a model that takes a discrete-time sine signal, processes it according to a dierence equation and plots the multiplexed input and output signals in a graph window. Complete this block diagram such that it implements the following dierence equation given in section 2.2 of the background exercises. y(n) = 0.9y(n 1) + 0.3x(n) + 0.24x(n 1)

Purdue University: ECE438 - Digital Signal Processing with Applications

11

Mux Sine Wave Mux Scope

1/z Unit Delay1

1 Gain1 Sum1

1/z Unit Delay 2

1 Gain2 Sum 2

1 Gain3

Figure 8: Incomplete Simulink setup for the discrete-time system analysis experiment. You are provided with the framework of the setup and the building blocks that you will need. You can change the values of the Gain blocks by double clicking on them. After you complete the setup, adjust the frequency of Sine Wave to the following frequencies: = /16, = /8, and = /4. For each frequency, make magnitude response measurements using the input and output sequences shown in the graph window. Compare your measurements with the values of the magnitude response |H(ej )| which you computed in the background exercises at these frequencies. An alternative way of nding the frequency response is taking the DTFT of the impulse response. Use your DTFT function to nd the frequency response of this system from its impulse response. The impulse response was calculated in section 2.2 of the background exercises. Plot the impulse response, and the magnitude and phase of the frequency response in the same gure using the subplot command. INLAB REPORT: Hand in the following: Printout of your completed block diagram Table of both the amplitude measurements you made and their theoretical values. Printout of the gure with the impulse response, and the magnitude and phase of the frequency response.