Open access peer-reviewed chapter

Fault Diagnosis in Analog Circuits via Symbolic Analysis Techniques

By Fawzi M Al-Naima and Bessam Z Al-Jewad

Submitted: April 12th 2012Reviewed: September 20th 2012Published: January 9th 2013

DOI: 10.5772/53643

Downloaded: 4565

1. Introduction

Fault diagnosis of analog circuits has been one of the most challenging topics for researchers and test engineers since the 1970s. Given the circuit topology and nominal circuit parameter values, fault diagnosis is to obtain the exact information about the faulty circuit based on the analysis of the limited measured circuit responses. Fault diagnosis of analog circuits is essential for analog and mixed-signal systems testing and maintenance both during the design process and the manufacturing process of VLSI ASICs.

There are three dominant and distinct stages in the process of fault diagnosis: fault detection to find out if the circuit under test (CUT) is faulty comparing with the fault-free circuit, or golden circuit (This stage is usually called test in industry), fault identification to locate where the faulty parameters are inside the faulty circuit, and parameter evaluation to obtain how much the faulty parameters deviated from their nominal values and to obtain values of other circuit parameters such as branch and nodal voltages. The bottlenecks of analog fault diagnosis primarily lie in the inherited features of analog circuits: nonlinearity, parameter tolerances, limited accessible nodes, and lack of efficient models. Multiple fault diagnosis techniques are even less developed than single fault diagnosis because it is more difficult to model and detect multiple faults, particularly in the presence of tolerance or measurement noise. In addition, in multiple fault situation, one fault’s effect on the circuit could be masked by the effects of the other faults. Generally speaking, there is no widely accepted paradigm for analog test or fault diagnosis even with the introduction of IEEE 1149.4 standard for mixed-signal test bus.

With recent sharp development of electronic design automation tools and widespread application of analog VLSI chips and mixed-signal systems in the area of wireless communication, networking, neural network and real-time control, the interests in analog test and fault diagnosis revives. System-on-chip solutions favored by modern electronics pose new challenges in this topic such as increased complexity and reduced die size and accessibility.

Several methods have been proposed for single fault diagnosis in linear analog circuit in the past. Multiple excitations are required and Woodbury formula in matrix theory is applied to locate the faulty parameters. This method is also applied to multiple fault diagnosis by decomposition technique assuming that each sub-circuit contains at most a single faulty parameter.

Among the different methods of fault diagnosis, the parametric fault diagnosis techniques hold an important part in the field of analog fault diagnosis. These techniques, starting from a series of measurements carried out on a previously selected test point set, given the circuit topology and the nominal values of the components, are aimed at determining the effective values of the circuit parameters by solving a set of equations generally nonlinear with respect to the component values. In this chapter the role of symbolic techniques in the automation of parametric fault diagnosis of analog circuits is investigated followed by a practical numerical procedure to evaluate the faults. Being in fact the actual component values that represent the unknown quantities, fault diagnosis aims also at finding the faults locations. Symbolic approach results are particularly suitable for the automation of parametric fault diagnosis techniques [1]. Obviously all this is applicable to linear analog circuits or to nonlinear circuits suitably linearized. On the other hand, present trend is moving as much as possible to design techniques that lead to linear analog circuits, so linearity is not a so serious restriction any more [2]. It is important to note that in the analog fault diagnosis two phases can be considered: the first one is the phase of testability analysis and ambiguity group determination, while the second one is the phase of fault location and fault value determination. Testability gives theoretical and rigorous upper limits to the degree of solvability of fault diagnosis problem once the test point set has been chosen, independently of the method effectively used in fault location phase. This becomes very important in the design stage of the linear circuit in which the designer can determine the list of accessible nodes for the operator and the fault detection ability that they can provide. Concerning ambiguity groups, they can be considered as sets of circuit components that, if considered as potentially faulty, yield an undetermined system of equations. For the testability evaluation problem symbolic approach is a natural choice, because a circuit description made by means of equations in which the component values are the unknowns is properly represented by symbolic relations. Also for ambiguity group determination the symbolic approach gives excellent results [3].

For the fault location phase several different approaches can be used and all of them can be considered as an optimization problem, because, starting from measurements carried out on the CUT, the component values better fitting them have to be determined. Generally symbolic techniques are suitable for optimization problems, because the relations required by specific optimization strategies are easily generated using symbolic methods.

The aim of this chapter is to present a unified treatment of the subject with emphasis on a generalized method for multiple fault diagnosis of linear analog circuits in frequency domain. In this approach, multiple excitations and Woodbury formula are also required for fault identification. However, a recently developed ambiguity group locating technique is applied for fault identification which reduces computational cost of the test method. Multiple faults can be located directly and efficiently, thus eliminating the requirement for decomposition and the corresponding restrictions. Moreover, the methodology developed in the proposed method, (i.e., constructing fault diagnosis equation on the basis of the analysis of the fault-free circuit and the measured responses of faulty circuit, then applying the ambiguity group locating technique to identify the faulty parameters, finally evaluating all parameter values of faulty circuit exactly), can be applied to other methods developed for multiple analog fault diagnosis.

The dominant differences among these methods are the distinct fault diagnosis equations resulting from distinct circuit analysis methods and distinct excitation and measurement methods. The methods proposed in this chapter can be classified as fault verification methods under the category of Simulation-after-Test (SAT), which can provide the exact solution to the circuit parameters and can be applied to detect large parameter changes when the number of independent measurements are greater than the number of faults in the CUT.

A major improvement of these techniques is achieved through the use of symbolic techniques in formulating the fault equations and in addressing the testability problem. Furthermore a developed method for minimum size ambiguity group locating technique is used based on QR factorization and is applied to detect and identify the multiple faults. Detailed procedures for a proposed fault diagnosis program are given to help practitioners and researchers as well to grasp the basic concepts of the topic and be able to contribute to this field.


2. Basic circuit formulations

Generally, the circuit topology as well as its parameters nominal values are known in advance. Consider for example a continuous-time, time-invariant, strongly connected, linear circuit with n+1 nodes and pparameters. The (n+1)stnode, denoted by zero, is assigned to be the grounded reference node while the remaining nnodes are ungrounded. All pparameters are divided into two categories: one is parameters which have admittance description such as conductance, capacitor and voltage- controlled-current source, another is parameters which have no admittance description such as impedance, inductor, current-controlled-source, operational amplifier, etc.

Of course the conventional method of analysis would be to apply the KCL to each circuit node to obtain nequations with variables being nodal voltages and parameter currents. Then constitutive equations in terms of nodal voltages and parameter currents, which define the characteristics of all parameters without admittance description, are appended to the above nKCL-based equations. The resulting system matrix from this approach would be


where Tgis a g×gcoefficient matrix consisting of circuit parameters, Xgis a g×1 solution vector of node voltage and parameter currents, and Wgis a g×1 excitation vector composed of independent current and voltage sources, and initial conditions of capacitors and inductors. The first nrows in Tg, Xgand Wgcorrespond to nnodes. The resulting system equation (1) is called the modified nodal analysis equation MNA[4]. Note that g=nfor normal nodal analysis of a circuit in which all parameters have admittance description, and g>nfor modified nodal analysis of a circuit in which some parameters have a non-admittance description.

Traditionally, the system matrix generated from the MNA method may still have many redundant variables for analysis purposes; for example: voltages of inaccessible nodes inside sub-circuits like op-amps or currents through nonphysical branches generated from sophisticated element models. A major development step to the MNA method is to eliminate all redundant variables to generate a compacted or reduced system matrix. The reduced system matrix is formulated by programming a lookup table for every element in the network. This table has conditioned link-lists that will test which variables of the element are actually needed in the final compacted matrix and introduce the element in a way so as to eliminate the redundant variables during the formulation process. This method is termed the compacted modified nodal analysis CMNA[5]. Provided that the circuit functions in a stable state, the parametric values of nodal voltages and parameter currents will be finite and unique. The coefficient matrix Tgis non-singular since the circuit is a strongly connected network.

Generally the system matrix described above cannot be formulated smoothly in a computerized solution without taking the circuit topology into consideration. One important fact about circuit topology is that each parameter, say hv( v= 1, 2,..., p), can be located by at most 4 circuit nodes [6]: 2 input nodes kvand lv,and 2 output nodes ivand jv. For 2-terminal parameters such as resistors and capacitors, the input nodes will be the same as the output nodes: kv= iνand lv= jν. Based on this fact, the circuit topology can be completely described by two g×pstructural matrices Pand Qwhich are defined as follows:

P=p1 p2  pp=[δi1-δj1 δi2-δj2  δip-δjp]Q=q1 q2  qp=[δk1-δl1 δk2-δl2  δkp-δlp]E2

where δνrepresents a g×1 vector of zeros except for the ventry, which is equal to one, and pνand qνrepresent g×1 vectors describing the locations of output nodes and input nodes, respectively. Matrices Pand Qare only determined by the locations, not the values of the circuit parameters. The columns of matrix Pcorrespond to the locations of the output nodes of circuit parameters while the columns of matrix Qcorrespond to the locations of the input nodes of circuit parameters.

Another important fact is that most parameters in linear circuits will enter the coefficient matrix Tgin the symbolic form

       kv     lvivjvhv-hv-hvhvE3

with the equivalent algebraic representation being


where superscript tdenotes transpose of matrix or vector. For any grounded node, the corresponding row or column in the symbolic form will be removed together with the δνin the algebraic form. Resistor, inductor, capacitor, dependent sources, and operational amplifier with its negative inverse gain being a parameter are examples of circuit devices described in this way. Thus the system matrix can be easily formulated using the equation

Tg=PdiaghQt E5

This topological formulation allows for the automatic direct translation of the Netlist (which is the list describing the values of the circuit elements and their connections to the corresponding nodes) into circuit equations. As an example consider the circuit shown in Figure 1 following [7]. This circuit will be used later in the analysis of fault equations.

Figure 1.

Example circuit

This circuit can be represented by the Netlist shown in Table 1. The first column of Table 1 is the element name and type (e.g. R for resistor, G for conductance, C for capacitor and so on). The second and third (possibly fourth and fifth depending on the element type) represent the connection nodes. Note that the ground is identified with node 0. The last column is the element value. It is possible for some elements (like a two port network or an active element model) to have more columns for its values. For the circuit in Figure 1 with unity current source J=1A applied to nodes {0,1}, it is assumed that the nominal values of the resistors are as follows (all resistors in Ω): R1=2.125, R2=3.6, R3=4.7, R4=11.5, R5=12.6, R6=21.2, R7=3.7, R8=0.54, R9=3.54, R10=3.125, R11=6.6, R12=5.7, R13=19.5, R14=12.8, R15=12.2, R16=3.2, R17=1.54, R18=8.7, R19=2.27, R20=3.16, R21=41.7, R22=31.5, R23=22.6, R24=51.2, R25=13.7, R26=3.44, R27=13.4, R28=31.9, R29=16.1, R30=11.7, R31=11.5, R32=17.8, R33=22.2, R34=23.2, R35=11.4, R36=18.7, R37=3.12, R38=33.2, R39=8.67.

ElementFrom NodeTo NodeValue

Table 1.

Netlist for the circuit shown in Figure 1

It can be shown easily [4] that the node numbers in the second column of Table 1 represent the values of iin equation (2) as long as they are not zero. Similarly column 3 gives the values of jin equation (2). Thus it is very straight forward to find the Pand Qmatrices by tracing through the Netlist. Additionally, the system matrix can be either formulated numerically or symbolically depending on the values used for hin the diagonal matrix of equation (5). This automatic formulation procedure can be easily extended to MNA and CMNA methods with proper rubber stamps for circuit elements.

3. Conventional fault analysis methods

As discussed earlier, the conventional method for multiple fault diagnosis can be divided into three steps: fault detection, fault location determination, and finding the faulty elements values. This conventional method is readily deemed to be a numerical method by its very own nature but it is presented here as it provides basic insight to the problem and the limitations facing all numerical methods. The problem is even more complicated for multiple faults due to ambiguity presented by element tolerances not to mention that different sets of certain faults may produce very similar measured values. Further complication is present owing to the fact that only a limited number of nodes are actually accessible for measurements and testing. The conventional method will be presented without derivation which can be found in [6-8]. Despite its effectiveness in dealing with ambiguity groups, the method has several limitations:

  1. The method requires multiple independent excitations among the accessible nodes. That is applying an independent source of excitation to a subset of the accessible nodes and measuring the circuit response for each source. By this, the method not only assumes that the circuit will remain linear and well-behaved under multiple excitations, it also destroyes the natural input-output relation of the circuit components and overlooks any form of signal isolation.

  2. The method needs a dictionary for the behavior of the fault-free circuit under multiple excitations. The dictionary must be extensive enough to enable detecting and locating a number of simultaneous faults. Yet even when the dictionary is extensive enoguh the method may still fail in differentiating certain ambiguity groups subject to rounding errors, inaccuracies and noise that may occur in the measurements.

  3. Depending on the set of accessible nodes the problem of testability and detectability of multiple faults immediately arises.

