Я бы с самого начала сделал чуть-чуть иначе:


где

- резонансная частота,

- затухание колебательного контура,

- добротность,

- нормированная частота.
Для АЧХ

В этом виде удобно строить график - мы получили семейство кривых с параметром

. Чтобы провести прикидочный анализ поведения кривых последнее выражение перепишем в виде:

Теперь смотрите: кроме анализа ситуации в нуле и на бесконечности, обращает на себя внимание, что в знаменателе под корнем у нас парабола относительно

, ветви этой пораболы направлены вверх, то есть где-то есть локальный минимум. Координата минимума пораболы по известной формуле

. Если

, этот минимум отбрасываем, так как

. Если

, то минимум учитываем, а при малых затуханиях

(больших добротностях

)

, то есть минимум приходится аккурат на резонансную частоту.
Между тем числитель АЧХ монотонен. Следовательно на частоте

, при

имеет место локальный максимум. При больших добротностях частота этого максимума близка к

.
Физически это означает существование в рассматриваемой цепи такого частотного режима, в котором в стационарном режиме при гармоническом воздействии наблюдается возрастание амплитуды гармонического сигнала на выходе цепи при приближении частоты к определённому значению. Такой частотный режим называется резонансом.
Вообще это хорошо, что Вы не доверяете маткаду. Но не доверяйте тогда как следует. А следовало бы вам рассмотреть и производную вашей функции и проанализировать её на предмет наличия корней, что указало бы вам на наличие локальных экстремумов. Иначе было бы очень просто: рассмотрел нуль и бесконечность и построил график. Но всё вокруг как будто бы против нас и приходится думать.