Open access peer-reviewed chapter

Low-Noise, Low-Sensitivity Active-RC Allpole Filters Using MATLAB Optimization

By Dražen Jurišić

Submitted: November 4th 2010Reviewed: April 22nd 2011Published: September 9th 2011

DOI: 10.5772/19745

Downloaded: 4428

1. Introduction

The application of Matlab, combining its symbolic and numeric calculation capabilities, to calculate noise and sensitivity properties of allpole active-RC filters is shown. Transfer function coefficients calculations, as well as plotting of amplitude-frequency and phase-frequency characteristics (Bode plots) have been performed using Matlab. Thus, using Matlab a comparison of different design strategies of active-RC filters is done. It is shown that active-RC filters can be designed to have low sensitivity to passive components and at the same time possess low output thermal noise. The classical methods were used to determine output noise of the filters. It was found that low-sensitivity filters with minimum noise have reduced resistance levels, low Q-factors, low-noise operational amplifiers (opamps) and use impedance tapering design. The design procedure of low-noise and low-sensitivity, positive- and negative-feedback, second- and third-order low-pass (LP), high-pass (HP) and band-pass (BP) allpole filters, using impedance tapering, is presented. The optimum designs, regarding both performances of most useful filter sections are summarized (as a cookbook programmed in Matlab) and demonstrated on examples. The relationship between the low sensitivity and low output noise, that are the most important performance of active-RC filters, is investigated, and optimum designs that reduce both performances are presented.

A considerable improvement in sensitivity of single-amplifier active-RC allpole filters to passive circuit components is achieved using the design technique called 'impedance tapering' (Moschytz, 1999), and as shown in (Jurisic et al., 2010 a)at the same time they will have low output thermal noise. The improvement in noise and sensitivity comes free of charge, in that it requires simply the selection of appropriate component values. Preliminary results of the investigation of the relation between low sensitivity and low thermal noise performances using impedance tapering on the numeric basis using Matlab have been presented in (Jurisic & Moschytz, 2000; Jurisic, 2002).

For LP filters of second- and third-order the complete analytical proofs for noise properties of the desensitized filters are given in (Jurisic et al., 2010 a). By means of classical methods as in (Jurisic et al., 2010 a) closed-form expressions are derived in (Jurisic et al., 2010 c), providing insight into noise characteristics of the LP, HP and BP active-RC filters using different designs. LP, HP and BP, low-sensitivity and low-noise filter sections using positive and negative feedback, that have been considered in (Jurisic et al., 2010 c) are presented here. These filters are of low power because they use only one opamp per circuit. The design of optimal second- and third-order sections referred to as 'Biquads' and 'Bitriplets', regarding low passive and active sensitivities has been summarized in the table form as a cookbook in (Jurisic et al., 2010 b). For common filter types, such as Butterworth and Chebyshev, design tables with normalized component values for designing single-amplifier LP filters up to the sixth-order with low passive sensitivity to component tolerances have been presented in (Jurisic et al., 2008). The filter sections considered in (Jurisic et al., 2010 c) and repeated here have been recommended in (Moschytz & Horn, 1981; Jurisic et al., 2010 b) as high-quality filter sections. It was shown in (Jurisic & Moschytz, 2000; Jurisic, 2002; Jurisic et al., 2008, 2010 a, 2010b , 2010c), that both noise and sensitivity are directly proportional to the pole Q’s and, therefore, to the pass band ripple specified by the filter requirements. The smaller the required ripple, the lower the pole Q’s. Besides, it is wise to keep the filter order nas low as the specifications will permit, because the lower the filter order, the lower the pole Q’s. Also, it was shown that positive-feedback filter blocks are useful for the realization of the LP and HP filters (belonging to class 4, according to the classification in (Moschytz & Horn, 1981), the representatives are SAK: Sallen and Key filters). Filters with negative feedback (class 3 SAB: Single-amplifier Biquad) are better for the BP filters, where the BP-C Biquad is preferable because it has lower noise than BP-R. A summary of figures and equations that investigate sensitivity and noise performance of active RC filters, and have been calculated in (Jurisic & Moschytz, 2000; Jurisic, 2002; Jurisic et al., 2008, 2010 a, 2010b , 2010c), by Matlab, will be presented here. Numeric and symbolic routines that were used in those calculations are shown here in details.In Section 2 a brief review of noise and sensitivity is given and the most important equations are defined. These equations will be used by Matlab in Section 3 to analyze a second-order LP filter as representative example. In Section 4 the results of analysis using Matlab of the LP, HP and BP sections of second- and third-order filters are summarized. Those results were obtained with the same Matlab algorithms as in Section 3 for the second-order LP filter, and are presented in the form of optimum-design procedures. The chapter ends with the conclusion in Section 5.


2. A brief review of noise and sensitivity of active-RC filters

2.1. Output noise and dynamic range

Thermal (or Johnson) noise is a result of random fluctuations of voltages or currents that seriously limit the processing of signals by analog circuits. Because this noise is caused by random motion of free charges and is proportional to temperature, it is referred to as thermal noise (Jurisic et al., 2010 a).The most important sources of noise in active-RC filters are resistors and opamps. For the purpose of noise analysis, appropriate noise models for resistors and opamps must be used. Resistors are represented by the well-known Nyquist voltage or current noise models shown in Figure 1(a) and (b), consisting of noiseless resistors and noise sources whose values are defined by the squared noise voltage density within the narrow frequency band f, i.e.,


or the squared noise current density given by


where k=1.381023[J/K] is Boltzmann's constant, Tis the absolute temperature of a conductor in Kelvin [K]. All examples are calculated for 22C (T=295K), i.e., room temperature.

Figure 1.

a) Voltage noise model of a resistor. (b) Current noise model of a resistor. (c) An opamp noise model.