Symbolic analysis techniques aim at resolving or at least reducing some or all of these limitations thus proving to be very vital to this subject. Having that said let us begin by showing the applications of symbolic analysis techniques for multiple fault diagnosis in linear circuits.

3.1. Symbolic analysis in fault diagnosis problem

As mentioned earlier, the symbolic circuit matrix Tcan be easily described by the multiplication of a row operator Pand a column operator Qwith a diagonal matrix of symbols, where Pand Qare the matrix operators (or topological matrices) that indicate location of the matrix element value in the symbolic matrix Tas was given in equation (5). One can immediately notice that the diagonal matrix can be represented by a simple array while both Pand Qare sparse numerical matrices. One benefit of this representation is that linear operations on rows such as addition or subtraction can be simply implemented on Pwhile linear operations on columns can be implemented on Qwithout altering the diagonal symbolic matrix. For most circuit analysis applications, system equations are rarely fully dense or fully symbolic. In that respect, some matrix elements may contain not only a single symbolic element value but constant values as well. These constants do not affect the structures of matrices Pand Q. In fact even if the element values are complete polynomials the representation is still not altered. Hence in solving a matrix equation like equation (1) only arrays of symbolic polynomial data structures need to be stored to represent the diagonal matrix while numerical matrices P, Qcan be manipulated to solve this system using topological methods like determinant decision diagrams DDD [9] or matrix reduction methods [10].

Traditionally symbolic techniques have been used along two separate paths in fault analysis: To introduce comprehensive fault models in small and moderate sized circuits and to find the optimum set of testable components in a faulty circuit. We will describe here both techniques and show ways of utilizing them later.

3.1.1. Symbolic techniques in comprehensive fault modeling

The method described before for tracing through the Netlist to generate the Pand Qtopological matrices is not the only way to generate the system matrix. In fact, most of the modified and compacted methods like the tableau, MNA and CMNA methods focus on tracing through the Netlist in an element-by-element fashion, increasing the size of the generated system matrix iteratively. This approach has the advantage of providing a way to reduce the system matrix during the formulation step as will be shown shortly [5].

Consider a general admittance yconnected between nodes iand jas shown in Figure 2. Assuming that the system matrix Tis already generated for the other branches, the impact of this admittance (following the tableau formulation) on the system matrix is to add an additional row and column corresponding to the new system variable iy


Now, if iyis not a solution variable then it can be eliminated from the system matrix to generate a compacted matrix with respect to the axis iy. Applying Kron’s reduction to eliminate axis-3 of this matrix we get [5]

yii+y    yij-yyji-y    yjj+y×ViVj=wiwj.E7

Equations (6) and (7) are the conditioned stamps for the admittance yand they can be programmed into a lookup table easily.

The occurrence of faults in circuit elements generally leads to a deviation in node voltages and branch currents from nominal values. The purpose of fault diagnosis is to use voltage measurements on a limited number of nodes to verify the presence of a fault in the circuit then identify the fault location and value through simulation. However additional care must be taken since the deviations in the measured values may very well result from normal parameter tolerances.

In general a fault is generated if the nominal element value is changed from yto y+dybeyond its tolerance. However instead of changing the value of yin the system matrix, owing to the symbolic approach we can simply introduce an extra faulty element dy and an extra faulty variable fas shown in Figure 3 for the passive admittance case. This deviation from nominal value will result in deviation in node voltages from Vto V+dVand branch currents from Ito I+Ifwhere Ifis the set of fault currents. As an example consider a linear admittance yconnected between nodes iand jas shown in Figure 3.

Figure 2.

A general admittance example

Figure 3.

Fault element model for a linear admittance

Due to a fault the value of the admittance is changed to y+dy. Now following [5] the fault rubber stamp of (7) is changed to


where we have appended the faulty element equation


in which ζ = (dy)-1. We must emphasize here that this equation is appended to the original system matrix after the last nodal equation so that the faulty ifvariable will appear after the last solution variable in the fault-free system equation. The impact of this on the solution will be apparent shortly. Assuming that ifis a solution variable, we can proceed by eliminating any non-solution variables in the stamp just like we did before to generate the faulty compacted system matrix stamp while leaving all the faulty currents as solution variables.

Clearly other fault models for the different circuit elements can also be developed by inspection. The stamping procedure for the faulty circuit elements generates the fault analysis equations. However, additional circuit elements and thus additional symbolic variables were introduced in the circuit to simulate the fault thus deeming this method suitable only for small and moderate sized networks. Assuming the original fault-free system equations were given by equation identical to (1) where Tis the compacted modified system matrix, Xis the chosen solution vector and Wis the vector of excitation sources which might be a combination of currents and voltages. With the introduction of the faulty elements the size of the system matrix has increased. Careful consideration of element stamps like (8) show that the new faulty system equations can be written as


where Xis the original solution vector while Xfis the solution vector of the fault currents and voltages. This formulation resulted from the fact that the fault variable was added after the last solution variable of the fault-free circuit. Expanding this equation, it can be shown that [5]


where Xis now the solution vector of the faulty system. Applying Woodbury formula [5] on (11) we get


Expanding (12) using (10) we get




This gives the variation in the solution vector in terms of the nominal fault-free solution vector, the topological matrices and the original system matrix. The benefit of having this variation solved symbolically is that it gives direct relationship between shifts in element values and the corresponding variation in circuit response. Once those variables are obtained symbolically it is very easy to carry out an analysis like Monte-Carlo analysis [5] to help solve the fault/tolerance ambiguity and verify the presence of a fault. Not only fault verification is possible with this equation but also locating the faulty element(s) can be done even with measurements taken from a limited set of accessible nodes using the k-fault method [11] or a linear combination matrix which will be explained later on where only a small set of the solution variables are measured to estimate the fault location.

Despite the usefulness of this approach in finding the fault model, it is highly restricted to small and moderate scale circuits. In addition, due to ambiguities, it is customary to find that the inner matrix R+QT-1Phas become singular therefore limiting the practical use of equation (14). Nevertheless the approach is still needed to model the faults symbolically and to tackle the testability problem.

3.1.2. Symbolic solution of the testability problem

