import matplotlib.pyplot as plt
import numpy as np
import sk_dsp_comm.sigsys as ss
fs = 100 # sampling rate in Hz
tau = 1
t = np.arange(-5,5,1/fs)
x0 = ss.rect(t-.5,tau)
plt.figure(figsize=(6,5))
plt.plot(t,x0)
plt.grid()
plt.ylim([-0.1,1.1])
plt.xlim([-2,2])
plt.title(r'Exact Waveform')
plt.xlabel(r'Time (s)')
plt.ylabel(r'$x_0(t)$')
plt.show()
#
import matplotlib.pyplot as plt
import numpy as np
import sk_dsp_comm.sigsys as ss
fs = 100 # sampling rate in Hz
tau = 1
t = np.arange(-5,5,1/fs)
x0 = ss.rect(t-.5,tau)
fe = np.arange(-10,10,.01)
X0e = tau*np.sinc(fe*tau)
plt.plot(fe,abs(X0e))
plt.grid()
plt.xlim([-10,10])
plt.title(r'Exact (Theory) Spectrum Magnitude')
plt.xlabel(r'Frequency (Hz)')
plt.ylabel(r'$|X_0e(f)|$')
plt.show()
#
import matplotlib.pyplot as plt
import numpy as np
import sk_dsp_comm.sigsys as ss
fs = 100 # sampling rate in Hz
tau = 1
t = np.arange(-5,5,1/fs)
x0 = ss.rect(t-.5,tau)
f,X0 = ss.ft_approx(x0,t,4096)
plt.plot(f,abs(X0))
plt.grid()
plt.xlim([-10,10])
plt.title(r'Approximation Spectrum Magnitude')
plt.xlabel(r'Frequency (Hz)')
plt.ylabel(r'$|X_0(f)|$');
plt.tight_layout()
plt.show()