using TySymbolicMath
@variables n,z
y = (1//2)^n;
Zy = ztrans(y) #这里返回y的z变换，结果里有z的收敛域，如何提取收敛域内关于z的分式？
Y = 1 / (1 + (-1//2) / z)
iZY = iztrans(Y) #这里为什么会返回NaN？如何让有理分式的z逆变换返回 指数+阶跃+斜坡等等 的形式？
# iZZy = iztrans(Zy) #这一句会报错，为什么？
Zy, iZY

(
  TySymbolicMath.PieceWise(TySymbolicMath.ExprCondPair(1 / (1 + (-1//2) / z), (1 / abs(z)) < 2), TySymbolicMath.ExprCondPair(TySymbolicMath.SymSum((z^(-n))*(2^(-n)), n, 0, Inf), true)),
  nan + TySymbolicMath.kroneckerDelta(n, 0)
)