如何在Scilab上获得简化分数？

Question

我正在使用函数F（s）= 1 / s5 + 5 * s4 + 9 * s3 + 7 * s2 + 2 * s，我想得到各个分数来应用拉普拉斯。

s= %s
n2=1
d2=s^5+5*s^4+9*s^3+7*s^2+2*s
sist2=syslin('c',n2/d2)
sist2desc=pfss(sist2)


sist2desc  =

   sist2desc(1)

0.5   
---   
 s    

   sist2desc(2)

 0.5    
-----   
2 + s   

   sist2desc(3)

            2     
- 2 - 2s - s      
--------------    
           2   3  
1 + 3s + 3s + s   

这是我用该代码得到的输出，但我希望得到0.5 / s + 0.5 /（s + 2）-1 / s + 1 -1 /（s + 1）

Answer 1

对于Scilab来说，你的分数很难，当处理多项式和分数时，它使用浮点算法，因为多个根降低了经典根寻找方法（伴随矩阵的特征值）的精度。在这里你的分母有-1作为三重根，三重多样性已经太多，如图所示

--> roots((s+1)^3*(s+2)*s)
 ans  =

  -2.  
  -1.0000109  
  -0.9999945 + 0.0000095i
  -0.9999945 - 0.0000095i
   0.  

有关替代根查找方法，请参阅讨论@ http://bugzilla.scilab.org/show_bug.cgi?id=15349。这种行为表明，使用Scilab这样的学术内容是一个坏主意。它可以用于线性代数，但不适用于符号微积分。