Там ничего хитрого, все через комплексный коэффициент передачи и делается. У Вас есть все ровно то же самое, только график построить остается. Строить те графики что в логарифмических единицах нужно тоже в логарифмическом масштабе, соответственно.
w = x / (rs + x) !Комплексный коэффициент передачи
afr = abs (w) !АЧХ
pfr = atan2 (aimag(w),real(w)) !ФЧХ
lafr=20*log10(afr) !Логарифмическая АЧХ
3-я строчка вычисляет аргумент комплексного числа (угол) с учетом правильного выбора квадранта, по сути это тангенс от мнимой и вещественной частей числа, по ней строится ФЧХ. Ну а АФЧХ строится так - берем модуль и аргумент (3-я и 4-я строчки) и строим в полярных координатах. Правильно по линии конечно наносить и подписывать точки с указанием частот, что бы понятно было. Тогда по кривой будет сразу видно какая частота например у самой правой точки (с максимальным радиусом круга), амплитуда максимальна, сдвига угла нет, это резонанс. Кривая начинается из нуля (минимальная частота), идет вверх и обходит по кругу, показывая угол и амплитуду одновременно.