Fault diagnosis and fault location in analog circuits are of fundamental importance for design validation and prototype characterization in order to improve yield through design modification. In the analog fault diagnosis field, an essential point is constituted by the concept of testability which, independently of the method that will be effectively used in fault location, gives theoretical and rigorous upper limits to the degree of solvability of the problem, once the test point set has been chosen by the circuit designer. A well-defined quantitative measure of testability can be deduced by referring to fault diagnosis techniques of the parametric kind [12]. These techniques, starting from a series of measurements carried out on previously selected test points, are aimed at determining the upper limit of solvable circuit parameters by solving a set of equations (the fault diagnosis equations as will be shown later) which are nonlinear with respect to the component values.

The solvability degree of these nonlinear equations constitutes the most used definition of testability measure [12]. This measure can be also interpreted as an indication of the ambiguity resulting from any attempt to solve the fault equations in a neighborhood of almost any failure. In addition to being valuable for the circuit designer in determining the number of accessible nodes, it is also very important for the circuit operator since attempting to address the fault-diagnosis equation without having an estimate on the maximum number of faults that can be detected from the available test set is highly prohibited. In other words, the testability measure provides information about the number of testable components with the selected test point set. When the testability value is not at its maximum, that is when it is less than the total number of potentially faulty circuit components, the problem is not uniquely solvable and it is necessary to consider further measurements, i.e., other test points. Alternatively we can accept a reduced number of potentially faulty components in order to locate the elements which have caused the incorrect behavior of the CUT.

Generally, the second alternative is used for two reasons. First, not all the possible test points can actually be considered because of practical and economic measurement problems strictly tied with the used technology and with the application field of the circuit under consideration. Second, the number of faulty components is generally smaller than the total number of circuit components. The single fault case is the most frequent while double or triple cases are less frequent, and the case of all faulty components is almost impossible. Therefore, as the testability is normally not at its maximum, the fault diagnosis problem is dealt with by assuming the quite realistic hypothesis that the number of faulty components is bounded; that is, the k-fault hypothesis is made. Under this hypothesis, in order to locate the faulty elements with as low ambiguity as possible, it is of fundamental importance to determine a set of components that is representative of all the circuit elements. This helps reducing the solution time by providing a stopping criterion instead of wasting computer resources seeking unattainable solutions. In this section a procedure for the determination of the optimum set of testable components in the k-fault diagnosis of analog linear circuits is presented, where by the optimum set we mean a set of components representing all the circuit elements and giving a unique solution.

The procedure is based on the testability evaluation of the circuit and on the determination of the canonical ambiguity groups. Referring again, for the sake of simplicity, to parametric fault diagnosis techniques we need to make some definitions first. An ambiguity group can be defined as a set of components that, if used as unknowns (i.e., if considered as potentially faulty), gives infinite solutions during the phase of fault location determination. A canonical ambiguity group is simply an ambiguity group that does not contain other ambiguity groups. It is worth pointing out that the proposed procedure gives information independently of the method that will be effectively used in the fault location phase (both simulation after test and simulation before test methods), even if it has been developed by referring to parametric fault diagnosis techniques. Furthermore, in the automation of the procedure the use of symbolic techniques is of fundamental importance because symbolic analysis, due to the fact that it gives symbolic rather than numerical results, is particularly suitable for applications such as testability and canonical ambiguity group determination, as will be shown later.

It is necessary in this procedure to determine a set of equations describing the circuit under test and solve it with respect to the component values. In the case of analog linear time-invariant circuits, the fault diagnosis equations can be constituted by the network functions relevant to the selected test points [12] which are nonlinear with respect to the potentially faulty circuit parameters. By assuming that the faults can be expressed as parameter variations without influencing the circuit topology (as was done in the previous section where faults like short and open are not considered), the testability measure τis given by the maximum number of linearly independent columns of the Jacobian matrix associated with the fault diagnosis equations, and it represents a measure of the solvability degree of the nonlinear fault diagnosis equations. The entries of the Jacobian matrix are rational functions depending on the complex frequency and the potentially faulty parameters. Thus, in order to evaluate the testability it is necessary to select fixed values for the potentially faulty parameter and the complex frequency. It can be shown that, once the frequency values are fixed, the rank of the obtained Jacobian matrix is constant almost everywhere, i.e., for all the potentially faulty parameter values except those lying in an algebraic variety [13]. Using this approach, the testability value, although independent of component values, is very difficult to handle and subject to round off errors if a numerical approach is used in its automation.

Generally the Jacobian matrix is very costly to find in fully symbolic form. It has been shown in [14] that starting from the network symbolic transfer functions expressed in the following way:

tlh,s=Nlh,sDh,s=i=0nlai(l),      l=1,, KE15

where h= [h1, h2, …, hp]t is the vector of the potentially faulty parameters and Kis the total number of equations, the testability is equal to the rank of a matrix Egiven by


This matrix is independent of the complex frequency whose entries are constituted by the derivatives of the coefficients of the fault diagnosis equations with respect to the potentially faulty circuit parameters. If the fault diagnosis equations are generated in a completely symbolic form, the testability evaluation becomes easy to perform. In this case, the entries of the matrix Ecan be simply led back to derivatives of sums of products and the computational errors are drastically reduced in the automation phase. Once the matrix Ehas been determined, testability evaluation can be performed by triangularizing Eand assigning arbitrary values to the components (since as was previously mentioned, testability does not depend on component values). Yet selecting the matrix Einstead of the Jacobian matrix as the testability matrix results in a different testability measure not directly related to the desired measure. However, this limitation can be overcome by splitting the fault diagnosis equation solution into two phases. In the first phase, starting from the measurements carried out on the selected test points at different frequencies, the coefficients of the fault diagnosis equations are evaluated, eventually exploiting a least-squares procedure in order to minimize the error due to measurement inaccuracy. In the second phase, the component values are obtained by solving the nonlinear system constituted by the equations expressing the previously determined coefficients as functions of the circuit parameters. In this way the following nonlinear system has to be solved:


where Ai(l)and Bj(i=0,…, nl, j=0, …, m-1) are the coefficients of the fault diagnosis equations in (15) which have been calculated in the previous phase. The Jacobian matrix of this system coincides with the matrix Ein (16), hence, all the information provided by a Jacobian matrix with respect to its corresponding nonlinear system can be obtained from the matrix E. In particular, if rank(E) is equal to the number of unknown parameters, the component values can be uniquely determined by solving the equations in (17) through the consideration of a set of measurements carried out on the test points. If the testability τ= rank(E) is less than the number of unknown parameters R, a locally unique solution can be determined only if R-τcomponents are considered not faulty.

