%求Z變換和反變換
syms a n z t
Fz1=ztrans(sym('Heaviside(t)'),n,z)  %求階躍函數的Z變換
Fz2=ztrans(sym('Dirac(t)'),n,z)  %求脈沖函數的Z變換
Fz3=ztrans(exp(-a*t),n,z)  %求e-at的Z變換
syms n z t
f1=iztrans(Fz1,z,n)
f2=iztrans(Fz2,z,n)
f3=iztrans(Fz3,z,n)