Part B Sliding-band Dynamic Range Compression ( N. Tiwari & P. C. Pandey, NCC 2014 )

1. P. C. Pandey, "Signal processing for improving speech perception by persons with sensorineural hearing loss: Challenges and some solutions", IEEE Workshop on Intelligent Computing, IIIT Allahabad, 13-15 Oct. 2014 Part B Sliding-band Dynamic Range Compression (N. Tiwari & P. C. Pandey, NCC 2014)

2. Overview • Introduction • Sliding-band Dynamic Range Compression • Offline & Real-time Implementations • Test Results • Summary & Conclusion

3. Introduction • Dynamic range compression • To present sounds comfortably within the limited dynamic range of the listener by amplifying the low level sounds without making the high level sounds uncomfortably loud. • Processing steps • Input level estimation • Gain calculation based on input level • Multiplication of input with gain function • Output resynthesis • Classification • On the basis of signal level calculation: single-band or multiband • On the basis of gain control method: feedback or feed-forward

4. Single-band dynamic range compression • Processing • Gain dependent on the dynamically varying signal level. • Parameters • Compression threshold (Th) • Compression ratio (CR) • Attack & release time • Problems • Does not account for frequency dependent loudness growth function • Power mostly contributed by low-frequency components • → amplification of high-frequency components depends low-frequency components • → Inaudibility of high frequency components, distortions in temporal envelope

5. Multiband dynamic range compression General scheme of processing Spectral components of the input signal divided in multiple bands and the gain for each band calculated on the basis of signal power in that band. Parameters (band specific): compression threshold Th, compression ratio CR, attack & release time for detection.

6. Lippmann et al. (1980): 16-channel compression • 9% improvement in recognition score over linear amplification. • Asano et al.(1991): Multiband dynamic range compression realized as a single time-varying FIR filter & implemented on a 32-bit DSP fixed-point processor • Less spectral distortion due to smoothened frequency response of FIR filter. • Stone et al. (1999): Comparison of single and four-channel compression schemes & effect of varying CR, Th, and attack & release times • Intelligibility & quality tests showed no specific preference for schemes. • Li et al. (2000): Wavelet-based compression (7 octave sub-band analysis using wavelet filter bank & resynthesis after applying a logarithmic compression on the wavelet coefficients) • Increase in intelligibility without introducing noticeable distortions. • Magotraet al. (2000): Multiband dynamic range compression using a 16-bit fixed-point processor • Taylor's series approximation used for the compression function to reduce computations in gain calculation.

7. Disadvantages of multiband dynamic range compression • Spurious spectral distortions • Reduction in spectral contrasts and modulation depth • Distortion in spectral shape of formants lying across the band boundaries • Distortion of formant transitions across the adjacent bands • Time-varying magnitude response without corresponding variation in the phase response leading to quality degradation • → Audible distortions, perceptible discontinuities, adverse effect on the perception of certain speech cues

8. Example of distortion due to multiband dynamic range compression during spectral transition Swept sinusoidal input: constant amplitude, 125 –250 Hz linearly swept frequency, 200 ms sweep duration Time (s) Processed output: multiband compression with 18 auditory critical bands, CR = 30, Ta = 6.4 ms, Tr = 192 ms Time (s)

9. Investigation objective • Real-time dynamic range compression to compensate for frequency-dependent loudness recruitment associated with sensorineural hearing loss for use in hearing aids with a low-power processor. • Low distortions • Low computational complexity & memory requirement • Low signal delay (algorithmic + computational)

10. Sliding-band compression • Proposed for significantly reducing the temporal and spectral distortions associated with the currently used single-band and multiband compressions in hearing aids. • Realized with computational complexity acceptable for implementation on a 16-bit fixed-point DSP processor and signal delay acceptable for real-time application. • Investigations using offline & real-time implementations • Selection of processing parameters • Evaluation of the implementations • Informal listening, PESQ measure

