# P wave def p_wave(t_center, t): return 0.25 * np.exp(-((t - t_center)**2) / (0.04))
beat_interval = 60 / bpm for i in range(int(duration / beat_interval)): center = i * beat_interval + 0.2 # offset to place QRS ecg += qrs(center, t) ecg += p_wave(center - 0.15, t) ecg += t_wave(center + 0.3, t) heartbeat simulator
# Add baseline drift and noise drift = 0.05 * np.sin(2 * np.pi * 0.2 * t) noise = 0.02 * np.random.randn(len(t)) return t, ecg + drift + noise t, signal = generate_ecg(bpm=85) plt.plot(t, signal) plt.title("Simulated ECG for Heartbeat Simulator") plt.xlabel("Time (s)") plt.ylabel("Amplitude (mV)") plt.grid(True) plt.show() # P wave def p_wave(t_center, t): return 0
# T wave def t_wave(t_center, t): return 0.35 * np.exp(-((t - t_center)**2) / (0.1)) # P wave def p_wave(t_center