def __init__(self, x, h_func, t, delay_inc=.5): self.t_s = t[1] - t[0] self.t = t self.x = x self.h_func = h_func self.h = h_func(t) self.y = conv(x, self.h, t) self.product = np.zeros_like(self.t) self.t0 = np.argmin(np.abs(self.t)) self.t0_mirror = np.argmin(np.abs(np.flip(self.t))) self.direction = 'stop' self.delay = 0 self.delay_inc = delay_inc
def pixaxis(ax): ticks_loc = ax.get_xticks().tolist() ax.set_xticks(ax.get_xticks().tolist()) ax.set_xticklabels(["{:.0f}".format(x / np.pi) + 'π' for x in ticks_loc]) # Questão 1a sp = plt.subplots(1, 2) ax_x = sp[0].axes[0] ax_y = sp[0].axes[1] t = np.arange(-3, 16, .01) xtick = np.arange(-3, 16) axis = [-1, 15, -.5, 10.2] x = u(t) h = 2*u(t) - u(t-3) - u(t-6) y = conv(x, h, t) plotit(x=t, y=x, ax=ax_x, axis=axis, xticks=xtick) plotit(x=t, y=y, ax=ax_y, axis=axis, wait=False, save=False, xticks=xtick) sp = plt.subplots(1, 1) ax = sp[0].axes[0] plotit(x=t, y=y, ax=ax, axis=[-1, 15, -1, 10], xlabel='t [s]', ylabel='y(t)', xticks=xtick, linewidth=3) # Questão 1b x = impulse(t)+impulse(t-1)+impulse(t-2) y = conv(x, h, t) sp = plt.subplots(1, 2) ax_x = sp[0].axes[0] ax_y = sp[0].axes[1] plotit(x=t, y=x, ax=ax_x, axis=axis, xticks=xtick) plotit(x=t, y=y, ax=ax_y, axis=axis, wait=False, save=False, xticks=xtick)
import matplotlib.pyplot as plt import numpy as np from funcoes import conv n = np.arange(-1, 7) x = np.zeros_like(n, dtype=float) x[1:5] = [2, 1, 2, 1] h = np.zeros_like(x) h[1:4] = [1, 2, 1] axis = [n[0], n[-1], -.1, 2.1] y = conv(x, h, n) y[0:-1] = y[1:] y[-1] = 0 plt.subplot(2, 1, 1) plt.title('Entrada e saída') plt.stem(n, x) plt.grid() plt.axis(axis) plt.subplot(2, 1, 2) plt.stem(n, h) plt.grid() plt.axis(axis) axis = [n[0], n[-1], -.1, 2.6] plt.figure() plt.suptitle('LIT') plt.subplot(4, 2, 1) x_ = np.zeros_like(x) x_[1] = 1 plt.stem(n, x_)
from graficos import myplot as plotit from funcoes import u, window, integrate, conv, movavgc, movavgpast, derivate, exp_u, sinc import numpy as np import matplotlib.pyplot as plt t = np.arange(-1, 5, .01) h = u(t) - u(t - 1) x = np.cos(3 * np.pi * t) * u(t) ya = conv(x, h, t) x = np.cos(2 * np.pi * t) * u(t) yb = conv(x, h, t) plt.subplot(3, 1, 1) plt.plot(t, h) plt.subplot(3, 1, 2) plt.plot(t, x) plt.subplot(3, 1, 3) plt.plot(t, ya) plt.figure() plt.subplot(2, 1, 1) plt.plot(t, ya) plt.grid() plt.xlabel('t [s]') plt.ylabel('Questão (1a)') plt.legend(['y(t)']) plt.subplot(2, 1, 2) plt.plot(t, yb) plt.grid() plt.xlabel('t [s]') plt.ylabel('Questão (1b)')
x_hat = x_hat + r_(t - i * width, width) * x[t0 + i * width_s] return x_hat def x_RLi(t, R, L): return np.exp(-t * R / L) * u(t) / L sp = plt.subplots(1, 1) ax_x = sp[0].axes[0] t = np.arange(-3, 10, .01) xtick = np.arange(-3, 10) axis = [-1, 5, -.2, 1.6] x = gaussian(t, 2, 1.5) h = x_RLi(t, 2, 1) y = conv(x, h, t) plotit(x=t, y=x, ax=ax_x, axis=axis, xticks=xtick, linewidth=2, save=False, wait=True) x_ = approx(x, t, 1) plotit(x=t, y=x, ax=ax_x, axis=axis, xticks=xtick, linewidth=2) plotit(x=t, y=x_, ax=ax_x, axis=axis, xticks=xtick,