A device for determining a sequence of heart rate HR values of a user over a time interval comprising: a PPG HR sensor configured for providing PPG HR signals over a sequence of PPG HR measurements; an accelerometer coupled to the user's body; and processing circuitry (32, 34, 36, 40) configured for: transforming the PPG HR signals into the frequency domain by utilization of a Fourier transform to provide frequency domain heart rate FDHR signal magnitudes, which correspond to the absolute values of the Fourier coefficients as calculated; transforming signals from the accelerometer into the frequency domain to provide frequency domain accelerometer FDA signal magnitudes; utilizing the FDHR signal magnitudes and the FDA signal magnitudes, providing compensated FDHR component magnitudes, each compensated FDHR component magnitude corresponding to an HR value, the compensated FDHR component magnitude compensating the FDHR signal magnitudes for artifacts from acceleration as measured by the accelerometer; and selecting a path being a time sequence of compensated FDHR component magnitudes from the sequence of PPG HR measurements wherein the path includes one compensated FDHR component magnitude selected among a multitude of compensated FDHR component magnitudes representing HR probabilities at each measurement and the path selection determines an optimum path score using a reward algorithm based on the magnitude of all compensated FDHR component magnitudes along the path and a penalty algorithm based on differences in consecutive HR values corresponding to consecutive compensated FDHR component magnitudes along the path. a PPG HR sensor configured for providing PPG HR signals over a sequence of PPG HR measurements; an accelerometer coupled to the user's body; and processing circuitry (32, 34, 36, 40) configured for: transforming the PPG HR signals into the frequency domain by utilization of a Fourier transform to provide frequency domain heart rate FDHR signal magnitudes, which correspond to the absolute values of the Fourier coefficients as calculated; transforming signals from the accelerometer into the frequency domain to provide frequency domain accelerometer FDA signal magnitudes; utilizing the FDHR signal magnitudes and the FDA signal magnitudes, providing compensated FDHR component magnitudes, each compensated FDHR component magnitude corresponding to an HR value, the compensated FDHR component magnitude compensating the FDHR signal magnitudes for artifacts from acceleration as measured by the accelerometer; and selecting a path being a time sequence of compensated FDHR component magnitudes from the sequence of PPG HR measurements wherein the path includes one compensated FDHR component magnitude selected among a multitude of compensated FDHR component magnitudes representing HR probabilities at each measurement and the path selection determines an optimum path score using a reward algorithm based on the magnitude of all compensated FDHR component magnitudes along the path and a penalty algorithm based on differences in consecutive HR values corresponding to consecutive compensated FDHR component magnitudes along the path. transforming the PPG HR signals into the frequency domain by utilization of a Fourier transform to provide frequency domain heart rate FDHR signal magnitudes, which correspond to the absolute values of the Fourier coefficients as calculated; transforming signals from the accelerometer into the frequency domain to provide frequency domain accelerometer FDA signal magnitudes; utilizing the FDHR signal magnitudes and the FDA signal magnitudes, providing compensated FDHR component magnitudes, each compensated FDHR component magnitude corresponding to an HR value, the compensated FDHR component magnitude compensating the FDHR signal magnitudes for artifacts from acceleration as measured by the accelerometer; and selecting a path being a time sequence of compensated FDHR component magnitudes from the sequence of PPG HR measurements wherein the path includes one compensated FDHR component magnitude selected among a multitude of compensated FDHR component magnitudes representing HR probabilities at each measurement and the path selection determines an optimum path score using a reward algorithm based on the magnitude of all compensated FDHR component magnitudes along the path and a penalty algorithm based on differences in consecutive HR values corresponding to consecutive compensated FDHR component magnitudes along the path.