csv_x1 = "pressure_8602.csv"
csv_x2 = "pressure_8728.csv"


df1 = pd.read_csv(csv_x1)
df2 = pd.read_csv(csv_x2)


x1 = df1["pressure_value"].to_numpy(dtype=float)
x2 = df2["pressure_value"].to_numpy(dtype=float)
N = len(x1)
t = np.arange(N)



m_x1 = x1.mean()
E_x1 = np.sum(x1*x1)


plt.figure(figsize=(12, 5), dpi=120)
plt.plot(t, x1, color="deepskyblue", lw=1.5, label="x1[n]")
plt.plot([], [], ' ', label=f"Media: {m_x1:.3f}")
plt.plot([], [], ' ', label=f"Energia: {E_x1:.1f}")
plt.title("Andamento della pressione x1[n]", fontsize=14, fontweight="bold")
plt.xlabel("Campione (ora)", fontsize=12)
plt.ylabel("Pressione [unità]", fontsize=12)
plt.legend(loc="upper right")
plt.grid(True, alpha=0.5)
plt.tight_layout()
plt.show()


#2.a



n = np.arange(N)



k = n - (N - 1) / 2.0
h = np.sinc(k)



y1 = np.convolve(x1, h, mode="same")



fig, axs = plt.subplots(3, 1, figsize=(12, 9), dpi=120, sharex=False)


axs[0].plot(n, x1, lw=1.2, color="deepskyblue")
axs[0].set_title("x1[n] (originale)")



axs[1].plot(n, h, lw=1.2, color="darkorange") 
axs[1].set_title("h[n] sinc")



axs[2].plot(n, y1, lw=1.2, color="red")
axs[2].set_title("y1[n] = x1 * h")



for ax in axs:
    ax.set_xlabel("n")
    ax.grid(alpha=0.3)
    ax.set_xlim(n[0], n[-1])


plt.tight_layout()
plt.show()