Generation of Ultrasound based on the frequency response characteristics of the ‘Koss Pro Headphone’ with R. David Case Sound wave files’ – A Case study

Generation of Ultrasound based on the frequency response characteristics of the ‘Koss Pro Headphone’ with R. David Case Sound wave files’ – A Case study

1R. David Case and 2D. Chuckravanen

1Inventor of the R. David Case Sound device and Research Collaborator, USA.
2 Research Scientist, British Association of Cognitive Neuroscience, Department of Psychology, University of Essex,
Colchester CO4 3SQ, United Kingdom

Abstract: R. David Case has been able to generate unique sounds that are able to treat people suffering from tinnitus successfully. There are currently many positive feedbacks from all those people who have listened to those unique sounds using specifically headsets or headphones of the Koss models (ktx Pro, ksc-75) and these discoveries have encouraged R. David Case to pursue positively his journey into finding out what is special about those sounds when being listened specifically with Koss models. In this research, we focus mainly on the technical aspects of R. David Case sound signals which he has recorded, and we will analyse this sound in the time domain, frequency domain as well as the effect of using the Koss Pro headphone frequency response characteristics applied to the sound. Results obtained from the analysis demonstrated that there is the generation of ultrasound waves which could be the underlying reason for the treatment of the tinnitus.

1. Introduction

Generally speaking, hearing depends on a sequence of complicated stages that change sound waves in the air into electrical signals. Then, the auditory nerve transports these signals to our brain. The sound waves passes through the outer ear and then it travels through the ear canal leading to the eardrum. Subsequently, the eardrum vibrates owing to the incoming sound waves and then sends these vibrations to tiny bones that are found in the middle ear. These bones amplify the sound vibrations which are then sent to the cochlea (snail-shaped structure) filled with fluid in the inner ear. An elastic partition runs from the beginning to the end of the cochlea, splitting it into an upper and lower part. This partition is called the basilar membrane because it serves as the base, or ground floor, on which key hearing structures sit. Once the vibrations cause the fluid inside the cochlea to ripple, a traveling wave forms along the basilar membrane. Hair cells—sensory cells sitting on top of the basilar membrane—ride the wave. Hair cells near the wide end of the snail- shaped cochlea detect higher-pitched sounds, such as an infant crying. Those closer to the center detect lower-pitched sounds, such as a large dog barking. As the hair cells move up and down, microscopic hair-like projections (known as stereocilia) that perch on top of the hair cells bump against an overlying structure and bend. Bending causes pore-like channels, which are at the tips of the stereocilia, to open up. When that happens, chemicals rush into the cells, creating an electrical signal. The auditory nerve carries this electrical signal to the brain, which turns it into a sound that we recognize and understand.

The normal range of human hearing is 20 Hz to 20,000 Hz. Under ideal lab conditions, human can hear sound as low as 12 Hz and as high as 28 kHz, though the threshold increases sharply at 15kHz in adults, corresponding to the last auditory channel of the cochlea. The next sub-section explains what is ultrasound and why it is important for this research work.

1.1 What is UltraSound?

Ultrasound is acoustic (sound) energy in the form of waves having a frequency above the human hearing range. The highest frequency that the human ear can detect is approximately 20 thousand cycles per second (20,000 Hz). This is where the sonic range ends, and where the ultrasonic range begins. In physics the term “ultrasound” applies to all acoustic energy with a frequency above human hearing (20,000 hertz or 20 kilohertz). Typical diagnostic sonographic scanners operate in the frequency range of 2 to 18 megahertz, hundreds of times greater than the limit of human hearing. Ultrasound is used in electronic, navigational, industrial, and security applications. It is also used in medicine to view internal organs of the body.

Since the inception of ultrasound as a therapeutic tool, its potential to treat disorders throughout the brain has been explored. The objective has been to utilize ultrasound’s focusing ability to target precisely within deep tissues, affecting only the interested volume while leaving all other structures unaltered (Hynynen and Clement, 2007). Varying the ultrasound parameters not only allows ablation of pathological tissue, such as brain tumours and metastases, or silencing of dysfunctional neuronal circuits, but also opens up the blood–brain barrier for targeted drug delivery and modulation of neural function (Martin and Werner, 2013). Before we go directly into analysing those unique sound waves generated by a ‘device’ built by R. David Case, it is important also to define what is tinnitus and how it affects a lot of people worldwide in the next sub-section.

1.2 Tinnitus

