700 likes | 712 Views
Learn about Trigonometric Fourier Series, Discrete Fourier Transform, DSP benefits, signal manipulation, coherence, and array processing in mixed-signal testing.
E N D
Outline • Trigonometric Fourier Series (FS) • Discrete-Time Fourier Series (DTFS) • Relationship to FS • Working directly with samples • Complex form • Discrete Fourier Transform (DFT) • Fast Fourier Transform (FFT) • Applications • Equivalence of Time and Frequency Domains • Frequency Domain Filtering • Summary
Advantages of DSP Based Testing • Reduced Test Time • DSP in this class will be limited to discrete (i.e. sampled) waveforms of finite length. • Advantages of coherent DSP based testing • reduced test time since we can create signals with multiple frequencies at the same time. • Once the output response of the DUT has been captured using a digitizer or capture memory, DSP allows the separation of test tones to give individual gain and phase measurements. • Also, by removing the input test tones, we can measure noise and distortion without running many separate tests.
Advantages of DSP Based Testing • Separation of Signal Components • By using coherent test tones, we are guaranteed that the harmonic distortion components will fall neatly into separate Fourier spectral bins rather than being smeared across many bins. • DSP based testing also has the major advantage in the elimination of errors and poor repeatability. • Advanced Signal Manipulations • DSP allows us to manipulate digitized output waveforms to achieve a variety of results • We can apply mathematical filters to remove noise thereby achieving better accuracy.
Digital Signal Processing • DSP and Array Processing • There is a slight difference between array processing and Digital Signal Processing. • An array (or vector) is a series of numbers (i.e. height of students in class) • A Digital Signal is also a series of numbers (i.e. voltages), yet the series is time stamped • Thus digital signal processing is a subset of array processing using time-ordered samples. • All DSP is accomplished on a special computer called the array processor (so much for the difference)
Digital Signal Processing • DSP and Array Processing - cont. • Array processing functions that are useful in mixed-signal testing: • averaging • To measure the RMS of a signal we must first remove the DC offset - this is accomplished by averaging the signal and subtracting the result from the original • Many functions like averaging are built into the ATE tester code set to allow easy use. • Built in functions are set up to maximally utilize the available computational resources to reduce test time.
Digital Signal Processing • DSP and Array Processing - cont. • Other built in functions include: • vector average - average value of an array • vector RMS - root mean square of the array values • max/min - maximum and minimum values in an array • vector add - add two arrays • add scalar to vector - add constant to each array value • subtract scalar from vector - subtract constant from each array value • vector multiply - multiply two arrays • multiply vector by scalar - multiply each array element by a constant • divide vector by scalar - divide each array element by a constant
Discrete Fourier Analysis • Fourier Transform • Jean Baptiste Joseph Fourier • French mathematician that found that any periodic waveform can be described as the sum of a series of sine and cosine waves at various frequencies plus a DC offset. • Developed for the study of heat transfer in solid bodies • A sequence is assumed to be periodic with a period T such that x(t) = x(t-T) for all values of t from minus infinity to plus infinity. • x(t) = a0+a1*cos(w0t)+b1*sin(w0t)+a2*cos(2 w0t)… … + to infinity
Discrete Fourier Analysis • Discrete Fourier Transform • Mathematical operation that allows us to split a composite signal into its individual frequency components. • A DFT operation is equivalent to a series of very narrow band pass filters followed by peak-responding voltmeters. The filters are not only frequency selective but also phase selective to determine the sine and cosine contributions individually. • x(n) = a0+a1*cos(2n/N)+b1*sin(2n/N)+a2*cos(2n/N)… … +a(N/2)*cos(2(N/2)n/N) + b(N/2)*sin (2(N/2)n/N)
Discrete Fourier Analysis • Discrete Fourier Transform - cont. • Digitizing spectrum analyzers and mixed-signal testers accomplish the filter and peak measurements using the DFT. The DFT uses a frequency sensitive correlation calculation for each value of a and b. • Functions that have zero correlation are called orthogonal • Superposition and orthogonality of coherent sine and cosine components allows us to extract the value of all a’s and b’s, even in the presence of other coherent test tones. The cosine correlation function is equivalent to a filter and peak measurement. Therefore we can measure many signals simultaneously, reducing test time.
Discrete Fourier Analysis • Complex form of the DFT • Most traditional DSP books use the Euler’s transform to convert sinusoids into exponentials. • e-j w t = cos(wt) – j*sin(wt)
Discrete Fourier Analysis • Complex form of the DFT • Notice that the complex form of the DFT correlates with a negative sine wave instead of a positive sine wave in the sine/cosine version. • This causes problems in the phase shift calculations!!! • Some testers will give the straight imaginary value, while others multiple by minus one to compensate for the difference. • The test engineer will need to find out whether the tester is reporting sine amplitudes or imaginary components before phase measurements can be made!!!
T k f DC f 2 f 3 f Fourier Analysis Of Periodic SignalsTrigonometric Form • For any periodic signal with a finite number of discontinuities, the signal can be represented by a Fourier Series:
Computing Fourier Coefficients • Coefficients are found from the following integral equations:
Fourier Series RepresentationMagnitude & Phase Form Rectangular Form: Magnitude&Phase Form: where
Amplitude (or RMS Value) c1 Time/Frequency c2 co c3 0 f fo 2fo 3fo 4fo 5fo 0 Spectral Plot = Phase f1 f3 f 0 f2 fo 2fo 3fo 4fo 5fo 0
Actual Vs. FS Representation • Increasing the number of terms in the FS increases the accuracy of the representation. • Gibbs phenomenon (overshoot at discontinuity) is a result of the finite sum of terms.
Outline • Trigonometric Fourier Series (FS) • Discrete-Time Fourier Series (DTFS) • Relationship to FS • Working directly with samples • Complex form • Discrete Fourier Transform (DFT) • Fast Fourier Transform (FFT) • Applications • Equivalence of Time and Frequency Domains • Frequency Domain Filtering • Summary
TS Discrete-Time Fourier SeriesFirst Principles Consider sampling x(t): But, FS=1/TS, allowing us to write
T TS Discrete-Time Fourier SeriesCoherent Sampling • Generally, we are interested in only those sample sets that are derived from a signal that satisfies T=NTSor fo=FS/N:
T TS Discrete-Time Fourier SeriesPeriodic Sample Sets • The fact that we are using coherent sample sets, implies periodicity in n. However, due to the symmetry of the formulation, x[n] is also periodic with respect to k with period N:
Discrete-Time Fourier SeriesRe-Grouping Formulation Split into 2 parts: To simplify further, use trig. substitutions:
Discrete-Time Fourier SeriesRe-Grouping Formulation Replace infinite summations with single parameter: DTFS:
Discrete-Time Fourier SeriesMagnitude & Phase Notation Rectangular Form: Magnitude&Phase Form: Used for spectral plot purposes where
10-5 s DTFS ExampleClock Signal Evaluate Infinite Summations: After 100 terms:
DTFS interpolation samples clock signal DTFS ExampleClock Signal DTFS:
DTFS interpolation samples clock signal FS Versus DTFSClock Signal Example • Unlike a FS that attempts to represent the periodic function over all time, a DTFS only attempts to represent the N periodic samples • Hence, a much simpler mathematical expression.
T TS Working Directly With DTFS Strategy to solve for unknown parameters: -Each sample must satisfy the DTFS for x[n] • A DTFS has N unknown parameters corresponding to N degrees of freedom. • A DTFS is a representation for a coherent sample set consisting of N samples.
Solving N Equations In N Unknowns 1st sample: (n=0) 2nd sample: (n=1) Nth sample: (n=N-1)
Matrix Formulation & Solution Compact notation: Unknown parameters:
Method of Orthogonal Basis • Even before Fourier’s development in the 1800’s , the famous mathematician, Euler had developed a closed-form solution for finding the unknown coefficients of the DTFS. • involves projections onto a set of orthogonal basis functions (harmonically-related sinusoids). • his efforts were dropped in the direction of Fourier analysis because of the conceptual difficulties that occurred with the step discontinuities in the signal. • The importance of this method is that it forms the basis of all modern methods related to Fourier Analysis, Wavelets, etc.
Method of Orthogonal Basis DTFS Coefficients: • The above formulae are found by multiplying the DTFS by (i) cos[k(2p/N)n] (ii) sin[k(2p/N)n], then summing n from 0 to N-1.
10-5 s DTFS ExampleClock Signal 10 samples bk coefficients ak coefficients All other coefficients are zero.
Frequency Denormalization FS = 100 kHz N =10 • DTFS is expressed in normalized time and frequency. • To return to proper time scale: • To return to proper frequency scale:
Complex Form of the DTFS • Through the application of Euler’s identity, we can convert the DTFS in trigonometric form to the complex form of the DTFS, where
Complex Form of the DTFSSeveral Examples Example 1: Example 2:
Outline • Trigonometric Fourier Series (FS) • Discrete-Time Fourier Series (DTFS) • Relationship to FS • Working directly with samples • Complex form • Discrete Fourier Transform (DFT) • Fast Fourier Transform (FFT) • Applications • Equivalence of Time and Frequency Domains • Frequency Domain Filtering • Summary
Discrete-Time Fourier Transform • Fourier greatest invention was the Fourier Transform (FT). • provides a frequency description (known as a Fourier transform) of an aperiodic signal (transient signal) • If y[n] exists for only finite time, then we can represent it by the following periodic function x[n] with period N (periodic extension of y[n]): y[n]
Discrete-Time Fourier TransformAperiodic Signal Description y[n] • Given some aperiodic signal y[n] that can be described in terms of a periodic signal x[n], then we can write • As x[n] is a periodic function, we can write y[n] as
Discrete-Time Fourier TransformInvestigating Impact of N-> add zeros • As the period N is made larger, a better match is made between y[n] & x[n]. As N->, y[n]=x[n] for all finite values of n. • Due to limiting argument, the infinite sum eqn. changes into an integral eqn: • The term Y(ejw) is called the D.T. Fourier Transform of y[n], given by
y[n] 1 F.T. 0 n 0 4 Discrete-Time Fourier TransformExample • Consider a set of samples from a unit-height rectangular pulse signal, the F.T. would be computed as follows: Note: Spectrum is continuous. |Y(w)| 5 w -2p/5 0 2p/5 4p/5 -4p/5 (w) p 4p/5 2p/5 w 0 -2p/5 -4p/5 -p
Relationship Between DTFS & FT |Y(w)| 5 w -2p/5 0 2p/5 4p/5 -4p/5 • The spectral coefficients of an N-point DTFS are samples of the FT: • Substituting the appropriate values for Y(ejw) gives