The invention belongs to the technical field of
biomedical signal processing, and particularly relates to an automatic
respiratory rate calculation method for respiratory sound
processing of an
electronic stethoscope, which comprises the following steps of: S1, acquiring a period of respiratory sound
signal with the duration of Ts (Ts is greater than or equal to 15 seconds) and the sampling rate of fs by using the
electronic stethoscope; s2, performing short-time
Fourier transform on the obtained
breath sound signal vector x; s3, calculating the amplitude S of the time-
frequency spectrum matrix S; s4, carrying out homomorphic envelope st on the
time domain vector # imgabs0 # of the
intermediate frequency component; and S5, searching all peak values greater than a preset threshold Eth on the
time domain homomorphic envelope st of the
intermediate frequency component. By calculating the
time domain homomorphic envelope of the
intermediate frequency component, main interference components including heart sound can be automatically eliminated, the detection specificity of the
breath sound signal component is greatly improved, and then the calculation precision of the
breath rate is improved;
respiratory cycle detection and
respiratory rate calculation can be realized only by searching a
peak value meeting a threshold condition on the time domain homomorphic envelope of the intermediate frequency component without any
pattern recognition method with relatively
high complexity.