Km=19.2;Tm=1.32;Kp=11.52;                                %參數來自圖4.2.6
Ti=0.66;Td=0.1650;                                        %例4.2.2計算得出的參數
K=0.6667;L=0.293;T=1.947;                               %參數來自圖4.2.3
theta=L/T;lambda=Km*K;beta=36/(27+5*lambda);          %(4.2.1)式和(4.2.10)式
num=10;den=conv([1,1],conv([1,3],[1,5]));G=tf(num,den);        %被控對象
num1=Kp*[Td,(1-beta)];den1=1;G1=tf(num1,den1);        %圖4.2.8并連反饋環節
GC=feedback(G,G1);                                      %圖4.2.8并連反饋傳遞函數
num2=Kp*[Ti*beta,1];den2=[Ti,0];G2=tf(num2,den2);     %圖4.2.8串連反饋環節
GG=G2*GC;GCG=feedback(GG,1);                     %GG系統開環；GCG閉環傳遞函數
theta,lambda,Kp,Ti,Td,                                     %參數顯示
step(GCG)                                               %系統閉環單位階躍響應