Tinnitus usually occurs when one experiences ringing or other noises in one or both ears (See Figure 1a). The noise you hear when you have tinnitus is not caused by an external sound, and other people usually cannot hear it. In fact, tinnitus is a common problem. It affects about 15% to 20% of people. Tinnitus is usually caused by an underlying condition, such as age-related hearing loss, an ear injury or a problem with the circulatory system. For many people, tinnitus improves with treatment of the underlying cause or with other treatments that reduce or mask the noise, making tinnitus less noticeable. In this research, we have investigated how the sound files when listening from Koss Pro headphones can really pose as Anti-tinnitus. Next section describes the methodology employed to analyse the sound wave signals as well as the effect of the frequency response characteristics of the Koss Pro headphones on R. David Case sounds.

Figure 1a: Tinnitus and the types of ringing 2. Research Methodology

All analyses are conducted on Scientific Matlab R2009a Platform running on a 64-bit computer system architecture with Intel Core I7 and processing speed is 2.90GHz. Sound files were recorded by Mr. R. David Case and they are in *.wav format and the total sound file is about 800 megabytes. All sound files (stereo sounds) were imported to the Matlab environment for analyses purposes. Spectral analyses are conducted on those sound waves, and also the effect of the frequency responses of the specific Koss pro headphone on those sound files are investigated too. Experiment is being conducted on Koss Pro headphone device with the following specifications: 15 Hz to 25000Hz, 60 Ohms, 101 dB SPL, 4 ft/1.2m. 3.5 mm (See figure below). In addition, follow-ups were conducted on all the people who utilised such sounds to improve their health conditions such as tinnitus, or vibrations in certain parts of the body or to reduce ‘voices’ that only them can hear.

Figure 1b: Koss Pro Headphone under analysis for frequency response using Matlab Software R2009a.

3. Results

3.1 Time-domain presentation of the R. David Case sound file

The sound wav file was imported to the Matlab environment so as the sound wav data can be used for further analyses as the sound data is in a matrix format. Once the sound wav was in the Matlab workspace, the following figure was generated in order to visualise in the time domain the sound wav. The sampling frequency of the sound wave is 44.1 kHz.

Figure 2: Sound 1 component of the device stereo sound wav (Amplitude vs. Time)

Figure 3: 2nd sound component in the stereo sound wav data (Amplitude vs. Time)

In contrast to Figure 2 which shows that the sound is varying randomly (where we need to see what is happening in frequency domain), in Figure 3, the amplitude of the sound changes abruptly at approximately half of the time period.

3.2 Frequency domain presentation of R. David Case sound file.

As shown in Figure 4, when the spectral analysis is applied to the component 1 of the Sound wav, there were peaks at around 0, 0.4, and 0.8 x 104 Hz. On the other hand, when the spectral analysis is applied to the component 2 of the sound wav, there were peaks at around 0 Hz, 0.20, 0.25, 0.4 and 0.9 x 104 Hz (Figure 5).

Figure 4: Spectral analysis with Frequency centred at f=0Hz for sound component 1

