Fourier series and gibbs phenomenon

That funny phase

You probably noticed in the last problem that even though the wave forms looked fairly different, the sound was similar. Let's look into this a bit deeper with a simpler sound.

• Create a script file called phasefun.m to put your code in for this problem.
• Pick two harmonic frequencies and generate a signal from two cosines at these frequencies added together and call it sig1 . Use Fs = 8000 (remember that you can reproduce only frequencies that are less than Fs/2 ).
• Now generate a second signal called sig2 exactly the same as the first one, except time delay the second cosine by a half cycle (half of its period).
• Use subplot to show a few periods of both signals, do they look different? Save the plot as phasesigs.tif . What did the time delay do to the phase?
• Play each signal with soundsc , do they sound different?
• Redo sig2 with a few different delays and compare the sound to the first signal.
• Create a sig3 that is one cosine at some frequency. Now add sig3 with a timed delayed version of itself and call it sig4 . Use a quarter cycle delay.
• Use subplot and plot a few periods of sig3 and sig4 . Play them with soundsc , do they sound different to you?
• What is suggested about our hearing capabilities from this experiment?
• You will need to show the TA the following files: phasefun.m phasesigs.tif

Gibbs phenomena

• Create a script file called gibbs.m to put your code in for this problem.
• Click here to download the MATLAB function Ck.m . Take a look at the contents of the function. This function takes one argument $k$ , and creates the $k$ th Fourier series coefficient for the squarewave above:
  $C_k=\{\begin{array}{ll}0\hfill & \text{if }k=0,k\text{ even}\hfill \\ \frac{1}{jk\pi }\left[\cos \left(\frac{2k\pi }{3}\right)-\cos \left(\frac{k\pi }{3}\right)\right]\hfill & \text{if }k\text{ odd}\hfill \end{array}$
  $C_k(1)=\frac{-1}{j\pi }=0+j0.3183$ . Plot the magnitude and phase of thecoefficients $C_k$ for $k\in \left\{-10,-9,\dots ,9,10\right\}$ . The magnitude and phase should be plotted separately using the subplot command, with the magnitude plotted in the top half of the window and the phase in the bottom half. Place frequency $w$ on the x axis. Use the MATLAB command stem instead of plot to emphasize that the coefficients are a function of integer-valued (not continuous) $k$ . Label your plots.
• Save the graph as Coeff.tif .
• Write whatever script/function files you need to implement the calculation of the signal $x(t)$ with a truncated Fourier series:
  $x(t)=\sum_{k=-K_{\max }}^{K_{\max }} C_{k}e^{(\mathrm{jk}{\omega }_{0}t)}=\sum_{k=0}^{K_{\max }} 2\left|C_k\right|cos(k{\omega }_{0}t+{\angle C}_k)$
  for a given $K_{max}$
  You can avoid numerical problems and ensure a real answer if you use the cosine form. For this example, $w_0=1$ .
• Produce plots of $x(t)$ for $t\in \left[-5 , 5\right]$ for each of the following cases: $K_{max}$ = 5; 15; and 30. For all the plots, use as your time values the MATLAB vector t=-5:.01:5 . Stack the three plots in a single figure using the subplot command and include your name in the title of the figure. Save the figure as FourTrunc.tif
• Add clear comments describing what the files do. You will need to show the TA the following files: gibbs.m Coeff.tifFourTrunc.tif