Certainly, <br>Try as in following :<br><br>num = [0 0.38 0.31]<br>den = [1 -0.86 0.55]<br><br>Hfilt      = poly(num,'z','c')/poly(den,'z','c');<br><br>[hfilt,fr] = frmag(Hfilt,1024);<br><br>
......... and then the story you want to build .................<br><br>Let me know if it helps.<br>Regards,<br><br>Sumit<br><br><br><br><br><div class="gmail_quote">On Fri, Feb 26, 2010 at 4:04 PM, David PICARD <span dir="ltr"><<a href="mailto:d.picard@opgc.univ-bpclermont.fr">d.picard@opgc.univ-bpclermont.fr</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Ooops, sorry, you're right. I mistyped the denominator (0.86z^-1 instead of -0.86z^-1) . Of course, I was confused by the weird step response...<br>

<br>
Is it possible to specify the transfer function with polynoms ? In this case, it would be something like :<br>
num = [0 0.38 0.31]<br>
den = [1 -0.86 0.55]<br>
<br>
Regards.<br>
<br>
<br>
Collette yann a écrit :<div><div></div><div class="h5"><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
For me it corresponds to the good transfert function:<br>
<br>
    0.31 + 0.38z                    z^2* (0.31*z^-2 + 0.38*z^-1)<br>
  ----------------         =       ------------------------------------<br>
  0.55 + 0.86z + z ^2            z^2*(0.55*z^-2 + 0.86z^-1 + 1)<br>
<br>
Scilab seems to have simplified in a way your transfert function<br>
<br>
YC<br>
<br>
<br>
David PICARD a écrit :<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
I need some help to define the following transfer function :<br>
<br>
 0.381417 z^-1 + 0.310306 z^-2<br>
---------------------------------<br>
1 - 0.857089 z^-1 + 0.548812 z^-2<br>
<br>
This is a 2nd order system : w0=1, Damp=1.3, Gain=1.0, Te=1s<br>
<br>
I tried :<br>
-->z=%z;<br>
-->sys=syslin(1.0,(0.38*z^-1+0.31*z^-2)/(1+0.86*z^-1+0.55*z^-2))<br>
sys  =<br>
<br>
     0.31 + 0.38z<br>
   ----------------<br>
                   2<br>
   0.55 + 0.86z + z<br>
<br>
<br>
But it NOT what I asked !!! Help !<br>
<br>
It is so straightforward with GNU Octave, that I am really frustrated !<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
sys=tf([0.000000 0.381417 0.310306],[1.000000 -0.857089 0.548812],1.0);<br>
step(sys)   // shows the oscillating, stabilizing step response<br>
</blockquote>
<br>
<br>
<br>
</blockquote>
<br>
<br>
<br>
</blockquote>
<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><br>--------------------------------------------<br>Sumit Adhikari<br>System Design Engineer<br>austriamicrosystems AG<br>Business Unit : Automotive<br>Mob : 00-91-9885271710/00-91-9000161710<br>