# Improving Accuracy in the Readout of Resistive Sensor Arrays

Authors: José A. Hidalgo-López; Raquel Fernández-Ramos; Jorge Romero-Sánchez; José F. Martín-Canales; Francisco J. Ríos-Gómez
Journal: Journal of Sensors
 (2018)

---

## Abstract

The most efficient way to carry out the reading of a set of resistive sensors is to organize them in an array form. This reduces the number of wires to the sum of the number of rows,M, and columns, N, and reading can be carried out using just M+N multiplexers and a single operational amplifier. The drawback in this procedure is the appearance of crosstalk (the influence of the resistance of some sensors on the measurements of others). Although different proposals have been presented in the literature to reduce this phenomenon, errors in determining the resistance values of each sensor still exist in all the proposals. This article presents a new method to determine these values, which considerably reduces errors without the need for any hardware other than the simplest reading circuits. The method consists of a very fast recursive convergence procedure. The results show that the new proposal outperforms other solutions in the first or second step.


---

## Body

## 1. Introduction

A wide range of applications uses resistive sensor arrays, such as gas detectors [1–3], nuclear electronics [4], temperature sensors [5, 6], thermal anemometry [7], foot plantar application [8], and tactile sensors or artificial electronic skins [9–14]. These arrays are organized in M rows and N columns, so each sensor is connected to the other sensors of the same row through one of its terminals and to the sensors of the same column through the other one. The sensor structure therefore only uses a total of M+N wires, which also allows the multiplexing circuits necessary to read each sensor’s information individually to be reduced. This reduction in reading electronics comes at the cost of the appearance of crosstalk from parasitic parallel paths. Errors caused by this phenomenon in each sensor reading have been widely studied in the literature [15–20] and can reach values as high as 30% for a 4 × 4 array [21].The main cause of the appearance of crosstalk in these arrays is the resistance of the different multiplexers and/or buffers used to select rows and columns in the reading electronics, which cause parasitic parallel paths to appear. A second group of causes is related to the nonidealities of the OAs (operational amplifiers) that are usually part of the reading circuits [19, 20, 22]. Naturally, the magnitude of the errors also depends on the type of reading circuit designed. Two excellent reviews of scanning approaches of two-dimensional resistive sensor arrays can be found in [23, 24]. There are two main ways to proceed when reading an array of resistive sensors: a reading circuit can be used for each array column, or a single circuit for the whole array. The first procedure allows faster reading of the array, but requires N OAs and M multiplexers. The second procedure uses a single OA and M+N multiplexers. This second procedure is slower but reduces both power consumption and total hardware. In [25], crosstalk is reduced using the first procedure one additional OA. However, crosstalk can be eliminated using the first reading procedure by adding calibration resistors and an extra OA to the hardware and carrying out additional array measurements, as shown in [19, 20].The second reading procedure with a single OA leads to two fundamental types of circuit. In the first one, the OA is used in its inverting configuration, only one column wire is connected to the inverter input of the OA, and only one-row wire is connected toVDD (the supply voltage of the array). On the other hand, the rest of the wires are connected to ground. This is the zero potential method (ZPM). In the second type of circuit, the OA output (in noninverting configuration) is used as a feedback for the array; this is known as the voltage feedback method (VFM). However, the appearance of crosstalk is inevitable in both types of circuit. References [17, 18] evaluate errors in the array reading using the VFM structure, and Liu et al. [16] compare errors made by circuits that use different variations of the ZPM and VFM, reaching the conclusion that, in general and under the same operating conditions, errors in circuits based on the ZPM are lower than those based on the VFM.In order to reduce crosstalk in the circuits based on the ZPM and VFM, Wu et al. present two designs: one based on the ZPM [26] using two OAs and an additional resistor (errors in this method are evaluated in [27]) and another based on the VFM [28, 29], which requires two additional resistors ([30] uses this approach for a piezoresistive composite sensor array). Although the results obtained improve measurement accuracy, they still show errors due to crosstalk. The increase in hardware in these proposals is offset by the fact that the value of the M·N resistances of the sensors can be found using uncomplicated algebraic equations, maintaining limited error.This paper will present a method to obtain theM·N values of the resistances of the sensors which, using the simplest versions of the ZPM and VFM, allows errors due to crosstalk to be almost completely eliminated without having to add further elements.The paper is structured in sections as follows: Section2 presents the resistive sensor array reading circuits with a single OA that appears in the literature. Section 3 presents the new method proposed to analyse these arrays, which uses the simplest versions of the ZPM and VFM circuits. Section 4 sets out the simulation results and their discussion. The final section provides a summary of the conclusions.


## 2. Reading an Array of Resistive Sensors Using ZPM and VFM Circuits