Figure 5: Power vs. Frequency spectrum (for sound component 2 of the recorded R. David Case sound wave. Frequency centred at f=0Hz for Sound component 2

3.3 Characteristic frequency response of Koss Pro headphone

From Figure 7, it is clearly observed that the power-frequency spectrum is not constant and varies throughout the normal hearing frequency spectrum. It decreases from 0 Hz to 4000 Hz and then the power increases at around 5000 Hz and again the power decreases to a minimum at about 7000 Hz and then increases to a plateau effect at around 10000 Hz and decreases till to 15000 Hz.

Figure 7: The Power (W)-frequency (Hz) response of the Koss Pro headphone

3.4 Convolution of the Sound wave signals with the frequency characteristics of the Koss
pro Headphone 3.4.1 Application of Convolution

One of the most important concepts in Fourier theory is that of a convolution. Mathematically, a convolution is defined as the integral over all space of one function at x times another function at u-x. The integration is taken over the variable x (which may be a 1D variable), typically from minus infinity to infinity over all the dimensions. So the convolution is a function of a new variable u, as shown in the following equations. The cross in a circle is used to indicate the convolution operation (Urynbassarova et al, 2017; McGillem et al., 1984; Proakis et al, 1996).

Note that it does not matter which function to take first, i.e. the convolution operation is commutative. The function f(x) represents the function of the sound wave and the function g(x) represents the Headphone frequency response characteristics.

The positive working of the sound wave coupled with the Koss headphones produce a frequency which is healing the inner organ of the human body and act as an anti-tinnitus, and therefore convolution algorithm is applied to the sound wav and also the frequency response characteristics of the Koss headphone to produce a new convolution spectrum elaborated in the next sub-section.

3.4.2 Analysed work

It was interesting to find that the convolved signal of the frequency response of the Koss headphone and the Sound component 1 did not show anything visually significant. However, when the convolution method was applied to the second sound wave component, the following spectrum figure is obtained.

Figure 8: Convolved signal (Sound wave signals and the headphone frequency response characteristics)
It is clearly observed that there are maximum peaks at the frequency 2.822 x 106 Hz is created from the machinery of the sound device and the headphone and this ultrasonic wave frequency needs to be further investigated.

3.5 Experiment repeated with a 80 minute sound wave file

The sound wave signals were brought to Matlab environment and as it was large size (memory), my program loops through every 2 minutes of the sounds in order to depict the various frequency components (frequency spectrum), the maximum frequency component also as well as the sounds through time. There was decrease and increase in the signal throughout the duration of the sound recordings. We have found the following peaked frequencies for each sound or rather ultra-sound component.

(i) Ultra-Sound component 1: 2.9949 x 106 Hz
(ii) Ultra-Sound Component 2: 2.8219 x 106 Hz

Figure 9: Frequency spectrum analysis demonstrated a maximum peak at 2.9949 MHz for the sound component 1 and this is the average for every 2 minutes of the sound wave file.

Figure 10: Frequency spectrum analysis demonstrated a maximum peak at 2.8219 MHz for the sound component 2 and this is the average for every 2 minutes for the sound wav file.


This research analysed the sound waves generated from a sound device by R. David Case and the sound is called after him which is R. David Case sound wave file. The sound has the capability to cure or accelerate the healing process or reverse the process of tinnitus. It is posited that there may be some frequencies that are found in the generated sound waves which could do this.

Based on various feedbacks from people who were administered this sound (at their own consents) described the sounds as follows. One of the people (subject 1) who listened to the long sound wave at low volume on and off state that he does not need Valium to circumvent the ringing in the ear and he states that the amplitude of the tinnitus seems to fluctuate high and low with improved cognitive concentration. Another person (subject 2) who listened to the sound continuously for couple of days explained that he is able to do his day to day activities in peace and he states that he is cured from the tinnitus. Another participant (subject 3) posits that the sound waves have the capability to eliminate vibrations that were occurring on chest area, and then torso and also he is looking forward to remove occurrences from his lower body.

Another participant (initial B.H, subject 4) claimed that the sound therapy he received from this peculiar sound really topped his hyperacusis in minimum duration (couple of weeks) and this sound has reduced his tinnitus greatly. Currently, he strongly wants to get rid of the tinnitus by listening to the sounds more consistently than he did before. The next participant (Subject 5) realises that after 7 consecutive days, after listening 6 hours of the sounds daily, he has observed that his right ear which normally can hear up to 10,400 Hz can now hear up to 14,600 Hz. His specialist advised him to continue using such sound device as his tinnitus has greatly improved as well as his hearing capabilities. The last but not the least, a participant (Subject 6) found that he can sleep better at night after listening to the sounds and at low volume.

Based on the digital signal analysis of the sound wav and the headphone frequency characteristics, it was found that there is a much higher ultrasound frequency which is being generated. This, in turn is causing a positive effect on those people suffering from tinnitus.


In this case study research, it is demonstrated how the unique sound wave comprising of both stereo and mono-waves generated by R. David Case coupled with the headphone characteristics (based on its Power-frequency response) can in fact alleviate the tinnitus problem. The ultrasound frequencies, as obtained in this research, can reverse the process of creating tinnitus and therefore it heals the human beings suffering from it. Based on the different interviews and feedbacks, it is clear that the sound wav and the headphone can be treated as a medical sound treatment or medical ultrasound treatment. Further research will look at analysing a larger sample of subjects/participants with their consent and observe if there are no long-term side effects. Couple with that, we need to devise an anti-tinnitus administration protocol and state clearly in that protocol the duration of listening this sound, and at what ideal amplitude, should the sound be played.


Martin, E. and Werner, B. (2013), ‘Focused Ultrasound Surgery of the Brain’, Current Radiolology Reports, 1, pp. 126-135.
K. Hynynen & G. Clement (2007) Clinical applications of focused ultrasound—The brain, International Journal of Hyperthermia, 23:2, 193-202, DOI: 10.1080/02656730701200094.
Didar Urynbassarova, Bing-Zhao Li, Zhi-Chao Zhang (2017), A Convolution Theorem for the Polynomial Fourier Transform’ , IAENG International Journal of Applied Mathematics 47(4):381- 387
McGillem, Clare D.; Cooper, George R. (1984). Continuous and Discrete Signal and System Analysis (2 ed.). Holt, Rinehart and Winston. p. 118 (3-102). ISBN 0-03-061703-0.
Proakis, John G.; Manolakis, Dimitri G. (1996), Digital Signal Processing: Principles, Algorithms and Applications (3 ed.), New Jersey: Prentice-Hall International, p. 297.
Pienkowski M, Tyler RS, Roncancio ER, Jun HJ, Brozoski T, Dauman N, Coelho CB, Andersson G, Keiner AJ, Cacace AT, Martin N, Moore BC (2014). “A review of hyperacusis and future directions: part II. Measurement, mechanisms, and treatment”. American Journal of Audiology. 23 (4): 420–36. doi:10.1044/2014_AJA-13-0037. PMID 25478787. S2CID 449625.
McCombe A, Baguley D, Coles R, McKenna L, McKinney C, Windle-Taylor P (2001). “Guidelines for the grading of tinnitus severity: the results of a working group commissioned by the British Association of Otolaryngologists, Head and Neck Surgeons, 1999”. Clinical Otolaryngology and Allied Sciences. 26 (5): 388–93. doi:10.1046/j.1365-2273.2001.00490.x. PMID 11678946. Archived (PDF) from the original on 2017-09-24.
Langguth B, Goodey R, Azevedo A, et al. (2007). “Consensus for tinnitus patient assessment and treatment outcome measurement: Tinnitus Research Initiative meeting, Regensburg, July 2006”. Tinnitus: Pathophysiology and Treatment. Progress in Brain Research. 166. pp. 525–36. doi:10.1016/S0079-6123(07)66050-6. ISBN 978-0444531674. PMC 4283806. PMID 17956816.
Meikle MB, Stewart BJ, Griest SE, et al. (2007). “Assessment of tinnitus: Measurement of treatment outcomes” (PDF). Tinnitus: Pathophysiology and Treatment. Progress in Brain Research. 166. pp. 511–21. doi:10.1016/S0079-6123(07)66049-X. ISBN 978-0444531674. PMID 17956815. Archived (PDF) from the original on 2017-09-25.


%copyright by Dr. D. Chuckravanen and Mr R. David Case %uiopen(‘C:\Users\User\Documents\MATLAB\device.wav’,1); figure(1), plot(data(:,1)), title(‘sound1′);
figure(2), plot(data(:,2),’g’), title(‘sound2’);
y1 = fft(data(:,1));
y2 = fft(data(:,2));
n = length(y1);
f = (0:n-1)*(fs/n);
power1 = abs(y1).^2/n;
axis([-100 50000 0 40]);
y01 = fftshift(y1);
f0 = (-n/2:n/2-1)*(fs/n); %
power01 = abs(y01).^2/n;
axis([-100 50000 0 40]);
power2 = abs(y2).^2/n;
axis([-100 50000 0 40]);
% number of samples
% frequency range
% power of the DFT
% shift y values
0-centered frequency range
% 0-centered power
% power of the DFT

y02 = fftshift(y2); % shift y values
f0 = (-n/2:n/2-1)*(fs/n); % 0-centered frequency range
power02 = abs(y02).^2/n;
axis([-100 50000 0 40]);
% 0-centered power
% creation of the frequency response of the headphones
dataheadphone=[20 35
30 38
40 42
50 43.5
60 44
70 44.2
80 44.5
90 44.2
100 44
200 42.5
300 41
400 40
500 39
600 37.5
700 37.0
800 36.5
900 36.3
1000 36
2000 35.5
3000 37
4000 35.1
5000 40
6000 30
7000 26
8000 30
9000 34.5
10000 34.5
15000 22.5];
%convert db to power
data11=dataheadphone(:,2)*0.001; dataheadphone1= [data11 dataheadphone(:,1)];
plot(dataheadphone1(:,2),dataheadphone1(:,1)) title(‘frequency response Power vs Frequency for headphone’); xlabel(‘frequency’);

%convolving the frequency response of both sound components with the %headphone characteristical specifications
C1= conv(power01,dataheadphone1(:,1)); figure(8), plot(C1)
C2= conv(power02,dataheadphone1(:,1)); figure(9), plot(C2)