US20160327418A1  Sensor system  Google Patents
Sensor system Download PDFInfo
 Publication number
 US20160327418A1 US20160327418A1 US15/097,356 US201615097356A US2016327418A1 US 20160327418 A1 US20160327418 A1 US 20160327418A1 US 201615097356 A US201615097356 A US 201615097356A US 2016327418 A1 US2016327418 A1 US 2016327418A1
 Authority
 US
 United States
 Prior art keywords
 sensor
 temperature
 signal
 compensation filter
 order
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Abandoned
Links
 230000001419 dependent Effects 0.000 claims abstract description 18
 238000009792 diffusion process Methods 0.000 claims description 22
 230000003595 spectral Effects 0.000 claims description 3
 230000001808 coupling Effects 0.000 claims description 2
 238000010168 coupling process Methods 0.000 claims description 2
 238000005859 coupling reaction Methods 0.000 claims description 2
 230000000875 corresponding Effects 0.000 abstract description 15
 238000010586 diagram Methods 0.000 description 11
 238000005070 sampling Methods 0.000 description 10
 238000004590 computer program Methods 0.000 description 4
 238000000034 method Methods 0.000 description 4
 230000003247 decreasing Effects 0.000 description 3
 239000011159 matrix material Substances 0.000 description 3
 238000005259 measurement Methods 0.000 description 3
 238000000354 decomposition reaction Methods 0.000 description 2
 238000009795 derivation Methods 0.000 description 2
 230000000694 effects Effects 0.000 description 2
 230000001747 exhibiting Effects 0.000 description 2
 239000012528 membrane Substances 0.000 description 2
 241001646071 Prioneris Species 0.000 description 1
 238000006243 chemical reaction Methods 0.000 description 1
 238000001514 detection method Methods 0.000 description 1
 230000002708 enhancing Effects 0.000 description 1
 230000003068 static Effects 0.000 description 1
 230000002195 synergetic Effects 0.000 description 1
 230000036962 time dependent Effects 0.000 description 1
Images
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
 G01D18/00—Testing or calibrating of apparatus or arrangements provided for in groups G01D1/00  G01D15/00
 G01D18/008—Testing or calibrating of apparatus or arrangements provided for in groups G01D1/00  G01D15/00 with calibration coefficients stored in memory

 G—PHYSICS
 G01—MEASURING; TESTING
 G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
 G01D3/00—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups
 G01D3/028—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure
 G01D3/036—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure on measuring arrangements themselves
 G01D3/0365—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure on measuring arrangements themselves the undesired influence being measured using a separate sensor, which produces an influence related signal
Abstract
Disclosed is a method for adjusting a sensor signal and a corresponding sensor system comprising a sensor for providing a sensor signal representative of a measure other than temperature, dynamic components of the sensor signal being dependent on temperature. In addition there is provided a temperature sensor for measuring the temperature. Dynamic components in the sensor signal are adjusted subject to the temperature sensed, and a compensated sensor signal is supplied. Such sensor system helps compensating for long response times of sensors.
Description
 This application claims the priority of European Patent Application 10005804.9, filed Jun. 4, 2010, the disclosure of which is incorporated herein by reference in its entirety.
 It is the nature of sensors that they react on real impacts with a certain time lag. Especially, this is true in fast varying environments in which the quantity to be measured may change in form of a step function, for example. However, the corresponding sensor signal may not step up to the new real measure value but rather gets there with a certain response time.
 There are many applications that only work properly with a sensor supplying a fast response to fast varying measures. However, for many sensors there are limitations in modifying the hardware in order to improve the response time.
 The problem to be solved by the present invention is therefore to provide a sensor and a method for providing a sensor signal with improved response time. It is also desired to provide a method for building such a sensor system or pieces of it respectively.
 The problem is solved by a sensor system according to the features of claim 1, by a method for adjusting a sensor signal according to the features of claim 7, and by a method for building a compensation filter for use in a sensor system according to the features of claim 12.
 According to a first aspect of the present invention, there is provided a sensor system comprising a sensor providing a sensor signal representative of a measure other than temperature, wherein dynamic components of the sensor signal are dependent on temperature. The system further includes a temperature sensor for providing a temperature signal. A compensation filter receives the sensor signal and the temperature signal. The compensation filter is designed for adjusting the dynamic components in the sensor signal subject to the temperature signal, and for providing a compensated sensor signal.
 According to another aspect of the present invention, there is provided a method for adjusting a sensor signal. According to this method a temperature is sensed and a sensor signal representative of a measure other than temperature is provided. Dynamic components of this sensor signal typically are dependent on temperature. The dynamic components in the sensor signal are adjusted subject to the temperature sensed. A compensated sensor signal is provided as a result of this adjustment.
 Sensors may not immediately react to changes in a measure but only react with a certain delay, also called response time. Such time lag may be owed to the appearance of diffusion processes, which may include a diffusion of the measure into the sensor—e.g. into a housing of the sensor—and, in addition, possibly a diffusion of the measure into a sensor element of the sensor. For some sensors and applications, chemical reactions of the measure with the sensor element may increase the response time, too. While a sensor may perfectly map a static measure into its sensor signal, dynamic components in the measure, such as swift changes, steps, or other highfrequent changes may be followed only with a delay in time.
 For such sensors/applications it may be beneficial if the sensor response is dynamically compensated. This means, that the sensor signal, and in particular the dynamic components of the sensor signal are adjusted such that the response time of the sensor is decreased. Such effort in decreasing the response time of a sensor is also understood as compensation of the sensor signal and especially its dynamic components.
 By using a dynamic compensation filter the response of a sensor, i.e. its output in response to a change in the measure, will be accelerated. If the dynamics of the sensor including its housing are known, an observer, i.e. a compensation filter, can be implemented to estimate the true physical value of the measure. This observer compensates for the sensor dynamics such that the response can be considerably accelerated in time, which means that the response time can be considerably decreased. This is why the proposed method and system enhance the system dynamics of sensors.
 It has been observed that the dynamic components of a sensor response—e.g. the gradient in the sensor signal—may be dependent on the ambient temperature, such that for example the response time of the sensor may be shorter at higher temperatures, and may take longer at lower temperatures. Consequently, in the present embodiments, it is envisaged to apply a compensating filter to the sensor output which compensating filter takes the measured temperature into account. As a result the compensated sensor signal supplied by the compensating filter is even more enhanced in that its deviation from the real measure is improved for the reason that temperature dependency of the sensor dynamics is taken into account of the compensation filter modelling.
 Another aspect of the invention provides a method for building a compensation filter for use in a sensor system. In this method, a sensor model of the sensor is built, the sensor model being characterized by a transfer function. A compensation filter is modelled based on an inverse of the transfer function of the sensor model. Temperature dependent terms are applied to the compensation filter.
 In another aspect of the present invention, there is disclosed a computer program element for adjusting a sensor signal, the computer program element comprising computer program instructions executable by a computer to receive a sensor signal representative of a measure other than temperature, dynamic components of the sensor signal being dependent on temperature, to adjust the dynamic components in the sensor signal subject to the temperature sensed, and to provide a compensated sensor signal.
 Other advantageous embodiments are listed in the dependent claims as well as in the description below. The described embodiments similarly pertain to the system, the methods, and the computer program element. Synergetic effects may arise from different combinations of the embodiments although they might not be described in detail.
 Further on it shall be noted that all embodiments of the present invention concerning a method might be carried out with the order of the steps as described, nevertheless this has not to be the only essential order of the steps of the method all different orders of orders and combinations of the method steps are herewith described.
 The aspects defined above and further aspects, features and advantages of the present invention can also be derived from the examples of embodiments to be described hereinafter and are explained with reference to examples of embodiments. The invention will be described in more detail hereinafter with reference of examples of embodiments but to which the invention is not limited.