The matrix Edoes not give only information about the global solvability degree of the fault diagnosis problem. In fact, by noting that each column is relevant to a specific element or parameter of the circuit and by considering the linearly dependent columns of E, other information can also be obtained. For example, if a column is linearly dependent with respect to another one, this means that a variation of the corresponding component provides a variation on the fault-equation coefficients, indistinguishable with respect to that produced by the variation of the component corresponding to the other column. This means that the two components are not testable and they constitute an ambiguity group of the second order. As an example two parallel connected resistors in a circuit where we cannot distinguish which one caused the fault. By extending this reasoning to groups of linearly dependent columns of E, ambiguity groups of a higher order can be found. Then, in summary, the following definition can be formulated.

Definition1:A set of components constitutes an ambiguity group of order jif the corresponding columns of the testability matrix Eare linearly dependent. In other words, the ambiguity groups of a circuit in which a certain test point set has been chosen can be determined by locating the linearly dependent columns of the testability matrix E. Furthermore, as was mentioned, an ambiguity group that does not contain other ambiguity groups is called canonical. Therefore, a canonical ambiguity group can be defined as follows.

Definition2:A set of kcomponents constitutes a canonical ambiguity group of order kif the corresponding kcolumns of the testability matrix Eare linearly dependent and every subset of this group of columns is constituted by linearly independent columns. It is important to notice that with this definition, the order of the canonical ambiguity groups cannot be greater than the testability value plus one τ+1.

In most cases the canonical ambiguity groups have some components in common. By unifying these types of groups, another ambiguity group, corresponding again to linearly dependent columns of the matrix Eis obtained. We define as global an ambiguity group of the following type.

Definition3:A set of mcomponents constitutes a global ambiguity group of order mif it is obtained by unifying canonical ambiguity groups having at least one element in common.

Obviously, a canonical ambiguity group which does not have components in common with any other canonical ambiguity group can be considered as a global ambiguity group. Finally, the columns of the matrix Ethat do not belong to any ambiguity group are linearly independent. We define these as surely testable a group of components of the following kind.

Definition 4:A set of ncomponents whose corresponding columns of the testability matrix Edo not belong to any ambiguity group constitutes a surely testable group of order n.

Obviously, the number of surely testable components cannot be greater than the testability value τ, that is, the rank of the matrix E.

With these definitions in mind, the optimum set of testable components can be determined as in [12].

3.2. Formulation of fault equations

Applying equation (1) to fault-free and faulty circuits, respectively, with the same excitation sources we get




It can be easily shown that


where the parameter variation can be found from the measured values of the faulty circuit, the original system matrix and an estimate of the change of the system matrix due to fault presence


It is customary to solve equation (23) as a constrained linear optimization problem. However such an approach is limited by the solution time and ambiguity leading to local minimum convergence. Suppose that the first fof pparameters are faulty and are changed from their nominal values h10, h20, …., hf0to new values h1=h10+d1, h2=h20+d2, …., hf=hf0+df, where d1, d2, …., dfare the parameter deviations and the deviation vector dis an f×1vector:

d=[d1 d2 df]tE24

Define Fas the faulty parameter set, and assume that each faulty parameter Fv( v= 1, 2,..., f) is located on intersection of the corresponding rows ivand jvand columns kvand lvof the coefficient matrix T. The deviation of the coefficient matrices now has the following form:


where diag(d) is an f×fdiagonal matrix and Pfand Qfare g×fmatrices which contain 0 and ±1 entries:

Pf=p1 p2  pf=[δi1-δj1 δi2-δj2  δif-δjf]Qf=q1 q2  qf=[δk1-δl1 δk2-δl2  δkf-δlf]E26

Note that Pfand Qfare sub-matrices of Pand Qrespectively. They can be constructed from Pand Qby selecting all columns in Pand Qcorresponding to faulty parameters. As an example assume that there are two faulty parameters: R9 is changed from 3.54Ω to 7.9Ω and R37 is changed from 3.12Ω to 2.8 Ω. The corresponding admittance deviations are ΔG9= 1/7.9 - 1/3.54 = -0.1559 / Ω and ΔG37=1/2.8 ̶ 2/3.12 = 0.03663 /Ω. The corresponding faulty parameter set F=[9,37] and the faulty nodes will be [1,14,16]. It can be easily verified that Pfand Qfare 20×2-matrices which can be obtained from the 9thand 37th columns of the matrices Pand Qrespectively. It can also be verified that ΔTwill have entries only at locations {1,1}, {14,14}, {14,16}, {16,14}, {16,16}.

Substituting (25) in (20) we get


and to obtain the solution vector for the faulty circuit we use


It can be shown using Woodbury formula that the value of dv(v=1,2,..., f) cannot be zero or infinity to meet with the requirement of inverting [6]. Since dvbeing zero means fault-free parameter and only faulty parameters will be identified by following fault diagnosis algorithm, we will have only one restriction: dvcannot be infinite, which corresponds to the case of open admittance or short impedance. But open or short faults can be dealt with by ideal switches introduced in modified nodal analysis [4]. Therefore, the proposed method can handle open and short faults as well but only if combined with a procedure that repeats the analysis of the circuit after introduction of ideal switches.

The solution vector for fault-free circuit is

X0=[x1,0 x2,0 xg,0]tE29

where subscript 0 indicates that the denoted parameters are for fault-free circuit. Hence the product of Qftand Xo can be written as

QftX0=δk1-δl1 δk2-δl2  δkf-δlftX0=xk1,0-xl1,0 xk2,0-xl2,0  xkf,0-xlf,0t=xk1l1,0 xk2l2,0  xkflf,0tE30

and it has the physical interpretation of controlling nominal signal values (e.g. voltages) on faulty parameter input terminals.

Let us define

β=[β1 β2  βn]t=T0-1Pfγ=QftT0-1Pf=QftβE31

It can be shown that the deviation vector ΔXcan be obtained by [6]



α=-βdiagd-1+γ-1=α11α12α21α22α1fα2fαg1αg2αgf=α1α2αg  E33

Usually voltage measurement is easier to carry out and is less invasive to analog circuit properties than current measurement. Therefore, we only consider the use of nodal voltage measurement in this formulation. As an example γfor our example circuit of Figure 1 assuming the aforementioned faults will be given by


Once we have the fault equations formulated and the faults simulated we can proceed to fault diagnosis.

3.3. Fault diagnosis

During the fault diagnosis we have the CUT with only a limited set of accessible nodes for measurement and excitation. Suppose the ithnode is accessible for measurement, then by equation (32)

Xi=αi1αi2αifxk1l1,0 xk2l2,0  xkflf,0tE35

