CN103776873B - A kind of method mapped by voltage-to-current map construction current-voltage - Google Patents
A kind of method mapped by voltage-to-current map construction current-voltage Download PDFInfo
- Publication number
- CN103776873B CN103776873B CN201410054956.XA CN201410054956A CN103776873B CN 103776873 B CN103776873 B CN 103776873B CN 201410054956 A CN201410054956 A CN 201410054956A CN 103776873 B CN103776873 B CN 103776873B
- Authority
- CN
- China
- Prior art keywords
- voltage
- lambda
- current
- matrix
- epsiv
- 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.)
- Active
Links
Landscapes
- Complex Calculations (AREA)
Abstract
The present invention relates to a kind of building method mapped by voltage-to-current map construction current-voltage, for the sensor in electricity tomography with N number of electrode, by the parameter such as eigenwert and proper vector of calculating voltage-current matrix, by derived formula, unique calculates current-voltage matrix, constructs current-voltage and maps.The present invention gives the direct building method that current-voltage maps, can be applicable in the direct reconstruction algorithm of electricity chromatography imaging field, its explicit physical meaning, and simple.
Description
Technical field
The present invention relates to electricity chromatography imaging field, particularly relate to a kind of method mapped by voltage-to-current map construction current-voltage.
Background technology
Electricity tomography (ElectricalTomography is called for short ET) technology is the one of chromatography imaging technique.By applying excitation to testee, and detecting the change of its boundary value, utilizing specific reconstruction algorithm to reconstruct the distribution of measurand internal electrical characterisitic parameter, thus obtaining the distribution situation of interior of articles.Compared with other chromatography imaging techniques, electricity tomography has radiationless, Noninvasive, portability, fast response time, the advantage such as cheap.
Electricity tomographic reconstruction algorithm generally can be divided into two classes: based on reconstruction algorithm and the direct reconstruction algorithm of sensitivity matrix.Use the former usually to need to separate ill linear equation, this just means all pixel values simultaneously must rebuilding measured zone.The latter is mapped by calculating current-voltage or voltage-to-current mapping realizes, and the gray-scale value of each pixel can through directly, independently calculating acquisition.
The building method that current-voltage maps is the important component part of the direct reconstruction algorithm of electricity tomography.For the sensor having N number of electrode, Law of Inner Product can be applied and map to construct current-voltage.But also do not map according to voltage-to-current the direct building method constructing current-voltage and map at present.
Summary of the invention
The object of the invention is to propose a kind of direct building method constructing current-voltage mapping according to voltage-to-current mapping.
Technical scheme of the present invention is:
Step one, have N number of electrode sensor electrode on current density have general expression
write as matrix form to have:
Wherein J is current density vector, and A is the internal surface area of each electrode, and C is capacitance matrix, and V is voltage vector, and expression formula is:
J
sthe current density on s (1≤s≤N) individual electrode inside surface, V
tthe magnitude of voltage on t electrode, C
s,tit is the capacitance between electrode s and electrode t.Self-capacitance C
s,sbe defined as the electric capacity summation between s electrode and other N-1 electrode, namely
So have
being number of poles is N, and specific inductive capacity is distributed as voltage-to-current mapping matrix during ε (z).
According to Kirchhoff's law, loop current sum is zero, and namely the order of current density vector J is N-1, and the order of Matrix C is also N-1, can show that the order of voltage vector is also N-1 thus, therefore make
So
Wherein
being number of poles is N, and specific inductive capacity is distributed as current-voltage mapping matrix during ε (z).
Can prove that front N-1 the eigenwert of V is-1, and last eigenwert is 0.By calculating eigenwert and the proper vector of V, can obtain:
V=P{diag([-1-1…-10])
N×N}P
T(7)
Wherein diag ()
n × Nrepresent the diagonal matrix on N rank, P
tp=I, P=[p
1p
2p
n-1p
n].Defined feature vector p
ii-th (1≤i≤N-1) row of P, i.e. Vp
i=(-1) p
i, p
neigenwert 0 characteristic of correspondence vector, i.e. Vp
n=(0) p
n.Because the mean value of often going in equation (5) is all 0, therefore
Wherein I
n × Nit is the unit matrix on N rank.
Similar, by calculating eigenwert and proper vector,
order be N-1, can be written as following form:
Wherein Σ
Λ=diag ([λ
1λ
2λ
n-10]
n × N) diagonal matrix that is made up of N number of eigenwert of R, the matrix that Q is made up of corresponding proper vector.Q
tq=I, Q=[q
1q
2q
n-1q
n], proper vector q
ii-th (1≤i≤N-1) row of Q, i.e. Rq
i=λ
iq
i, q
nproper vector 0 characteristic of correspondence vector, i.e. Rq
n=(0) q
n,
Therefore have
Step 2, by formula (7) substitute into formula (6) have:
(11)
Both members is premultiplication P simultaneously
t, the right side takes advantage of P to obtain:
(12)
The matrix that a Matrix Multiplication obtains with diagonal matrix if easily know is still for diagonal matrix, then this matrix is also diagonal matrix.So known:
Again at formula (13) premultiplication P, P is taken advantage of on the right side
t:
Formula (7) and formula (9) are substituted into formula (14) obtain
At formula (15) premultiplication Q
t, the right side takes advantage of Q to obtain:
Formula (10) is substituted into formula (16) and get final product:
So have:
At formula (18) premultiplication Q, Q is taken advantage of on the right side
t:
Formula is known thus, matrix
can be now uniquely determined.
Further, with reference to this method, any orthogonal set of current excitation pattern may be used to calculating current-voltage and maps.
Accompanying drawing explanation
Fig. 1 is implementing procedure figure.
Fig. 2 is embodiment isoboles.
Embodiment
See Fig. 1, a kind of building method algorithm block diagram constructing current-voltage mapping according to voltage-to-current mapping.For 16 end to end ring resistance networks shown in Fig. 2, the embodiment of this method is described.
Said method comprising the steps of:
Step one, have N=16 electrode sensor electrode on current density have general expression
write as matrix form to have:
Wherein J is current density vector, and A is the internal surface area of each electrode, and C is capacitance matrix, and V is voltage vector, and expression formula is:
J
sthe current density on the individual electrode inside surface of s (1≤s≤16), V
tthe magnitude of voltage on t electrode, C
s,tit is the capacitance between electrode s and electrode t.Self-capacitance C
s,sbe defined as the electric capacity summation between s electrode and other 15 electrodes, namely
So have
According to Kirchhoff's law, loop current sum is zero, and namely the order of current density vector J is 15, and the order of Matrix C is also 15, can show that the order of voltage vector is also 15 thus, therefore make
So
Can prove that front 15 eigenwerts of V are-1, and last eigenwert is 0.By calculating eigenwert and the proper vector of V, can obtain:
V=P{diag([-1-1…-10])
16×16}P
T(26)
Wherein diag ()
16 × 16represent the diagonal matrix on 16 rank, P
tp=I, P=[p
1p
2p
15p
16].Defined feature vector p
ii-th (1≤i≤15) row of P, i.e. Vp
i=(-1) p
i, p
16eigenwert 0 characteristic of correspondence vector, i.e. Vp
16=(0) p
16.Because the mean value of a line every in equation (24) is all 0, therefore
So V can be written as:
Wherein I
16 × 16it is the unit matrix on 16 rank.
Similar, by calculating eigenwert and proper vector,
order be 15, can be written as following form:
Wherein Σ
Λ=diag ([λ
1λ
2λ
150]
16 × 16) diagonal matrix that is made up of 16 eigenwerts of R, the matrix that Q is made up of corresponding proper vector.Q
tq=I, Q=[q
1q
2q
15q
16], proper vector q
ii-th (1≤i≤15) row of Q, i.e. Rq
i=λ
iq
i, q
16proper vector 0 characteristic of correspondence vector, i.e. Rq
16=(0) q
16,
From formula (24), in ring resistance network, can be calculated current matrix:
It can thus be appreciated that voltage-to-current mapping matrix is
Calculate according to formula (28)
(31)
Step 2, by formula (26) substitute into formula (25) have:
Obtain
Formula is known thus, matrix
can be now uniquely determined.
For proving this conclusion, doing following checking: under adjacent actuators pattern, calculating node potential matrix
With setting V be consistent, can this method be verified thus.
Described a kind ofly construct by voltage-to-current the method that current-voltage maps, give a kind of computing method that current-voltage maps, the method explicit physical meaning, simple.With reference to this method, any orthogonal set of excitation measurement pattern may be used to calculating current-voltage and maps.
Above to the description of the present invention and embodiment thereof, being not limited thereto, is only one of embodiments of the present invention shown in accompanying drawing.When not departing from the invention aim, designing the structure similar with this technical scheme or embodiment without creation, all belonging to scope.
Claims (1)
1. map the direct building method constructing current-voltage and map according to voltage-to-current, it is characterized in that, the method comprises the steps:
Step one, the current density had on each electrode of sensor of N number of electrode have general expression
wherein A is the internal surface area of each electrode, and C is capacitance matrix, and V is voltage vector, so voltage-to-current mapping matrix when having specific inductive capacity to be distributed as ε (z)
According to Kirchhoff's law, loop current sum is zero, and namely the order of current density vector J is N-1, and the order of Matrix C is also N-1, can show that the order of voltage vector is also N-1 thus, therefore make
So
Calculate eigenwert and the proper vector of V, can obtain:
V=P{diag([-1-1…-10])
N×N}P
T(4)
Wherein diag ()
n × Nrepresent the diagonal matrix on N rank, the matrix that P is made up of corresponding proper vector, P
tp=I, P=[p
1p
2p
n-1p
n], proper vector p
ii-th row of P, wherein 1≤i≤N-1, i.e. Vp
i=(-1) p
i, p
neigenwert 0 characteristic of correspondence vector, i.e. Vp
n=(0) p
n,
Namely V can be written as:
Wherein I
n × Nthe unit matrix on N rank,
Similar, by calculating eigenwert and proper vector,
order be N-1, following form can be written as:
Wherein Σ
Λ=diag ([λ
1λ
2λ
n-10]
n × N) diagonal matrix that is made up of N number of eigenwert of R, the matrix that Q is made up of corresponding proper vector, Q
tq=I, Q=[q
1q
2q
n-1q
n], proper vector q
ii-th row of Q, wherein 1≤i≤N-1, i.e. Rq
i=λ
iq
i, q
nproper vector 0 characteristic of correspondence vector, i.e. Rq
n=(0) q
n,
So have
Step 2, by formula (4) substitute into formula (3) have:
Then through deriving, the matrix that current-voltage mapping pair is answered
can be written as:
Formula is known thus, matrix
can be now uniquely determined.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410054956.XA CN103776873B (en) | 2014-02-18 | 2014-02-18 | A kind of method mapped by voltage-to-current map construction current-voltage |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410054956.XA CN103776873B (en) | 2014-02-18 | 2014-02-18 | A kind of method mapped by voltage-to-current map construction current-voltage |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103776873A CN103776873A (en) | 2014-05-07 |
CN103776873B true CN103776873B (en) | 2016-04-06 |
Family
ID=50569372
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410054956.XA Active CN103776873B (en) | 2014-02-18 | 2014-02-18 | A kind of method mapped by voltage-to-current map construction current-voltage |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103776873B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107091858A (en) * | 2017-03-27 | 2017-08-25 | 北京航空航天大学 | A kind of method mapped by Current Voltage map construction voltage x current |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1274044A1 (en) * | 1997-04-17 | 2003-01-08 | Ge Medical Systems Israel, Ltd. | Direct tomographic reconstruction |
WO2006102388A1 (en) * | 2005-03-22 | 2006-09-28 | The Ohio State University | 3d and real-time electrical capacitance volume-tomography: sensor design and image reconstruction |
CN101520478B (en) * | 2009-03-13 | 2011-08-17 | 北京航空航天大学 | Direct image reconstruction method based on capacitance tomography of round sensor |
CN103163404B (en) * | 2013-02-01 | 2015-07-15 | 北京航空航天大学 | Current-voltage mapping construction method based on adjacent stimulus measurement model |
CN103412189B (en) * | 2013-07-30 | 2015-07-01 | 北京航空航天大学 | Information filtering demodulation method for electrical tomography system |
-
2014
- 2014-02-18 CN CN201410054956.XA patent/CN103776873B/en active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107091858A (en) * | 2017-03-27 | 2017-08-25 | 北京航空航天大学 | A kind of method mapped by Current Voltage map construction voltage x current |
Also Published As
Publication number | Publication date |
---|---|
CN103776873A (en) | 2014-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schmidt et al. | The distribution of relaxation times as basis for generalized time-domain models for Li-ion batteries | |
CN106530367B (en) | A kind of electricity tomography sparse reconstruction method based on Firm threshold value iteration | |
Ren et al. | Inclusion boundary reconstruction and sensitivity analysis in electrical impedance tomography | |
CN111281385B (en) | Electrical impedance imaging method based on tissue space distribution characteristics and impedance variation characteristics along with frequency | |
CN101520478A (en) | Direct image reconstruction method based on capacitance tomography of round sensor | |
CN107369187A (en) | The electricity tomography regularization reconstruction method for the sum that is deteriorated based on adjoint point | |
Bera et al. | A MATLAB‐Based Boundary Data Simulator for Studying the Resistivity Reconstruction Using Neighbouring Current Pattern | |
Wang et al. | Reconstruction of electrical impedance tomography (EIT) images based on the expectation maximum (EM) method | |
de Castro Martins et al. | Electrical impedance tomography reconstruction through simulated annealing with incomplete evaluation of the objective function | |
CN108376124B (en) | Multi-conductor system admittance matrix fast calculation method for electrical imaging | |
CN104965134B (en) | Grounding net of transformer substation surface potential distribution calculation method based on echo state network | |
CN104535294B (en) | Corrected L-curve electrical tomography reconstruction method based on second-order differential | |
CN103163404B (en) | Current-voltage mapping construction method based on adjacent stimulus measurement model | |
CN103776873B (en) | A kind of method mapped by voltage-to-current map construction current-voltage | |
CN103149472B (en) | Voltage-current mapping constructing method based on double-terminal excitation measuring model | |
CN105701824B (en) | Electricity based on arc-shaped mesh generation chromatographs image rebuilding method | |
CN103116101B (en) | Adjacent excitation measurement mode based voltage-current mapping construction method | |
CN114155354B (en) | Method and device for reconstructing capacitance tomography based on graph convolution network | |
CN107845119A (en) | A kind of electricity tomography mixed method | |
CN103054577A (en) | Sparse reconstruction method for electrical impedance tomography | |
CN110470743B (en) | Electrical/ultrasonic information fusion bimodal tomography method | |
Ding et al. | Fuzzy optimal solution of electric tomography imaging: Modelling and application | |
Yan et al. | A Detection Method for Fast Electrical Impedance Imaging of Grounding Grid Based on Optimized Differential-Multigrid-Homotopy Algorithm | |
Kim et al. | Image reconstruction using genetic algorithm in electrical impedance tomography | |
Lim et al. | An efficient forward solver in electrical impedance tomography by spectral element method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |