Page 86 - 363_

87

om0 = le-2; omk = le2.
Тепер треба задатися кількістю точок на майбутньому графіку n = 41, і

сформувати масив точок за частотою:
ОМ = logspace(-2, 2, 41),

де значення -2 і +2 відповідають десятковим порядкам початкового om0 і
кінцевого omk значень частоти.

Користуючись функцією polyval, можна обчислити спочатку вектор

комплексних значень “ch” чисельника частотної функції передачі, яка
відповідає заданій функції передачі, якщо як аргумент функції polyval

використати сформований вектор частот ОМ, елементи якого помножені на

уявну одиницю i. Аналогічно обчислюється комплексно значний вектор “zn”
знаменника ЧФП.

Вектор значень АЧХ (амплітудно-частотної характеристики) можна

знайти, обчислюючи модулі векторів чисельника і знаменника ЧФП і ділячи
поелементно одержані вектори. Щоб знайти вектор значень ФЧХ (фазо-

частотноі характеристики) треба поділити поелементно коммплексно-значні
вектори чисельника і знаменника ЧФП і визначити вектор аргументів елементів

одержаного вектора. Для того, фазу подати у градусах, одержані результати

слід помножити на 180 і поділити на .
Нарешті, для побудови графіка АЧХ у логарифмічному масштабі,

достатньо застосувати функцію loglog, а для побудови ФЧХ зручніше

користуватися функцією semilogx.

У цілому послідовність дій може бути такою:
>> ОМ = logspace(-2 ,2 ,40)
>> ch = polyval(Pl, i*OM);
>> zn = polyval(P2, i*OM);
>> ACH = abs(ch). /abs(zn);
>> loglog(OM,ACH); grid; title('Гpaфік Амплітудно-Частотної Хар-кі')
>> xlabel('Частота (рад/с'); уlabel('Відношення амплітуд')
>> FCH = angle(ch. /zn)*180 / pi;
>> semilogx(OM,FCH); grid, title('Фазо-Частотна Характеристика')
>> xlabel('Частота (рад/с) '), уlabel('Фаза (градуси) ')

В результаті одержуються графіки, наведені на рисунках 8.10 і 8.11.