The noise defined by (Eq. 1) and (Eq. 2) has a constant spectrum over the frequency band, and is referred to as 'white noise'. The squared noise spectral density in (Eq. 1) has the dimension [V2/Hz], unless writtenenR2(ω)=(2kT/π)R; in which case it has the dimension [V2/rad/s]. The dimension of the spectrum in (Eq. 2) is [A2/Hz], unless writteninR2(ω)=2kT/(πR); in which case it has the dimension [A2/rad/s]. The noise in real capacitors is also of thermal origin and is negligible. The noise in opams is caused by the built-in semiconductors and resistors. The equivalent schematic of a noisy opamp is shown in Figure 1(c), i.e., a noiseless opamp combined with voltage and current noise sources. For the TL081/TI (Texas instruments) FET input opamp, typical values found in the data-sheets are ena(f)=17nV/Hz and ina1(f)ina2(f)=0.01pA/Hz. These values are considered constant within the frequency interval up to about 50 kHz and have been used in the noise analysis here.The noise is additive and the spectral power density of the noise voltage at the output terminal is obtained by adding the contributions from each source. Thus, the squared output noise spectral density, derived from all the noise sources and their corresponding noise transfer functions, is given by (Schaumann et al., 1990):


where Ti,k(j) is the transfer impedance, i.e. the ratio of the output voltage and input current of the kth current noise source (in)k, and T,l(j) is the corresponding voltage transfer function, i.e. the ratio of the output voltage and the input voltage of the lth voltage noise source (en)l.The total output noise power is obtained by the integration of the mean-square noise spectral density e2no() in (Eq. 3) over the total frequency band from 0 to ∞; thus:


The dynamic range is defined by:

DR=20log(Vsorms)max(Eno)rms   [dB]E5

where (Vsorms)maxrepresents the maximum undistorted rms voltage at the output, and the denominator is the noise floor defined by the square root of (Eq. 4). (Vsorms)maxis determined by the opamp slew rate, power supply voltage, and the corresponding THD factor of the filter. In our examples we use a 10Vpp signal which yields


2.2. Sensitivity to passive component variations

Sensitivity analysis provides information on network changes caused by small deviations of passive component values. Given the network function F(s, x1,, xN), where s is a complex variable and xk(k=1,, N) are real parameters of the filter, the relative deviation of F, F/F, due to the relative deviations xk/xk(k=1,, N) is given to the first approximation by:


where SxkFrepresents the relative sensitivity of the function F to variations of a single parameter (component) xk, namely:


If several components deviate from the nominal value, a criterion for assessing the deviation of the function F due to the change of several parameters must be used. With xk/xkconsidered to be an independent random variable with zero mean and identical standard deviation x, the squared standard deviation 2Fof the relative change F/F is given by:


Fis therefore dependent on the component sensitivitiesSxkF, but also on the number of passive components N. The more components the circuit has, the larger the sensitivity. Equation 9 defines multi-parametric measure of sensitivity (Schoeffler, 1964; Laker & Gaussi, 1975; Schaumann et al., 1990).In the following Section, all Matlab calculations regarding noise and sensitivity performance will be demonstrated on the second-order LP filter circuit with positive feedback (class-4 or Sallen and Key). All Matlab commands and variables will appear in the text using Courier New font.

3. Application to second-order LP filter

3.1. Calculating transfer function coefficients and parameters using 'symbolic toolbox' in Matlab

Consider the second-order low-pass active-RC allpole filter circuit (Biquad) shown in Figure 2(a). This circuit belongs to the positive feedback or class-4 (Sallen and Key) filters (Moschytz & Horn, 1981). In Figure 2(b)there is a simplified version of the same circuit with the voltage-controlled voltage source (VCVS) having voltage gain . For an ideal opamp in the non-inverting mode it is given by


Note that the voltage gain of the class-4 circuit is positive and larger than or equal to unity. Voltage transfer function for the filters in Figure 2 expressed in terms of the coefficients ai(i=0, 1, 2) is given by


and in terms of the pole frequency p, the pole Q, qpand the gain factor K by:




Figure 2.

Second-order Sallen and Key LP active-RC filter. (a) With ideal opamp having feedback resistors RFand RG, and nodes for transfer-function calculus. (b) Simplified circuit with the gain element replaced by VCVS .

To calculate the voltage transfer function T(s)=Vout(s)/Vg(s) of the Biquad in Figure 2(a), consider the following system of nodal equations (note that the last equation represents the opamp):


The system of Equation 12 can be solved using 'Symbolic toolbox' in Matlab. The following Matlab code solves the system of equations:Matlab command syms defines symbolic variables in Matlab's workspace:

syms A R1 R2 C1 C2 RF RG s Vg V1 V2 V3 V4 V5;

Matlab command solve is used to solve analytically above system of five Equation 12 for the five voltages V1 to V5 as unknowns. The unknowns are defined in the last row of command solve. Note that all variables used in solve are defined as symbolic.



'-V1*1/R1 + V2*(1/R1+1/R2+s*C1)-V3*1/R2 - V5*s*C1=0',...

'-V2*1/R2 + V3*(1/R2+s*C2)=0',...


'(V3-V4)*A =V5',...


Once all variables are known simple symbolic division of V5/V1 yields the desired transfer function (limit value for A∞ has to be applied, as well):



Another way of presentation polynomials is by collecting all coefficients that multiply 's':


Transfer function coefficients and parameters readily follow.To obtain coefficients, it is useful to separate numerator and denominator using the following command:


syms a2 a1 a0 wp qp k;



Now coefficients follow




And parameters




Typing command whos we obtain the following answer about variables in Matlab workspace:

>> whos

Name Size Bytes Class

A 1x1 126 sym object

C1 1x1 128 sym object