According to definition of g×fmatrix αin equations (33) and (31), matrix αdoes not depend on the location of excitation sources. Thus matrix αis invariant when applying the multiple excitation method, i.e., the same coefficients αijlinks deviation of measurements ΔXiand nominal signal values on faulty parameter xkjljindependent of the excitation vector applied. After measuring the corresponding nodal voltages on the ithnode with mindependent excitation vectors We(e= 1, 2,..., m), we then obtain

Xi(1)=αi1αi2αifxk1l1,0(1) xk2l2,0(1)  xkflf,0(1)tXi(2)=αi1αi2αifxk1l1,0(2) xk2l2,0(2)  xkflf,0(2)tXi(m)=αi1αi2αifxk1l1,0(m) xk2l2,0(m)  xkflf,0(m)tE36

or in matrix form

XiM=Xi1Xi2Xim=xk1l1,01xk2l2,01  xkflf,01xk1l1,02xk2l2,02  xkflf,02xk1l1,0mxk2l2,0m  xkflf,0mαi1αi2αif=XbMFαiE37

where superscript Mdenotes the set of multiple excitations and mis the number of these excitations. The single measurement node can be one of the nodes used for multiple excitation method, and then the total number of accessible excitation nodes should be m. Assume that fm-1 ≤ p, then the coefficient matrix XbMFhas more rows than columns thus guaranteeing the uniqueness of the solution to equation (37) with verification. Equation (37) establishes the linear relationship between the measured responses of the faulty circuit ΔXiMand the faulty parameter deviations dsince vector αiis a linear functions of daccording to equation (33). Therefore equation (37) is called the fault diagnosis equation, and the coefficient matrix XbMFis called the fault diagnosis matrix[6].

As said earlier, with only a limited number of accessible nodes the issue of testability and consistency of the selected set of accessible nodes to detect fnumber of simultaneous faults immediately arises. However, testability is not the focus of this chapter. We assume that the given measurement set can give at least one finite solution to circuit parameters.

As the first stage of fault diagnosis, fault detection is easily implemented. If the measurement deviation vector ΔXiMin the fault diagnosis equation is a zero vector, obviously the CUT is judged as fault-free for the given excitation and measurement sets. Otherwise, at least one fault is judged detected by the given measurement set. To identify the faulty parameters, first let us analyze the fault diagnosis equation. The left-side of equation (37) is a known vector from measurements; the right side is the product of an unknown coefficient matrix XbMFand an unknown solution vector αi. According to equation (30), matrix XbMFis determined by faulty parameter locations and X0, solution vector for the fault-free circuit. Hence the columns in XbMFrepresent the differences between the nominal values of nodal voltages or parameter currents across the 2 input nodes of the faulty parameters. Although we do not know matrix XbMFinitially for the CUT since we do not know initially the location or number of faults, but we really know all of the nodal voltages and parameter currents in the fault-free circuit!

Similarly as in equation (30), we can construct a new m×pmatrix XbMPas follows

QtX0=δk1-δl1 δk2-δl2  δkp-δlptX0=xk1,0-xl1,0 xk2,0-xl2,0  xkp,0-xlp,0t=xk1l1,0 xk2l2,0  xkplp,0tE38
XbMP=xk1l1,01xk2l2,01  xkplp,01xk1l1,02xk2l2,02  xkplp,02xk1l1,0mxk2l2,0m  xkplp,0mE39

where superscript Pdenotes the set of all pcircuit parameters. Each column of XbMPcorresponds to one circuit parameter. Evidently, the fault diagnosis matrix XbMFis a sub-matrix of XbMPand can be constructed by collecting all columns in XbMPcorresponding to the faulty parameters. Apparently matrix XbMFhas more rows than columns whereas XbMPhas less rows than columns due to the restriction fm-1 ≤ p.

For the purpose of fault identification, we need to find out which set or sets of columns in XbMPcan satisfy the fault diagnosis equation, i.e. the dependency between ΔXiMand the desired coefficient matrix in fault diagnosis matrix.

Basically ΔXiMvector for all pparameters has to be generated from the fault-free circuit and stored as a dictionary of fault-free response to mmultiple excitations over the designated maccessible nodes. This dictionary will be used later to determine whether the CUT is faulty and will be used in locating the faults. It must be emphasized that only one node for voltage measurement is sufficient for this method although multiple linearly-independent excitations are required across all maccessible nodes for successful fault location. It is thus possible to use only one of the accessible mnodes to carry out the measurements while using the rest to carry out the excitations. As an example node {2} in Figure1 is selected as the only measurement node, while nodes {2, 4, 15, 16, 17} are selected as accessible nodes for the multiple excitations. That is the unit current source is applied to these nodes respectively and the corresponding nodal voltage at node {2} is measured. Thus the measured changes of nodal voltage will be


One obvious way is to have a combinatorial search through all columns in XbMP, which is the traditional way in the fault verification method [15] and requires a number of operations of the order Oi=1fpifor flimited faults among pparameters. This is equivalent to assuming that any number of faults up to fsimultaneous faults have occurred randomly in any subset of the pparameters then evaluate the response to such faults and compare it to the measured response. However, the method being described here is more efficient than that and involves locating the minimum size ambiguity group which satisfies the fault diagnosis equation. An ambiguity group is defined as a set of parameters corresponding to linearly dependent columns of XbMPwhich in general does not give a unique solution in fault identification. Minimum size ambiguity groups (called canonical ambiguity groups) can be found using a linear combination matrix with minimum number of non-zero entries as will be shown shortly. But to generate this we need to perform a Gaussian elimination step.

3.3.1. Gaussian eliminaion step

Let us first denote an augmented m×(p+1)matrix BSas the concatenation of the stored dictionary vector ΔXiMand the matrix XbMP:


Then we will normalize the first column of matrix BSto have a unity in its first row,

B^Si,1=BSi,1BS1,1,   i=1, 2, , m.E42

If the first entry of matrix BS, BS(1,1)happens to be zero, just permute or swap the rows of BSso that the first entry BS(1,1)is non-zero. Such a nonzero entry must exist since ΔXiMis a non-zero vector for faulty circuit. Eliminate the remaining entries in the first row of matrix BSby performing a similar operation to Gaussian elimination as follows:

B^Si,j=BSi,j-BSi,1BS1,1BS1,j,   i=1, 2, , m;j=2,3,, p+1.E43

Finally we obtain m×(p+1)matrix B^Sin the following form:


where the superscript represents the size of a vector or a matrix. Matrix Bis obtained from XbMPafter elimination of dependence on ΔXiMand is called the verification matrix[6]. The dependency of the desired columns of matrix Bsurely indicates the dependency between ΔXiMand the desired columns of matrix XbMP. Thus we can only concentrate on the dependency among the columns of the verification matrix B.

3.3.2. QR factorization

The rank rof the matrix Bdetermines a maximum number of faults that can be uniquely identified by solving the fault diagnosis equation. Because m-1<p, Bcan be permuted column wise and decomposed into two linearly dependent sub-matrices as follows


where permrefers to column-wise permutation, (m-1)×rmatrix B1has the full column rank equal to the rank rof the matrix B, and r×( pr) matrix Cis called linear combination matrix whose columns expand a set of basis columns from B1into the corresponding columns of B2. It can be easily shown that B1is a sub-matrix of B with all the rows and only a subset of the columns (called the basis set) while B2is a sub-matrix of B with all the rows and the remaining set of the columns (called the co-basis set).Note that the selection of independent columns of B1is not unique and is an important issue in solving the fault diagnosis equation in the presence of ambiguities. Different partitions define different linear combination matrices C.

Since an ambiguity group is a set of circuit parameters corresponding to linearly dependent columns of B, we define a canonical ambiguity group as a minimal set of parameters corresponding to linearly dependent columns of B. This means that if any single parameter is removed from the canonical ambiguity group, then the remaining set corresponds to independent columns of Band can be uniquely solvable. A combination of canonical ambiguity groups with at least one common element was defined as ambiguity cluster.

To efficiently deal with fault verification problem, we will look for a partition (45) with the matrix Cin a minimum form, which is defined as such a matrix that one or several of its columns have the maximum number of entries equal to zero. Thus, we can get the minimum number of columns in XbMPsatisfying the fault diagnosis equation (37). The corresponding partition (45) is called a canonical form of the fault diagnosis equation. Notice that according to fault verification principles [15] it is enough to find a single entry in one column of Cequal to zero to solve the fault diagnosis equation. Yet since many such solutions exist we will select the column with the maximum number of zeros assuming that the faulty response was caused by the smallest number of faults. This column and all rows with non-zero entries will correspond to the faulty parameters as indicated by the element of co-basis B2and elements of basis B1, respectively.

One way to find these matrices from the matrix Bwith high numerical stability is based on QRfactorization [8], which can find a solution of over determined system of linear equations that minimizes the least square error. As a result of the QRfactorization of (m-1)×pverification matrix B, we obtain:


where Eis p×pcolumn selection matrix, Qis (m-1)×(m-1)orthogonal matrix, and Ris (m-1)×pupper triangular matrix. Each column of matrix Ehas only one nonzero entry, which is equal to one. Matrix product BErepresents the permutation of the original columns of the verification matrix Brequested in equation (45). Matrix Rhas its rank equal to the rank of matrix B. Since Ris an upper triangular matrix and m-1<p, Rcan be written as


where R1is r×rupper triangular and has its rank equal to the rank of the verification matrix B. Having this factorization computed, it can be shown that [8]


Furthermore the basis set will be the row values of the non-zero elements in the first rcolumns of Ewhile the co-basis will be the row values of the remaining p-rcolumns of E. As an example, the values of R1 for our example circuit of Figure 1 is

R1=10.5475-2.944403.54350.00282.3965-0.0014-2.80110              00              02.4161-0.00050   1.7202E52

The column permutation is {39, 15, 2, 35, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 3, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 4, 36, 37, 38, 1}. Thus the basis is {39, 15, 2, 35} and co-basis is {5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 3, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 4, 36, 37, 38, 1}.

3.3.3. Swapping performance

A single QRrun cannot guarantee that the matrix Cwill be obtained with one or several of its columns having the maximum number of zero entries if the proper basis is not selected. To find the minimum form partition, we have to swap one parameter of the basis with one parameter of the co-basis in the ambiguity cluster in order to increase number of nonzero entries in C. Note that swapping parameters of the basis and the co-basis can be performed independently in different ambiguity clusters, since different clusters have mutually disjoint sets of parameters. There are simply two conditions to consider in swapping performance:

  1. The necessary condition for swapping to increase the number of zero entries in Cis that the columns of basis and co-basis to be swapped have a singular 2×2 sub-matrix of nonzero entries.

Let us consider a linear combination matrix Cwith a 2×2 singular sub-matrix

R1=10.5475-2.944403.54350.00282.3965-0.0014-2.80110              00              02.4161-0.00050   1.7202E53

with all nonzero entries. If we swap the jthelement of the basis with kthelement of the co-basis, then after swapping, the kthcolumn of Cchanges to


In addition, all other columns of matrix Cwill be equal to


Such that all zero locations in the kthcolumn of Cwill be zero as they were in the original C. However, as can be deducted from (53), a nonzero location cimin row iand column mwill become zero. It is understood that if one element in the current basis has been swapped into the basis by the previous swapping performance, then this element will not be considered during the later swapping.

Any columns of Cwith zero entries form an ambiguity group Fand has to be considered for further processing. Since ambiguities may exist in the original matrix XbMPthen Fcontains all faults in the CUT only if the corresponding columns in XbMPare independent. Hence we must consider the following condition

  1. The necessary condition for an ambiguity group Fof the linear combination matrix Cto contain the set of all faults in the tested circuit is that the rank of the corresponding columns in matrix XbMPis equal to the cardinality of F

rankcolumns in XbMPcorresponding to F=cardFE56

Thus according to this condition any ambiguity group of the verification matrix which do satisfy (55) needs to be further analyzed. The stopping criterion for the above procedure can simply be τthe testability measure found from the symbolic analysis.

As an example for our circuit of Figure 1, careful study of the generated Cmatrix reveals multiple zero entry at columns 5,12, 32, 33, and 34 (corresponding to nodes {9}, {16}, {36}, {37}, and {38} from the co-basis). The non-zero row entries will either be on rows 1, 2, 4 (corresponding to nodes {39}, {15} and {35} from the basis) or on rows 2, 3, and 4 (corresponding to nodes {15}, {2} and {35}). Thus the corresponding ambiguity clusters include {39, 15, 9, 35}, {16, 15, 2, 35}, {39, 15, 36, 35}, {39, 15, 37, 35}, and {39, 15, 38, 35}. Yet none of these ambiguity clusters satisfies condition (b) except for the first one. Accordingly only one suspicious faulty group F={39, 15, 35, 9} is qualified with parameter {9} from the co-basis and parameters {39, 15, 35} from the basis. The current minimum size of qualified F is 4.

