Plot(real(W),imag(W))

Далее с помощью пункта меню Insert →X Label → Label Y→ Title окна Figure № вводятся обозначения осей и заголовок рисунка. Общий вид изображенияамплитудно-фазовой характеристики показан на рис. 10.

Рис. 10. Изображение амплитудно-фазовой характеристики

Щелкнув правой кнопкой мыши на поле рисунка, можно вызвать всплывающее меню, с помощью которого вводится координатная сетка (пункт Grid). С помощью пункта меню View→Plot Browser можно изменить тип, форму и цвет линий графика.

Построение расчетных логарифмических амплитудной LР(w) и фазовой arР(w) частотных характеристик производится в соответствии со следующими ниже правилами.

1. Характеристики строятся в полулогарифмических координатах (ось Х – логарифмическая, ось У – естественная) с использованием функции semilogx, у которой правила ввода аргументов и значений спецификатора кривой аналогичны функции plot.

2. В силу особенностей функции MATHLAB atan(x) = arctg(x) (выдаются значения только в диапазоне от -π/2 до π/2), функция ar(w), набирается в виде двух выражений (см. текст последовательности вычислений ниже).

3. Каждая характеристика набирается отдельным графиком, и ввод заканчивается командой hold on, обеспечивающей вывод последующего графика в то же текущее окно рисунка.

4. Из набранных осевых надписей MATHLAB выводит только последние (ar2 и log(w2)). Для вывода всех необходимых надписей (как показано на рис.1.4) следует:

–в окне Figure No.1 командой меню Edit – Axes Properties вызвать окно редактора заголовка и осевых надписей (Property Editor) и ввести необходимые тексты; кнопки Edit рядом с каждым оконцем позволяют менять у текста тип шрифта и его параметры;

–нажатием кнопки «А» и щелчком левой кнопки мышки ввести текстовое поле (с последующим вводом текста) для пометки выбранной кривой (стрелка вставляется нажатие кнопки «Стрелка» с последующим вычерчиванием линии u1083 левой кнопкой мышки).

Характеристики набираются в командном окне согласно формул (10) и (11) в соответствии с цепочкой команд:

>> k=12;

z=0.1;

T=1;

w=0:0.01:100;

W=k./(T.^2.*(w.*j).^2+2.*T.*z.*w.*j+1);

L=20.*log10(k./sqrt((2.*T.*z.*w).^2+(1-T.^2.*w.^2).^2));

semilogx(w,L);

hold on;

w1=0.1:0.01:1;

ar1=-(180./pi).*atan(2.*T.*z.*w1./(1-T.^2.*w1.^2));

semilogx(w1,ar1);