C2 1x1 128 sym object

CircuitEquations 1x1 2828 struct array

Tofs 1x1 496 sym object

Tofsa 1x1 252 sym object

Tofsc 1x1 248 sym object

R1 1x1 128 sym object

R2 1x1 128 sym object

RF 1x1 128 sym object

RG 1x1 128 sym object

V1 1x1 128 sym object

V2 1x1 128 sym object

V3 1x1 128 sym object

V4 1x1 128 sym object

V5 1x1 128 sym object

Vg 1x1 128 sym object

a0 1x1 150 sym object

a1 1x1 210 sym object

a2 1x1 126 sym object

denTa 1x1 232 sym object

denLP2 1x3 330 sym object

k 1x1 144 sym object

numTa 1x1 134 sym object

numLP2 1x1 144 sym object

qp 1x1 254 sym object

s 1x1 126 sym object

wp 1x1 166 sym object

Grand total is 1436 elements using 7502 bytes

It can be seen that all variables that are defined and calculated so far are of symbolic type.We can now check the values of the variables. For example we are interested in voltage transfer function Tofsa. Matlab gives the following answer, when we invoke the variable:

>> Tofsa

Tofsa =


The command pretty presents the results in a more beautiful way.

>> pretty(Tofsa)




s C2 R2 RG + R2 s C1 R1 C2 RG - s C1 R1 RF + RG + R1 s C2 RG

Or we could invoke variable Tofsc (see above that Tofsc is the same as Tofsa, but with collected coefficients that multiply powers of 's').

>> pretty(Tofsc)




R2 s C1 R1 C2 RG + (C2 R2 RG - C1 R1 RF + R1 C2 RG) s + RG

Other variables follow using pretty command.

>> pretty(a0)



R2 C1 R1 C2

>> pretty(a1)

C2 R2 RG - C1 R1 RF + R1 C2 RG


RG R2 C1 R1 C2

>> pretty(a2)


>> pretty(wp)

/ 1 \1/2


\R2 C1 R1 C2/

>> pretty(qp)

/ 1 \1/2

|-----------| RG R2 C1 R1 C2

\R2 C1 R1 C2/


C2 R2 RG - C1 R1 RF + R1 C2 RG

>> pretty(k)




Next, according to simplified circuit in Figure 2(b) having the replacement of the gain element by defined in (Eq. 10), we can substitute values for RFand RGusing the command subs and obtain simpler results [in the following example we perform substitution RFRG(–1)]. New symbolic variable is beta

>> syms beta

>> a1=subs(a1,RF,'(beta-1)*RG');

>> pretty(a1)

C2 R2 RG - C1 R1 (beta - 1) RG + R1 C2 RG


RG R2 C1 R1 C2

Note that we have obtained RGboth in the numerator and denominator, and it can be abbreviated. To simplify equations it is possible to use several Matlab commands for simplifications. For example, to rewrite the coefficient a1 in several other forms, we can use commands for simplification, such as:

>> pretty(simple(a1))

1 beta 1 1

----- - ----- + ----- + -----

C1 R1 R2 C2 R2 C2 R2 C1

>> pretty(simplify(a1))

-C2 R2 + C1 R1 beta - C1 R1 - R1 C2

- -----------------------------------

R2 C1 R1 C2

The final form of the coefficient a1 is the simplest one, and is the same as in (11c) above. Using the same Matlab procedures as presented above, we have calculated all coefficients and parameters of the different filters' transfer functions in this Chapter.If we want to calculate the numerical values of coefficients ai(i=0, 1, 2) when component values are given, we simply use subs command. First we define the (e.g. normalized) numerical values of components in the Matlab’s workspace, and then we invoke subs:

>> R1=1;R2=1;C1=0.5;C2=2;

>> a0val=subs(a0)

a0val =


>> whos a0 a0val

Name Size Bytes Class

a0 1x1 150 sym object

a0val 1x1 8 double array

Grand total is 15 elements using 158 bytes

Note that the new variable a0val is of the double type and has numerical value equal to 1, whereas the symbolic variable a0 did not change its type. Numerical variables are of type double.

3.2. Drawing amplitude- and phase-frequency characteristics of transfer function using symbolic and numeric calculations in Matlab

Suppose we now want to plot Bode diagram of the transfer function, e.g. of the Tofsa, using the symbolic solutions already available (see above). We present the usage of the Matlab in numeric way, as well. Suppose we already have symbolic values in the Workspace such as:

>> pretty(Tofsa)




RG - C1 R1 RF s + C2 R1 RG s + C2 R2 RG s + C1 C2 R1 R2 RG s

Define set of element values (normalized):

>> R1=1;R2=1;C1=1;C2=1;RG=1;RF=1.8;

Now the variables representing elements R1, R2, C1, C2, RG, and RFchanged in the workspace to double and have values; they become numeric. Substitute those elements into transfer function Tofsa using the command subs.

>> Tofsa1=subs(Tofsa);

>> pretty(Tofsa1)



/ 2 s \

5 | s + - + 1 |

\ 5 /

Note that in new transfer function Tofsa1 an independent variable is symbolic variable s. To calculate the amplitude-frequency characteristic, i.e., the magnitude of the filter's voltage transfer function we first have to define frequency range of , as a vector of discrete values in wd, make substitution s=j into T(s) (in Matlab represented by Tofsa1) to obtain T(j), and finally calculate absolute value of the magnitude in dB by α(ω)=20 log |T(jω)|. The phase-frequency characteristic is φ(ω)=arg T(jω) and is calculated using atan2(). This can be performed in following sequence of commands:

wd = logspace(-1,1,200);

ad1 = subs(Tofsa1,s,i*wd);


semilogx(wd, Alphad, 'g-');

axis([wd(1) wd(end) -40 30]);

