??
? ?? ????/????? ??? ???? ?? ??? ?? ??? HR ???? ??? ECG ?????? ??? ????.

?? ??:

• ???
• ??? ?????
• 1/4?? ???/?? ???/?? ???(??? ?????? ?? ???? ????? ?? ?)

?? ????

?? ??? ?? ??? ?????. ??? ?? ? ??? ?? ???? ??? ? ????.

?? ??? ???? ??? ??? ??? ? ????. ??? ??? ??? ???? ??? ?? ??? ???? ????. ???? ??? ?? ?? ???? ??? 1~40Hz ????? ???? ????.

????

?? ??? ? 1/4?? ???? ?? ?? ??? ????. ?? ?? USB ??? ?????? ???? ???? ??? ???? ???? ?????. . ????? Python?? ????? ?????.

??

1??: 1/4?? ???? ???? ?? ? ???? ?????. ? ?? ?? ??? ??? ???. ??? ???? ?? ???? ?? ???. ??/?? ??? ??(?? ????? ? ?? ??? ?? ? ????. ????? ?? ???? ??? ?????.) ??? ?????? ??? ?????(?? ??? ????).

2??: ?? ??? ?????. input_device_index ??? ??? ?????? ????? ?????. ??? ?? ?? ?? ???? ???? ??? ??? fft? ???? ??? ???? ???? ?? ???? ???? 0?? ??? ?? ?? ?? ???? ???? ????. ???? HR? ???? ?? ??? ?? ?? ????? ???? ???? ?????.

import numpy as np 
import pyaudio as pa 
import struct 
import matplotlib.pyplot as plt 
from scipy.signal import decimate, find_peaks

CHUNK = 4410 #.1 second
FORMAT = pa.paInt16
CHANNELS = 1
RATE = 44100 # in Hz
fstep = RATE/CHUNK
p = pa.PyAudio()

values = []
dsf=44 #down sample factor
rds=RATE/dsf #down sampled rate

stream = p.open(
    format = FORMAT,
    channels = CHANNELS,
    rate = RATE,
    input_device_index=3, #adjust based on input
    input=True,
    frames_per_buffer=CHUNK
)

#set up graph
fig,ax = plt.subplots(1)
x = np.arange(0,2*CHUNK,2)
line, = ax.plot(x, np.random.rand(CHUNK))
ax.set_ylim(-100,100) 
ax.set_xlim(0,2500) 
text = ax.text(0.05, 0.95, str(0), transform=ax.transAxes, fontsize=14,
       verticalalignment='top')
fig.show()

def getFiltered(x,hp=1,lp=41): #this sets the unneeded freqs to 0

    fft=np.fft.fft(x)
    hptrim=len(fft)/RATE*hp
    lptrim=len(fft)/RATE*lp
    fft[int(lptrim):-int(lptrim)]=0 
    fft[0:int(hptrim)]=0 
    return np.real(np.fft.ifft(fft))

def getHR(x): 

    pdis = int(0.6 * rds) #minimum distance between peaks. stops rapid triggering. also caps max hr, so adjust
    peaks, _ = find_peaks(x, distance=pdis, height=0.1)
    intervals = np.diff(peaks)/rds # in seconds
    hr = 60 / intervals # in BPM
    return peaks,round(np.mean(hr),0) #peaks,avg hr

while 1:
    data = stream.read(CHUNK)
    dataInt = struct.unpack(str(CHUNK) + 'h', data)

    filtered=getFiltered(dataInt) #filter (working with full chunk)
    dsed=decimate(filtered, 44) #down sample (turns chunk into ds chunk)
    values=np.concatenate((values,dsed)) #puts the chunks into an array
    peaks,hr = getHR(values*-1) # gets the peaks and determins avg HR. 

    text.set_text(str(hr))
    line.set_xdata(np.arange(len(values))) 
    line.set_ydata(values*-10) #the negative is bc it comes in upside down with my set up. the *10 is just for fun
    ax.set_xlim(max(0,len(values)-2500),len(values)) #keep the graph scrolling
    vlines = ax.vlines(peaks,ymin=-100,ymax=100,colors='red', linestyles='dashed') # pop some lines at the peaks 

    fig.canvas.draw()
    fig.canvas.flush_events()
    vlines.remove()

    if len(values)>10000: #keeps the array managably sized, and graph scrolling pretty
        values=values[5000:] #5 seconds @ ~1000 sr.

??
???? ??? ?????. ??? ???? ???? ??? ? ? ? ?? ???? ? ?? ????. Garmin ??? ???? ??? ?? ????? ??? ?? ???????.

??

????
????? ?? ??? ??? ??? ??? ?? ?????. ???? ?? ?? ???? ??? ???? ??? ?????? ???? ????... ? ??? ?? ??? ????? ????. ?? ???? ????. ?? ??? ??? ?? ?? ???? ???? ??? ????.

?? ??
? ??? ECG ??? ?? ?? ??? ???? ?? ?? ??? ? ???? ????. ??? ?? ??? ???? ???? ????.

EMG? ?? ?? ??? ???? ?? ????? ????.

???? ????? ??? ? ??? ???? ?? ??? ??? ??? ?? ??? ??? ?? ??? ??? ? ????. ??? ??? ?? ??? Amazon?? ? ?????(?????. ???? ???????). ??? ?? ? ?? ?? ??? ??????. ??? ?? ???? ?? ?????? ??? ?? JFET opamp(?????? ?? ??)? ???? ??? ?? ??? ??????. ?? 3?, ??? ???? ?????? ?????. ADC? ??? ?????? ???? ?? ??? ?? ??? 3?? ????? ???? ???? ???(??? ???? ? ?? ??)

?
? ?? ????? ??? ?? ???? ?? DAW? EQ ????? ??? ??? ?????.

? ??? ??? DIY HR ??? ECG ?????? ?? ?????. ??? ??? PHP ??? ????? ?? ?? ??? ?????!