二元系统的Barkatou-LeRoux-Levelt秩约简算法(ISSAC’2006) 

Levet算法(1990年)的Barkatou-LeRoux泛化(ISSAC’2006)计算了与输入系统等效的系统,其子系统分别和同时具有真正的Poincare秩。输出由等效系统的矩阵、变换和计算时间组成。
函数调用:
 

 

LeveltFaff(输入矩阵,其中极点位于两个变量或两个变量中的任何一个,第一个变量、第二个变量、第一个变量的截断(如果是幂级数条目)、第二变量的截断)。 

 

LeveltPfaffM(输入矩阵,两个变量或两个变量中的任何一个都有极点,第一个变量、第二个变量、第一个变量的截断(如果是幂级数条目)、第二变量的截断)。LeveltFaffM使用Moser约简准则来避免不必要的计算。 

 

在以下示例中,输入系统为: 

dF/dx=A(x,y)F 

dF/dy=B(x,y)F 

 

> A:=矩阵(2,2,[`/`(`*`(`+`(`**`(`^`(x,3)),y)),`*` y)),`*`(`^`(x,4))]);1; ...
A:=矩阵(2,2,[`/`(`*`(`+`(`**`(`^`(x,3)),y)),`*` y)),`*`(`^`(x,4))]);1; ...
 

 

 

A:=矩阵(%id=18446744078265552638)
B:=矩阵(%id=18446744078265553134)
等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078265554574),矩阵 (2.1)
 

> A:=矩阵(%id=18446744078265556614)-1; B:=矩阵(%id=18446744078265556734)-1; 等价_A,等价-B,转换,计算时间:=级别Faff(A,B,x,y,10,10);1
A:=矩阵(%id=18446744078265556614)-1; B:=矩阵(%id=18446744078265556734)-1; 等价_A,等价-B,转换,计算时间:=级别Faff(A,B,x,y,10,10);1
 

等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078265566382),矩阵 (2.2)
 

>
 

> A:=矩阵(%id=18446744078265568062)-1; B:=矩阵(%id=18446744078265568182)-1; 等价_A,等价-B,转换,计算时间:=级别Faff(A,B,x,y,10,10);1
A:=矩阵(%id=18446744078265568062)-1; B:=矩阵(%id=18446744078265568182)-1; 等价_A,等价-B,转换,计算时间:=级别Faff(A,B,x,y,10,10);1
 

等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078265569022),矩阵(%id=18446744078 265569142),矩阵 (2.3)
 

> A:=矩阵(%id=18446744078272776182)-1; B:=矩阵(%id=18446744078272776302)-1; 等价_A,等价-B,转换,计算时间:=级别Faff(A,B,x,y,10,10);1
A:=矩阵(%id=1844674407822776182)-1; B:=矩阵(%id=18446744078272776302)-1; 等价_A,等价-B,转换,计算时间:=级别Faff(A,B,x,y,10,10);1
 

等值_A,等值_B,变换,Computation_time:=矩阵(%id=1844674407822777022),矩阵(%id=1844674407822777142),矩阵(%id=18446744078227772762),0.29e-1 (2.4)
 

> A:=矩阵(%id=18446744078272778958)-1; B:=矩阵(%id=18446744078272779078)-1; 等价_A,等价-B,转换,计算时间:=级别Faff(A,B,x,y,10,10);1
A:=矩阵(%id=18446744078272778958)-1; B:=矩阵(%id=18446744078272779078)-1; 等价_A,等价-B,转换,计算时间:=级别Faff(A,B,x,y,10,10);1
 

等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272780398),矩阵 (2.5)
 

> A:=矩阵(%id=18446744078272783174)-1; B:=矩阵(%id=18446744078272783294)-1; 等价_A,等价-B,转换,计算时间:=级别Faff(A,B,x,y,10,10);1; 等效_。。。
A:=矩阵(%id=18446744078272783174)-1; B:=矩阵(%id=18446744078227782394)-1; 等价_A,等价-B,转换,计算时间:=级别Faff(A,B,x,y,10,10);1; 等效_。。。
 

 

等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272785814),矩阵
等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272785814),矩阵
等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272796782
等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272796782
(2.6)
 

>