IP Library Granted Patent US 7,552,154
Granted Patent B2
US 7,552,154 · App. 11/109,061 · Granted Jun 23, 2009

System and method for statistically separating and characterizing noise which is added to a signal of a machine or a system

View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 7,552,154
App. No.
11/109,061
Granted
Jun 23, 2009
Kind
B2
Abstract

Method for finding the probability density function type and the variance properties of the noise component N of a raw signal S of a machine or a system, said raw signal S being combined of a pure signal component P and said noise component N, the method comprising: (a) defining a window within said raw signal; (b) recording the raw signal S; (c) numerically differentiating the raw signal S within the range of said window at least a number of times m to obtain an m order differentiated signal; (d) finding a histogram that best fits the m order differentiated signal; (e) finding a probability density function type that fits the distribution of the histogram; (f) determining the variance of the histogram, said histogram variance being essentially the m order variance σ 2 (m) of the noise component N; and (g) knowing the histogram distribution type, and the m order variance σ 2 (m) of the histogram, transforming the m order variance σ 2 (m) to the zero order variance σ 2 (0) , said σ 2 (0) being the variance of the pdf of the noise component N, and wherein the histogram type as found in step (e) being the probability density function type of the noise component N.

Claims (105)

1. A method for an apparatus to determine the probability density function type and the variance properties of the noise component N of a raw signal S of a machine or a system, said raw signal S being combined of a pure signal component P and said noise component N, the method comprising:

a. defining a window within said raw signal;

b. recording the raw signal S;

c. numerically differentiating the raw signal S within the range of said window at least a number of times m to obtain an m order differentiated signal;

d. finding a histogram that best fits the m order differentiated signal;

e. finding a probability density function type that fits the distribution of the histogram;

f. determining the variance of the histogram, said histogram variance being essentially the m order variance σ 2 (m) of the noise component N;

g. knowing the histogram distribution type, and the m order variance σ 2 (m) of the histogram, transforming the m order variance σ 2 (m) to the zero order variance σ 2 (0) , said σ 2 (0) being the variance of the probability density function of the noise component N, and wherein the histogram type as found in step (e) being the probability density function type of the noise component N; and

h. the apparatus outputting at least one of the zero order variance σ 2 (0) or the probability density function type.

2. The method according to claim 1 , which is repeatedly performed as the raw signal S progresses.

3. The method according to claim 2 which is performed in real-time.

4. The method according to claim 3 , wherein the one probability density function type is the Gaussian type.

5. The method according to claim 1 , wherein the probability density function type that fits the distribution of the histogram is the one from a list that found to be best fitting the distribution of the histogram.

6. The method according to claim 5 , wherein the list comprises only one probability density function type.

7. The method according to claim 1 , wherein the transform is performed by means of a specific expression suitable for the said fitted probability density function, wherein said specific expression is derived from the following general expression

F

n

i

(

m

)

=

D

s

{

j

=

1

m

S

j

m

f

(

ξ

)

}

ξ

(

m

)

.

8. The method according to claim 1 , wherein when the fitted probability distribution function is Gaussian, the transform is performed by means of the following specific expression

m

N

(

0

,

σ

0

2

)

i

m

=

α

(

m

)

N

(

0

,

β

(

m

)

σ

0

2

)

.

9. An apparatus for determining the probability density function type and the variance properties of the noise component N of a raw signal S of a machine or a system, said raw signal S being combined of a pure signal component P and said noise component N, the apparatus comprising:

a. a differentiating module, for receiving and numerically differentiating the raw signal S within the range of a predefined window at least a number of times m to obtain an m order differentiated signal;

b. a module for finding a histogram that best fits the m order differentiated signal;

c. a list containing at least one type of predefined probability density function;

d. a module for finding one probability density function type from said list that best fits the distribution of the histogram;

e. a module for determining the variance of the histogram, said histogram variance being essentially the m order variance σ (m) of the noise component N; and

f. a module for, given the histogram distribution type and them m order variance σ 2 (m) of the histogram, transforming the m order variance σ 2 (m) to the zero order variance σ 2 (0) , said σ 2 (0) being the variance of the probability density function of the noise component N, wherein the histogram type as found in step (d) being the probability density function type of the noise component N.

10. The apparatus according to claim 9 , wherein the apparatus components operate repeatedly to find the updated probability density function type, and the variance properties of the noise component as the signal S progresses.

11. A system for receiving a raw signal S which is combined from a pure signal P and a noise component N, and for outputting a signal which is essentially said pure signal, comprising:

a. an apparatus according to claim 9 for receiving said raw signal and outputting the probability density function and distribution type of said noise component into a filter; and

b. a filter receiving said raw signal and also receiving said probability density function and distribution type of said noise component from said apparatus, and given said received data, processing and outputting a signal which is essentially said pure signal.

12. The system according to claim 11 , wherein the filter is an adaptive filter.

13. The system according to claim 12 , wherein the filter is a Kalman Filter.

14. The system according to claim 11 , for operating continuously in real time.

Assignments (3)
MERGER Recorded Nov 11, 2015
From: NEMOR PROPERTIES LLC
To: MINERAL LASSEN LLC
Reel/Frame 037017/0001 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded May 9, 2008
From: SIONET LTD.
To: NEMOR PROPERTIES LLC
Reel/Frame 020931/0215 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Aug 15, 2005
From: MORIYA, NETZER
To: SIONET LTD.
Reel/Frame 016894/0190 →
Priority Claims (1)
IL 166837 · Feb 10, 2005 · national
Continuity (1)
Related Publication 20060179089A1 · Aug 10, 2006