FIG. 1 shows a diagram illustrating a change in a measure and associated sensor signals for different ambient temperatures, 
FIG. 2 shows a schematic illustration of a sensor system according to an embodiment of the present invention in a block diagram, 
FIG. 3A illustrates diagrams with sample measures, associated sensor signals, and associated compensated sensor signals, 
FIG. 3B illustrates diagrams with sample measures, associated sensor signals, and associated compensated sensor signals, 
FIG. 4A shows a graph of signals in the discrete time domain which signals appear in a sensor system of an embodiment according to the present invention, 
FIG. 4B shows graphs of signals in the discrete time domain which signals appear in a sensor system of an embodiment according to the present invention, 
FIG. 4C shows a graph of a signal in the discrete time domain which signals appear in a sensor system of an embodiment according to the present invention, 
FIG. 5A illustrates graphs for supporting the understanding of how to build a compensation filter according to an embodiment of the present invention, 
FIG. 5B illustrates graphs for supporting the understanding of how to build a compensation filter according to an embodiment of the present invention, 
FIG. 6A illustrates graphs for supporting the understanding of how to build a compensation filter according to another embodiment of the present invention, 
FIG. 6B illustrates graphs for supporting the understanding of how to build a compensation filter according to another embodiment of the present invention, 
FIG. 6C illustrates graphs for supporting the understanding of how to build a compensation filter according to another embodiment of the present invention, and 
FIG. 7 illustrates on a temperature scale with different temperature subranges a measure may be exposed to. 
FIG. 1 shows a diagram with a step function st representing an immediate change in the quantity of a measure to be measured. An xaxis of the diagram indicates time in s, the y axis indicates the quantity of the measure to be measured, in this particular case the relative humidity RH(t) in %. Accordingly, at t=0 the relative humidity RH steps up from 20% to 80%. The other graphs in the diagram indicate sensor signals over time provided by a relative humidity sensor in response to the step function st of the relative humidity. It can be derived from these graphs that the humidity sensor cannot immediately follow the change in the measure and approaches the new relative humidity of 80% only with a delay which delay is also called response time. The various graphs represent sensor signals subject to different ambient temperatures in Celsius degrees. Basically, the higher the ambient temperature is the faster the response time of the sensor signal is. The lower the ambient temperature is the longer the response time of the sensor signal is. 
FIG. 2 shows a schematic illustration of a sensor system according to an embodiment of the present invention. A relative humidity sensor 1 senses the ambient relative humidity RH. In the following, the term humidity is used as a synonym for relative humidity. The physics of the humidity sensor 1 is such that dynamics of the sensor signal RH_{sensor }at its output is temperature dependent. Insofar, the humidity sensor 1 ofFIG. 2 may behave according to a humidity sensor providing a temperature dependent output according to the graphs inFIG. 1 . This means that the response time of the humidity sensor is dependent from the ambient temperature T. Accordingly, inFIG. 2 the ambient temperature T is illustrated as input to the humidity sensor 1 although the humidity sensor 1 is not meant to measure the ambient temperature 1 but rather the dynamics of the sensor signal RH_{sensor }at its output is varying subject to the ambient temperature T.  A second sensor is provided, which is a temperature sensor 2 for measuring the ambient temperature T. The output of the temperature sensor 2 provides a temperature signal T representing the ambient temperature T.
 A compensator 3 receives the sensor signal RH_{sensor}(t) over time—which sensor signal is also denoted as u(t)—and the temperature signal T(t) over time, and adjusts the dynamics in the sensor signal RH_{sensor}(t) subject to the temperature signal T(t). As a result of this adjustment, the compensator 3 provides at its output a compensated sensor signal RH_{compensated}(t) over time. In qualitative terms, the compensated sensor signal RH_{compensated}(t) shall compensate for the response time of the humidity sensor 1 at its best and take a gradient that is more close to the gradient of the measure to be measured. As such, the compensator 3 adjusts for the dynamics of the sensor signal RH_{sensor}(t) and consequently for the physics of the sensor 1 not allowing for better response times.
 In order to build a compensator 3 that actually compensates for the dynamics of the sensor 1 in the desired way, the behaviour of the sensor 1 needs to be understood. While the following sections are described in connection with humidity sensing, it is understood that the principles can be generalized to any other sensor which shows a response time not satisfying for an application the sensor is used in and in particular a temperature dependent response time.
 Modelling the Sensor:
 The following desrcibes the modelling of a humidity sensor out of which model a suitable compensator may be derived. Such compensator may be applied downstream to the real sensor and compensate for dynamics in the sensor signal output by such sensor.
 For a humidity sensor, the humidity as the relevant measure needs to reach a sensing element of the humidity sensor. Such sensing element of a humidity sensor preferably is a membrane. For further information on humidity sensing background it is referred to EP 1 700 724. Consequently, the humidity to be measured needs to diffuse into the membrane of the humidity sensor. Prior to this, the humidity needs to diffuse into a housing of the humidity sensor provided the humidity sensor has such housing.
 For humidity sensors for which both of these diffusion processes—i.e. from the outside into the housing and from the housing into the sensor element—are relevant both processes are preferred to be respected in a corresponding sensor model.
 Both of the processes may be described with a differential diffusion equation with two independent diffusion time constants. In a first approximation, the sensor can be described by a transfer function in the frequency domain with two poles and one zero such that the generic design of a corresponding transfer function G_{2}(s) may look like, specifically in a second order model:

$\begin{array}{cc}{G}_{2}\ue8a0\left(s\right)=\frac{\mathrm{Ks}+1}{\left({T}_{1}\ue89es+1\right)\ue89e\left({T}_{2}\ue89es+1\right)},& \left(\mathrm{Eq}.\phantom{\rule{0.8em}{0.8ex}}\ue89e1\right)\end{array}$  where s denotes the complex Laplace variable and K, T_{1}, T_{2 }are constants to be identified. T_{1 }and T_{2 }are time constants of the respective diffusion processes and K defines a coupling between the two processes. The transfer function G(s) generally describes the characteristics of the sensor in the frequency domain by

RH _{sensor}(s)=G(s)*RH(s)  If the housing of the humidity sensor is very complex, an additional pole and zero may be added, and the transfer function of the humidity sensor may be amended accordingly.
 Some other applications may require a less complex sensor representation. In this case, a first order model may approximate the sensor, which first order model is characterized by a transfer function in the frequency domain with one pole such that the generic design of such transfer function G_{1}(s) may look like, specifically in a first order model:

$\begin{array}{cc}{G}_{1}\ue8a0\left(s\right)=\frac{1}{\left({T}_{1}\ue89es+1\right)}.& \left(\mathrm{Eq}.\phantom{\rule{0.8em}{0.8ex}}\ue89e2\right)\end{array}$  Such model may be sufficient, for example, if the humidity sensor does not include a housing such that the diffusion process into the housing can be neglected, or, if one of the two diffusion processes—either from the outside into the housing or from the housing into the sensor element—is dominant over the other, such that the transfer function may be approximated by the dominant diffusion process only.
 There are different ways for identifying the parameters T_{1}, T_{2 }and K. One approach is to find the parameters by trial and error. In a first trial, the sensor output is simulated by the sensor model wherein the sensor model makes use of a first estimation of the parameters. The output of the sensor model then is compared with the sensor signal supplied by the real sensor. Afterwards, the parameters are adjusted until a deviation between the simulated sensor output and the real sensor signal is acceptably small.
 For implementing the sensor model in a digital system such as a microcontroller the sensor model preferably is implemented in the discrete time domain rather than in the continuous frequency domain as described by equations 1 or 2. Therefore, the sensor model needs to be digitised, i.e. transformed into a set of difference equations.
 First, the differential equations in the frequency domain, e.g. the equations 1 and 2, can be reverse transformed into the time continuous domain by the Laplace backtransform.
 For the second order model of the sensor according to equation 1 the equivalent time continuous state space description in the control canonical form is:

$\left[\begin{array}{c}\frac{\uf74c{x}_{1}\ue8a0\left(t\right)}{\uf74ct}\\ \frac{\uf74c{x}_{2}\ue8a0\left(t\right)}{\uf74ct}\end{array}\right]={A}_{c}^{M}\xb7\left[\begin{array}{c}{x}_{1}\ue8a0\left(t\right)\\ {x}_{2}\ue8a0\left(t\right)\end{array}\right]+{B}_{c}^{M}\ue89ew\ue8a0\left(t\right)$ $v\ue8a0\left(t\right)={C}_{c}^{M}\ue89ex\ue8a0\left(t\right)+{D}_{c}^{M}\ue89ew\ue8a0\left(t\right)$  with
 w(t) denoting the real relative humidity RH at time t in the time domain,
 x1(t) and x2(t) denoting internal states of the second order sensor model, and
 v(t) denoting the sensor model output over time t.
 In this time continuous state space representation, the coefficients A, B, C and D are determined by:

${A}_{c}^{M}=\left[\begin{array}{cc}\frac{{T}_{1}+{T}_{2}}{{T}_{1}\ue89e{T}_{2}}& \frac{1}{{T}_{1}\ue89e{T}_{2}}\\ 1& 0\end{array}\right],{B}_{c}^{M}\ue8a0\left[\begin{array}{c}1\\ 0\end{array}\right]$ ${C}_{c}^{M}=\left[\begin{array}{cc}\frac{K}{{T}_{1}\ue89e{T}_{2}}& \frac{1}{{T}_{1}\ue89e{T}_{2}}\end{array}\right],{D}_{c}^{M}=0$  For a first order sensor model the time continuous state space description in the control canonical form is:

$\frac{\uf74cx\ue8a0\left(t\right)}{\uf74ct}={A}_{c}^{M}\xb7x\ue8a0\left(t\right)+{B}_{c}^{M}\ue89ew\ue8a0\left(t\right)$ $v\ue8a0\left(t\right)={C}_{c}^{M}\ue89ex\ue8a0\left(t\right)+{D}_{c}^{M}\ue89ew\ue8a0\left(t\right)$  again, with
 w(t) denoting the real relative humidity RH at time t in the time domain,
 x(t) denoting an internal state of the first order sensor model, and
 v(t) denoting the sensor model output over time t.
 In this time continuous state space representation, the coefficients A, B, C and D are determined by:

${A}_{c}^{M}=\frac{1}{{T}_{1}},{B}_{c}^{M}=1,{C}_{c}^{M}=\frac{1}{{T}_{1}},{D}_{c}^{M}=0$  Second, the representation in the time continuous domain may be transformed into the time discrete domain.
 A time discrete state space representation equivalent to the time continuous state space representation for the second order sensor model may be:

$\left[\begin{array}{c}{x}_{1}\ue8a0\left(k+1\right)\\ {x}_{2}\ue8a0\left(k+1\right)\end{array}\right]={A}_{d}^{M}\xb7\left[\begin{array}{c}{x}_{1}\ue8a0\left(k\right)\\ {x}_{2}\ue8a0\left(k\right)\end{array}\right]+{B}_{d}^{M}\ue89ew\ue8a0\left(k\right)$ $v\ue8a0\left(t\right)={C}_{d}^{M}\ue8a0\left[\begin{array}{c}{x}_{1}\ue8a0\left(k\right)\\ {x}_{2}\ue8a0\left(k\right)\end{array}\right]+{D}_{d}^{M}\ue89ew\ue8a0\left(k\right)\ue89e\text{}\left[\begin{array}{c}{x}_{1}\ue8a0\left(0\right)\\ {x}_{2}\ue8a0\left(0\right)\end{array}\right]={\left(I{A}_{d}^{M}\right)}^{1}\ue89e{B}_{d}^{M}\xb7w\ue8a0\left(0\right)$  where w denotes the measured relative humidity at time step k with the sampling time t(k+1)−t(k)=T_{s}. There are two internal states x_{1}(k), x_{2}(k) of the second order sensor model. This means x(k)=(x_{1}(k), x_{2}(k))^{T}. v(k) denotes the sensor model output in the discrete time domain.
 In this time discrete state space representation, the coefficients A, B, C and D are determined by:

A _{d} ^{M} =e ^{A} ^{ c } ^{ M } ^{T} ^{ s } , B _{d} ^{M}=(A _{c} ^{M})^{−1}(e ^{A} ^{ c } ^{ M } ^{T} ^{ s } −I)·B _{c} ^{M } 
C _{d} ^{M} =C _{c} ^{M} , D _{d} ^{M} =D _{c} ^{M }  In addition,

$I=\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]$  and T_{s }is the sampling time.
 A time discrete state space representation of the time continuous state space representation for the first order sensor model may be:

$x\ue8a0\left(k+1\right)={A}_{d}^{M}\xb7x\ue8a0\left(k\right)+{B}_{d}^{M}\xb7w\ue8a0\left(k\right)$ $v\ue8a0\left(k\right)={C}_{d}^{M}\xb7x\ue8a0\left(k\right)+{D}_{d}^{M}\xb7w\ue8a0\left(k\right)$ $x\ue8a0\left(0\right)=\frac{{B}_{d}^{M}}{1{A}_{d}^{M}}\ue89ew\ue8a0\left(0\right)$  where w(k) denotes the measured relative humidity RH at time step k with the sampling time t(k+1)−t(k)=T_{s}. There is an internal state x(k) of the first order sensor model. v(k) denotes the modelled sensor output in the discrete time domain.
 In this time discrete state space representation, the coefficients A, B, C and D are determined by:

${A}_{d}^{M}={\uf74d}^{{T}_{1}\ue89e{T}_{s}},{B}_{d}^{M}=\frac{1}{{T}_{1}}\ue89e\left({\uf74d}^{{T}_{1}\ue89e{T}_{s}}1\right)$ ${C}_{d}^{M}=1,{D}_{d}^{M}=0$  where T_{s }is the sampling time.
 The derivation of the matrices for the first and the second order sensor model requires Schur decomposition or series expansion and matrix inversion. Computer software may help to calculate the coefficients.
 The time discrete state space representations of the first or the second order sensor model may be run on a microprocessor, and the parameters T_{1 }or T_{1}, T_{2 }and K respectively may be varied until the sensor model output v(k) is close enough to the real sensor signal u(t) which may be present in digitized form u(k), too, or may be digitized for comparing purposes.
 More sophisticated methodologies for determining the parameters of the sensor model may use system identification tools that automatically build dynamical models from measured data.
 Modelling the Compensator:
 In a next step, the compensator 3 out of
FIG. 2 is determined and implemented. When the humidity sensor 1 is modelled by a second order model then, advantageously, the compensator is modelled by a second order model, too. The transfer function C_{2}(s) of such a second order compensation filter 3—also denoted as compensator—in the frequency domain may be described by: 
$\begin{array}{cc}{C}_{2}\ue8a0\left(s\right)=\frac{\left({T}_{1}\ue89es+1\right)\ue89e\left({T}_{2}\ue89es+1\right)}{\left(\mathrm{Ks}+1\right)\ue89e\left(\mathrm{Ps}+1\right)}& \left(\mathrm{Eq}.\phantom{\rule{0.8em}{0.8ex}}\ue89e3\right)\end{array}$  If a first order sensor model is applied, the following first order compensation filter is proposed, which may be described by a transfer function C_{1}(s) in the frequency domain by:

$\begin{array}{cc}{C}_{1}\ue8a0\left(s\right)=\frac{\left({T}_{1}\ue89es+1\right)}{\left(\mathrm{Ps}+1\right)}.& \left(\mathrm{Eq}.\phantom{\rule{0.8em}{0.8ex}}\ue89e4\right)\end{array}$  For both models, s denotes the complex Laplace variable and K, T_{1}, T_{2 }are the constants that were identified when determining the sensor model.
 C(s) generally denotes the transfer function of the compensation filter in the frequency domain, wherein

RH _{compensated}(s)=C(s)*RH _{sensor}(s)  Preferably, the compensation filter transfer function C(s) is the inverse to the sensor model transfer function G(s), i.e. the diffusion function, such that

C(s)=1/G(s)  The term (Ps+1) is introduced in the compensator transfer function to make the function physically applicable. Parameter P is kept small in order to keep impact on filter function low, but can be used to filter measurement noise.
 An important feature of the compensation filter transfer function C(s) is that the final value of C(s) converges to 1, i.e.

$\underset{s\to 0}{\mathrm{lim}}\ue89eC\ue8a0\left(s\right)=1$  This means that the compensation filter 3 only changes the sensor output characteristic during transition. When the compensation filter 3 is in steady state, it does not affect the sensor output, even if the modelling of the sensor and its housing is inaccurate. Please note that overshoots may occur if the real system response RH(t) is faster than modelled.
 Typically, the compensation filter 3 is implemented in a digital system such as a microcontroller which operates on samples of the measured humidity rather than on the continuous signal. As a consequence, microcontrollers cannot integrate and not implement differential equations like those in equation 3 or 4. Therefore, the compensator needs to be digitised, i.e. transformed into a set of difference equations.
 First, the differential equations in the spectral domain, i.e. the equations 3 and 4 in the present example, can be reverse transformed into the time continuous domain by the Laplace backtransform.
 For the second order compensation filter the time continuous state space description in the control canonical form is:

$\left[\begin{array}{c}\frac{\uf74c{x}_{1}\ue8a0\left(t\right)}{\uf74ct}\\ \frac{\uf74c{x}_{2}\ue8a0\left(t\right)}{\uf74ct}\end{array}\right]={A}_{c}^{C}\xb7\left[\begin{array}{c}{x}_{1}\ue8a0\left(t\right)\\ {x}_{2}\ue8a0\left(t\right)\end{array}\right]+{B}_{c}^{C}\ue89eu\ue8a0\left(t\right)$ $y\ue8a0\left(t\right)={C}_{c}^{C}\ue89ex\ue8a0\left(t\right)+{D}_{c}^{C}\ue89eu\ue8a0\left(t\right)$  with u(t) denoting the measured humidity at time t in the time domain, i.e. the sensor signal RH_{sensor}(t),
 x1(t) and x2(t) denoting internal states of the second order compensation filter, and
 y(t) denoting the compensated sensor signal, i.e. RH_{compensated}(t) according to
