%-----------------------------------------------------------------
% exa020806, for example 2.8.6
% to test tf2zp,zp2tf,sort,and poly;
%-----------------------------------------------------------------
clear; 

B=[1 -3.3 7.25 -6.7 3 -0.8];
L=length(B);
A=zeros(1,L);
A(1)=1;   % 得到  H(z)=B(z)/A(z),A(z)=1;
%
[Z,P,K]=tf2zp(B,A)                 
% 求H(z)的極零點，即實現由有理多項式到極零點的轉換;
% 
sort(Z)                            
%將零點按模的大小，以遞減的方式排序；

[b,a]=zp2tf(Z,P,K)   
% 實現由極零點到有理多項式的轉換; 
%
Z1=roots(B)
% 求多項式的根；
poly(Z1)
% 由多項式的根重構該多項式