28

Does the Nyquist frequency of the Cochlear nerve impose the fundamental limit on human hearing? No. A quick run-through the human auditory system: The outer ear (pinnae, ear canal), spatially "encodes" the sound direction of incidence and funnel the sound pressure towards the ear drum, which converts sound into physical motions, i.e. mechanical ...


21

It is actually not distorted, it is sampled at high enough rate. What fools you is the straight lines drawn between sample points, it gives you a false impression of the waveform. It shows you a linear interpolation of the signal. It does not represent how the signal would actually look like. A sampled signal exists only at the sample points, and to convert ...


15

Figure 1.(c) shows the Test image reconstructed from MAGNITUDE spectrum only. We can say that the intensity values of LOW frequency pixels are comparatively more than HIGH frequency pixels. Actually, this is not correct. The phase values determine the shift in the sinusoid components of the image. With zero phase, all the sinusoids are centred at the same ...


13

The FFT can only be performed over a limited chunk of data. The basic math is based on the assumption that the time domain signal is periodic, i.e. your chunk of data is repeated in time. That typically results in a major discontinuity at the edges of the chunk. Let's look at a quick example: FFT size = 1000 points, Sample Rate = 1000 Hz, Frequency ...


11

Phase Noise and Frequency Noise are not two different noise sources, they are artifacts of the same noise, it is just a matter of what units you want to use. Frequency and Phase are directly related as frequency is phase changing with time, so if you have one you will always have the other; frequency and phase are related by derivatives and integrals: the ...


10

One trick, for even-length signals, is what to do with the "middle" sample. Here, I've split it half and half between each side of the FFT. The other trick is to ensure that you have the right amplitudes in the resampled signal. Here's it's a factor of 2. Try this in scilab: x = rand(1,100,'normal'); X = fft(x); XX = 2*[X(1:50) X(51)/2 zeros(1,99) X(51)/...


10

Well, first of all the Sound Level Pressure decreases by $6 \; \mathtt{dB}$ when doubling the distance - this plays a big role. We do also have sound attenuation coming from our medium - air. Let's take a closer look onto sound absorption coefficient for different frequencies: Knowing that human speech is mostly concentrated at the range of $300\;\mathtt{Hz}...


9

This comes from music terminology. The name "octave" comes from the fact that in the heptatonic musical scales (which are the prevalent scales in western music), the note with a 2:1 frequency ratio is the eighth note in the scale. For example, in the C major scale (C D E F G A B C) the eighth note is one octave above / has a 2:1 frequency ratio with the ...


9

The term Doppler Shift is actually a bit of a misnomer. The frequencies are not actually shifted but they are scaled (see http://fourier.eng.hmc.edu/e101/lectures/handout3/node2.html for definition of shifting vs. scaling). It's a relative change not an absolute one. Both time and frequency domains are scaled: when the source is moving towards you, the ...


9

The actual requirement is to sample at GREATER then twice the bandwidth, not at a rate equal to it... So only your 80Hz same set actually meets the requirement, because the 60Hz case is ambiguous in general, consider if you were sampling sin (2PiFt) instead then you would get a flat line at zero amplitude.... And changing the angle between sin and cos would ...


8

The OP's opening statement is incorrect: $f_s > f_{max}/2$ prevents frequency aliasing for a bandlimited signal, but not amplitude aliasing $f_s > 2 f_{max}$ prevents aliasing. It's as simple as that. There is no such distinction as "amplitude aliasing". Since the OP has stated the signal is band-limited; as long as we can assume that means ...


7

The HUP follows directly from the properties of the Fourier Transform, because time and frequency are orthogonal bases in which we can expand the co-efficient sequence of our signal. In fact all pairs of orthonormal bases will have some kind of Uncertainty Principle associated with them. In traditional Fourier analysis, the either the time axis or the ...


7

For the source, go to end of the answer Suppose one day you got one note which has some thing written to it, say "Major frequency components are 10 Hz, 25Hz, 50 Hz and 100 Hz". Somehow, you understood that its time-series representation is a very important thing (may be master-piece work of a great musician, or some national security matter, anything). So ...


7

Note: I originally posted this answer for the Stack Overflow copy of this question, before realizing that it had also been asked here. It somewhat duplicates pichenettes' answer, but I felt it still worth (re)posting here, since it includes some extra details. (Whether those details are useful or not, I'll leave for you and the OP to judge.) If you know ...


7

From the ones I've been using I can recommend: YAAFE - very pleasant to work with in Python ESSENTIA - another one I like particularly due to Python integration aubio FEAPI Aquila - friend of mine used it extensively and he likes it a lot Recently I came across this paper and I believe that this should perfectly answer your question. Moffat D. et al - ...


7

