Signals and Systems/Second Order Transfer Function

The second order transfer function is the simplest one having complex poles.
Its analysis allows to recapitulate the information gathered about analog filter design and serves as a good starting point for the realization of chain of second order sections filters.
All-Pole Second Order Transfer Function.
Transfer Function.
The transfer function of a continuous-time all-pole second order system is:
Note that the coefficient of formula_2 has been set to 1.
This simplifies the writing without any loss of generality, as numerator and denominator can be multiplied or divided by the same factor.
The frequency response, taken for formula_3, has a DC amplitude of: 
For very high frequencies, the most important term of the denominator is formula_5 and the frequency response gets closer and closer to:
At high frequencies, the amplitude response looks like a (squared) hyperbol in a linear plot and like a straight line with a negative slope in a log-log plot.
Plotting the frequencies in decades and the amplitude in decibels reveals a slope of -40 [dB/decade].
Having a given amplitude at DC and an amplitude nearing zero at high frequencies indicates that the transfer function is of lowpass type.
Poles.
The poles of the system are given by the roots of the denominator polynomial:
If the term inside the square root is negative, then the poles are complex conjugates.
This is the general case in filter design: there is poor interest in a second order transfer function having two real poles.
From the location of the poles, the transfer function can be rewritten as:
The amplitude of the poles gives the corner frequency of the filter.
The corner frequency is defined as the abscissa of the point where the horizontal and the -40 [dB/decade] lines meet in the log-log magnitude response plot.
Note that this is not necessarily the -3 [dB] attenuation frequency of the filter.
The ratio between the real part of the poles and the corner frequency is proportional to the damping, or inversely proportional to the quality factor of the system.
The closer the poles are to the imaginary axis, the more a resonance will appear at a frequency smaller but close to the corner frequency of the system.
Transfer Function Rewritten.
With this, the transfer function with unity gain at DC can be rewritten as a function of the corner frequency and the damping in the form:
Both formula_13 and formula_14 have a unit of [s-1].
formula_15 has a unit of [1] and so does the total transfer function.
Biquadratic Second Order Transfer Function.
Transfer Function.
The name biquadratic stems from the fact that the functions has two second order polynomials:
Poles.
The poles are analysed in the same way as for an all-pole second order transfer function.
They determine the corner frequency and the quality factor of the system.
Zeroes.
The zeroes are used to affect the shape of the amplitude response:
Example.
6th Order Normalized Butterworth Filter.
 This example analyses the 3 second order sections of a 6th order normalized Butterworth filter.
The poles of the Butterworth filter are regularly spaced on the left half of a circle centered at the origin of the complex plane.
Each complex conjugate pole pair builds a second order all-pole transfer function.
The product of these second order functions gives the 6th order Butterworth transfer function.
Pole Pairs.
In the figure on the side, the pole formula_21 and its complex conjugate are close to the imaginary axis.
Their amplitude response will show an overshoot at the corner frequency.
This corresponds to an underdamped case and the second order section will show some resonance at frequencies close to the corner frequency.
The pole formula_22 and its complex conjugate are at 45° in respect to the imaginary axis.
Their amplitude response will show 3 dB loss at the corner frequency.
It is the limiting case where the amplitude response shows no overshoot.
Nevertheless, this doesn't correspond to a critically damped case: the step response will have overshoots before stabilization.
The pole formula_23 and its complex conjugate are far away from the imaginary axis.
Their amplitude response will show a large attenuation at the corner frequency.
This corresponds to an overdamped case.
Amplitude Responses.
The following Octave code allows to plot the amplitude responses of the individual second order sections and of the global Butterworth amplitude response:
 # Specifications
 order = 6;
 pointNb = 100;
 wLog = logspace(-1, 1, pointNb);
 AdbMax = 10;
 AdbMin = -40;
 #-------------------------------------------------------------------------------
 # Butterworth filter
 [num, den] = butter(order, 1, 's');
 while ( length(num) &lt; length(den) )
 num = [0, num];
 endwhile;
 #-------------------------------------------------------------------------------
 # Normalized filter, transfer functions plot
 Adb = 20*log10(abs(freqs(num, den, wLog)));
 figure(2);
 semilogx(wLog, Adb);
 hold on;
 for index = 1:order
 if imag(poles(index)) &gt; 0
 num = [0, 0, 1];
 den = [1, -2*real(poles(index)), abs(poles(index))^2];
 Adb = 20*log10(abs(freqs(num, den, wLog)));
 semilogx(wLog, Adb, 'g');
 endif
 endfor
 hold off;
 axis([wLog(1), wLog(length(wLog)), AdbMin, AdbMax]);
 xlabel('Frequency');
 ylabel('Amplitude responses');
 grid;
The blue curve on the side shows the global amplitude response.
The green curves are the responses of the individual second order sections.
They all have a hozizontal asymptote towards DC.
They also all have a -40 dB/decade asymptote for high frequencies.
Both asymptotes cross at the point (formula_24 [s-1], formula_25 [dB]).
The corner frequency is found at formula_24 [s-1] or formula_27 [Hz]
High-Q Second Order Function.
The top green amplitude response shows what a response with a high quality factor looks like.
It has a maximum of more than 0 dB (here 6.02 dB) at a frequency a little below the corner frequency.
At the corner frequency, the amplitude has already fallen down (here to 5.68 dB).
Maximally Flat Second Order Function.
The middle green amplitude response shows what a maximally flat response looks like.
It has an amplitude of -3.02 dB at the corner frequency.
Whatever its order, a Butterworth function shows the same -3.02 dB loss at the corner frequency.
Low-Q Second Order Function.
The bottom green amplitude response shows what a response with a low quality factor looks like.
It has an amplitude of less than -3 dB (here -5.72 dB) at the corner frequency.