[Scilab-users] Limit command

[duce]

z = %z;
syms n z1;
z = z1
X =z *(3*z -(5/6) ) /((z -(1/4) )*(z -(1/3) ))
X1 = denom (X);
zp = roots (X1);
X1 = z1 *(3* z1 -(5/6) ) /(( z1 -(1/4) )*(z1 -(1/3) ))
F1 = X1 *( z1 ^(n -1) )*(z1 -zp (1) );
F2 = X1 *( z1 ^(n -1) )*(z1 -zp (2) );
h1 = limit (F1 ,z1 ,zp (1) );
disp (h1 , ’ h1 [ n]= ’ )
h2 = limit (F2 ,z1 ,zp (2) );
disp (h2 , ’ h2 [ n]= ’ )
h = h1+h2;
disp (h, ’ h [ n]= ’ )
// R e s u l t
//h [ n]= ( 1 / 4 ) ˆn +(2/3) ˆn

This is the original script. It is used to find the inverse z-transform of
the expression given as *X* . The example is from a textbook by Alan V.
Oppenheim, Alan V. Willsky called "Scilab Code for Signals and Systems" so
I'm guessing it's not a self-defined function. However I'm not sure if it's
from an atom module.



--
View this message in context: http://mailinglists.scilab.org/Limit-command-tp4032725p4032732.html
Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com.