11. Sliding-band Dynamic Range Compression Processing Applying a frequency-dependent gain function, with the gain for each spectral sample determined by the short-time power in auditory critical bandwidth centered at it & in accordance with the specified hearing thresholds, compression ratios, and attack and release times. Processing steps • Short-time spectral analysis: windowing, zero-padding, DFT calculation • Spectral modification: gain calculation, output spectrum calculation • Resynthesis: IDFT calculation, windowing, overlap-add

12. Short-time spectral analysis: windowing (length L, shift S), zero-padding, N-point DFT Spectral modification Pmc(k): Power at upper comfortable listening level CR(k): Compression ratio Resynthesis: N-point IDFT, overlap-add

13. Gain calculation • Auditory critical bandwidth • BW(k) = 25 + 75(1 + 1.4f 2)0.69, freq. sample = k, freq. = f • Target gain calculation • Power at upper comfortable listening level: Pmc(k) • Compression ratio: CR(k) • Input power: Pic(k), Output power: Poc(k) • Target gain: Gt(k) = Poc(k) / Pic(k) • Compression relation • dB scale: [Poc(k) / Pmc(k)]dB = [Pic(k) / Pmc(k)]dB / CR(k) • linear scale: Poc(k) / Pmc(k) = [Pic(k) / Pmc(k)]1/ CR(k) • Target gain for kth spectral sample • [Gt(k)]dB = [1 − 1 / CR(k)] [Pmc(k) / Pic(k)]dB

14. Gain calculation (contd.) Gain changed in steps from the previous value towards the target value with settable attack and release times Number of steps during attack phase = sa Number of steps during release phase = sr Target gain corresponding to min. input level = Gmax Target gain corresponding to max. input level = Gmin Gain ratio for attack phase γa = (Gmax/ Gmin)1/sa Gain ratio for release phase γr = (Gmax/ Gmin)1/sr Gain for ith window & kth spectral sample G(i,k) = max[G(i − 1,k) / γa, Gt(i,k)] for Gt(i,k) < G(i − 1,k) min[G(i − 1,k) γr, Gt(i,k)] for Gt(i,k) > G(i − 1,k) Attack time Ta = saS/ fs, Release time Tr = srS/ fs [fs= sampling freq., S = window shift] Fast attack: to avoid the output level from exceeding UCL during transients Slow release: to avoid the pumping effect or amplification of breathing

15. Implementation related challenges • Modifications in the short-time magnitude spectrum without corresponding changes in the phase spectrum can cause audible distortions. • Computational complexity: log or series approximation based gain calculations not suitable for use in sliding-band compression. Solution • Analysis-synthesis using least-square error based signal estimation from modified STFT (Griffin & Lim, 1984): Processing artifacts reduced by masking the effect of phase discontinuities in the modified short-time complex spectrum. • Look-up table based gain calculation: Two-dimensional look-up table relating the input power with gain as a function of frequency. • Reduces computations for real-time implementation. • Permits compression function most suited to compensate for the abnormal loudness growth.

16. Offline & Real-time Implementations • Implementation for offline processing • Implementation using Matlab 7.10 for evaluating the performance of the proposed technique and the effect of processing parameters. • Processing parameters • ◦ fs= 10 kHz ◦ Frame length = 25.6 ms (L = 256) • ◦ Overlap = 75% (S = 64) ◦ FFT size N =512 • 2D look-up table for frequency-dependent compression based on a linear relation between input-dB and output-dB, with settable CR(k) and Pmc(k). • ◦ Input range: 20 log intervals (trade-off: small gain increments, look-up table size). • ◦ Look-up table with 256×20 entries • Attack and release times • ◦ sa=1, Ta = 6.4 ms: Fast attack to avoid uncomfortable level during transients • ◦ sr=30, Tr = 192 ms: Slow release to avoid pumping & amplification of breathing