title('Amplitude Characteristic');

legend('Circuit 1 (normalized)');

xlabel('Frequency /rad/s');ylabel('Magnitude / dB');



semilogx(wd, Phid, 'g-');

axis([wd(1) wd(end) -180 0]);

title('Phase Characteristic');

legend('Circuit 1 (normalized)');

xlabel('Frequency /rad/s');ylabel('Phase / deg');


Commands are self-explanatory. The amplitude- and phase-frequency characteristics thus obtained are shown in Figure 3. Note that we have generated vectors of values wd, Alphad and Phid to be plotted in logarithmic scale by the command semilogx (instead, we could have used command plot to generate linear axis).The next example defines new set of second-order LP filter element values (those are obtained when above normalized elements are denormalized to the frequency 0=286103 rad/s and impedance R0=37k; see in (Jurisic et al., 2008) how):

>> R1=37e3;R2=37e3;C1=50e-12;C2=50e-12;RG=1e4;RF=1.8e4;

Those element values were calculated starting from transfer function parameters p= 286103 rad/s and qp=5 and are represented as example 1) non-tapered filter (=1, and r=1) (see Equation 18 and Table 3in Section 4 below). We refer to those values as 'Circuit 1'.

>> Tofsa2=subs(Tofsa);

>> pretty(Tofsa2)




800318296602402496323046008438980478515625 s 4473025532574128109375 s

-------------------------------------------------- + ------------------------- + 10000

23384026197294446691258957323460528314494920687616 1208925819614629174706176

Figure 3.

Transfer-function (a) magnitude and (b) phase for Circuit 1 (normalized).

It is seen that the denormalized-transfer-function presentation in symbolic way is not very useful. It is possible rather to use numeric and vector presentation of the Tofsa2. First we have to separate numerator and denominator of Tofsa2 by typing:

>> [num2, den2]=numden(Tofsa2);

then we have to convert obtained symbolic data of num2 and den2 into vectors n2 and d2:

>> n2=sym2poly(num2)

n2 =


>> d2=sym2poly(den2)

d2 =

1.0e+053 *

0.0000 0.0000 2.3384

and finally use command tf to write transfer function which uses vectors with numeric values:

>> tf(n2,d2)

Transfer function:



8.003e041 s^2 + 8.652e046 s + 2.338e053

If we divide numerator and denominator by the coefficient of s2 in the denominator, i.e., d2(1), we have a more appropriate form:

>> tf(n2/d2(1),d2/d2(1))

Transfer function:



s^2 + 1.081e005 s + 2.922e011

Obviously, the use of Matlab (numeric) vectors provides a more compact and useful representation of the denormalized transfer function.Finally, note that when several (N) filter sections are connected in a cascade, the overall transfer function of that cascade can be very simply calculated by symbolic multiplication of sections' transfer functions Ti(s) (i=1, , N), i.e. T=T1**TN, if Ti(s) are defined in a symbolic way. On the other hand, if numerator and denominator polynomials of Ti(s) are defined numerically (i.e. in a vector form), a more complicated procedure of multiplying vectors using (convolution) command conv should be used.

3.3. Calculating noise transfer function using symbolic calculations in Matlab

Using the noise models for the resistors and opamps from Figure 1, we obtain noise spot sources shown in Figure 4(a).

Figure 4.

a) Noise sources for second-order LP filter. (b) Noise transfer function for contribution of R1.

The noise transfer functions as in (Eq. 3) Tx(s)=Vout/Nxfrom each equivalent voltage or current noise source to the output of the filter in Figure 4(a) has to be evaluated.As a first example we find the contribution of noise produced by resistor R1 at the filter's output. We have to calculate the transfer resistance Ti,R1(s)=Vout(s)/InR1(s). According to Figure 4(b) we write the following system of nodal equations:


The system of Equation 13 can be solved using Matlab Symbolic toolbox in the same way as the system of Equation 12 presented above. The following Matlab code solves the system of Equation 13:



'-V1*1/R1 + V2*(1/R1+1/R2+s*C1)-V3*1/R2 - V5*s*C1=InR1',...

'-V2*1/R2 + V3*(1/R2+s*C2)=0',...


'(V3-V4)*A =V5',...





syms a2 a1 a0 b0



%Coefficients of the transfer function





In Matlab workspace we can check the value of each coefficient calculated by above program, simply, by typing the corresponding variable. For example, we present the value of the coefficient b0 in the numerator by typing:

>> pretty(b0)

-- RF + RG --

| ----------- |

-- C1 C2 R2 RG --

The coefficients a0, a1 and a2 are the same as those of the voltage transfer function calculated in Section 3.1 above, which means that two transfer functions have the same denominator, i.e., D(s). Thus, the only useful data is the coefficient b0. The transfer resistance Ti,R1(s) is obtained.The noise transfer functions of all noise spot sources in Figure 4(a)have been calculated and presented in Table 1 in the same way as Ti,R1(s) above. We use current sources in the resistor noise model. Nxis either the voltage or current noise source of the element denoted by x.

3.4. Drawing output noise spectral density of active-RC filters using numeric calculations in Matlab

Noise transfer functions for second-order LP filter, generated using Matlab in Section 3.3, are shown in Table 1. We can retype them and use Matlab in only numerical mode to calculate noise spectral density curves at the output, that are defined as a square root of (Eq. 3). Define set of element values (Circuit 1)

>> R1=37e3;R2=37e3;C1=50e-12;C2=50e-12;RG=1e4;RF=1.8e4;

inR1,inR11, inR12β1R2C1C2/D(s)
ina2, inRG, inRF, ena*RF[s2+R2C2+R1C2+R1C1R1R2C1C2s+1R1R2C1C2]/D(s)

Table 1.

Noise transfer functions for second-order LP filter (*enahas instead –RF).

