基于Moser-Barkatou算法推广的二元系统秩约简算法(ISSAC’2014) 

MoserPfaff计算输入系统的等效系统,其子系统分别和同时具有真实的Poincare秩。输出由等效系统的矩阵、变换和计算时间组成。 


函数调用:
 

 

MoserPfaff(输入矩阵,两个变量或两个变量中的任何一个都有极点,第一个变量,第二个变量,在第一个变量中截断(如果是幂级数条目),在第二个参数中截断)。 

 

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

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=18446744078272799318)
B:=矩阵(%id=18446744078272799798)
等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272801238),矩阵 (3.1)
 

> A:=矩阵(%id=18446744078272803294)-1; B:=矩阵(%id=18446744078272803414)-1; 等价_A,等价-B,变换,计算时间:=MoserPfaff(A,B,x,y,10,10);1
A:=矩阵(%id=18446744078272803294)-1; B:=矩阵(%id=18446744078272803414)-1; 等价_A,等价-B,变换,计算时间:=MoserPfaff(A,B,x,y,10,10);1
 

等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272804854),矩阵(%id=18446744078 272804974),矩阵 (3.2)
 

>
 

> A:=矩阵(%id=18446744078272806534)-1; B:=矩阵(%id=18446744078272806654)-1; 等价_A,等价-B,变换,计算时间:=MoserPfaff(A,B,x,y,10,10);1
A:=矩阵(%id=18446744078272806534)-1; B:=矩阵(%id=18446744078272806654)-1; 等值_A,等值_B,变换,计算时间:=MoserPfaff(A,B,x,y,10,10);1
 

等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272807510),矩阵(%id=18446744078 272807630),矩阵 (3.3)
 

> A:=矩阵(%id=18446744078272809790)-1; B:=矩阵(%id=18446744078272809910)-1; 等值_A,等值_B,变换,计算时间:=MoserPfaff(A,B,x,y,10,10);1
A:=矩阵(%id=18446744078272809790)-1; B:=矩阵(%id=18446744078272809910)-1; 等值_A,等值_B,变换,计算时间:=MoserPfaff(A,B,x,y,10,10);1
 

等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272810630),矩阵(%id=18446744078 272810750),矩阵 (3.4)
 

> A:=矩阵(%id=18446744078272833046)-1; B:=矩阵(%id=18446744078272833166)-1; 等价_A,等价-B,变换,计算时间:=MoserPfaff(A,B,x,y,10,10);1
A:=矩阵(%id=18446744078272833046)-1; B:=矩阵(%id=18446744078272833166)-1; 等价_A,等价-B,变换,计算时间:=MoserPfaff(A,B,x,y,10,10);1
 

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

> A:=矩阵(%id=184467440782728366662)-1; B:=矩阵(%id=18446744078272836782)-1; 等价_A,等价-B,变换,计算时间:=MoserPfaff(A,B,x,y,10,10);1
A:=矩阵(%id=184467440782728366662)-1; B:=矩阵(%id=18446744078272836782)-1; 等价_A,等价-B,变换,计算时间:=MoserPfaff(A,B,x,y,10,10);1
 

等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272839302),矩阵(%id=18446744078 272839422),矩阵
等价_A,等价-B,转换,计算时间:=矩阵(%id=18446744078272839302),矩阵(%id=18446744078 272839422),矩阵
(3.6)
 

>
 

>
 

>
 

>
 

>