Consider a liner discrete-time system. Assume we can define it in terms of an input-output relation as follows (you can assume a more general model but it is enough for our purpose): $$a_0y[n]+a_{1}y[n-1]+\cdots+a_{N}y[n-N]=b_0x[n]+b_{1}x[n-1]+\cdots+b_{M}x[n-M]\tag{1}$$ When the coefficients $\{a_i\}$ and $\{b_i\}$ are constant, we call it a finite-order ...


7

A zero coefficient at DC simply means the mean of the waveform is 0; bin 0 of the DFT is identical to calculating the mean value scaled by $N$. Consider the general expression for the DFT: $$X(k) = \sum_{n=0}^{N-1}x(n)e^{-j k \omega_n n}$$ with $\omega_n = 2\pi/N$ and $k = 0,1 \ldots N-1$ (This is a correlation of the waveform x(n) with every integer ...


6

A signal which is the sum of sinusoids with different frequencies $f_i$ is only periodic if $$f_i=k_if_0,\quad i=1,2,\ldots$$ is satisfied for some $f_0$ and integer values $k_i$. Or, in other words, the frequencies $f_i$ must be rational multiples of each other. In your example, the signal is periodic but its fundamental frequency is not $f_1=800\,\text{...


6

The theory behind sweep-sine measurements of LTI systems requires a signal with constantly changing the frequency. You cannot simply playback few tones - the whole frequency range is necessary. So that if you want to identify your system with the impulse response $h[n]$, you feed the sweep sine signal $s[n]$ into it and record the output. Obviously output ...


6

The article Handling Spectral Inversion in Baseband Processing that you refer to in your question is not about simple inversion of the frequency axis, but it is about inversion of the frequency axis and conjugation in the frequency domain. This is necessary because due to different mixing conventions in several conversion stages, the upper and lower side-...


6

You need to build a time varying delay, where you can modulate the delay amount over time. The peak delay modulation is a function of your maximum desired frequency shift and the modulation frequency. This is not trivial since it will require fractional sample delays with some kind of interpolation algorithm. You can't round to the nearest integer delay ...


6

In digital signal processing, one almost never deals directly with high-frequency signals. The reason is that the frequency band that a signal occupies is completely irrelevant to the information that the signal conveys; the information exists in the envelope of the signal, which usually can be handled either by a generic CPU, a DSP, FPGA, or ASIC. So, let'...


6

I think it means the "apparent" frequency of oscillation. It's poorly worded. What's happening is that the frequency is increased to the Nyquist rate in the first half and then above it, causing aliasing, in the second. What looks like the frequency increases in the first half and decreases in the second because of aliasing.


6

$\Omega$ is the usual angular frequency in radians per second, and is equal to $2 \pi f$. It is the way to measure frequency for continuous-time signals. In discrete-time, frequency is measured in radians per sample, and is denoted as $\omega$. Here, a frequency component with $\omega = 2\pi$ is indistinguishable from $\omega = 0$. $\Omega$ and $\omega$ is ...


6

There is no aliasing as 𝑓 = 30 Hz is less than or equal to the folding frequency, 30 Hz and 40 Hz, respectively. Yes and no. There isn't significant aliasing when you're sampling at 80Hz, because the resulting signal has frequency components at 30Hz and 50Hz. The result is thus unambiguous as long as you take that 50Hz signal into account. There is ...


5

The usual way to estimate the amplitude of a particular frequency is to use the Goertzel algorithm. There is a good write-up by Rick Lyons here. Even though Rick's writeup is about single tone detection, it can be applied when multiple tones are present, too.


5

The correct way to group multiple bins together is to multiply each complex fft bin output by its complex conjugate (which gives the bin power) then add all the bin powers together and divide by the number of bins in the group. If you want to display in db (which is the conventional approach) then take 10*log10() of the result. Negative db values are normal ...


5

Nyqvist-shannon sampling theorem only says that your sampling frequency should be greater than twice the bandwidth of the signal and not twice the carrier frequency of the signal. Hence you can modulate your 20 MHz signal at any carrier frequency you need and still get back your original signal by sampling at greater than 40 MHz. Hence you need not sample ...


5

In your line mag1=abs(f1); you are leaving the total intensity of the image unchanged (test this by summing up the intensities over all pixels). Rejecting the phase information in Fourier space just leads to a spatial redistribution of the intensity in real space, such that r1 will have the same total insity as i1. In your line phase1=angle(f1); you are ...


5

If the DFT is the Uniform sampling from $ 0 $ to $ 2 \pi $ then the first bin is given by: $$ x[k] = \sum_{n = 0}^{N - 1} x[n] $$ Namely it is the sum of all the samples. Hence in order to remove the DC (Mean) all you need is a filter which has zero in its DC bin. Since, the filtered signal, which is a convolution (Circular) of the DFT of the input signal ...


Only top voted, non community-wiki answers of a minimum length are eligible