We draw the curve:



Fstart=1e4; %Hz

Fstop=1e6; %Hz

fd =logspace(log10(Fstart),log10(Fstop),Nfreq);

% NOISE SOURCES at temperature T=295K (22 deg C)









D=1/(R1*R2*C1*C2) - (fd*2*pi).^2 +...








TIM=-RF*(1/(R1*R2*C1*C2)-(fd*2*pi).^2 +...







TR1A =(abs(TR1)).^2;

TR2A =(abs(TR2)).^2;

TIPA =(abs(TIP)).^2;

TIMA =(abs(TIM)).^2;



TEPA =(abs(TEP)).^2;


UR1 =TR1A*IR1^2;

UR2 =TR2A*IR2^2;









titletext=sprintf('Output Noise');title(titletext);

xlabel('Frequency / kHz');

ylabel('Noise Spectral Density / \muV/\surdHz');

axis ([fd(1) fd(Nfreq) 0 3e-6]); grid;

% Numerical integration of Total Noise Power at the Output (RMS)

Eno = sqrt(sum(U22(1:Nfreq))/(Nfreq-1)*(fd(Nfreq)-fd(1)));

To draw the second curve, apply the following method. Define the second set of element values, that are represented as example 4) ideally tapered filter (=4, and r=4), (see Equation 18and Table 3 in Section 4 below). We refer to those values as 'Circuit 2'.

>> R1=23.1e3;R2=92.4e3;C1=80e-12;C2=20e-12;RG=1e4;RF=1.05e4;

>> hold on;

>> redo all above equations; use 'r--' for the second curve shape

>> hold off;

>> legend('Circuit 1', 'Circuit 2');

