Line Spectra and their Applications

Transcription

1 In [ ]: cd matlab pwd Line Spectra and their Applications Scope and Background Reading This session concludes our introduction to Fourier Series. Last time ( Resources/blob/master/week5/exp_fs.ipynb) we saw that we could exploit the complex exponential e jωt to redefine trigonometric Fourier Series into the Exponential Fourier Series and in so doing we eliminate one integration and at the same time simplify the calculation of the coefficients of the Fourier series. In this session we show how the Exponential form of the Fourier Series leads us to the ability to present wavefoms as line spectra, simplifies the calculation of power for systems with harmonics and leads in the limit as approaches infinity to the Fourier Transform. T The material in this presentation and notes is based on Chapter 7 (Starting at Section 7.0) of Steven T. Karris, Signals and Systems: with Matlab Computation and Simulink Modelling, 5th Edition. ( ppg=3&docid=338497&tm= ) from the Required Reading List. Some clarificattion was needed and I used Chapter 4 of Benoit Boulet, Fundamentals of Signals and Systems ( ppg=50&docid=33597&tm= ) from the Recommended Reading List for this. Agenda Last time Exponents and Euler's Equation The Exponential Fourier series Symmetry in Exponential Fourier Series Example Page of 3

2 This Time Line spectra Power in periodic signals Steady-state response of an LTI system to a periodic signal Line Spectra When the Exponential Fourier series are known is is useful to plot the amplitude and phase of the harmonics on a frequency scale. This is the spectrum of the Exponential Fourier Series calculated last time Page of 3

3 Line Spectra for Trig. FS If we take the results for the Exponential Fourier Series and gather terms, the amplitudes for the Trig. Fourier Series are given by: = a 0 C 0 = ( + ) a C k C k = j( ) b k C k C k Applying this to the previous result we get Page 3 of 3

4 Example 3 Compute the exponential Fourier series for the waveform shown below and plot its line spectra. Solution The recurrent rectangular pulse is used extensively in digital communication systems. To determine how faithfully such pulses will be transmitted, it is necessary to know the frequency components. What do we know? T/w The pulse duration is. The recurrence interval T is w times the pulse duration. is the ratio of pulse repetition time to the pulse duration normally called the duty cycle. w Coefficients of the Exponential Fourier Series? Given jk( t) C k = f( Ω t) d( t) π 0 e Ω 0 Ω 0 Is the function even or odd? Does the signal have half-wave symmetry? What are the cosequencies of symmetry on the form of the coefficients C k? What function do we actually need to integrate to compute C k? π π Page 4 of 3

5 DC Component? Let k = 0 then perform the integral Page 5 of 3

6 Harmonic coefficients? Integrate for k 0 Exponential Fourier Series? Page 6 of 3

7 Effect of pulse width on frequency spectra Recall pulse width = T/w w = Ω 0 = w = T = π T/w = π rad/s; ; s; s. Page 7 of 3

8 w = 5 Ω 0 = w = 5 T = π T/w = π/5 rad/s; ; s; s. Page 8 of 3

9 w = 0 Ω 0 = w = 0 T = π T/w = π/5 rad/s; ; s; s. Implications As the width of the pulse reduces the width of the freqency spectra needed to fully describe the signal increases more bandwidth is needed to transmit the pulse. Note Text book seems to get the wrong results. Karris plots rather than in producing the diagrams shown in Figs However, if you view as in indication of the bandwidth needed to transmit a pulse of width the plots Karris gives make more sense. T/w sin(wx)/wx sin(wx)/(wx) sin(x/w)/(x/w) Page 9 of 3

10 Example 4 Use the result of Example to compute the exponential Fourier series of the impulse train shown below δ(t ± πk) Solution To solve this we take the previous result and choose amplitude (height) so that area of pulse is unity. Then we let width go to zero while maintaining the area of unity. This creates a train of impulses. δ(t ± πk) and, therefore C n = π A f(t) = π k= e jk Ω 0 t Try it! Page 0 of 3

11 Proof! From the previous result, and the pulse width was defined as T/w C n, that is A =. w sin(kπ/w) kπ/w T w = π w Let us take the previous impulse train as a recurrent pulse with amplitude w A = = =. T/w π/w π Pulse with unit area The area of each pulse is then π w w = π and the pulse train is as shown below: Page of 3

12 New coefficents The coefficients of the Exponential Fourier Series are now: π/w 0 and as impulse train. C n w/π sin(kπ/w) = = w kπ/w π sin(kπ/w) kπ/w each recurrent pulse becomes a unit impulse, and the pulse train reduces to a unit Also, recalling that sinx lim x 0 x = the coefficents reduce to C n = π That is all coefficients have the same amplitude and thus f(t) = π n= e jk Ω 0 t Spectrum of Unit Impulse Train The line spectrum of a sequence of unit impulses δ(t ± kt) is shown below: Page of 3

13 Another Interesting Result Consider the pulse train agin: What happens when the pulses to the left and right of the centre pulse become less and less frequent? That is what happens when? T Well? T Ω 0 0 t = 0 As the fundamental frequency We are then left with just one pulse centred around. The frequency difference between harmonics also becomes smaller. Line spectrum becomes a continous function. This result is the basis of the Fourier Transform which is coming next. Power in Periodic Signals Page 3 of 3