The simplest circuit for reading an array of resistive sensors based on the ZPM is shown in Figure1. The reading procedure for this circuit consists of selecting, using the row selection switches, a positive voltage for one of the wires of row i1≤i≤M, while the rest are connected to ground. Likewise, all column wires except one, j1≤j≤N, are connected to ground. Meanwhile, the column switch j connects the wire of the jth array column to the inverting input of the OA. Using this switch layout (if the switches were ideal), there is only current flowing through the sensor connected to row i and column j, since voltages Vri of the row wires and Vcj of the column wires of the array of Figure 1 would be VDD or ground, depending on the selection in the switches. This current could be calculated by the following expression:
(1)Rij=−VDDVoi,jRf,where Voi,j is the output voltage of the OA when the switch of row i selects VDD and the switch of column j connects this column’s wire to the inverting input terminal of the OA. Unfortunately, (1) is merely an approximation to the real value of Rij, since voltages Vri and Vcj do not match VDD or ground due to the resistances of the switches, Rs.Figure 1
Circuit based on the ZPM for reading the array of resistive sensors. The position of the switches allows the reading ofR12.On the other hand, Figure2 shows a circuit that uses the VFM for the array reading. This circuit is a slightly modified version of the classical VFM circuits that uses a single resistor, Rg (as shown in Figure 2), rather than the one with the same value for each of the columns of the array, such as in [16]. The readout procedure consists of using the row selection switches to select a positive voltage, VDD, for one of the row wires, i, while the rest are connected to the output of the OA. Moreover, all column switches except one, j, are connected to the output voltage of the OA. The column switch, j, connects the wire of this column to the resistor, Rg, which in turn is connected to the noninverting input terminal of the OA. Once again, if the row and column switches were ideal, only the current would flow through the resistor Rij, the value of which could be calculated by the following expression:
(2)Rij=VDDVoi,j−1⋅Rg,where Voi,j is the output voltage of the OA when the switch of row i selects VDD and the switch of column j connects this column’s wire to the noninverting input terminal of the OA.Figure 2
Circuit based on the VFM for reading the array of resistive sensors. The position of the switches allows the reading ofR12.As in the case of the ZPM circuit, (2) would only be valid if the resistors of the column and row switches, Rs, are 0. The error that occurs in estimating Rij using these equations will depend on Rs, on the size of the array (M and N), and the range of values of the resistances of the sensors [16, 26].To reduce the errors due toRs in (1), Wu et al. [31] propose modifying the circuit based on the ZPM of Figure 1 by adding a second OA and a resistor (Figure 3). For this circuit, Wu et al. propose calculating Rij, with the value of Rs known, according to the following equation:
(3)Rij=−VDD+Vcg⋅RsRcgRlVLij−2Rs,Figure 3
Modified circuit based on the ZPM proposed by Wu et al. [31] to reduce the influence of Rs in calculating Rij.Vcg and VLij are the output voltages of the OAs in Figure 3 and Rl and Rcg are the feedback resistors of these OAs (as one example in Figure 3, multiplexers are selected to the readout of R11).Although the circuit in Figure3 reduces the number of parasitic parallel paths due to Rs, these are not completely eliminated and, therefore, expression (3) still shows errors in estimating Rij. These errors again depend on Rs, on the size of the array, and on the resistors to be measured.Based on the VFM, in [26], Wu et al. propose a new design (Figure 4) in order to reduce crosstalk. Wu et al. propose evaluating Rij according to the following equation:
(4)Rij=VDDVoi,j−1⋅Rg+VDDVoi,j⋅Rs.Figure 4
Modified circuit based on the VFM proposed by Wu et al. [26] to reduce the effects of crosstalk.The second term of the member on the right is the only difference from (2). Although, for some ranges of Rij and Rs, the errors using (4) for the circuit of Figure 4 are smaller than those obtained using (2) for the circuit (Figure 2), these errors do not disappear.


## 3. Method for Reducing Crosstalk Errors in Reading an Array of Resistive Sensors

To reduce errors, due to the resistances of the switches that appear in (1), (2), (3), and (4), a new method to calculate Rij is presented in this section. The method is the same, with minor modifications, for the simplest version of the ZPM- and VFM-based circuits (Figures 1 and 2). In essence, the method uses the estimates values of Rij by using (1) and (2) as the starting point to obtain the voltages in the row and column wires of the array (Vri and Vcj) that will be used to, in successive steps, improve the estimation of Rij. The method is described below in its two variants.


### 3.1. Method for Reducing Crosstalk Errors in Reading Circuits Based on ZPM

For the circuit of Figure1, the proposed method begins by obtaining an initial estimation of all Rij, Rij0, through (1):
(5)Rij0=−VDDVoi,jRf,where the superscript will indicate the current step of the process. Next, if resistors of multiplexers, Rs, are considered for improving accuracy in the readout of Rij, it is necessary to analyse the equivalent circuit in Figure 5. In this circuit, the current flowing towards the Rf resistor when connecting the wire of row i to VDD and the wire of column j to the OA, Ioi,j, is
(6)Ioi,j=−Voi,jRf=Iji,where Iji is the current flowing through the multiplexer of the jth column wire (see Figure 5). It should be noted that Iji does not depend on the column selected for reading since, in any selection state of the column switches, these connect the columns of the array at a voltage of 0 V. This current is used to calculate the voltage values of the column wires in Figure 5,
(7)Vcji=Iji⋅Rs=−Voi,jRsRf,where Vcji is the voltage of the wire of column j when connecting row i to VDD.Figure 5
Equivalent circuit, based on the ZPM, including the resistor of multiplexersRs for improving accuracy in the readout of Rij.The recursive process begins by calculating the voltages of any row wire in the array,k, when VDD is connected row wire i: Vrk1i (note that Vrk do not depend on which column is selected). This voltage can be calculated using Millman’s theorem [32] and (7) as
(8)Vrk1i=VDD⋅Ys⋅δik+∑p=1NVcpi⋅Ykp0∑p=1NYkp0+Ysand replacing Vcpi from (7):
(9)Vrk1i=VDD⋅Ys⋅δik−Rs/Rf⋅∑p=1NVoi,pYkp0∑p=1NYkp0+Ys,Ykp0 and Ys are the inverse of Rkp0 and Rs, respectively, while δik is the Kronecker delta. On the other hand, the current passing through Rij in this first step, IRij1 (see Figure 5), is
(10)IRij1=Ioi,j−∑k=1,k≠iMIRkj1=Ioi,j−∑k=1,k≠iMVrk1i−Vcji⋅Ykj0=−Voi,jRf1+Rs⋅∑k=1,k≠iMYkj0−∑k=1,k≠iMVrk1i⋅Ykj0.Now,Rij1 can be found using the results of (5), (6), (7), (8), and (9):
(11)Rij1=Vri1i−VcjiIRij1=Vri1i+Voi,jRs/Rf−Voi,j/Rf1+Rs⋅∑k=1,k≠iMYkj0−∑k=1,k≠iMVrk1i⋅Ykj0.Equations (9) and (11)are the recursive equations of the proposed method. Thus, in the qth step, new values of Vrkqi would be calculated using (9), with Rijq−1 from the q−1 step (instead of Rij0). The step finishes using (11) to calculate Rijq with Vrkqi instead of Vrk1i.


