%首先來構造單位函數并對其性能進行驗證
syms n
Delta = sym('charfcn[0](n)');            %定義單位函數δ(n)
D0 = subs(Delta,n,0);                  %計算δ(0)
D15 = subs(Delta,n,15);                %計算δ(15)
disp('[D0,D15]');
disp([D0,D15])
syms z 
fn = 2*Delta+6*(1-(1/2)^n);
FZ = ztrans(fn,n,z);
disp('FZ = ');
pretty(FZ)
FZ_n = iztrans(FZ,z,n)