Searching for the 2×2 singular matrix with non-zero entries in Creveals that parameter {9} from the co-basis should be swapped with the parameter {39} from the basis according to the swapping procedure in condition (a), and a new matrix Cresults. Re-applying condition (b) to the new matrix C, 5 qualified suspicious faulty groups are obtained: F={9, 2, 35, 5}, F={9, 15, 35, 39}, F={9, 15, 35, 36}, F={9, 15, 38} and F={9, 37}. Obviously, F={9, 37} is the unique solution with the minimum size equal to 2. Since no smaller size of faulty set F can be found by swapping, thus F={9, 37} is the only solution located by the procedure of fault diagnosis which is the exact solution for the given CUT.

Once the fault locations are determined the fault values must be evaluated and compared with the element tolerances before giving a final judgment on the circuit whether being faulty or not.

3.3.4. Parameter evaluation

After locating the faulty parameters, the matrix XbMFcan be found from the matrix XbMPby taking only the columns corresponding to the fault locations. Then the invariant vector αican be uniquely solved from equation (37)


where this form is used since the system is over-determined with XbMFbeing non-square. Finally the deviation vector dcan be exactly computed by

d=αi ./ β-αiγE58

where. / is an element-by-element division of two vectors. What remains after evaluating the deviations is to compare them to the element allowable tolerances to decide finally whether the measured CUT performance is still considered acceptable or deemed faulty. Basically, when the fault locations and parameter deviations are found, all the circuit can be re-solved to get all the node voltages and element currents of the CUT.

3.4. Mixed symbolic numerical algorithm for fault diagnosis

A computer program which implements the fault diagnosis discussed above can be easily advised. In Phase 1, a topological description of the circuit is obtained and the circuit is solved numerically. Since nominal values of circuit parameters are known, all nodal voltages in fault-free circuit can be solved by (18). In phase 2, an upper limit for the testability τneeds to be determined for the provided set of accessible nodes. It is not generally required to obtain a fully symbolic solution of the circuit and only a partial symbolic solution would be sufficient. In phase 3 we need to measure the nodal voltages of the ithnode in the CUT under multiple excitation method to obtain measurement deviation vector ΔXiM. In phase 4 we need to generate the fault locator matrix XbMPfrom equations (38) and (39) then use it to find the linear combination matrix Cafter the Gaussian elimination step and the QR factorization. In Phase 5, analysis of the combination matrix is done where Fdenotes one suspicious fault set and min(size(F)) represents a scalar which is equal to the minimum size of all suspicious fault sets. In Phase 6, if several suspicious fault sets have the same minimum size, min(size(F)), select one of them arbitrarily for analysis. Only one parameter in the selected Fis from the co-basis and the remaining parameters are from the basis. Swap that co-basis parameter which corresponds to column kin matrix Cwith one of basis parameters which corresponds to row jin the matrix C. By (53) and (54), all zero entries in the column k of matrix Cwill be held after swapping while new zero-entry will appear in another column of new matrix C, thus the new value of min(size(F)) will be equal to, or less than the old value before swapping.

There are two rules for swapping. One is that row jis selected with nonzero cjkon the intersection of row jand column kof matrix C. Another rule is that if one parameter in the current basis has been swapped into the basis by the previous swapping operation, then this element will not be considered during the later swapping operation. Usually m-1 is far less than p, and the rank of r×(p-r) matrix C, ris not greater than m-1, thus there are far less basis parameters than co-basis parameters. The comprehensive swapping between the co-basis parameter kand the basis parameters are very limited, as a result of the two swapping principles.

In Phase 7 equivalent adjoint suspicious fault sets are recorded. In Phase 8 the corresponding fault diagnosis matrices XbMFare found from the fault locator matrix. In Phase 9 the invariant vector αiis evaluated. Phase 10 is used for verification. One or several suspicious fault sets with minimum size are used to compute the deviation vector ΔX. If a computed vector matches the real measured vector ΔXiM, the corresponding fault set Fis our final solution to faulty parameters. Otherwise, we discard this set, and turn to the adjoint suspicious fault sets recorded in Phase 7. Verification in this phase continues until at least one qualified solution to faulty parameters is found. Otherwise, the CUT is concluded as un-solvable because the restriction f≤ m - 1 is not satisfied. In the final Phase the parameter variations are compared to the element tolerances to decide if the circuit response is indeed faulty or just shifted within the accepted tolerance.

4. Conclusion

In this chapter, a generalized fault diagnosis and verification approach for linear analog circuits was discussed. Fault verification methods intend to obtain the information about the faulty circuit based on the limited measured responses of the faulty circuit. There are two easily implemented prerequisites: one is that the circuit topology and nominal values of circuit parameters should be known, another is that the number of measurements minus one is not less than the number of faulty parameters. A symbolic method is proposed to solve the testability problem during the detection, and location of the multiple faults in a linear analog circuit in frequency domain, then to exactly evaluate the faulty parameter deviations.

Applying the Woodbury formula in the matrix theory to the modified nodal analysis, fault diagnosis equation is constructed to establish the relationship between the measured responses and the faulty parameter deviations in a linear way. A numerically robust approach has been modified to fit the condition stated in this chapter in order to implement fault location, i.e., location of the minimum size ambiguity group in the fault diagnosis equation based on QR factorization. Parameter evaluation is then performed from results of the analysis of fault diagnosis equation.

One node for voltage measurement is sufficient for the proposed method although multiple excitations are required for fault location. Although the faulty parameter deviation cannot be infinity, open or short condition can be dealt with well by switches in modified nodal analysis.

Therefore, the faults can be parametric or catastrophic. The proposed method is extremely effective for large parameter deviations and a very limited number of accessible nodes used for excitations and measurements. The computation cost for the fault location is on the order of O(p3), and compares favorably with the combinatorial search traditionally used in fault verification methods which requires the number of operations Oi=1fpi.

© 2013 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Fawzi M Al-Naima and Bessam Z Al-Jewad (January 9th 2013). Fault Diagnosis in Analog Circuits via Symbolic Analysis Techniques, Analog Circuits, Yuping Wu, IntechOpen, DOI: 10.5772/53643. Available from:

chapter statistics

4565total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

A Successive Approximation ADC using PWM Technique for Bio-Medical Applications

By Tales Cleber Pimenta, Gustavo Della Colletta, Odilon Dutra, Paulo C. Crepaldi, Leonardo B. Zocal and Luis Henrique de C. Ferreira

Related Book

First chapter

Analog CMOS Design Automation Methodologies for Low-Power Applications

By Alessandro Girardi and Lucas C. Severo

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us