Patent ID: 8554501

Claim:
A method for non-intrusive determination of electrical power consumed by a subscriber installation, in which periodic sampling of values of supply voltage and of electric current intensity delivered to said subscriber installation is performed in order to obtain a consumption signal, wherein the periodic sampling comprises the steps of, based on the sampled values of the consumption signal: grouping by a processor said values sampled in successive observation time windows; discriminating by a processor in each successive current observation time window a set of poles and of residues of the consumption signal by a Pencil method, said set of poles and residues being representative, in a current observation window, of a combination of signatures of separate electrical loads including at least defined poles and residues associated with at least one singular value obtained by the Pencil method in operation in said subscriber installation; calculating by a processor, for at least a duration of said current observation window, at least a subset active power consumed by at least one subset of the electrical loads in operation in the subscriber installation, expressed as a sum of separate active powers consumed by each separate electrical load of said subset complying with a consumed power function; wherein, for said subset active power, said consumed power function complies with a relation: P = ∑ i = 1 M ⁢ ∑ i ′ = 1 M ′ ⁢ R ci · ⁢ R vi ′ · ⅇ ( S ci + S vi ′ ) ⁢ ( t kd + Δ ⁢ ⁢ t 2 ) · sinh ⁢ ⁢ c ⁡ ( ( S ci + S vi ′ ) ⁢ Δ ⁢ ⁢ t 2 ) where P: denotes the subset active power consumed by the subset of separate electrical loads C i of signature (S i ,R i ,SG i ); S ci and R ci : denote the poles and respectively the residues of the sampled current delivered to the subscriber installation; S vi′ and R vi′ : denote the poles and respectively the residues of the sampled voltage delivered to the subscriber installation, SGi denoting the singular values; t kd : denotes the start time of the observation window of rank k; Δt: denotes the duration of the observation window; sinh c(X): a hyperbolic sine function of a cardinal of a value X, wherein: sinh ⁢ ⁢ c ⁡ ( ( S ci + S vi ′ ) ⁢ Δ ⁢ ⁢ t 2 ) = sinh ⁡ ( ( S ci + S vi ′ ) ⁢ Δ ⁢ ⁢ t 2 ) / ( ( S ci + S vi ′ ) ⁢ Δ ⁢ ⁢ t 2 ) ; M: denotes the number of poles and residues of the sampled current delivered to the subscriber installation; M′: denotes the number of poles and residues of the sampled voltage delivered to the subscriber installation.