%-----------------------------------------------------------------
% exa020807, for example 2.8.7
% to test tf2sos,sos2tf,zp2sos and sos2zp;
%-----------------------------------------------------------------
clear; 

B=[0.0201 0 -0.0402 0 0.0201];
A=[1 -1.637 2.237 -1.307 0.641];
%
[sos,G]=tf2sos(B,A)
% 將有理分式（轉移函數）分解為二階子系統的級聯；
%
[B,A]=sos2tf(sos,G)
%由級聯的二階子系統重構出轉移函數；
%
[Z,P,K]=tf2zp(B,A)
%
%
[sos,G]=zp2sos(Z,P,K)
% 將一轉移函數的極零點轉換為二階子系統的級聯；
%
[Z,P,K]=sos2zp(sos,G)
%由級聯的二階子系統重構出轉移函數的極零點