### 3.2. Method for Reducing Crosstalk Errors in Reading Circuits Based on VFM

For the circuit in Figure2, the method is very similar to the one presented in the previous subsection, beginning with the calculation of all the Rij0 from (2). Again, if Rs is considered, it is necessary to analyse the circuits of Figure 6 for improving accuracy. In this circuit, the current flowing towards the Rg resistor when switching row i to VDD and connecting column j to the OA, Ioi,j, is calculated:
(12)Ioi,j=Voi,jRg.Figure 6
Equivalent circuit, based on the VFM, including the resistor of multiplexersRs for improving accuracy in the readout of Rij.For this circuit (if any rowi is selected), it is possible to find the current flowing through the switch of a selected column t, Iti,t as a function of the current flowing through the switch of the same column, when another column j was selected:
(13)Iti,jIti,t=VDD−Voi,jVDD−Voi,t.This equation shows the effect of the voltage scaling that has taken place in the array when columnt or column j was selected (keeping i as the selected row). Having in mind that Iti,t=Ioi,t, (12) and (13) can be used to calculate the voltage in the wire of the column t when selecting row i and column j, Vcti,j:
(14)Vcti,j=Iti,j⋅Rs+Voi,j=Voi,t⋅RsRg⋅VDD−Voi,jVDD−Voi,t+Voi,j.The recursive process for this circuit starts by using Millman’s theorem and (14) to recalculate the voltages of kth row of the array when selecting for reading row i and column j, Vrk1i,j:
(15)Vrk1i,j=∑t=1NVcti,j⋅Ykt0∑t=1NYkt0+Ys+VDD⋅δik+Voi,j⋅1−δik⋅Ys∑t=1NYkt0+Ys or
(16)Vrk1i,j=Voi,j⋅1+Rs/Rg⋅∑t=1NVDD−Voi,j/VDD−Voi,t⋅Voi,t/Voi,j⋅Ykt0∑t=1NYkt0+Ys+VDD−Voi,j⋅δik⋅Ys∑t=1NYkt0+Ys.The current passing throughRij when selected row i and column j (in the first step of the recursive method) IRij1i,j is (see Figure 6):
(17)IRij1i,j=Ioi,j−∑k=1,k≠iMIRkj1i,j=Ioi,j−∑k=1,k≠iMVcji,j−Vrk1i,j⋅Ykj0=Voi,jRg−Voi,j⋅1+RsRg⋅∑k=1,k≠iMYkj0+∑k=1,k≠iMVrk1i,j⋅Ykj0.Thus,Rij1 can be calculated as
(18)Rij1=Vri1i,j−Vcji,jIRij1i,j=Vri1i,j−Voi,j⋅1+Rs/RgVoi,j/Rg−Voi,j⋅1+Rs/Rg⋅∑k=1,k≠iMYkj0+∑k=1,k≠iMVrk1i,j⋅Ykj1.For VFM-based circuits, (16) and (18) are the recursive equations of the proposed method. Thus, in the qth step, new values of Vrkqi,j would be calculated using (16), with Rijq−1 founded in the q−1 step (instead of Rij0). The step finishes using (18) to calculate Rijq with Vrkqi,j instead of Vrk1i,j.


## 3.1. Method for Reducing Crosstalk Errors in Reading Circuits Based on ZPM

For the circuit of Figure1, the proposed method begins by obtaining an initial estimation of all Rij, Rij0, through (1):
(5)Rij0=−VDDVoi,jRf,where the superscript will indicate the current step of the process. Next, if resistors of multiplexers, Rs, are considered for improving accuracy in the readout of Rij, it is necessary to analyse the equivalent circuit in Figure 5. In this circuit, the current flowing towards the Rf resistor when connecting the wire of row i to VDD and the wire of column j to the OA, Ioi,j, is
(6)Ioi,j=−Voi,jRf=Iji,where Iji is the current flowing through the multiplexer of the jth column wire (see Figure 5). It should be noted that Iji does not depend on the column selected for reading since, in any selection state of the column switches, these connect the columns of the array at a voltage of 0 V. This current is used to calculate the voltage values of the column wires in Figure 5,
(7)Vcji=Iji⋅Rs=−Voi,jRsRf,where Vcji is the voltage of the wire of column j when connecting row i to VDD.Figure 5
Equivalent circuit, based on the ZPM, including the resistor of multiplexersRs for improving accuracy in the readout of Rij.The recursive process begins by calculating the voltages of any row wire in the array,k, when VDD is connected row wire i: Vrk1i (note that Vrk do not depend on which column is selected). This voltage can be calculated using Millman’s theorem [32] and (7) as
(8)Vrk1i=VDD⋅Ys⋅δik+∑p=1NVcpi⋅Ykp0∑p=1NYkp0+Ysand replacing Vcpi from (7):
(9)Vrk1i=VDD⋅Ys⋅δik−Rs/Rf⋅∑p=1NVoi,pYkp0∑p=1NYkp0+Ys,Ykp0 and Ys are the inverse of Rkp0 and Rs, respectively, while δik is the Kronecker delta. On the other hand, the current passing through Rij in this first step, IRij1 (see Figure 5), is
(10)IRij1=Ioi,j−∑k=1,k≠iMIRkj1=Ioi,j−∑k=1,k≠iMVrk1i−Vcji⋅Ykj0=−Voi,jRf1+Rs⋅∑k=1,k≠iMYkj0−∑k=1,k≠iMVrk1i⋅Ykj0.Now,Rij1 can be found using the results of (5), (6), (7), (8), and (9):
(11)Rij1=Vri1i−VcjiIRij1=Vri1i+Voi,jRs/Rf−Voi,j/Rf1+Rs⋅∑k=1,k≠iMYkj0−∑k=1,k≠iMVrk1i⋅Ykj0.Equations (9) and (11)are the recursive equations of the proposed method. Thus, in the qth step, new values of Vrkqi would be calculated using (9), with Rijq−1 from the q−1 step (instead of Rij0). The step finishes using (11) to calculate Rijq with Vrkqi instead of Vrk1i.