Output noise spectral density is shown in Figure 5.Furthermore, two values of rms voltages Eno(representing total noise power at the output or the noise floor) as defined by the square root of (Eq. 4), have been calculated as a result of numerical integration in Matlab code given above, and they are as follows: Eno1=176.0 V (Circuit 1 or example #1 in Table 3) and Eno2=127.7 V (Circuit 2 or example #4 in Table 3). They are shown in the last column of Table 3, in Section 4. For all filter examples the rms total output noise Enowas calculated numerically using Matlab and presented in the last column of Tables.To plot output noise spectral density and calculate total output noise voltage it was easy to retype the noise transfer function expressions from Table 1 in Matlab code. In the following Section 3.5 it is shown that retyping of long expressions is sometimes unacceptable (e.g. to calculate the sensitivity). Then we have another option to use Matlab in symbolic mode.

Figure 5.

Output noise spectral density of Circuit 1 and Circuit 2 (denormalized).

3.5. Sensitivity characteristic of active-RC filter using both symbolic and numeric calculations in Matlab

To efficiently calculate multi-parametric sensitivity in (Eq. 9), we use a mixture of symbolic and numeric capabilities of Matlab.Suppose F in (Eq. 7)–( Eq. 9) is our transfer function T(s)=N(s)/D(s) defined by (Eq. 11), where xkare elements R1, R2, C1, C2, RF and RG. We will use previous symbolic results of transfer functions numerator numLP2 and denominator denLP2, and Matlab operation of symbolic differentiation diff to produce relative sensitivity in (Eq. 8). To calculate the transfer function sensitivity as defined by (Eq. 8) we will also apply the following rule:


To construct (Eq. 14), we proceed as follows. The following code reveal numerator and denominator as function of components. (Division of both numerator and denominator by RGis just to have nicer presentation.) First we make the substitution s=j into N(s) and D(s). Then we have to produce absolute values of N(j) and D(j). In the subsequent step we perform symbolic differentiation using Matlab command diff or the operator D.

>> den=simplify(denTa/RG);

>> pretty(den)

2 C1 R1 RF s

C2 R1 s + C2 R2 s + C1 C2 R1 R2 s - ---------- + 1


>> denofw = subs(den,s,i*wd)

denofw =

C2*R2*wd*i - C1*C2*R1*R2*wd^2 + 1 - (C1*R1*RF*wd*i)/RG + C2*R1*wd*i

(To calculate all components and frequency values as real variables we have to retype real and imaginary parts of denofw.)

>> syms wd;

>> redenofw= - C1*C2*R1*R2*wd^2 + 1;

>> imdenofw= C2*R2*wd - (C1*R1*RF*wd)/RG + C2*R1*wd;

>> absden=sqrt(redenofw^2+imdenofw^2);

>> pretty(absden)

/ / C1 R1 RF wd \2 2 2 \1/2

| | C2 R1 wd + C2 R2 wd - ----------- | + (C1 C2 R1 R2 wd - 1) |

\ \ RG / /

>> SDR1=diff(absden,R1)*R1/absden;

>> pretty(SDR1)

/ / C1 RF wd \ / C1 R1 RF wd \ 2 2 \

R1 | 2 | C2 wd - -------- | | C2 R1 wd + C2 R2 wd - ----------- | + 2 C1 C2 R2 wd (C1 C2 R1 R2 wd - 1) |

\ \ RG / \ RG / /


/ / C1 R1 RF wd \2 2 2 \

2 | | C2 R1 wd + C2 R2 wd - ----------- | + (C1 C2 R1 R2 wd - 1) |

\ \ RG / /

The same calculus (with simpler results) can be done for the numerator:

>> num=simplify(numTa/RG);

>> pretty(num)


-- + 1


>> numofw = subs(num,s,i*wd)

numofw =

RF/RG + 1

>> renumofw= RF/RG + 1;

>> imnumofw= 0;

>> absnum=sqrt(renumofw^2+imnumofw^2);

>> pretty(absnum)

/ / RF \2 \1/2

| | -- + 1 | |

\ \ RG / /

>> SNR1=diff(absnum,R1)*R1/absnum;

>> pretty(SNR1)


Sensitivity of the numerator to R1 is zero. We have obviously obtained too long result to be analyzed by observation. We continue to form sensitivities to all remaining components in symbolic form.

>> SDR2=diff(absden,R2)*R2/absden;

>> SDC1=diff(absden,C1)*C1/absden;

>> SDC2=diff(absden,C2)*C2/absden;

>> SDRF=diff(absden,RF)*RF/absden;

>> SDRG=diff(absden,RG)*RG/absden;

>> SNR2=diff(absnum,R2)*R2/absnum;

>> SNC1=diff(absnum,C1)*C1/absnum;

>> SNC2=diff(absnum,C2)*C2/absnum;

>> SNRF=diff(absnum,RF)*RF/absnum;

>> SNRG=diff(absnum,RG)*RG/absnum;

By application of rule (Eq. 14), we form sensitivities to each component, whose squares we finally have to sum, and form (Eq. 9).



The resulting analytical form of multi-parametric sensitivity is as follows:

>> SigmaAlpha=sqrt(SCH)*0.01*8.68588964;

The multiplication by 0.01 defines the standard deviation of all passive elements xin (Eq. 9) tobe 1%. The multiplication by 8.68588965 converts the standard deviation Fin (Eq. 9) into decibels.When typing SigmaAlpha in Matlab's workspace, a very large symbolic expression is obtained. We do not present it here (it is not recommended to try!). Because it is too large neither is it useful for an analytical investigation, nor can it be retyped, nor presented in table form. Instead we will substitute in this large analytical expression for SigmaAlpha component values and draw it numerically. This has more sense.Define first set of element values (Circuit 1 with equal capacitors and equal resistors):

>> R1=37e3;R2=37e3;C1=50e-12;C2=50e-12;RG=1e4;RF=1.8e4;

By equating to values, elements changed in the workspace to double and they have become numeric. Substitute those elements into SigmaAlpha.

>> Schoefler1=subs(SigmaAlpha);

Note that in new variable Schoefler1 independent variable is symbolic wd. To calculate its magnitude, we have to define first the frequency range of , as a vector of discrete values in wd. When the frequency in Hz is defined, we have to multiply it by 2. The frequency assumed ranges from 10kHz to 1MHz.

>> fd = logspace(4,6,200);

>> wd = 2*pi*fd;

>> Sch1 = subs(Schoefler1,wd);

>> semilogx(fd, Sch1, 'g-.');

>> title('Multi-Parametric Sensitivity');

>> xlabel('Frequency / kHz'); ylabel('\sigma_{\alpha} / dB');

>> legend('Circuit 1');

>> axis([fd(1) fd(end) 0 2.5])

>> grid;

This is all needed to plot the sensitivity curve of Circuit 1.To add the second example, we set the element values of Circuit 2 in the Matlab workspace:

>> R1=23.1e3;R2=92.4e3;C1=80e-12;C2=20e-12;RG=1e4;RF=1.05e4;

Then we substitute symbolic elements (components) in the SigmaAlpha with the numeric values of components in the workspace to obtain new numeric vales for sensitivity

>> Schoefler2=subs(SigmaAlpha);

>> Sch2 = subs(Schoefler2,wd);

Finally, to draw both curves we type

>> semilogx(fd, Sch1, 'k-', fd, Sch2, 'r--');

>> title('Multi-Parametric Sensitivity');

>> xlabel('Frequency / kHz'); ylabel('\sigma_{\alpha} / dB');

>> legend('Circuit 1', 'Circuit 2');

>> axis([fd(1) fd(end) 0 2.5])

>> grid;

Sensitivity curves of Circuit 1 and Circuit 2 are shown in Figure 6. Recall that both circuits realize the same transfer-function magnitude which is shown in Figure 3a above. Note that only several lines of Matlab instructions have to be repeated, and none of large analytical expressions have to be retyped.In the following Chapter 4, we will use Matlab routines presented so far to construct examples of different filter designs. According to the results obtained from noise and sensitivity analyses we prove the optimum design.

Figure 6.

Standard deviation of magnitudes of Circuit 1 and Circuit 2 (sensitivity).

4. Application to second- and third-order LP, BP, and HP filters

4.1. Second-order Biquads

Consider the second-order Biquads that realize LP, HP and BP transfer functions, shown in Figure 7. Those are the Biquads that are recommended as high-quality building blocks; see (Moschytz & Horn, 1981; Jurisic et al., 2010 b,2010c). In (Moschytz & Horn, 1981) only the design procedure for min. GSP is given (and by that providing the minimum active sensitivity design). On the basis of component ratios in the passive, frequency-dependent feedback network of the Biquads in Figure 7, defined by:


the detailed step-by-step design of those filters, in the form of cookbook, for optimum passive and active sensitivities as well as low noise is considered in (Jurisic et al., 2010 b2010c). The optimum design is presented in Table 1 in (Jurisic et al., 2010c) and is programmed using Matlab.Note that the Biquads in Figure 7 shown vertically are related by the complementary transformation, whereas those shown horizontally are RC–CR duals of each other. Thus, complementary circuits are LP (class-4: positive feedback) and BP-C (class-3: negative feedback), as well as HP (class-4) and BP-R (class-3). In class-4 case there is , whereas in class-3 there isβ¯, that are related by:


Dual Biquads in Figure 7 are LP and HP (class-4), as well as BP-C and BP-R (class-3); they belong to the same class.

Figure 7.

Second-order LP, HP and BP active-RC filters with impedance scaling factors r and .

Voltage transfer functions for all the filters shown in Figure 7 in terms of the pole frequency p, the pole Q, qpand the gain factor K, are defined by:


where numerators n(s) are given by:


Parameters p, qpand K, as functions of filter components, are given in Table 2. They are calculated using Matlab procedures presented in Section 3.1. Referring to Figure 7, the voltage attenuation factor (0<1), which decouples gains K and , see (Moschytz, 1999), is defined by the voltage divider at the input of the filter circuits. Note that all filters in Figure 7 have the same expressions for p, and that the expressions for pole Q, qpare identical only for complementary circuits. This is the reason why complementary circuits have identical sensitivity properties and share the same optimum design, see (Jurisic et al., 2010 c).

(a) LP and (c) BP-C(b) HP and (d) BP-R
ωp=1R1R2C1C2,qp=R1R2C1C2R1(C1+C2)+R2C2βR1C1, K=((for LPandK=αβ¯qpR1C1/(R2C2)for BP-C.ωp=1R1R2C1C2,qp=R1R2C1C2(R1+R2)C2+R1C1βR2C2, K=((for HP andK=αβ¯qpR2C2/(R1C1)for BP-R.

Table 2.

Transfer function parameters of second-order active-RC filters in Figure 7.

No.Filter\Design Parameterr((C1C2CTOTR1R2RTOTEno
1Non Tapered110.3332.85050100373774176.0
2Capacitively Tapered140.3331.4802010046.346.392.5102.5
3Resistively Tapered410.3335.6505010018.57492.5360.9
4Ideally Tapered440.4442.05802010023.192.5115.6127.7
5Cap-Taper and min. GSP1.8540.3971.58802010034.0262.996.94103.9

Table 3.

Component values and rms output noise Enoof design examples of second-order LP and BP-C filters as in Figure 7(a) and (c) with p=286krad/s and qp=5 (resistors in [k], capacitors in [pF], noise in [V]).

No.Filter\Design Parameterr((C1C2CTOTR1R2RTOTEno
1Non Tapered110.3332.85050100373774201.6
2Capacitively Tapered140.3335.6802010046.346.392.5460.1
3Resistively Tapered410.3331.4505010018.57492.596.73
4Ideally Tapered440.4442.05802010023.192.5115.6137.0
5Res-Taper and min. GSP41.850.3971.58653510019.477.697.0100.3

Table 4.

Component values and rms output noise Enoof design examples of second-order HP and BP-R filters as in Figure 7 (b) and (d)with p=286krad/s and qp=5 (resistors in [k], capacitors in [pF], noise in [V]).

On the other hand, two 'dual' circuits will have dual sensitivities and dual optimum designs. Dual means that the roles of resistor ratios are interchanged by the corresponding capacitor ratios, and vice versa. It is shown in (Jurisic et al., 2010 c) that complementary Biquads have identical noise transfer functions and, therefore, the same output noise. An optimization of both sensitivity and noise performance is possible by varying the general impedance tapering factors (Eq. 15) of the resistors and capacitors in the passive-RC network of the filters in Figure 7, see (Moschytz, 1999; Jurisic et al., 2010 b). By increasing r>1 and/or >1, the R2 and C2 impedances are increased. High-impedance sections are surrounded by dashed rectangles in Figure 7.For illustration, let us consider the following practical design example as one in (Moschytz, 1999):

ωp=2π86 kHz;   qp=5;    CTOT=100 pF.E21

As is shown in (Moschytz, 1999), there are various ways of impedance tapering a circuit. By application of various impedance scaling factors in (Eq. 15), the resulting component values of the different types of tapered LP (and BP-C) circuits are listed in Table 3, and the components of HP (and BP-R) filters are listed in Table 4. The corresponding transfer function magnitudes are shown in Figure 8 using Matlab (see Section 3.2). In order to compare the different circuits with regard to their noise performance, the total capacitance for each is held constant, i.e. CTOT=100pF.A multi-parametric sensitivity analysis was performed using Matlab (see Section 3.5) on the filter examples in Table 3 and 4 with the resistor and capacitor values assumed to be uncorrelated random variables, with zero-mean and 1% standard deviation. The standard deviation ()[dB] of the variation of the logarithmic gain Δα=8.68588Δ|T(ω)|/|T(ω)| [dB] was calculated, with respect to all passive components, and plotted for the cases in Table 3 and 4in Figure 9. There exist four different plots for all four Biquads in Figure 7.In Figure 9(a) and (c) it is shown that the LP and BP-C filters no. 2, i.e. the capacitively-tapered filters with equal resistors (=4 and r=1) have the minimum sensitivity to passive component variations (Moschytz, 1999). The next best result is obtained with filter no. 5, i.e. the capacitively-tapered filter with minimum Gain-Sensitivity-Product (GSP). It is shown in Figure 9(b) and (d)that the HP and BP-R filters no. 3, i.e. the resistively tapered filters with equal resistors (having component values in the third row in Table 4) have the minimum sensitivity to passive component variations. The next best result is the 'optimum' design no. 5.To conclude, the sensitivity curves in Figure 9 confirm that complementary Biquads have identical optimum design, whereas dual Biquads have dual optimum designs. All complementary and dual Biquads in Figure 7have identical sensitivity figure of merit (all corresponding Schoeffler sensitivity curves in Figure 9 are equally high).

Figure 8.

Transfer function magnitudes of LP, HP and BP second-order filter examples [with (18) and K=1].

The output noise spectral density enodefined by square roof of (Eq. 3) has been calculated using Matlab (see Sections 3.3 and 3.4) and for these filters is shown in Figure 10. Note that there are only two figures; one for both the (complementary) LP and BP-C filters, i.e. Figure 10 (a), because they have identical noise properties, and the other for HP and BP-R filters, i.e. Figure 10(b). The total rms output noise voltage Enodefined by square root of (Eq. 4) are presented in the last columns of Table 3 and Table 4 (Jurisic et al., 2010 c).Considering the noise spectral density in Figure 10(a) and the Enocolumn in Table 3, we conclude that the LP and BP-C filters, with the lowest output noise and maximum dynamic range, are again filters no. 2. The second best results are obtained with filters no. 5, and these results are the same as those for minimum sensitivity shown above (see Figures 9a and c).

Figure 9.

Schoeffler sensitivities of second-order (a) LP, (c)BP-C filter examples inTable 3and (b) HP, (d) BP-R filter examples given inTable 4. (e) Legend.

Figure 10.

Output noise spectral densities of second-order (a) LP/BP-C and (b) HP/BP-R filter examples given inTables 3and4.

Analysis of the results in Figure 10(b) and the Enocolumn in Table 4 leads to conclusion that designs no. 3 and no. 5 of the HP and BP-R filters have best noise performance, as well as minimum sensitivity (see Figures 9b and d). The noise analysis above confirms that complementary circuits have identical noise properties and, on the other hand, those related by the RC–CR duality have different noise properties. Thus, there is a difference between LP and its dual counterpart HP filter in an output noise value. From inspection of Figure 10it results that the noise of the HP filter is larger than that of the LP filter, for all design examples.Consequently, we propose to use the LP and BP-C Biquads in Figure 7(a) and (c)as recommended second-order active filter building blocks, because they have better noise figure-of-merit, and the HP Biquad in Figure 7(b)as a second-order active filter building block for high-pass filters, if low noise and sensitivity properties are wanted. Unfortunately, it is unavoidable, that HP realizations will have a little bit worse noise performance.

4.2. Third-order Bitriplets

The extension to third-order filter sections follows precisely the same principles as those above. Unlike with second-order filters, third-order filters cannot be ideally tapered; instead only capacitive or resistive tapering is possible (Moschytz, 1999).Let us consider the third-order filter sections (Bitriplets) that realize LP and HP transfer functions, shown in Figure 11. Optimum design of those filters for low passive and active sensitivities, as well as low noise, is given in (Jurisic et al., 2010 b, Jurisic et al., 2010c). The optimum design is presented in Table 6 in (Jurisic et al., 2010 c)and is programmed using Matlab. In (Jurisic et al., 2010 a2010c), the detailed noise analysis on the analytical basis is given for the third-order LP and the (dual) HP circuits in Figure 11. Both sensitivity and noise analysis are performed using Matlab routines in Section 3. Voltage transfer functions for the filters in Figure 11are given by:


where numerators n(s) are given by:


Coefficients ai(i=0, 1, 2), and gain K as functions of filter components are given in Table 5.

Figure 11.

Third-order LP and HP active-RC filters with impedance scaling factors riand i(i=2, 3).

Coefficient(a) LP
Coefficient(b) HP

Table 5.

Transfer function coefficients of third-order active-RC filters with positive feedback in Figure 11.An optimization of both sensitivity and noise performance is possible by varying the general impedance scaling factors of the resistors and capacitors in the passive network of the filters in Figure 11, see (Moschytz, 1999):


The quantity referred to as 'design frequency' is defined by 0=1/(RC) (Moschytz, 1999).The third-order LP and HP filters with the minimum sensitivity to component tolerances as well as the lowest output noise and maximum dynamic range are the circuits designed in the optimum way as presented in Table 6 in (Jurisic et al., 2010 c). The LP filter circuit was designed by capacitive impedance tapering with 2=, 3=2; >1 and 0 chosen to provide r2r3. In the case of the third-order HP filter, the optimum design is dual: circuit has to be designed by resistive impedance tapering with r2=r, r3=r2; r>1 and 0 chosen to provide 23. Thus, the minimum-noise and minimum-sensitivity designs coincide. Comparing the output noise of two third-order dual circuits we see again that HP filter has larger noise than LP filter, c)although their sensitivities are identical, see (Jurisic et al., 2010 .

5. Conclusion

In this paper the application of Matlab analysis of active-RC filters performed regarding noise and sensitivity to component tolerances performance is demonstrated. All Matlab routines used in the analysis are presented. It is shown in (Jurisic et al., 2010 c) and repeated here that LP, BP and HP allpole active-RC filters of second- and third-order that are designed in (Jurisic et al., 2010 b) for minimum sensitivity to component tolerances, are also superior in terms of low output thermal noise when compared with standard designs. The filters are of low power because they use only one opamp. What is shown here is that the second-order, allpole, single-amplifier LP/HP filters with positive feedback, designed using capacitive/resistive impedance tapering in order to minimize sensitivity to component tolerances, also posses the minimum output thermal noise. The second-order BP-C filter with negative feedback is recommended filter block when the low noise is required. The same is shown for low-sensitivity, third-order, LP and HP filters of the same topology. Using low-noise opamps and metal-film small-valued resistors together with the proposed design method, low-sensitivity and low-noise filters result simultaneously. The mechanism by which the sensitivity to component tolerances of the LP, HP and BP allpole active-RC filters is reduced, also efficiently reduces the total noise at the filter output. Designs are presented in the form of optimum design tables programmed in Matlab [see Tables 1 and 6 in (Jurisic et al., 2010 c)].All curves are constructed by the presented Matlab code, and all calculations have been performed using Matlab.

© 2011 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike-3.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited and derivative works building on this content are distributed under the same license.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Dražen Jurišić (September 9th 2011). Low-Noise, Low-Sensitivity Active-RC Allpole Filters Using MATLAB Optimization, Applications of MATLAB in Science and Engineering, Tadeusz Micha&#322;owski, IntechOpen, DOI: 10.5772/19745. Available from:

chapter statistics

4428total 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

On Design of CIC Decimators

By Gordana Jovanovic Dolecek and Javier Diaz-Carmona

Related Book

First chapter

PID Control Design

By A.B. Campo

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