EE513 Audio Signals and Systems

1. EE513Audio Signals and Systems Digital Signal Processing (Systems) Kevin D. DonohueElectrical and Computer EngineeringUniversity of Kentucky

2. Laplace to begin with The Laplace transform is used to characterize and analyze signal and system interactions. If x(t) is the time domain signal, its Laplace transform is defined as: where s belongs to the set of complex numbers over which the integral converges. What is the Laplace transform of a system impulse response? What does the Laplace transform become if it is evaluated along the imaginary axis?

3. Z-Transform If the Laplace transform is made discrete by sampling the time axis with interval Ts, it becomes: Now let:

4. Z-Transform Substitute z to obtain: Normalize sampling rate and define Z-transform as: where z belongs to the set of complex numbers for which summation converges. Describe mapping from S-plane to Z-plane.

5. j Axis in s to z Mapping j axis in s corresponds to a zero real part,  = 0 for all imaginary values, . In z this results in: S-plane Z-plane IM IM RE RE

6. Z-plane S-plane Relationships Imaginary axis (j) in S-plane, maps into the unit circle in Z-plane, where segments on the j axis described by: map into a complete unit circle for every integer k. Since the mapping from S to Z is not unique, let the range for k=0 be referred to as the root frequency range and all other values of k with the aliased frequency ranges.

7. Negative Real Axis s to z Mapping Negative real axis in s corresponds to  < 0. In z this results in: S-plane Z-plane IM IM RE RE

8. Z-plane S-plane Relationships The negative real axis of the S-plane ( < 0 ) maps into the area inside the unit circle of Z (| z | < 1) . Therefore the stable region of the S-plane (left-half plane) corresponds to the area inside the unit circle of the Z-plane. Similarly, the unstable region of the S-plane (right-half plane) corresponds to the area outside the unit circle of the Z-plane. Because of aliasing from the sampling and scaling from the exponential transformation, there is no simple (linear) scaling between the Z and S planes. A warping or distortion occurs when matching domain points: Warping Possible Aliasing

9. Z-Transform One-Sided Many applications assume the input starts at t = 0 (n=0 for discrete) and no response exists before t = 0. So the Z-transform is often written as: Examples: Find the z-transforms of x[n] = u[n] and x[n] = anu[n]; Assuming the z-transform of x[n] is X(z), find the Z-transform of x[n-k] for k>0.

10. Homework(1) Use definition Find the z-Transforms of: a) b) c) for k > 0. Use definition Use ZT properties (delay)

11. Convolution Given the impulse response of a discrete linear system h[n] the input-output relationship is described by discrete convolution: For x[n] and h=[n] below, graphically demonstrate their convolution.

27. Sinusoidal Response Consider sinusoidal input: Note this input is always on (steady-state). Show that for impulse response h(n), the convolution sum for evaluating the output becomes: Important Concepts: • The response to a sinusoidal input is a phasor multiplication between input phasor and transfer function value at the excitation frequency. • The frequency response (Transfer Function) of a discrete system is the z-transform of its impulse response evaluated on the unit circle! • Convolution in time domain is equivalent to multiplication in the frequency domain. Input Complex Coefficient

28. Sinusoidal Response Example For sinusoidal input: And system described by: Derive an expression for and plot the frequency response (phase and magnitude) of the system output.

29. Homework(2) For sinusoidal inputs: And system described by: Plot the frequency response (phase and magnitude) and compute the corresponding outputs.

30. Sampling Sampling rate determines the highest signal frequency that can be reconstructed from the signal samples without error. At least 3 samples (2 complete sampling intervals) must fall within a period for digitization without aliasing. In other words the sampling rate must be greater than twice the highest signal frequency for a band limited signal. Fs =200 Hz Ts = 5 ms

31. Sampling Aliasing I – The Movie (FS=200, Range 50-150 Hz) Run following mfile at: http://www.engr.uky.edu/~donohue/ee513/mfiles/aliasexm.m

32. Sampling Aliasing II – The Sequel (FS=200, Range 350-450 Hz). Change beginning and ending frequency parameters in following mfile to run: http://www.engr.uky.edu/~donohue/ee513/mfiles/aliasexm.m

33. The Sampling Theorem A band-limited continuous signal s(t) can be reconstructed without error from its samples provided: where fs is the sampling frequency in samples per second, and fb is the frequency above which s(t) has no energy.

34. Aliased Signal Spectra Sampling in time with sampling frequency fs creates an infinite pattern of a shifted analog spectra so that frequency domain has a periodicity of fs. Let Sa(f) be the spectra of the original analog signal, the spectrum of the sampled signal becomes:

35. Aliased Signal Spectra Spectral periodicity of a low-pass signal (not really band-limited) resulting from an 8 kHz sampling

36. Restoring Sampled Signals A sampled signal is reconstructed by low-pass filtering the samples with a cut-off near the folding frequency

37. Aliased Signal Example Before sampling at a given rate, signals are often low-pass filtered (anti-aliasing filter) to limit distortions from aliasing. • Original Sound • Limited Bandwidth (LPF with 900 Hz cutoff) and sampled at 2 kHz • Original Sound sampled at 2 kHz (aliasing)

38. Homework(3) Determine the aliased frequencies in the range of For the following sampling frequency and signal pairs:

39. Linear Difference Equations Discrete systems are described by Z-transforms in the frequency domain and difference equations in the time domain. Digital filters can be designed in either domain. Find the impulse response of the following filters. (FIR) (IIR) • Compute impulse response directly by hand • Use Matlab function “filter” • Take inverse of Z transform • Examine poles and zeros of filters

40. Linear Difference Equations (FIR)

41. Linear Difference Equations (IIR)

42. Homework (4) Find the impulse response of the following filters. a) b) c) • Use Matlab function “filter” • Take inverse of Z transform • Examine poles of filter and comment on expected stability