FIG. 2 .  In this time continuous state space representation, the coefficients A, B, C and D are determined by:

${A}_{c}^{C}=\left[\begin{array}{cc}\frac{K+P}{\mathrm{KP}}& \frac{1}{\mathrm{KP}}\\ 1& 0\end{array}\right],{B}_{c}^{C}=\left[\begin{array}{c}1\\ 0\end{array}\right],{D}_{c}^{C}=\frac{{T}_{1}\ue89e{T}_{2}}{\mathrm{KP}}$ ${C}_{c}^{C}=\left[\begin{array}{cc}\frac{{T}_{1}+{T}_{2}}{\mathrm{KP}}\frac{{T}_{1}\ue89e{T}_{2}\ue8a0\left(K+P\right)}{{\left(\mathrm{KP}\right)}^{2}}& \frac{1}{\mathrm{KP}}\frac{{T}_{1}\ue89e{T}_{2}}{{\left(\mathrm{KP}\right)}^{2}}\end{array}\right]$  For the first order compensation filter the time continuous state space description in the control canonical form is:

$\frac{\uf74cx\ue8a0\left(t\right)}{\uf74ct}={A}_{c}^{C}\xb7x\ue8a0\left(t\right)+{B}_{c}^{C}\ue89eu\ue8a0\left(t\right)$ $y\ue8a0\left(t\right)={C}_{c}^{C}\ue89ex\ue8a0\left(t\right)+{D}_{c}^{C}\ue89eu\ue8a0\left(t\right)$  with u(t) denoting the measured humidity at time t in the time domain, i.e. the sensor signal RH_{sensor}(t) in
FIG. 2 ,  x(t) denoting an internal state of the second order compensation filter, and
 y(t) denoting the compensated sensor signal, i.e. RH_{compensated}(t) in
FIG. 2 .  In this time continuous state space representation, the coefficients A, B, C and D are determined by:

${A}_{c}^{C}=\frac{1}{P},{B}_{c}^{C}=1,{C}_{c}^{C}=\frac{P{T}_{1}}{{P}^{2}},{D}_{c}^{C}=\frac{{T}_{1}}{P}$  Second, the representation in the time continuous domain may be transformed into the time discrete domain.
 A time discrete state space representation of the time continuous state space representation for the second order compensation filter may be:

$\begin{array}{cc}\left[\begin{array}{c}{x}_{1}\ue8a0\left(k+1\right)\\ {x}_{2}\ue8a0\left(k+1\right)\end{array}\right]={A}_{d}^{C}\xb7\left[\begin{array}{c}{x}_{1}\ue8a0\left(k\right)\\ {x}_{2}\ue8a0\left(k\right)\end{array}\right]+{B}_{d}^{C}\ue89eu\ue8a0\left(k\right)\ue89e\text{}\ue89ey\ue8a0\left(t\right)={C}_{d}^{C}\ue8a0\left[\begin{array}{c}{x}_{1}\ue8a0\left(k\right)\\ {x}_{2}\ue8a0\left(k\right)\end{array}\right]+{D}_{d}^{C}\ue89eu\ue8a0\left(k\right)\ue89e\text{}\left[\begin{array}{c}{x}_{1}\ue8a0\left(0\right)\\ {x}_{2}\ue8a0\left(0\right)\end{array}\right]={\left(I{A}_{d}^{C}\right)}^{1}\ue89e{B}_{d}^{C}\xb7u\ue8a0\left(0\right)& \mathrm{Equ}.\phantom{\rule{0.8em}{0.8ex}}\ue89e5\end{array}$  where u(k) denotes the measured humidity at time step k with the sampling time t(k+1)−t(k)=T_{s}. There are two internal states of the second order compensator, i.e. x(k)=(x_{1}(k), x_{2}(k))^{T}. y(k) denotes the compensated sensor signal in the discrete time domain.
 In this time discrete state space representation, the coefficients A, B, C and D are determined by:

A _{d} ^{C} =e ^{A} ^{ c } ^{ C } ^{T} ^{ s } 
B _{d} ^{C}=(A _{c} ^{C})^{−1}(e ^{A} ^{ c } ^{ C } ^{T} ^{ s } −I)·B _{c} ^{C } 
C _{d} ^{C} =C _{c} ^{C } 
D _{d} ^{C} =D _{c} ^{C }  A time discrete state space representation of the time continuous state space representation for the first order compensation filter may be:

$x\ue8a0\left(k+1\right)={A}_{d}^{C}\xb7x\ue8a0\left(k\right)+{B}_{d}^{C}\xb7u\ue8a0\left(k\right)$ $y\ue8a0\left(k\right)={C}_{d}^{C}\xb7x\ue8a0\left(k\right)+{D}_{d}^{C}\xb7u\ue8a0\left(k\right)$ $x\ue8a0\left(0\right)=\frac{{B}_{d}^{C}}{1{A}_{d}^{C}}\ue89eu\ue8a0\left(0\right)$  where u(k) denotes the measured humidity at time step k with the sampling time t(k+1)−t(k)=T_{s}. There is an internal state x(k) of the first order compensator. y(k) denotes the compensated sensor signal in the discrete time domain.
 In this time discrete state space representation, the coefficients A, B, C and D are determined by:

${A}_{d}^{C}={\uf74d}^{{\mathrm{PT}}_{s}},{B}_{d}^{C}=\frac{1}{P}\ue89e\left({\uf74d}^{{\mathrm{PT}}_{s}}1\right)$ ${C}_{d}^{C}=\frac{P{T}_{1}}{{P}^{2}},{D}_{d}^{C}=\frac{{T}_{1}}{P}$  The derivation of the matrices for the first and the second order compensation filter requires Schur decomposition or series expansion and matrix inversion. Computer software may help to calculate the coefficients.

FIGS. 4A, 4B, and 4C illustrate sample signals, primarily in the discrete time domain.FIG. 4A illustrates a sample ambient humidity RH over time t. The measured humidity as continuous signal over time is denoted by u(t), i.e. the sensor signal u(t) supplied by the humidity sensor. The sensor signal u(t) is sampled at points k=0, k=1, . . . with the sampling time T_{s}=t(k+1)−t(k). u(k) denotes the sampled sensor signal which is a time discrete signal.  For a second order compensation filter the internal states x1 and x2 over sampling points k are illustrated in the graph in
FIG. 4B . In the last graph inFIG. 4C , the compensated sensor signal is illustrated as step function in the discrete time domain. For illustration purposes, an equivalent Y in the continuous time domain is depicted as dot and dash line.  The compensation effect with respect to a sample humidity characteristic over time is also illustrated in the diagrams of
FIGS. 3A, 3B and 3C . InFIG. 3A , the relative humidity in the environment, denoted as “Ambient condition” is shown as a function of time and includes three step wise changes and a gradient in form of a ramp. The humidity sensor supplies a sensor signal—denoted as “Sensor output”—which exhibits significant response times to the fast changes in the Ambient condition. As mentioned before, the Ambient condition is filtered within the sensor because of the internal diffusion dynamics and the dynamics of the housing. A compensation filter following the sensor supplies a compensated sensor signal—denoted as “Compensated output”. In the present example, the compensation filter is implemented as a second order model. As can be derived from the diagram, the Compensated output much more resembles the Ambient condition than the sensor output does without a compensator applied to the sensor output. Consequently, the compensator helps in improving the quality of the measuring system such that the overall output of the sensor system including the compensator provides for a better quality output especially for fast changes in the Ambient condition.FIG. 3B illustrates the Sensor output and the Compensated Output for another Ambient condition characteristic in form of a saw tooth. In this example, the compensation filter is implemented as a first order model.  The diagrams in
FIGS. 5A and 5B support the illustration of a method for implementing a sensor system according to an embodiment of the present invention. In a first step, the sensor itself is modelled for the purpose of deriving a compensator model from the data that describes the sensor model. In this first step, the real sensor again is exposed to a step wise increase of the measure, i.e. the relative humidity denoted as “Ambient condition” again. The relative humidity RH is increased in a single step from about 15% to about 65%. The corresponding Sensor output is measured, and preferably stored. The Sensor output exhibiting a response time is plotted in bothFIG. 5A andFIG. 5B . Upfront, it is determined, that the humidity sensor used may best be represented by a first order model. Consequently, in a second step, the sensor model is built by means of the transfer function: 
${G}_{1}\ue8a0\left(s\right)=\frac{1}{\left({T}_{1}\ue89es+1\right)}$  and the parameter T_{1 }needs to be determined. In order to determine the parameter T_{1 }such that the output of the sensor model fits the real Sensor output as depicted in
FIGS. 5A and 5B best, the first order type sensor model is implemented on a microcontroller in the time discrete state space according to the equations that were explained previously: 
$x\ue8a0\left(k+1\right)={A}_{d}^{M}\xb7x\ue8a0\left(k\right)+{B}_{d}^{M}\xb7w\ue8a0\left(k\right)$ $v\ue8a0\left(k\right)={C}_{d}^{M}\xb7x\ue8a0\left(k\right)+{D}_{d}^{M}\xb7w\ue8a0\left(k\right)$ $x\ue8a0\left(0\right)=\frac{{B}_{d}^{M}}{1{A}_{d}^{M}}\ue89ew\ue8a0\left(0\right)$  With coefficients

