#Scilab state space to transfer function code#

Here is a Scilab implementation of Algorithm 3: function =bal(Ain) In the algorithm the operations are quite simple (computing inf and 2-norms, swaping columns or rows of a matrix), maybe this could easily translated ?Ī more readable version of the algorithm can be found on page 3 (Algorithm 2) of the following document.

In Scilab balanc() is hard-coded and based on LAPACK's dgebal (see the Fortran source at Netlib). Thanks in advance for the help in doing this. Note: I've tried to make the calculations in wxMaxima to have Ab and X matrices as elements with variables but the system of equations remains with too many variables and couldn't be solved.

#Scilab state space to transfer function code#

Most of my code in wxMaxima to reach in the near future has been done by translating the Scilab code into wxMaxima, currently I'm writing the tf2ss() function an inside that function the balanc() function is called, the problem is that I couldn't find the code for that function in Scilab installation directory, I've searched info in books and papers but every example starts with the Ab matrix given as an input to the problem, Scilab instead has the option to have as an input only the A matrix and it calculates the Ab and X matrices, so, I need help to make this function exactly as Scilab has it programmed to been able to compare all the steps that I'm doing.įinally, wxMaxima has a function to calculate similarity transforms but it don't have the same output as Scilab what it means to me that they uses different criteria to calculate the similarity transform. XĪs approximately equal row and column norms. Right now I'm stranded with the function balanc() which finds a similarity transform such that Ab = X^(-1). Scilab to calculate the z-transform first converts the transfer function to the state space, after that the system must be discretized and after that converted to z transfer function, I need this because of some algebraic calculations that I need to do to analyze stability of a system in function of the sample period. I'm trying to program the z-transform in wxMaxima which doesn't have it programmed but not by definition but by using the Scilab approach.