Confirmation of presence of narrowband interference by harmonic analysis
One or more processing units confirm existence of narrow band interference in a signal by using an estimate f of the frequency, to check for one or more harmonics. In illustrative embodiments, the estimate f is automatically identified as a second harmonic if a predetermined criterion is satisfied by the signal (in the frequency domain) at either of two frequencies namely (A) frequency f/2 and (B) frequency (M−f)/2 and whichever of these two frequencies is stronger is identified as the fundamental frequency. In several such embodiments, the estimate f is automatically identified as a third harmonic if a predetermined criterion is satisfied by the signal (in the frequency domain) at any of three frequencies namely (C) frequency f/3 and (D) frequency (M−t)/3 and (E) frequency (M+f)/3. If the predetermined criteria are not met at all five frequencies (A)-(E) then f is identified as the fundamental frequency.
1. A method of processing a signal, the method comprising:
retrieving from a computer-readable memory, an estimate f of a frequency of narrowband interference in said signal:
using at least one processing unit coupled to the computer-readable memory to check if a first predetermined criterion is satisfied, by a first element in a frequency domain representation S[0]-S[M−1] of said signal, on at least one of (a) frequency f/q and (b) frequency (M−f)/q, wherein S is an absolute weighted average, and wherein q is an integer greater than 1 and M is a quantity of samples used in a frequency domain transformation; and
storing in said computer-readable memory coupled to said at least one processing unit, a flag indicative of a result of said checking.
2. The method of claim 1 further comprising: checking if a second predetermined criterion is satisfied by a second element in the frequency domain representation S[0]-S[M−1] of said signal on the other of (a) frequency f/q and (b) tency (M−f)/q.
3. The method of claim 2 wherein q is equal to 3 and the method further comprises: checking if a third. predetermined criterion is satisfied by a third. element in the frequency domain representation S[0]-S[M−1] of said signal on (c) frequency (M+f)/3.
4. The method of claim 2 wherein q is equal to 2 and the method further comprises:
repeating at least said using with q equal to 3 and additionally checking if a third predetermined criterion is satisfied by a third element in the frequency domain representation S[0]-S[M−1] of said signal at (c) frequency (M+f)/3; and
identifying f as a harmonic frequency if no values of said signal at respective frequencies f/2, (M−f)/2, f/3, (M−f)/3 and (M+f)/3 are found to satisfy corresponding predetermined criteria.
5. The method of claim 1 wherein said using comprises: computing a first approximation of an absolute value of each of elements S[A(f/q)] and S[A((M−f)/q)] in the frequency domain representation S[0]-S[M] of the signal; wherein A is a function to round an argument thereof to an integer.
6. The method of claim 5 wherein: said integer is closest to said argument of function A.
7. The method of claim 5 wherein said integer is smaller than said argument, and said function A is round-to-floor function, the method further comprising:
computing a second approximation of the absolute value of each of elements S[B(f/q)] and S[B((M−f)/q)] wherein B is a round-to-ceiling function; and
computing a weighted average of said first approximation and said second approximation at each of frequencies f/q and (M−f)/q based on a plurality of predetermined weights.
8. The method of claim 1 further comprising: identifying A(f/q) as a harmonic frequency depending on an outcome of the first predetermined criterion;
wherein A is a function to round an argument thereof to an integer.
9. A circuit for processing a signal, the circuit comprising at least a processor, the circuit comprising:
a memory coupled to the processor, the memory comprising an estimate f of a frequency of narrowband interference in said signal;
means, coupled to the memory, for checking if a first predetermined criterion is satisfied, by a first element in a frequency domain representation S[0]-S[M−1] of said signal, on at least one of (A) frequency f/q and (B) frequency (M−f)/q, wherein S is an absolute weighted average, and wherein q is an integer greater than 1 and M is a quantity of samples used in a frequency domain transformation; and
means for storing in said memory an indication of a result generated at least partially by said means for checking.
10. The circuit of claim 9 further comprising: means for checking if a second predetermined criterion is satisfied by a second element in the frequency domain representation S[0]-S[M−1] of said signal on the other of (A) frequency f/q and (B) frequency (M−f)/q.
11. The circuit of claim 10 wherein q is equal to 3 and the circuit further comprises: means for checking if a third predetermined criterion is satisfied by a third clement in the frequency domain representation S[ 0 ]-S[M−1] of said signal (C) frequency (M+t)/3.
12. The circuit of claim 10 wherein q is equal to 2 and the circuit further comprises:
additional means for checking with q equal to 3 and checking if a third predetermined criterion is satisfied by a third element in the frequency domain representation S[0]-S[M−1] of said signal at (C) frequency (M+t)/3; and
means for identifying f as a harmonic frequency if no values of said signal at respective frequencies f/2, (M−f)/2, f/3, (M−f)/3 and (M+f)/3 are found to satisfy corresponding predetermined criteria.
13. The circuit of claim 9 wherein said means for checking comprises:
means for computing a first approximation of an absolute value of each of elements S[A(f/q)] and S[A((M−f)/q)] in the frequency domain representation S[0]-S[M] of the signal; wherein A is a function to round an argument thereof to an integer.
14. The circuit of claim 13 wherein: said integer is closest to said argument of function A.
15. The circuit of claim 13 wherein said integer is smaller than said argument, and said function A is round-to-floor function, the circuit further comprising:
means for computing a second approximation of the absolute value of each of elements S[B(f/q)] and S[B((M−f)/q)] wherein B is a round-to-celling function; and
means for computing a weighted average of said first approximation and said second approximation at each of frequencies f/q and (M−f)/q based on a plurality of predetermined weights.
16. The circuit of claim 9 further comprising:
means for additionally checking if S[A(f/q)] is greater than S[A(M−f)/q)]; and
means for identifying A(f/q) as a harmonic frequency if a result generated by said means for additionally checking is true; wherein A is a function to round an argument thereof to an integer.