${A}_{d}^{M}={\uf74d}^{{T}_{1}\ue89e{T}_{s}},{B}_{d}^{M}=\frac{1}{{T}_{1}}\ue89e\left({\uf74d}^{{T}_{1}\ue89e{T}_{s}}1\right)$ ${C}_{d}^{M}=1,{D}_{d}^{M}=0$  In this time discrete state space representation of the sensor model, an initial value of the parameter T_{1 }is chosen, e.g. T_{1}=6s. The sensor model output with T_{1}=6s is illustrated in
FIG. 5A . The sensor model is also run with parameters T_{1}=8s and T_{1}=12s. The corresponding sensor model output is illustrated inFIG. 5A , too.  In a next step, the value of the parameter is chosen that makes the sensor model output come closest to the real sensor signal. From
FIG. 5A it can be derived, that this is the case for T_{1}=8s. Accordingly, the complete transfer function of the sensor model in the frequency domain is 
${G}_{1}\ue8a0\left(s\right)=\frac{1}{\left(8\ue89es+1\right)}$  While its representation in the time discrete domain is

$x\ue8a0\left(k+1\right)=0.8465\xb7x\ue8a0\left(k\right)+0.9211\xb7w\ue8a0\left(k\right)$ $v\ue8a0\left(k\right)=0.1667\ue89e\phantom{\rule{0.3em}{0.3ex}}\ue89ex\ue8a0\left(k\right)+0\xb7w\ue8a0\left(k\right)$ $x\ue8a0\left(0\right)=\frac{0.9211}{10.8465}\ue89ew\ue8a0\left(0\right)$  The method can be modified in that the sensor is modelled first with a first value of the parameter. If the deviation of the sensor output based on the model with the first parameter value is too large, another value for the parameter is chosen. Iteratively, as many parameter values are chosen as long as there is considered to be sufficient similarity between the sensor model output and the real sensor signal, i.e. the deviation between those two is below a threshold.
 In a next step, the compensator model is determined to be a first order model with a general representation in the frequency domain of

${C}_{2}\ue8a0\left(s\right)=\frac{{T}_{1}\ue89es+1}{\mathrm{Ps}+1}$  The representation in the time discrete state space according to the above is:

$x\ue8a0\left(k+1\right)={A}_{d}^{C}\xb7x\ue8a0\left(k\right)+{B}_{d}^{C}\xb7u\ue8a0\left(k\right)$ $y\ue8a0\left(k\right)={C}_{d}^{C}\xb7x\ue8a0\left(k\right)+{D}_{d}^{C}\xb7u\ue8a0\left(k\right)$ $x\ue8a0\left(0\right)=\frac{{B}_{d}^{C}}{1{A}_{d}^{C}}\ue89eu\ue8a0\left(0\right)$  with u(k) denoting the discrete input to the compensator, i.e. the discrete representation of output of the real sensor, and y(k) denoting the discrete output of the compensator, i.e. the compensated sensor signal in the time discrete domain. The corresponding coefficients are:

${A}_{d}^{C}={e}^{{\mathrm{PT}}_{s}},\text{}\ue89e{B}_{d}^{C}=\frac{1}{P}\ue89e\left({e}^{{\mathrm{PT}}_{s}}1\right)$ ${C}_{d}^{C}=\frac{P{T}_{1}}{{P}^{2}},\text{}\ue89e{D}_{d}^{C}=\frac{{T}_{1}}{P}$  With the determination of the parameter T_{1 }upon implementation of the sensor model, the coefficients AD of the compensator can now be determined. In this step, parameter P is chosen, too. Parameter P effects reducing the signal noise.
 The compensator can now be implemented. In its time discrete state space, the compensator is described by:

$x\ue8a0\left(k+1\right)=0.3679\xb7x\ue8a0\left(k\right)+0.6321\xb7u\ue8a0\left(k\right)$ $y\ue8a0\left(k\right)=8.0\xb7x\ue8a0\left(k\right)+7.0\xb7u\ue8a0\left(k\right)$ $x\ue8a0\left(0\right)=\frac{0.6321}{10.3679}\ue89eu\ue8a0\left(0\right)$  The compensator now can be validated on the data measured. P can be adjusted to best performance in the signal to noise ratio SNR.
 In the present embodiment, parameter P is chosen as P=1.
 The Compensated output is shown in
FIG. 5B .  The response time is limited compared to the uncompensated Sensor output. The compensator, of course works for all kinds of signals
 The diagrams in
FIGS. 6A, 6B and 6C support the illustration of a method for implementing a sensor system according to another embodiment of the present invention. In contrast to the embodiment ofFIGS. 5A and 5B the sensor model and the compensator now are described as a second order model. Besides, the steps of implementing the sensor system including the compensator are the same as described in connection withFIGS. 5A and 5B . Again, in a first step, the sensor itself is modelled for the purpose of deriving a compensator model from the data that describes the sensor model. In this first step, the real sensor again is exposed to a step wise increase of the measure, i.e. the relative humidity denoted as “Ambient condition” again. The relative humidity RH is increased in a single step from about 15% to about 65%. The corresponding Sensor output is measured, and preferably stored. The Sensor output exhibiting a response time is plotted in allFIGS. 5A, 5B and 5C . Consequently, in a second step, the sensor model is built by means of the transfer function: 
${G}_{2}\ue8a0\left(s\right)=\frac{\mathrm{Ks}+1}{\left({T}_{1}\ue89es+1\right)\ue89e\left({T}_{2}\ue89es+1\right)}$  and the parameter T_{1}, T2 and K need to be determined. In order to determine the parameter such that the output of the sensor model fits the real Sensor output as depicted in
FIG. 6A best, the sensor model is implemented on a microcontroller in the time discrete state space according to the equations as given above. Initial values of the parameters are chosen and in an iterative way the parameters are set such that the output of the sensor model fits the sensor signal of the real sensor best. InFIG. 6A , there are illustrated multiple outputs of sensor models implemented with different parameter settings. It can be derived, that the parameter setting of T_{1}=16, T_{2}=155, and K=90 fits best. With such a parameter setting, the description of the sensor model in the time discrete state space is finalized.  In a next step, the compensator model is built as a second order model with a representation in the frequency domain—including the parameters as determined while building the sensor model, and assuming parameter P to be set to P=20:

${C}_{2}\ue8a0\left(s\right)=\frac{\left(16\ue89es+1\right)\ue89e\left(155\ue89es+1\right)}{\left(90\ue89es+1\right)\ue89e\left(20\ue89es+1\right)}$  This compensator model can be described in the time discrete state space by:

$\left[\begin{array}{c}{x}_{1}\ue8a0\left(k+1\right)\\ {x}_{2}\ue8a0\left(k+1\right)\end{array}\right]=\left[\begin{array}{cc}9.78\xb7{10}^{1}& 2.20\xb7{10}^{4}\\ 1.98\xb7{10}^{1}& 1.00\xb7{10}^{1}\end{array}\right]\xb7\left[\begin{array}{c}{x}_{1}\ue8a0\left(k\right)\\ {x}_{2}\ue8a0\left(k\right)\end{array}\right]+\left[\begin{array}{c}0.198\\ 0.020\end{array}\right]\ue89eu\ue8a0\left(k\right)$ $y\ue8a0\left(k\right)=\left[\begin{array}{cc}0.116& 0.00195\end{array}\right]\ue8a0\left[\begin{array}{c}{x}_{1}\ue8a0\left(k\right)\\ {x}_{2}\ue8a0\left(k\right)\end{array}\right]+2.756\ue89eu\ue8a0\left(k\right)\ue89e\text{}\left[\begin{array}{c}{x}_{1}\ue8a0\left(0\right)\\ {x}_{2}\ue8a0\left(0\right)\end{array}\right]=\left[\begin{array}{c}0\\ 1800\end{array}\right]\xb7u\ue8a0\left(0\right)$  Again, the compensator now can be validated on the data measured. The Compensated output is shown in
FIG. 6B andFIG. 6C . Both Compensated outputs differ in that different parameters P are applied. While inFIG. 6B , P is set to P=10, the Compensated output is exposed to significant noise. InFIG. 6C P is set to P=20 which provides a much better signal to noise ratio SNR.  However, diffusion processes in general depend very much on temperature. Therefore, the model for the compensating filter preferably takes temperature dependency into account—especially if temperature varies by more than 1020° C. This can be easily achieved by making the constants T1, T2, . . . , K, . . . time dependent or by using a set of sensor models G_{i}(s) and corresponding compensator models C_{i}(s) for different temperature ranges and switching amongst them.
 In a preferred embodiment, a second order temperature dependent sensor model may be described in the frequency domain by:

${G}_{2}\ue8a0\left(s,T\right)=\frac{\mathrm{Ks}+1}{\left({T}_{1}\ue8a0\left(T\right)\ue89es+1\right)\ue89e\left({T}_{2}\ue8a0\left(T\right)\ue89es+1\right)}$  with T being the temperature.
 The corresponding compensation filter may be described by

${C}_{2}\ue8a0\left(s,T\right)=\frac{\left({T}_{1}\ue8a0\left(T\right)\ue89es+1\right)\ue89e\left({T}_{2}\ue8a0\left(T\right)\ue89es+1\right)}{\left(\mathrm{Ks}+1\right)\ue89e\left(\mathrm{Ps}+1\right)}$  In the time continuous state space, the following equations represent the temperature dependent compensator:

$\frac{{\mathrm{dx}}_{1}\ue8a0\left(t\right)}{\mathrm{dt}}={a}_{11}\ue8a0\left(T\right)\ue89e{x}_{1}+{a}_{12}\ue8a0\left(T\right)\ue89e{x}_{2}+{b}_{1}\ue8a0\left(T\right)\ue89eu\ue8a0\left(t\right)$ $\frac{{\mathrm{dx}}_{2}\ue8a0\left(t\right)}{\mathrm{dt}}={a}_{21}\ue8a0\left(T\right)\ue89e{x}_{1}+{a}_{22}\ue8a0\left(T\right)\ue89e{x}_{2}+{b}_{2}\ue8a0\left(T\right)\ue89eu\ue8a0\left(t\right)$ $y\ue8a0\left(t\right)={c}_{1}\ue89e{x}_{1}+{c}_{2}\ue89e{x}_{2}$  In the time discrete state space, the following equations represent the temperature dependent compensator:

x(k+1)=A _{d} ^{C}(T)·x(k)+B _{d} ^{C}(T)·u(k) 
y(k)=C _{c} ^{D}(T)·x(k)+D _{d} ^{C}(T)·u(k)  with coefficients

${A}_{d}^{C}\ue8a0\left(T\right)=\left[\begin{array}{cc}{a}_{11}^{1}+{a}_{12}^{2}\ue89eT& {a}_{12}^{1}+{a}_{12}^{2}\ue89eT\\ {a}_{21}^{1}+{a}_{21}^{2}\ue89eT& {a}_{22}^{1}+{a}_{22}^{2}\ue89eT\end{array}\right]$ ${B}_{d}^{C}\ue8a0\left(T\right)=\left[\begin{array}{c}{b}_{1}^{1}+{b}_{1}^{2}\ue89eT\\ {b}_{2}^{1}+{b}_{2}^{2}\ue89eT\end{array}\right]$ ${C}_{d}^{C}\ue8a0\left(T\right)=\left[\begin{array}{cc}1& 1\end{array}\right]$ ${D}_{d}^{C}\ue8a0\left(T\right)=\left[{d}^{1}+{d}^{2}\ue89eT\right]$  The determination of A_{d} ^{C}(T) may be complex. A linear interpolation of matrix elements may be an appropriate way to determine the coefficients.
 Alternatively, the temperature range is divided into n subranges i, with temperature subrange i=0 being the first one, and temperature subrange i=n being the last one spanning the temperature range covered. For each subrange i a different sensor model and consequently a different compensator model is determined. For example, the following compensation filter models are determined in the frequency domain for temperature subranges [i=0, . . . , i=n]:

 C_{2} ^{1}(s), C_{2} ^{2}(s), C_{2} ^{3}(s), C_{2} ^{4}(s), . . . , C_{2} ^{n}(s)
 And accordingly, the different transfer functions are transformed each in the time discrete sate space such that for each temperature range a corresponding compensator model, which may be different from the compensator models for the other temperature ranges, may be provided, stored, and applied whenever the ambient temperature is identified to fall within the corresponding temperature subrange.
 In execution, as illustrated in connection with
FIG. 2 , the temperature sensor 2 determines the temperature T and supplies such temperature signal to the compensator 3. The compensator 3 determines in which temperature subrange i the measured temperature T falls into, and applies the compensator model that is defined for such temperature subrange. The sensor signal u(t) is then compensated by such selected compensator model 3. As a result the entire sensor system comprising the sensor 1, the temperature sensor 2 and the compensator 3 not only compensates for response times, it compensates also for temperature dependencies in the sensor signal. By this means, the sensor system can be applied to environments with fast changing measures which environments are also characterized to be operated under varying temperatures.  Whenever temperature is determined to be in subrange i, i.e. if T is in T_{i}, then k=I and the corresponding compensator model is applied:

x(k+1)=A _{d} ^{k}(T)·x(k)+B _{d} ^{k}(T)·u(k) 
y(k)=C _{d} ^{k}(T)·x(k)+D _{d} ^{k}(T)·u(k)  In
FIG. 7 there is illustrated a sample temperature range T divided into i=4 subranges, with [i=0, . . . i=n=4]. The dashed line shows a temperature characteristic taken by the temperature sensor 2 out ofFIG. 2 . The temperature rises first from a temperature in subrange i=2 to a temperature in subrange i=3, then falls back to a temperature in subrange i=2, and further drops to a temperature in subrange i=1. Consequently, during this measurement cycle, three compensator models are applied in sequence to the humidity sensor signal, i.e. the models i=2, i=3, i=2 again and finally i=1 corresponding to the temperature ranges.  The compensator 3 including its models may preferably be implemented in software, in hardware, or in a combination of software and hardware.
 Preferably, the sensor system and the corresponding methods may be applied in the antifogging detection for vehicles. In such application, it is preferred that the sensor 1 is a humidity sensor. The humidity sensor may be arranged on or close to a pane such as a windscreen, for example. In addition, the temperature sensor may be arranged on or close to the pane, too, for measuring the ambient temperature, preferably at the same location the humidity sensor covers.
 The results of the measurements may be used to take action against a fogged windscreen, for, example, and start operating a blower.