14 In your previous courses you may have come across the definitions of Signal Energy, Average Signal Power and Root Mean Square Power: T E = f(t) dt 0 P av = f(t) dt T 0 T P RMS = f(t) dt T T 0 Parseval's Theorem Parseval's Theorem ( states that the total average power of a a periodic signal is equal to the sum of the average powers of all its harmonic components. f(t) k C k e jkω 0 t The power in the th harmonic is given by P k = dt = dt = T C k e jkω 0 t C T k P k = P k k P k Since, the total power of the th harmomic is. 0 T 0 T C k Parseval's theorem states that T P = f(t) dt =. T C k 0 k= RMS Power By a similar argument: P RMS T = f(t) dt =. T C 0 k k= Page 4 of 3

15 Example 4 Compute the average power of a pulse train for which the pulse width is previous result: T/ (duty cycle 50%). Use the C n A =. w sin(kπ/w) kπ/w as your starting point. Power Spectrum C k C k = Ck = C k The power spectrum of signal is the sequence of average powers in each complex harmonic: real perodic signals the power spectrum is a real even sequence as.. For Page 5 of 3

16 T/8 This is the power spectrum for a pulse with width. Note that most of the power is concentrated at DC and in the first seven harmonic components. That is in the frequency range rad/s. [ 4π/T, + 4π/T] Page 6 of 3

17 Total Harmonic Distortion Suppose that a signal that is supposed to be a pure sine wave of amplitude A is distorted as shown below This can occur in the line voltages of an industrial plant that makes heavy use of nonlineear loads such as electric arc furnaces, solid state relays, motor drives, etc (E.g. Tata Steel!) THD Defined k ± Clearly, some of the harmonics for are nonzero. One way to characterie the distortion is to compute the ratio of average power in all the harmonics that "should not be present", that is for k >, to the total average power of the distorted sine wave. The square-root of this ratio is called the total harmonic distortion (THD) of the signal. Page 7 of 3

18 If the signal is real and based on a sine wave (that is odd), then C k = 0 and and we can define the THD as the ratio of the RMS value for all the harmonics for to the RMS of the fundamental which is f RMS = C k k= C : C k k= THD = 00 % C K > (the distortion) Computation of THD Page 8 of 3

19 Steady-State Response of an LTI System to a Periodic Signal The response of an LTI system with impulse response to a complex exponential signal e st is the same complex exponential multiplied by a complex gain:, where: h (t) y(t) = H(s)e st H(s) = h (τ) dτ. e sτ s = jω In particular, for, the output is simply. H(s) H(jω) y(t) = H(jω)e jωt The complex functions and are called the system's transfer function and frequency response, respectively. Page 9 of 3

20 By superposition The output of an LTI system to a periodic function with period by: T represented by a Fourier series is given where Ω 0 = T/π is the fundamental frequency. y(t) = C k H(jk Ω 0 ) e jkω 0 t k= Thus y(t) is a Fourier series itself with coefficients D k : D k = C k H(jk Ω 0 ) Implications of this important result The effect of an LTI sustem on a periodic input signal is to modify its Fourier series through a multiplication by its frequency response evaluated at the harmonic frquencies. Illustration This picture below shows the effect of an LTI system on a periodic input in the frequency domain. Page 0 of 3

21 Filtering A consequence of the previous result is that we can design a system that has a desirable frequency spectrum H(jk Ω 0 ) that retains certain frequencies and cuts off others. We will return to this idea later. End of Second Hour Summary Line spectra Power in periodic signals Steady-state response of an LTI system to a periodic signal Next Time The Fourier Transform Answers Given π jk( t) C k = f(t) e d( t) π Ω 0 Ω 0 π f(t) = f( t) Is the function even or odd? even! Does the signal have half-wave symmetry? No! What are the cosequencies of symmetry on the form of the coefficients C k? C k will be real values. Trig. equivalent no sine terms. What function do we actually need to integrate to compute C k? We only need to integrate between the limits π/w π/w π/w jk( C k = Ae Ω t) A 0 jk( d( Ω t) = d( t) π 0 e Ω 0t) Ω π 0 π/w π/w π/w Page of 3

22 Solution: DC component! C 0 π/w A = ωt = π π/w A π π ( w + π w ) or C 0 = A w Harmonic coefficients! C k π/w π/w A = = jkπ e jk(ωt) C k = A e jkπ/w e jkπ/w kπ ( ) = A sin kπ j kπ ( w ) A w sin (kπ/w) kπ/w Exponential Fourier Series! f(t) = k= A w sin(kπ/w) kπ/w e k Ω 0 t Solution 4 w = so: Write down an expression for P C n A =. sin(kπ/) kπ/ using Parseval's Theorem P P = C k = A kπ sinc = k= k= A ( + kπ sinc 4 k= 4 ) (kπ/) = 0 k k = 0,, 4, 6, sinc for even ( ) so...? Page of 3

23 P for k odd P = ( + kπ sinc 4 ) = A + k=,3,5, 4 A k=,3,5, sin kπ kπ sin(kπ/) = k k =, 3, 5, 7, for odd ( ) so...? P after eliminating sine P = A ( + [ ] = 4 π 9 5 ) k=,3,5, A ( + π 4 π [ 8 ]) P = A Check P from f(t) π/ P = f(t) dωt = A dωt = π π A π ( + π π ) = A. π/ π/ π/ Page 3 of 3