%首先來(lái)構(gòu)造單位函數(shù)并對(duì)其性能進(jìn)行驗(yàn)證
syms n
Delta = sym('charfcn[0](n)');            %定義單位函數(shù)δ(n)
D0 = subs(Delta,n,0);                  %計(jì)算δ(0)
D15 = subs(Delta,n,15);                %計(jì)算δ(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)