数字信号处理--MATLAB.ppt

Convolution Sum The M-file conv implements the convolution sum of two finite-length sequences If a=[-2 0 1 -1 3] b=[1 2 0 -1] then conv(a,b) yields [-2 -4 1 3 1 5 1 -3] DTFT Computation Using MATLAB The function freqz can be used to compute the values of the DTFT of a sequence, described as a rational function in the form of DTFT Computation Using MATLAB For example, the statement H = freqz(num, den, w) returns the frequency response values as a vector H of a DTFT defined in terms of the vectors num and den containing the coefficients {pi} and {di}, respectively at a prescribed set of frequencies between 0 and 2p given by the vector w DTFT Computation Using MATLAB There are several other forms of the function freqz The Program 3_1.m in the text can be used to compute the values of the DTFT of a real sequence It computes the real and imaginary parts, and the magnitude and phase of the DTFT DTFT Computation Using MATLAB Example - Plots of the real and imaginary parts, and the magnitude and phase of the DTFT DTFT Computation Using MATLAB The Unwrapped Phase Function In numerical computation, when the computed phase function is outside the range [-π, π], the phase is computed modulo 2π, to bring the computed value to this range Thus, the phase functions of some sequences exhibit discontinuities of 2π radians in the plot The Unwrapped Phase Function For example, there is a discontinuity of 2π at ω=0.72 in the phase response below The Unwrapped Phase Function In such cases, often an alternate type of phase function that is continuous function of ω is derived from the original phase function by removing the discontinuities of 2π Process of discontinuity removal is called unwrapping the phase The unwrapped phase function will be denoted as θc(ω) The Unwrapped Phase Function Note: The phase spectrum displays a discontinuity of 2π at w = 0.72 This discontinuity can be removed using the function unwrap as indicated below Fre