Claims (19)
1. A sensor system, comprising
a sensor providing a sensor signal representative of a measure other than temperature, dynamic components of the sensor signal being dependent on temperature,
a temperature sensor for providing a temperature signal, and
a compensation filter receiving the sensor signal and the temperature signal,
wherein the compensation filter is designed for adjusting the dynamic components in the sensor signal subject to the temperature signal, and for providing a compensated sensor signal, and
wherein the compensation filter is modeled based on an inverse of a transfer function of a sensor model of the sensor, which compensation filter contains temperature dependent terms, said model being provided in a spectral domain.
2. The sensor system of claim 1 ,
wherein the sensor model modelled for building the compensation filter is a sensor model of second order.
3. The sensor system of claim 2 ,
wherein the sensor comprises a housing,
wherein the sensor comprises a sensor element sensitive to a component in a medium outside the housing,
wherein the sensor model of second order represents two diffusion processes, a first one of which diffusion processes is a diffusion process from the outside into the housing, a second one of which diffusion processes is a diffusion process from within the housing into the sensor element.
4. The sensor system of claim 2 ,
wherein the transfer function of the sensor model of second order is
where
s denotes the complex Laplace variable,
K, T_{1}, T_{2 }are constants to be identified,
T_{1 }and T_{2 }are time constants of respective diffusion processes, and
K defines a coupling between the two diffusion processes.
5. The sensor system of claim 2 ,
wherein the compensation filter is a compensation filter of second order.
6. The sensor system of claim 4 ,
wherein a transfer function C_{2}(s) of the compensation filter is:
7. The sensor system of claim 1 ,
wherein the sensor model modelled for building the compensation filter is a sensor model of first order.
8. The sensor system of claim 7 ,
wherein the sensor comprises a sensor element sensitive to a component in a medium outside the housing,
wherein the sensor model of first order represents a dominant diffusion process for the medium to reach the sensor element.
9. The sensor system of claim 7 ,
wherein the transfer function of the sensor model of first order is
where
s denotes the complex Laplace variable, and
T_{1 }is a time constant of the respective diffusion process.
10. The sensor system of claim 7 ,
wherein the compensation filter is a compensation filter of first order.
11. The sensor system of claim 9 ,
wherein a transfer function C_{1}(s) of the compensation filter is:
12. The sensor system of claim 1 ,
wherein the compensation filter comprises a temperature continuous compensation model.
13. The sensor system according to claim 1 ,
wherein the sensor is a humidity sensor.
14. The sensor system according to claim 13 ,
wherein the humidity sensor and the temperature sensor are arranged close to a pane of a vehicle for detecting fogging of the pane.
15. The sensor system of claim 1 ,
wherein the compensation filter is designed to decrease a response time of the sensor signal subject to the temperature signal.
16. The sensor system of claim 1 ,
wherein the same order selected for the transfer function of the sensor model is applied to the transfer function of the compensating filter.
17. A method for adjusting a sensor signal, comprising the steps of:
(a) sensing a temperature,
(b) providing a sensor signal representative of a measure other than temperature, dynamic components of the sensor signal being dependent on temperature, and
(c) by a computing device configured for doing so,
(i) adjusting the dynamic components in the sensor signal subject to the temperature sensed, and
(ii) providing a compensated sensor signal,
wherein the computing device comprises a compensation filter receiving the sensor signal and the temperature signal,
wherein the compensation filter is designed for adjusting the dynamic components in the sensor signal subject to the temperature signal, and for providing a compensated sensor signal, and wherein the compensation filter is modeled based on an inverse of a transfer function of a sensor model of the sensor, which compensation filter contains temperature dependent terms, said model being provided in a spectral domain, and wherein said steps are executed on a sensor system according to claim 1 .
18. The method according to claim 17 ,
wherein the dynamic components in the sensor signal are adjusted by applying a temperature continuous compensation model.
19. The method according to claim 17 ,
wherein the sensor signal represents measured humidity, and
wherein the compensated sensor signal is used in an antifogging application in a vehicle.
Priority Applications (4)
Application Number  Priority Date  Filing Date  Title 

EP10005804.9  20100604  
EP10005804.9A EP2392898B1 (en)  20100604  20100604  Sensor system 
US13/068,628 US9341499B2 (en)  20100604  20110517  Sensor system 
US15/097,356 US20160327418A1 (en)  20100604  20160413  Sensor system 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

US15/097,356 US20160327418A1 (en)  20100604  20160413  Sensor system 
Related Parent Applications (1)
Application Number  Title  Priority Date  Filing Date  

US13/068,628 Continuation US9341499B2 (en)  20100604  20110517  Sensor system 
Publications (1)
Publication Number  Publication Date 

US20160327418A1 true US20160327418A1 (en)  20161110 
Family
ID=43302981
Family Applications (2)
Application Number  Title  Priority Date  Filing Date 

US13/068,628 Active 20310928 US9341499B2 (en)  20100604  20110517  Sensor system 
US15/097,356 Abandoned US20160327418A1 (en)  20100604  20160413  Sensor system 
Family Applications Before (1)
Application Number  Title  Priority Date  Filing Date 

US13/068,628 Active 20310928 US9341499B2 (en)  20100604  20110517  Sensor system 
Country Status (4)
Country  Link 

US (2)  US9341499B2 (en) 
EP (1)  EP2392898B1 (en) 
JP (2)  JP5718163B2 (en) 
KR (2)  KR101900242B1 (en) 
Cited By (1)
Publication number  Priority date  Publication date  Assignee  Title 

US20180209929A1 (en) *  20150817  20180726  Denso Corporation  Humidity detector 
Families Citing this family (20)
Publication number  Priority date  Publication date  Assignee  Title 

EP2392898B1 (en) *  20100604  20171213  Sensirion AG  Sensor system 
EP2642289A1 (en)  20120320  20130925  Sensirion AG  Portable electronic device 
EP2682715B1 (en) *  20120702  20150311  Sensirion AG  Portable electronic device 
US9772317B2 (en)  20120726  20170926  Sensirion Ag  Method for operating a portable electronic device 
EP2728327B1 (en)  20121102  20200219  Sensirion AG  Portable electronic device 
EP2762880B1 (en)  20130131  20170315  Sensirion AG  Portable electronic device with breath analyzer 
US9404812B2 (en) *  20130314  20160802  Samsung Electronics Co., Ltd.  Method for detecting environmental value in electronic device and electronic device 
EP2793450B1 (en) *  20130418  20170315  Sensirion AG  System and method for determining sensor accuracy of a portable electronic device 
EP2796841B1 (en) *  20130425  20170201  Sensirion AG  Dynamic compensation of a sensor in a mobile device 
EP2797000A1 (en)  20130425  20141029  Sensirion AG  Portable electronic device 
EP2802128B1 (en)  20130506  20180711  Sensirion AG  Selfcalibrating temperature sensor within a mobile terminal device 
EP2801804B1 (en)  20130506  20180711  Sensirion AG  Selfcalibrating temperature sensor within a mobile terminal device 
EP2808651A1 (en)  20130531  20141203  Sensirion AG  Portable electronic device with integrated temperature sensor being compensated by other sensing data 
EP2808650B1 (en)  20130531  20170322  Sensirion AG  Portable electronic device 
EP2808652B1 (en)  20130531  20161116  Sensirion AG  Portable electronic device with integrated temperature sensor being compensated by other sensing data 
EP2846135B1 (en)  20130910  20200429  Sensirion AG  Portable Electronic Device with Environmental Sensor 
US9551599B2 (en) *  20130923  20170124  Rosemount Inc.  Normalized process dynamics 
EP2910938B1 (en)  20140225  20170906  Sensirion AG  Portable electronic device 
EP2930475B1 (en)  20141222  20171115  Sensirion AG  Flow sensor arrangement 
EP3324164B1 (en) *  20161121  20190123  Sensirion AG  Compensating a sensor signal 
Citations (9)
Publication number  Priority date  Publication date  Assignee  Title 

US20020101135A1 (en) *  20000728  20020801  Marco Giovanardi  Method and device for noise damping 
US20040190570A1 (en) *  20010412  20040930  Bernd Hubner  Method and device for regulating the average wavelength of a laser, especially a semiconductor laser 
US20050008170A1 (en) *  20030506  20050113  Gerhard Pfaffinger  Stereo audiosignal processing system 
US20050283330A1 (en) *  20040616  20051222  Laraia Jose M  Reactive sensor modules using pade' approximant based compensation and providing modulesourced excitation 
US20070288180A1 (en) *  20020719  20071213  Celerity, Inc.  Methods and apparatus for pressure compensation in a mass flow controller 
US20090080287A1 (en) *  20040827  20090326  Westerngeco, L.L.C.  Method for Estimating Absorption Parameter Q(T) 
US20110282542A9 (en) *  20090403  20111117  Certusview Technologies, Llc  Methods, apparatus, and systems for acquiring and analyzing vehicle data and generating an electronic representation of vehicle operations 
US20120255354A1 (en) *  20091223  20121011  Beijing Polymer Sensing Technology Co., Ltd.  Humidity measurement device and method 
US9341499B2 (en) *  20100604  20160517  Sensirion Ag  Sensor system 
Family Cites Families (35)
Publication number  Priority date  Publication date  Assignee  Title 

