Electronic – How to go about determining the component values of the described filter? I know what they are, but not how to get them using proper formulas

component-valuesfilteroperational-amplifier

The filter in question takes a mono signal of amplitude between 5 mVp and 200 mVp, and filters the signal based on the frequency. The potentiometer that leads into the inverting input of the op amp determines how this signal is ultimately affected. Depending on which extreme the potentiometer is dialed to, the circuit should amplify either low or high input frequencies, while attenuating the other end of the frequency extrema. When the potentiometer is dialed to the opposite extreme however, its behavior inverts; that is to say that if it amplified high frequencies by a voltage gain of A1 and attenuated low frequencies by A2, it should now amplify the low frequencies by A1 and attenuate the high frequencies by A2. When the potentiometer is set to exactly 50%, the gain should be -1, as no filtering is actually being done.

For this specific purpose, I need A1 to be 3 and A2 to be 1/3, i.e. A2 = 1/A1. I have values that I obtained somewhat coincidentally that achieve this behavior, which are shown in the circuit below. My question, finally arises from my lack of understanding of how to actually approach this circuit. It's almost like a band-pass filter, but not quite as one of the ends will amplify the voltage by the same amount it would have been attenuated. That's pretty much where my confusion lies. In this schematic, the circuit is designed to work towards the gains previously described, and is built to handle frequencies ranging from 20 Hz to 20 kHz.

schematic

simulate this circuit – Schematic created using CircuitLab

I know that:
C1 = C2
R1 = R7
R2 = R6
Due to the gain of the circuit being the impedance of the right hand divided by the impedance of the left, and at 50% of the potentiometer, this must equal 1. If the potentiometer is split 10k and 10k for both impedances, then in order for it to equal -1 the above relationships must hold true.

R3 is always a 20k potentiometer. However, as far as the other components go, they can vary so long as they fall in accordance with those relationships. My issue lies in that I'm not quite sure what to do next to get the answers I came up with. I think I maybe have to compare how the expression for gain changes at each extreme with the expression for a gain of -1, but I'm not sure. If anyone can walk me through this, as well as correct any incorrect assumptions I've made, it would be a great help.

Best Answer

I did not realize this circuit would be so complicated to analyse : ) when I started to look at it using the FACTs. First off, I redrawn the circuit in which the inputs were no longer floating but ground-referenced. I also considered the open-loop gain \$A_{OL}\$ to later push it to infinity. The newly-redrawn circuit is below. I have split the potentiometer in two distinct resistances linked by a ratio \$k\$ and reflected in the ground-referenced source driving the op amp:

enter image description here

Both dc and ac responses are similar. The next step is to find the dc transfer function for \$s=0\$: open the capacitors and determine the gain in this configuration:

enter image description here

If you do the maths ok, you should get \$H_0=-\frac{R_3(1-k_1)+R_7}{R_1+k_1R_3}\$

For the next steps, I will determine the time constants \$\tau\$ involving capacitors \$C_1\$ and \$C_3\$. You need to reduce the excitation voltage to 0 V and replace it by a short circuit. You obtain \$\tau_1\$ and \$\tau_3\$. Final step, determine the time constant involving \$C_3\$ while \$C_1\$ is replaced by a short circuit. You obtain \$\tau_{13}\$. Assemble these time constants to form the denominator:

\$D(s)=1+s(\tau_1+\tau_3)+s^2\tau_1\tau_{13}=1+\frac{s}{\omega_0Q}+(\frac{s}{\omega_0})^2\$

Then, calculate and check if the quality factor \$Q\$ is well below 1 to apply the low-\$Q\$ approximation. You can factor \$D\$ as two cascaded poles if \$Q<<1\$: \$D(s)=(1+\frac{s}{\omega_{p1}})(1+\frac{s}{\omega_{p2}})\$ in which \$\omega_{p1}=Q\omega_0\$ and \$\omega_{p2}=\frac{\omega_0}{Q}\$. Unfortunately here, Q is around 0.47 (\$k=0.6\$) so the poles are not far from each other and the low-\$Q\$ approximation does not hold. The resonant frequency in \$D(s)\$ changes from 492 Hz (\$k=0.1\$) to 1160 Hz (\$k=0.9\$).

The numerator can be determined using several approaches: inspect the transformed circuit and see if any impedance combination could lead to a null in \$V_{out}\$ despite the presence of \$V_{in}\$. You could for example check the condition for which \$\epsilon=0\$ and solve the equation. The other option is to apply an Null Double Injection (NDI) and determine the new time constants (as we did for \$D(s)\$) but with \$V_{in}\$ back in place and while \$V_{out}=0\$. These two paths give the simplest expression for \$N(s)\$. I have adopted the generalized form in which I will calculate three gains \$H\$ when the capacitors are alternatively replaced by a short circuit (for \$a_1\$) or when both are shorted (\$a_2\$). This gives a slightly more complex expression than the two other approaches and you will need to rework the final result to simplify it. For a dynamic response analysis like in here, the generalized formula is okay:

\$N(s)=H_0+s(H^1\tau_1+H^2\tau_2)+s^2H^{13}\tau_1\tau_{13}\$

You can also rearrange this expression in a second-order polynomial form as we did for the denominator. Then assemble the transfer function as \$H(s)=H_0\frac{N(s)}{D(s)}\$. The expressions are given in the below Mathcad files:

enter image description here

enter image description here

I have tweaked these equations so that they depend on \$k\$ as a variable to pass and obtained the following responses:

Magnitude:

enter image description here

Phase:

enter image description here

I did not do it but I believe you can obtain the frequency at which all the curves cross by differentiating the magnitude expression with respect to \$k\$ and solve the value of \$\omega\$ which cancels the expression. You actually calculate the function sensitivity to \$k\$ and check the frequency at which it is 0.

The cool thing with the FACTs is that you do not need to really bother about the circuitry itself: determine the time constants in certain conditions (\$V_{in}=0\$ for \$D\$ and (\$V_{out}=0\$ for \$N\$) and assemble them as recommended to obtain the final expression. If you are interested by the technique, you can check a seminar taught at APEC in 2016. It is an introduction to show the power of the FACTs and to apply them to simple passive and active circuits.