## 3.2. Method for Reducing Crosstalk Errors in Reading Circuits Based on VFM

For the circuit in Figure2, the method is very similar to the one presented in the previous subsection, beginning with the calculation of all the Rij0 from (2). Again, if Rs is considered, it is necessary to analyse the circuits of Figure 6 for improving accuracy. In this circuit, the current flowing towards the Rg resistor when switching row i to VDD and connecting column j to the OA, Ioi,j, is calculated:
(12)Ioi,j=Voi,jRg.Figure 6
Equivalent circuit, based on the VFM, including the resistor of multiplexersRs for improving accuracy in the readout of Rij.For this circuit (if any rowi is selected), it is possible to find the current flowing through the switch of a selected column t, Iti,t as a function of the current flowing through the switch of the same column, when another column j was selected:
(13)Iti,jIti,t=VDD−Voi,jVDD−Voi,t.This equation shows the effect of the voltage scaling that has taken place in the array when columnt or column j was selected (keeping i as the selected row). Having in mind that Iti,t=Ioi,t, (12) and (13) can be used to calculate the voltage in the wire of the column t when selecting row i and column j, Vcti,j:
(14)Vcti,j=Iti,j⋅Rs+Voi,j=Voi,t⋅RsRg⋅VDD−Voi,jVDD−Voi,t+Voi,j.The recursive process for this circuit starts by using Millman’s theorem and (14) to recalculate the voltages of kth row of the array when selecting for reading row i and column j, Vrk1i,j:
(15)Vrk1i,j=∑t=1NVcti,j⋅Ykt0∑t=1NYkt0+Ys+VDD⋅δik+Voi,j⋅1−δik⋅Ys∑t=1NYkt0+Ys or
(16)Vrk1i,j=Voi,j⋅1+Rs/Rg⋅∑t=1NVDD−Voi,j/VDD−Voi,t⋅Voi,t/Voi,j⋅Ykt0∑t=1NYkt0+Ys+VDD−Voi,j⋅δik⋅Ys∑t=1NYkt0+Ys.The current passing throughRij when selected row i and column j (in the first step of the recursive method) IRij1i,j is (see Figure 6):
(17)IRij1i,j=Ioi,j−∑k=1,k≠iMIRkj1i,j=Ioi,j−∑k=1,k≠iMVcji,j−Vrk1i,j⋅Ykj0=Voi,jRg−Voi,j⋅1+RsRg⋅∑k=1,k≠iMYkj0+∑k=1,k≠iMVrk1i,j⋅Ykj0.Thus,Rij1 can be calculated as
(18)Rij1=Vri1i,j−Vcji,jIRij1i,j=Vri1i,j−Voi,j⋅1+Rs/RgVoi,j/Rg−Voi,j⋅1+Rs/Rg⋅∑k=1,k≠iMYkj0+∑k=1,k≠iMVrk1i,j⋅Ykj1.For VFM-based circuits, (16) and (18) are the recursive equations of the proposed method. Thus, in the qth step, new values of Vrkqi,j would be calculated using (16), with Rijq−1 founded in the q−1 step (instead of Rij0). The step finishes using (18) to calculate Rijq with Vrkqi,j instead of Vrk1i,j.


## 4. Results and Discussion

