Using MATLAB

A Fourier series representation for a square wave, with frequency fo, is given by the following infinite summation: 4 4 4 11? Note, the above equation assumes time in units of seconds and frequency in units of Hertz. Write an m-file that computes the summation shown above, truncated after different numbers of terms. The truncated series provides an approximation of the square wave function. . Set the fundamental frequency to fo-500kHz. . Use a for loop to calculate and store the summation, truncated after 1, 2, 3, 4, and 5 terms. That is, compute five separate vectors, each the same length as the time vector, each including a different number of terms from the Fourier series. . Compute each of the truncated series over a time range of 0 t S 4usec .Store each resulting vector as a row in a matrix, f. This matrix should, therefore, have dimensions 5 x length(t). . So, for example, the second row in f should include the first two terms of the series and be computed using the following function: f(t)--sin(2??t) + 3-sin(21-3 *Plot rows 2-5 of the matrix f in subplots 1-4 of a 2 x 2 subplot grid. . Plot each function versus time in microseconds (not seconds), as shown below . The title of each subplot should indicate the number of terms used in each truncated series approximation, as shown below Your plot should look exactly like the following plot: N- 2 N-3 -2 -2 time usec) N-4 time [usec N-5 -2 -2 time lusec) time [usec