17. Implementation for real-time processing • Implementation on a 16-bit fixed-point DSP board to examine suitability of the technique for use in hearing aids. • DSP chip: TI/TMS320C5515 • 16 MB memory space (320 KB on-chip RAM with 64 KB dual access, 128 KB on-chip ROM) • Three 32-bit programmable timers • 4 DMA controllers each with 4 channels • FFT hardware accelerator (up to 1024-point FFT) • Max. clock speed: 120 MHz • DSP Board: eZdsp • 4 MB on-board NOR flash for user program • Stereo codec TLV320AIC3204: 16/20/24/32-bit ADC & DAC, 8 – 192 kHz sampling • Software development: C using TI's 'CCStudio ver. 4.0

18. Implementation • Input-output operations: DMA based I/O with cyclic buffers • ADC and DAC: one codec (left channel) with 16-bit quantization • Processing parameters (same as for offline processing): fs= 10 kHz, L = 256, S = 64, N =512 • Data representation (input samples, spectral values, processed samples): 16-bit real & 16-bit imaginary

19. Data transfers & buffering operations (S = L/4) • DMA cyclic buffers • 5-block S-sample input buffer • 2-block S-sample output buffer • Pointers • Current input block • Just-filled input block • Current output block • Write-to output block • (incremented cyclically on DMA interrupt) • Signal delay • Algorithmic: 1 frame (25.6 ms) • Computational ≤ frame shift (6.4 ms)

20. Test Results • Tests for verification and evaluation • Offline processing • Verification of the compression technique for speech input with a large level variation and examination of the effect of different set of processing parameters. • Assessment of output speech quality (using informal listening) for different input speech materials and time varying levels. • Comparison of distortions introduced by different compression techniques during spectral transitions. • Real-time processing • Comparison of the processed outputs from offline & real-time implementation: informal listening, PESQ measure (0 – 4.5). • Signal delay & computational requirement.

21. Results from offline processing Example: "you will mark ut please" concatenated with scaling factors for variation in the input level. CR = 2, Ta = 6.4 ms, Tr = 6.4 & 192 ms. Input waveform Scaling factor Unprocessed waveform Processed Tr= 6.4 ms, low Pmc Processed Tr= 192 ms, low Pmc Processed Tr= 6.4 ms, highPmc Processed Tr= 192 ms, highPmc Time (s) • Processing of different speech materials with varying levels: No audible roughness or distortion during informal listening.

22. Distortions during spectral transitions: Example of swept sinusoidal input. Input: constant amplitude, 125 –250 Hz linearly swept frequency, 200 ms sweep duration Single-band compression output Multiband compression (18 auditory critical bands) output Sliding band compression output CR = 30, Ta = 6.4 ms, Tr = 192 ms. Time (s)

23. Results from real-time processing Example: "you will mark ut please" concatenated with scaling factors for variation in the input level. CR = 2, Ta = 6.4 ms, Tr = 192 ms, low Pmc. Unprocessed waveform Offline processed waveform Real-time processed waveform Time (s) Informal listening: real-time output perceptually similar to the offline output PESQ for real-time w.r.t. offline : 3.5 Signal delay = 36 ms Use of processing capacity: 41% (lowest proc. clock for satisfactory operation = 50 MHz, max. clock = 120 MHz)

24. Summary & Conclusions • Sliding-band dynamic range compression presented to compensate for frequency-dependent loudness recruitment associated with sensorineural hearing loss without introducing the distortions associated with single-band & multiband compression. • Realized using modified fixed-frame analysis-synthesis for low computational complexity & without distortions associated with phase discontinuities. • Suitable for speech & non-speech audio & provision for settable attack time, release time, & compression ratios. • Implemented using 16-bit fixed-point DSP chip & tested for satisfactory operation: 36 ms signal delay, 41% use of processing capacity, indicating scope for combination with other processing techniques.

25. Thank You To be continued to Part C.