matlab ss2tf函数 C实现

继上一篇matlab zp2ss函数 C实现

这一篇实现ss2tf(将状态空间系统转换为传递函数形式)

代码

void ss2tf(double *a, double *b, double *c, double *d, int n, double *num, double *den)
{
   
   
    poly(a, n, den);

    double *temp_a = malloc(sizeof(double) * n * n);

    trmul(b, c, n, 1, n, temp_a);
    
    for (size_t i = 0; i < n; i++)
    {
   
   
        for (size_t j = 0; j < n; j++)
        {
   
   
            temp_a[i * n + j] = a[i * n + j] - temp_a[i * n + j];
        }
    }

    poly(temp_a, n, num);