DE2311103C2 (en)  19730306  19840112  Helmut Dipl.Chem. 8000 Muenchen De Ulrich  
DE2505669C3 (en)  19750211  19790802  Draegerwerk Ag, 2400 Luebeck  
JPS6011849B2 (en)  19760210  19850328  Nippon Electric Co  
US4403296A (en)  19801218  19830906  Electromedics, Inc.  Measuring and determination device for calculating an output determination based on a mathematical relationship between multiple different input responsive transducers 
JPH0319569B2 (en)  19831230  19910315  Hitachi Ltd  
US4607530A (en) *  19841101  19860826  Schlumberger Technology Corporation  Temperature compensation for pressure gauges 
JP2526872B2 (en) *  19860821  19960821  日本電装株式会社  Defroster for window glass 
CA1293568C (en)  19871001  19911224  Romet Limited  Electronic volume correctors 
US4909078A (en)  19871014  19900320  Rosemount Inc.  Fluid flow detector 
US4918995A (en)  19880104  19900424  Gas Research Institute  Electronic gas meter 
US4961348A (en)  19881216  19901009  Ulrich Bonne  Flowmeter fluid composition correction 
US4885938A (en)  19881216  19891212  Honeywell Inc.  Flowmeter fluid composition correction 
GB8903744D0 (en)  19890218  19890405  Endress & Hauser Ltd  Flowmeter 
US4986243A (en)  19900119  19910122  Siemens Automotive L.P.  Mass air flow engine control system with mass air event integrator 
US5237523A (en)  19900725  19930817  Honeywell Inc.  Flowmeter fluid composition and temperature correction 
DE69109236T2 (en)  19901109  19960111  Hewlett Packard Co  Methods and systems for identifying liquids and determining flow. 
US5196851A (en)  19910524  19930323  Samsung Electronics Co., Ltd.  Linearizing nonlinear analogtodigital process and circuit 
US5774378A (en) *  19930421  19980630  The Foxboro Company  Selfvalidating sensors 
JPH08110317A (en)  19941012  19960430  Hitachi Ltd  Integrated micro sensor 
US5815418A (en)  19960131  19980929  Analogic Corporation  Continuous selfcalibrating data acquistion system 
US5951625A (en)  19970630  19990914  Truevision, Inc.  Interpolated lookup table circuit 
JPH11190653A (en)  19971226  19990713  Tokyo Gas Co Ltd  Flow amount compensator 
US6389364B1 (en)  19990710  20020514  Mykrolis Corporation  System and method for a digital mass flow controller 
WO2001018500A1 (en)  19990909  20010315  Sensirion Ag  Method and device for precision mass flow measurement 
US6744389B2 (en)  20000204  20040601  Sensirion Ag  A/Dconverter with lookup table 
AT310986T (en)  20010424  20051215  Celerity Group Inc  Method for determining a valve opening for a mass flow controller 
JP2003098141A (en) *  20010925  20030403  Mitsuba Corp  Dew point measuring device and its manufacturing method 
JP3975443B2 (en) *  20030120  20070912  株式会社山武  Sensor device 
JP4074556B2 (en) *  20030604  20080409  株式会社堀場製作所  Method for improving response speed of analysis system and exhaust gas measurement system 
US7375723B2 (en) *  20030610  20080520  Samsung Electronics Co., Ltd.  Display device and method of compensating primary image data to increase a response speed of the display 
JP4246006B2 (en) *  20030704  20090402  大和製衡株式会社  Weight signal creep error compensation apparatus and compensation method 
JP4432723B2 (en) *  20041021  20100317  株式会社デンソー  Control device 
US7398173B2 (en) *  20050504  20080708  Ami Semiconductor, Inc.  Providing nonlinear temperature compensation for sensing means by use of Padé approximant function emulators 
DE602006008844D1 (en)  20060719  20091015  Sensirion Holding Ag  Moisture sensor to detect fogging on a window 
EP2479681A1 (en)  20110124  20120725  Sensirion AG  Approximating a nonlinear characteristic in a sensor system 

2010
 20100604 EP EP10005804.9A patent/EP2392898B1/en active Active

2011
 20110517 US US13/068,628 patent/US9341499B2/en active Active
 20110602 KR KR1020110053148A patent/KR101900242B1/en active IP Right Grant
 20110603 JP JP2011124966A patent/JP5718163B2/en active Active

2015
 20150318 JP JP2015054969A patent/JP6072112B2/en active Active

2016
 20160413 US US15/097,356 patent/US20160327418A1/en not_active Abandoned

2018
 20180228 KR KR1020180024173A patent/KR101936569B1/en active IP Right Grant
Patent Citations (10)
Publication number  Priority date  Publication date  Assignee  Title 

US20020101135A1 (en) *  20000728  20020801  Marco Giovanardi  Method and device for noise damping 
US20040190570A1 (en) *  20010412  20040930  Bernd Hubner  Method and device for regulating the average wavelength of a laser, especially a semiconductor laser 
US20070288180A1 (en) *  20020719  20071213  Celerity, Inc.  Methods and apparatus for pressure compensation in a mass flow controller 
US20050008170A1 (en) *  20030506  20050113  Gerhard Pfaffinger  Stereo audiosignal processing system 
US20050283330A1 (en) *  20040616  20051222  Laraia Jose M  Reactive sensor modules using pade' approximant based compensation and providing modulesourced excitation 
US7006938B2 (en) *  20040616  20060228  Ami Semiconductor, Inc.  Reactive sensor modules using Pade' Approximant based compensation and providing modulesourced excitation 
US20090080287A1 (en) *  20040827  20090326  Westerngeco, L.L.C.  Method for Estimating Absorption Parameter Q(T) 
US20110282542A9 (en) *  20090403  20111117  Certusview Technologies, Llc  Methods, apparatus, and systems for acquiring and analyzing vehicle data and generating an electronic representation of vehicle operations 
US20120255354A1 (en) *  20091223  20121011  Beijing Polymer Sensing Technology Co., Ltd.  Humidity measurement device and method 
US9341499B2 (en) *  20100604  20160517  Sensirion Ag  Sensor system 
NonPatent Citations (1)
Title 

LARRY M. MILOSHEVICH, "Development and Validation of a TimeLag Correction for Vaisala Radiosonde Humidity Measurements" JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY SEPTEMBER 2004 * 
Cited By (2)
Publication number  Priority date  Publication date  Assignee  Title 

US20180209929A1 (en) *  20150817  20180726  Denso Corporation  Humidity detector 
US10724980B2 (en) *  20150817  20200728  Denso Corporation  Humidity detector 
Also Published As
Publication number  Publication date 

KR101900242B1 (en)  20180920 
KR101936569B1 (en)  20190110 
US20110307208A1 (en)  20111215 
JP5718163B2 (en)  20150513 
EP2392898A1 (en)  20111207 
KR20110133436A (en)  20111212 
US9341499B2 (en)  20160517 
JP2015143702A (en)  20150806 
KR20180027462A (en)  20180314 
JP2012021975A (en)  20120202 
JP6072112B2 (en)  20170201 
EP2392898B1 (en)  20171213 
Similar Documents
Publication  Publication Date  Title 

Londoño et al.  Identification of backbone curves of nonlinear systems from resonance decay responses  
Wang et al.  Data filtering based recursive least squares algorithm for Hammerstein systems using the keyterm separation principle  
Hagenblad et al.  Maximum likelihood identification of Wiener models  
US6714878B2 (en)  System and method for a digital mass flow controller  
US8798281B2 (en)  Control of a loudspeaker output  
RU2328753C2 (en)  Device and method for evaluating state of charge of battery with use of neutral network  
Zheng  A bias correction method for identification of linear dynamic errorsinvariables models  
US7398173B2 (en)  Providing nonlinear temperature compensation for sensing means by use of Padé approximant function emulators  
Adams et al.  A frequency domain method for estimating the parameters of a nonlinear structural dynamic model through feedback  
US6760716B1 (en)  Adaptive predictive model in a process control system  
Bolzon et al.  Material model calibration by indentation, imprint mapping and inverse analysis  
US7499819B2 (en)  Flow measuring device  
TWI673631B (en)  Method of generating at least one filter of a device, touch sensitive device, and nontransitory computer readable medium  
EP3207388B1 (en)  Battery condition monitoring  
Terwilliger  Using primeandswitch phasing to reduce model bias in molecular replacement  
US8024164B1 (en)  Adaptive lookup table: a graphical simulation component for recursively updating numeric data storage in table form  
US7496414B2 (en)  Dynamic controller utilizing a hybrid model  
Gershgorin et al.  Test models for improving filtering with model errors through stochastic parameter estimation  
Devriendt et al.  Operational modal analysis in the presence of harmonic excitations by the use of transmissibility measurements  
JP2008511827A (en)  An effective mathematical model for determining and calibrating the output of linear sensors  
EP1775451B1 (en)  A method for estimating engine friction torque  
US20070271072A1 (en)  Sensor adaptors and methods  
US6041172A (en)  Multiple scale signal processing and control system  
US8898030B2 (en)  Method for adaptively compensating position error of resolver  
CN101346618B (en)  Calibration in a laboratory reference method 
Legal Events
Date  Code  Title  Description 

AS  Assignment 
Owner name: SENSIRION AG, SWITZERLAND Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:GRAF, MARKUS;HESS, VINCENT;NIEDERBERGER, DOMINIK;SIGNING DATES FROM 20110630 TO 20110704;REEL/FRAME:038264/0244 

STCB  Information on status: application discontinuation 
Free format text: ABANDONED  FAILURE TO RESPOND TO AN OFFICE ACTION 