%exp4_12.m
%margin函數通常用在bode函數之后，先由bode函數得到
%幅值、相角和頻率矢量，然后由margin繪制出幅值裕度
%和相角裕度的波特圖。
clear
clc
close all
%某系統的開環傳遞函數為：G(s)=k/s(s+1)(0.2s+1)
%求k分別為2和20時的幅值裕度與相角裕度
num1=2;num2=20;
den=conv([1 0],conv([1 1],[0.2 1]));
w=logspace(-1,2,100);
figure(1)
[mag1,pha1]=bode(num1,den,w);
margin(mag1,pha1,w)
figure(2)
[mag2,pha2]=bode(num2,den,w);
margin(mag2,pha2,w)