An electronic system for heart rate detection comprising: a random sampling sensor module comprising a first sensor circuit configured for providing nonuniform random samples below Nyquist rate of a photoplethysmographic signal; a second sensor circuit configured for providing nonuniform random samples below Nyquist rate of a motion signal, said motion signal being sampled with the same nonuniform pattern as the photoplethysmographic signal; and a heart rate detection module configured for receiving a plurality of said photoplethysmographic signal nonuniform random samples and calculating a heart rate value; and characterized in that, the heart rate detection module is further configured for calculating a power spectral density value set based on a Lomb-Scargle periodogram of said photoplethysmographic signal samples; receiving a plurality of said motion signal nonuniform random samples and calculating a power spectral density value set based on a Lomb-Scargle periodogram of said motion signal samples; normalizing the calculated photoplethysmographic and motion signal power spectral density value sets; subtracting the normalized motion signal power spectral density values from the normalized photoplethysmographic signal power spectral density values; renormalizing the motion subtracted photoplethysmographic signal power spectral density value set; and calculating a heart rate value based on the frequency corresponding to a highest power peak of the renormalized photoplethysmographic signal power spectral density value set in a frequency range of interest. a random sampling sensor module comprising a first sensor circuit configured for providing nonuniform random samples below Nyquist rate of a photoplethysmographic signal; a second sensor circuit configured for providing nonuniform random samples below Nyquist rate of a motion signal, said motion signal being sampled with the same nonuniform pattern as the photoplethysmographic signal; and a first sensor circuit configured for providing nonuniform random samples below Nyquist rate of a photoplethysmographic signal; a second sensor circuit configured for providing nonuniform random samples below Nyquist rate of a motion signal, said motion signal being sampled with the same nonuniform pattern as the photoplethysmographic signal; and a heart rate detection module configured for receiving a plurality of said photoplethysmographic signal nonuniform random samples and calculating a heart rate value; and characterized in that, the heart rate detection module is further configured for calculating a power spectral density value set based on a Lomb-Scargle periodogram of said photoplethysmographic signal samples; receiving a plurality of said motion signal nonuniform random samples and calculating a power spectral density value set based on a Lomb-Scargle periodogram of said motion signal samples; normalizing the calculated photoplethysmographic and motion signal power spectral density value sets; subtracting the normalized motion signal power spectral density values from the normalized photoplethysmographic signal power spectral density values; renormalizing the motion subtracted photoplethysmographic signal power spectral density value set; and calculating a heart rate value based on the frequency corresponding to a highest power peak of the renormalized photoplethysmographic signal power spectral density value set in a frequency range of interest. receiving a plurality of said photoplethysmographic signal nonuniform random samples and calculating a heart rate value; and characterized in that, the heart rate detection module is further configured for calculating a power spectral density value set based on a Lomb-Scargle periodogram of said photoplethysmographic signal samples; receiving a plurality of said motion signal nonuniform random samples and calculating a power spectral density value set based on a Lomb-Scargle periodogram of said motion signal samples; normalizing the calculated photoplethysmographic and motion signal power spectral density value sets; subtracting the normalized motion signal power spectral density values from the normalized photoplethysmographic signal power spectral density values; renormalizing the motion subtracted photoplethysmographic signal power spectral density value set; and calculating a heart rate value based on the frequency corresponding to a highest power peak of the renormalized photoplethysmographic signal power spectral density value set in a frequency range of interest.