To compare the results provided by the proposed methods to those provided by the classical equations of the ZPM and VFM circuits, (1) and (2), and also the results of the circuits proposed by Wu et al., (3) and (4), a batch of simulations has been carried out using Cadence Orcad-Pspice 16.6. A general purpose OA, the Texas Instrument OPA4188, was used for the readout circuit in every test with the power supplies adjusted to +5 V and −5 V. The OA is designed with autozeroing techniques to provide low offset voltage (25 μV, maximum) and near zero-drift over time and temperature. The DC gain is 136 dB, and the input bias current is 16 pA. Errors due to the nonidealities of the OA are minimized thanks to these characteristics. If, for example, the offset voltage of the OA was higher, the technique reported in [22] could be used to reduce its influence. The sensing resistor Rij range used for simulation was [100 Ω, 10 kΩ]. This range is typically for touch sensors [33]. Another reason to use this range of resistors is that it shows the difference in operation and limitations of circuits based on the ZPM and VFM, as will be shown in the results. For this range of sensor resistors, Rf=75 Ω, Rg=350 Ω, Rl=75 Ω, and Rcg=75 Ω have been selected in order not to saturate the output voltages of the OA and to ensure the difference between the maximum and the minimum Vo which is approximately equal in all circuits.In general, simulations will be presented to compare the operation of the circuits of Figures1 and 2 (evaluated, resp., by (1) and (2) and hereinafter classical approach), the two circuits of Figures 3 and 4 proposed by Wu et al. (evaluated, resp., by (3) and (4) and hereinafter ZPM-Wu and VFM-Wu), and the proposed method presented in the previous section.Table1 shows simulation results for the different circuits and methods presented based on the ZPM. The aim of this series of simulations is to analyse the influence of nonscanned resistors, Rns, on the measurement of the resistor, Rij. In all cases, it is an 8 × 8 array in which the resistance of the switches is set to 1 Ω (a resistance that is relatively easy to obtain for switches). Different values of Rns are used for each value of the resistor to be measured, Rij. Both Rns and Rij vary in the selected range [100 Ω, 10 kΩ].Table 1
Comparison of the results obtained in an 8 × 8 array with different values ofRij if Rs=1 Ω and the rest of the resistors of the array, Rns, take different values. (a) Classical approach. (b) ZPM circuit proposed by Wu et al. (c) Results for the proposed method in the first step. (d) Results for the proposed method in the second step.
(a)
ZPM classical approach
Rs = 1 ΩRns (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
100
16.00
12.39
5.48
5.68
13.30
19.78
200
9.10
7.19
5.95
4.07
2.59
1.16
300
6.73
5.02
4.42
3.77
3.27
2.79
400
5.55
3.88
3.44
3.11
2.88
2.65
600
4.36
2.73
2.38
2.22
2.12
2.03
800
3.77
2.15
1.82
1.70
1.64
1.58
1200
3.18
1.57
1.25
1.15
1.11
1.08
1800
2.79
1.18
0.86
0.78
0.74
0.72
3000
2.47
0.87
0.55
0.47
0.44
0.42
6000
2.24
0.63
0.32
0.24
0.21
0.19
10000
2.14
0.54
0.23
0.15
0.12
0.10
(b)
ZPM-Wu
Rs = 1 ΩRns (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
100
97.34
68.28
92.43
97.00
98.00
98.50
200
15.26
84.06
95.65
98.21
98.80
99.10
300
10.83
72.51
91.38
96.37
97.55
98.15
400
18.78
54.82
83.04
92.46
94.85
96.08
600
1.16
1.15
1.09
0.99
0.90
0.81
800
0.87
0.87
0.84
0.78
0.73
0.68
1200
0.58
0.58
0.56
0.53
0.50
0.48
1800
0.39
0.39
0.37
0.35
0.33
0.32
3000
0.23
0.23
0.22
0.20
0.19
0.17
6000
0.12
0.11
0.11
0.09
0.07
0.06
10000
0.07
0.07
0.06
0.04
0.03
0.01
(c)
ZPM proposed method, 1st step
Rs = 1 ΩRns (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
100
1.83
1.22
0.35
3.30
5.62
7.84
200
0.60
0.45
0.28
0.01
0.26
0.50
300
0.32
0.22
0.16
0.06
0.02
0.10
400
0.21
0.13
0.09
0.03
0.02
0.06
600
0.12
0.06
0.04
0.00
0.03
0.06
800
0.08
0.04
0.01
0.02
0.04
0.07
1200
0.06
0.02
0.00
0.03
0.05
0.07
1800
0.04
0.01
0.01
0.03
0.05
0.07
3000
0.03
0.00
0.01
0.03
0.05
0.06
6000
0.02
0.00
0.01
0.03
0.05
0.06
10000
0.02
0.00
0.01
0.03
0.04
0.06
(d)
ZPM proposed method, 2nd step
Rs = 1 ΩRns (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
100
0.24
0.14
0.13
0.64
1.06
1.48
200
0.04
0.02
0.02
0.10
0.17
0.23
300
0.02
0.00
0.02
0.07
0.12
0.16
400
0.01
0.00
0.02
0.06
0.10
0.13
600
0.00
0.00
0.02
0.05
0.08
0.11
800
0.00
0.00
0.02
0.05
0.07
0.09
1200
0.00
0.00
0.02
0.04
0.06
0.08
1800
0.00
0.00
0.01
0.04
0.06
0.07
3000
0.00
0.00
0.01
0.03
0.05
0.07
6000
0.00
0.00
0.01
0.03
0.05
0.06
10000
0.00
0.00
0.01
0.03
0.04
0.06In the first place, it should be noted that the errors forRns<600 Ω are very large in the case of the ZPM-Wu circuit, Table 1(b). This is due to the additional OA decreasing its output voltage as the current flowing through the parasitic parallel paths increases (which occurs for low values of Rns), reaching a point at which its output becomes saturated. On the other hand, the four tables of Table 1 show how the errors are greater as Rns decreases. This is justified again by the increase in the current in the parasitic parallel paths with lower values of Rns. With the exception of cases in which Rns is low, the errors in the ZPM-Wu circuit are always lower than those that appear in the ZPM evaluated using the classical approach (Table 1(a)). For its part, the proposed method outperforms the previous methods. It can be seen how the convergence of this method is very fast, since better results are obtained in the first step (Table 1(c)) than with the other two methods (except for errors below 0.03% in ZPM-Wu). In the second step, the maximum error obtained is only 1.48%, more than an order of magnitude lower than that obtained by the classical approach, and with better results in all situations. It also improves the results of the ZPM-Wu circuit in all situations for errors over 0.06%. The results provided in Table 1(a) quite closely match those reported in [10] for a 10 × 10 array, in which the aim is to measure the value of a sensor when the rest of the sensors has resistances between 3 kΩ and 19 kΩ.The aim of the simulations presented in Table2 is to analyse the influence of the resistances of the switches, Rs, on the reading of Rij. In the first place, it is observed how the errors increase for the classical approach and for ZPM-Wu when Rs is increased. The same happens when the value of Rij decreases. The proposed method also shows these trends, except when the errors are very low (below 0.06%). For any combination of Rns and Rs, the results obtained by the proposed method in the second step outperform the results of the other methods. Even in the first step, except for Rij=100 Ω and Rs>9 Ω, the results are better than in the other methods. It should be noted that the maximum error in the second step is only 2.12% when the resistances of the switches are a quarter of the value to be measured.Table 2
Comparison of the results obtained in an 8 × 8 array with different values ofRij if Rns=5 kΩ and the resistances of the switches take different values. (a) Classical approach. (b) ZPM circuit proposed by Wu e al. (c) Results for the proposed method in the first step. (d) Results for the proposed method in the second step.
(a)
ZPM classical approach
Rns = 5 kΩRs (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
1
2.28
0.68
0.37
0.29
0.26
0.24
2
4.57
1.36
0.75
0.61
0.57
0.54
4
9.17
2.73
1.51
1.25
1.18
1.13
6
13.79
4.10
2.27
1.89
1.79
1.73
9
20.76
6.17
3.42
2.85
2.70
2.62
15
34.86
10.35
5.73
4.77
4.52
4.38
20
46.77
13.87
7.67
6.36
6.03
5.84
25
58.82
17.42
9.61
7.96
7.54
7.29
(b)
ZPM-Wu
Rns = 5 kΩRs (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
1
0.14
0.14
0.13
0.11
0.10
0.08
2
0.28
0.28
0.27
0.25
0.24
0.22
4
0.56
0.56
0.55
0.53
0.51
0.49
6
0.84
0.84
0.82
0.80
0.78
0.77
9
1.26
1.25
1.24
1.21
1.19
1.17
15
2.09
2.09
2.06
2.02
1.98
1.95
20
2.79
2.78
2.75
2.68
2.63
2.57
25
3.48
3.47
3.42
3.33
3.26
3.18
(c)
ZPM proposed method, 1st step
Rns = 5 kΩRs (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
1
0.03
0.00
0.01
0.03
0.05
0.06
2
0.10
0.01
0.01
0.03
0.04
0.06
4
0.38
0.04
0.01
0.02
0.03
0.05
6
0.80
0.09
0.03
0.00
0.02
0.03
9
1.67
0.20
0.08
0.04
0.02
0.00
15
4.07
0.55
0.23
0.16
0.14
0.11
20
6.56
0.94
0.41
0.31
0.27
0.24
25
9.39
1.43
0.63
0.48
0.44
0.40
(d)
ZPM proposed method, 2nd step
Rns = 5 kΩRs (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
1
0.00
0.00
0.01
0.03
0.05
0.06
2
0.00
0.00
0.01
0.03
0.05
0.06
4
0.02
0.00
0.01
0.03
0.05
0.06
6
0.06
0.00
0.01
0.03
0.05
0.06
9
0.17
0.00
0.01
0.03
0.05
0.06
15
0.63
0.03
0.00
0.03
0.04
0.06
20
1.26
0.07
0.01
0.02
0.03
0.05
25
2.12
0.14
0.03
0.00
0.02
0.04Table3 shows simulation results for the different circuits and methods presented based on the VFM. The aim of this series of simulations is to analyse the influence of nonscanned resistors, Rns, on the measurement of the resistor, Rij. Again, it is an 8 × 8 array in which the resistance of the switches is set to 1 Ω. Different values of Rns are used for each value of the resistor to be measured, Rij. Both Rns and Rij vary in the selected range [100 Ω, 10 kΩ]. The first question to consider is that the errors for any method are very large if Rns<600 Ω. This comes about since the parasitic current flowing through the Rns of row i is excessive for the OA (even though the feedback resistor, Rg, has been chosen so there is no saturation of voltages of the OA in the selected range of individual resistors). The situation gets worse with increased Rij, since Voi,j decreases, and in consequence, the current flowing through all the Rns of row i increases. This situation does not occur in the ZPM circuits since the row selection switches are responsible for providing this current. The range limitation of Rns is therefore a drawback of VFM-based circuits compared to those based on the ZPM. It is worth to note that situations with all Rns very low are uncommon in arrays of sensors. For example, in [33], a 16 × 16 tactile sensor array designed with a continuous electroactive material using a piezoresistive sheet of capLINQ (code MVCF-40012BT50KS/2A) has 10 kΩ for no-pressed tactels; however, only the few pressed ones might have resistances around 300 Ω.Table 3
Results obtained for methods based on the VFM in an 8 × 8 array with different values ofRij if Rs=1 Ω and the rest of the resistors in the array, Rns, take different values. (a) Classical approach. (b) VFM circuit proposed by Wu et al. (c) Results for the proposed method in the first step. (d) Results for the proposed method in the second step.
(a)
VFM classical approach
Rs = 1 ΩRns (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
100
13.85
80.64
95.03
98.00
98.67
99.00
200
9.17
54.87
88.01
95.14
96.75
97.56
300
6.79
14.46
76.08
90.19
93.42
95.05
400
5.60
3.92
51.04
79.41
86.12
89.53
600
4.40
2.76
2.43
2.30
2.24
2.18
800
3.81
2.18
1.86
1.78
1.75
1.72
1200
3.21
1.59
1.28
1.22
1.21
1.20
1800
2.81
1.19
0.89
0.84
0.83
0.82
3000
2.50
0.88
0.58
0.52
0.52
0.51
6000
2.26
0.64
0.34
0.29
0.28
0.28
10000
2.16
0.55
0.25
0.19
0.18
0.18
(b)
VFM-Wu
Rs = 1 ΩRns (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
100
0.19
73.03
93.09
97.22
98.14
98.61
200
0.54
37.21
83.35
93.25
95.49
96.61
300
0.66
0.19
66.77
86.38
90.87
93.13
400
0.69
0.12
31.97
71.40
80.71
85.45
600
0.72
0.09
0.27
0.36
0.41
0.47
800
0.72
0.08
0.24
0.29
0.31
0.33
1200
0.73
0.07
0.22
0.26
0.26
0.27
1800
0.73
0.07
0.22
0.25
0.25
0.25
3000
0.74
0.07
0.22
0.25
0.25
0.25
6000
0.74
0.07
0.22
0.25
0.25
0.25
10000
0.74
0.07
0.22
0.25
0.25
0.25
(c)
VFM proposed method, 1st step
Rs = 1 ΩRns (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
100
27.48
83.53
95.76
98.30
98.86
99.15
200
0.31
58.05
88.85
95.48
96.98
97.73
300
0.04
18.45
77.19
90.64
93.72
95.28
400
0.03
0.17
52.71
80.11
86.59
89.88
600
0.12
0.09
0.05
0.02
0.07
0.08
800
0.12
0.06
0.05
0.06
0.08
0.02
1200
0.08
0.04
0.04
0.04
0.05
0.05
1800
0.07
0.02
0.03
0.03
0.03
0.04
3000
0.05
0.02
0.02
0.02
0.03
0.03
6000
0.05
0.01
0.01
0.02
0.02
0.02
10000
0.04
0.01
0.01
0.01
0.02
0.02
(d)
VFM proposed method, 2nd step
Rs = 1 ΩRns (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
100
29.45
83.93
95.87
98.34
98.89
99.17
200
1.00
58.27
88.91
95.50
96.99
97.74
300
0.38
18.63
77.24
90.66
93.74
95.29
400
0.19
0.04
52.77
80.14
86.60
89.89
600
0.00
0.03
0.03
0.03
0.04
0.04
800
0.04
0.02
0.01
0.03
0.04
0.04
1200
0.03
0.02
0.02
0.02
0.03
0.04
1800
0.03
0.01
0.02
0.01
0.03
0.03
3000
0.02
0.01
0.01
0.01
0.02
0.03
6000
0.02
0.01
0.01
0.01
0.02
0.02
10000
0.02
0.01
0.01
0.01
0.02
0.02If this saturation does not come about (due to the current to be provided by the OA), the proposed method outperforms the other methods in both the first and the second steps for all combinations ofRij and Rns. Finally, the results obtained by the proposed method are very similar in the VFM and ZPM circuits (outside of the saturation range of the VFM circuit).Table4 shows the influence of the resistances of the switches, Rs, on the error in estimating Rij for the methods based on the VFM. Once again, the error increases with the value of Rs and, except for VFM-Wu, also increases with decreasing Rij. The biggest errors occur, as in all previous cases, for the classical approach, while the proposed method presents the best results for all combinations, even in the first step. Comparing the errors obtained for ZPM-based methods to the errors of VFM-based methods, it is found that they are quite similar for the classical approach. The errors are similar (although slightly higher) in both steps for the proposed methods and a great deal higher in the case of VFM-Wu.Table 4
Results obtained for methods based on the VFM in an 8 × 8 array with different values ofRij if Rns=5 kΩ and the resistances of the switches take different values. (a) Classical approach. (b) VFM circuit proposed by Wu et al. (c) Results for the proposed method in the first step. (d) Results for the proposed method in the second step.
(a)
VFM classical approach
Rns = 5 kΩRs (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
1
2.31
0.69
0.39
0.34
0.33
0.32
2
4.59
1.37
0.77
0.66
0.63
0.62
4
9.19
2.74
1.53
1.30
1.24
1.22
6
13.81
4.12
2.30
1.94
1.86
1.82
9
20.78
6.19
3.45
2.90
2.77
2.71
15
34.88
10.36
5.76
4.82
4.59
4.48
20
46.79
13.88
7.69
6.42
6.11
5.94
25
58.85
17.43
9.64
8.02
7.62
7.39
(b)
VFM-Wu
Rns = 5 kΩRs (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
1
0.74
0.07
0.22
0.25
0.25
0.25
2
0.16
0.64
0.79
0.81
0.82
0.82
4
0.98
1.76
1.90
1.93
1.93
1.93
6
2.12
2.85
2.99
3.02
3.02
3.03
9
3.83
4.47
4.59
4.62
4.63
4.63
15
7.23
7.58
7.66
7.69
7.71
7.72
20
10.05
10.07
10.09
10.12
10.15
10.18
25
12.86
12.46
12.41
12.45
12.49
12.54
(c)
VFM proposed method, 1st step
Rns = 5 kΩRs (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
1
0.05
0.01
0.01
0.02
0.02
0.02
2
0.12
0.02
0.01
0.02
0.02
0.03
4
0.40
0.05
0.03
0.03
0.03
0.04
6
0.82
0.11
0.05
0.05
0.05
0.05
9
1.70
0.22
0.10
0.09
0.09
0.09
15
4.09
0.56
0.25
0.21
0.21
0.21
20
6.59
0.96
0.44
0.36
0.34
0.34
25
9.42
1.45
0.65
0.54
0.51
0.50
(d)
VFM proposed method, 2nd step
Rns = 5 kΩRs (Ω)
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
1
0.02
0.01
0.01
0.02
0.02
0.02
2
0.02
0.01
0.01
0.02
0.02
0.02
4
0.04
0.01
0.01
0.02
0.02
0.02
6
0.08
0.01
0.01
0.02
0.02
0.02
9
0.19
0.02
0.01
0.02
0.02
0.03
15
0.66
0.05
0.03
0.02
0.03
0.03
20
1.29
0.09
0.04
0.04
0.04
0.04
25
2.15
0.16
0.06
0.05
0.05
0.06Two sets of simulations have been carried out for the circuits of the ZPM of Table5 and the VFM of Table 6 in order to analyse the errors by changing the size of the array. In both simulations, Rs=5 Ω and Rns=2 kΩ. For all methods, errors increase for bigger arrays. The classical approach shows higher errors than Wu proposals, while the proposed methods give the lowest errors. On the other hand, very similar errors appear for the classical approach in both circuits while errors in ZPM-Wu are lower than those obtained in VFM-Wu. It is worth noting that the proposed method shows the best results in all cases and, when the size increases, the advantage grows. So when the size varies from 4 × 4 to 16 × 16, errors for the proposed methods in the biggest array are less than the error for the smallest array with the other approaches.Table 5
Results obtained for methods based on the ZPM by changing the size of the array and maintainingRns=2 kΩ and Rs=5 Ω. (a) Classical approach. (b) ZPM circuit proposed by Wu et al. (c) Results for the proposed method in the first step. (d) Results for the proposed method in the second step.
(a)
ZPM classical approach
Rns = 2 kΩSize
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
4 × 4
11.58
3.51
1.99
1.66
1.57
1.51
6 × 6
12.63
4.53
2.99
2.64
2.53
2.46
8 × 8
13.69
5.55
3.98
3.61
3.49
3.39
10 × 10
14.76
6.57
4.98
4.58
4.43
4.31
12 × 12
15.83
7.59
5.98
5.54
5.35
5.20
14 × 14
16.91
8.62
6.98
6.49
6.26
6.07
16 × 16
17.99
9.66
7.98
7.43
7.16
6.93
(b)
ZPM-Wu
Rns = 2 kΩSize
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
4 × 4
0.75
0.74
0.73
0.70
0.68
0.66
6 × 6
1.25
1.24
1.22
1.18
1.14
1.11
8 × 8
1.74
1.74
1.70
1.64
1.58
1.53
10 × 10
2.24
2.23
2.18
2.09
2.01
1.93
12 × 12
2.74
2.72
2.65
2.52
2.41
2.30
14 × 14
3.24
3.21
3.12
2.95
2.80
2.65
16 × 16
3.73
3.70
3.59
3.36
3.17
2.98
(c)
ZPM proposed method, 1st step
Rns = 2 kΩSize
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
4 × 4
0.58
0.07
0.02
0.00
0.02
0.04
6 × 6
0.70
0.13
0.06
0.03
0.01
0.01
8 × 8
0.84
0.20
0.12
0.09
0.06
0.04
10 × 10
0.99
0.29
0.20
0.16
0.13
0.10
12 × 12
1.15
0.40
0.29
0.24
0.21
0.17
14 × 14
1.33
0.52
0.40
0.33
0.29
0.25
16 × 16
1.52
0.66
0.52
0.44
0.39
0.34
(d)
ZPM proposed method, 2nd step
Rns = 2 kΩSize
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
4 × 4
0.03
0.00
0.01
0.03
0.05
0.06
6 × 6
0.05
0.00
0.01
0.03
0.05
0.07
8 × 8
0.06
0.00
0.01
0.03
0.05
0.07
10 × 10
0.08
0.01
0.01
0.03
0.05
0.07
12 × 12
0.10
0.02
0.00
0.03
0.05
0.07
14 × 14
0.13
0.03
0.01
0.02
0.05
0.07
16 × 16
0.16
0.04
0.02
0.01
0.04
0.06Table 6
Results obtained for methods based on the VFM by changing the size of the array and maintainingRns=2 kΩ and Rs=5 Ω. (a) Classical approach. (b) VFM circuit proposed by Wu et al. (c) Results for the proposed method in the first step. (d) Results for the proposed method in the second step.
(a)
VFM classical approach
Rns = 2 kΩSize
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
4 × 4
11.60
3.53
2.01
1.71
1.64
1.60
6 × 6
12.66
4.54
3.01
2.69
2.61
2.56
8 × 8
13.72
5.56
4.02
3.67
3.57
3.50
10 × 10
14.79
6.59
5.02
4.64
4.52
4.43
12 × 12
15.86
7.62
6.02
5.61
5.45
5.33
14 × 14
16.94
8.65
7.02
6.56
6.37
6.21
16 × 16
18.02
9.68
8.02
7.51
7.27
7.07
(b)
VFM-Wu
Rns = 2 kΩSize
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
4 × 4
10.77
2.76
1.25
0.95
0.88
0.84
6 × 6
11.27
3.26
1.75
1.43
1.35
1.30
8 × 8
11.77
3.75
2.23
1.89
1.79
1.73
10 × 10
12.27
4.25
2.71
2.35
2.22
2.13
12 × 12
12.77
4.74
3.19
2.79
2.63
2.52
14 × 14
13.27
5.24
3.66
3.22
3.03
2.88
16 × 16
13.77
5.73
4.13
3.63
3.41
3.21
(c)
VFM proposed method, 1st step
Rns = 2 kΩSize
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
4 × 4
0.61
0.08
0.05
0.04
0.05
0.05
6 × 6
0.73
0.14
0.09
0.09
0.09
0.09
8 × 8
0.87
0.22
0.15
0.15
0.15
0.15
10 × 10
1.02
0.31
0.24
0.22
0.22
0.22
12 × 12
1.18
0.42
0.33
0.31
0.30
0.30
14 × 14
1.15
0.55
0.44
0.41
0.39
0.38
16 × 16
0.07
0.68
0.56
0.52
0.50
0.48
(d)
VFM proposed method, 2nd step
Rns = 2 kΩSize
Error (%),Rij (Ω)
100
500
2000
5000
7500
10000
4 × 4
0.06
0.01
0.01
0.02
0.02
0.02
6 × 6
0.07
0.02
0.01
0.02
0.03
0.03
8 × 8
0.09
0.02
0.02
0.03
0.03
0.04
10 × 10
0.11
0.03
0.03
0.03
0.04
0.04
12 × 12
0.13
0.04
0.04
0.04
0.05
0.05
14 × 14
0.08
0.06
0.05
0.05
0.06
0.06
16 × 16
1.56
0.07
0.06
0.06
0.07
0.07


## 5. Conclusions

The main problem in determining sensor resistances in a resistive array is the appearance of crosstalk. This phenomenon is due mainly to the resistance of the sensor selection switches to be measured. Crosstalk appears in the two types of circuits known in the literature for reading the array: the zero potential method (ZPM) and the voltage feedback method (VFM). This article presents a recursive method to obtain the individual values of the resistors in an array of resistive sensors that reduces the influence of crosstalk in determining the resistances of the sensors. The new method proposed is similar for both types of circuits and is based on the recursive calculation of the row and column voltages of the array.The proposed method converges very quickly, outperforming, in its first or second iteration, the classical equations of the ZPM and VFM and also the circuits proposed by Wu to reduce errors. The method does not require any additional hardware and is applied directly to the classic ZPM and VFM circuits. The errors found in the calculations using the proposed method are similar for both types of circuit and increase slowly as the array size and selection switch resistance increase.As a further conclusion, it should be noted that a decrease has been observed in the range of resistance values that can be measured in VFM circuits compared to ZPM circuits for equal output voltage ranges.



---

*Source: 9735741-2018-07-16.xml*