SmithDecomposition

SmithDecomposition[m]

给出整数矩阵 m 的史密斯标准型分解.

更多信息

  • 结果以 {u,r,v} 形式给出,其中 uv 为幺模矩阵, r 是对角矩阵,每个对角线元素除以下一个元素,并且 u.m.vr.
  • 幺模矩阵 uv 是整数矩阵,并有 Abs[Det[u]]1,它们的逆矩阵也是整数矩阵.

范例

打开所有单元关闭所有单元

基本范例  (1)

m 分解成幺模矩阵 uv 及一个对角矩阵 r

uv 的行列式具有绝对值 1:

r 的对角线上的每个元素都整除它的下一个元素:

验证该分解:

范围  (12)

基础用法  (8)

求整数方阵的史密斯分解:

格式化结果:

奇异矩阵的史密斯分解:

r 中的非零行数等于 m 的秩:

矩形矩阵的史密斯分解:

方阵 uvm 带入对角矩形:

SmithDecomposition 用于有理矩阵:

高斯整数矩阵的史密斯分解:

高斯有理矩阵的史密斯分解:

史密斯分解仅可用于有理矩阵:

有效计算大型整数矩阵的史密斯分解:

特殊矩阵  (4)

求稀疏矩阵的史密斯分解:

结构化矩阵的史密斯分解:

IdentityMatrix 的史密斯分解由三个单位矩阵组成:

HilbertMatrix 的史密斯分解:

应用  (4)

史密斯分解可用于求解整数上的线性方程. 考虑 ,其中 如下:

由于 m=(TemplateBox[{u}, Inverse].r.TemplateBox[{v}, Inverse]),等式 m.x=(TemplateBox[{u}, Inverse].r.TemplateBox[{v}, Inverse]).x=b 表示 x=v.TemplateBox[{r}, Inverse].u.b

Solve 的结果相比:

对于可逆 ,当且仅当上述过程生成整数结果时有解:

史密斯分解给出了有限生成的阿贝尔群的规范表示. 考虑以关系 为模的整数两向量 的商群. 也就是说,如果两个点相差向量 的整数倍,则认为其为等价的. 计算矩阵 的史密斯分解,其行是

可以看出群的结构是 ,或者简单地说 因为 是平凡群:

因此,平面上的每个点都等价于九个点之一 . 考虑点 . 为了找到 的规范表示,使用右侧的 旋转到对角基并计算 的对角项的模数:

等价于 ,因此应该映射到相同的规范表示,结果也确实如此:

因此,平面中的每个整数点都可以简化为规范表示. 相反,矩阵 TemplateBox[{v}, Inverse] 可以被视为将规范基映射到起始基. 由于 生成的点阵只是与原点等效的所有点的集合,因此将该点阵旋转 TemplateBox[{v}, Inverse] 得到的点阵将与 生成的点阵相同:

使用 SmithDecomposition 来计算两个整数 的扩展最大公约数:

为项为 的列矩阵:

史密斯分解给出 ,其中 是唯一的正 幺模矩阵:

因此 First[u].mFirst[r]

但是由于 的第一行是 并且等式的左侧是 ,因此上述等式是 Bézout 恒等式. 这意味着广义最大公约数是:

ExtendedGCD 比较:

使用 SmithDecomposition 计算两个矩阵 的右侧广义最大公约数:

为由 的行组成的 矩阵:

史密斯分解中的矩阵 将是一个 矩阵:

的左上角和右上角:

矩阵 和矩阵 r.TemplateBox[{v}, Inverse] 一样是 矩阵:

表示 r.TemplateBox[{v}, Inverse] 的上半部分:

矩阵 m_a=a.TemplateBox[{g}, Inverse]m_b=b.TemplateBox[{g}, Inverse] 有整数项:

因此 的右除数:

由于 的行列式服从 gcd(TemplateBox[{a}, Det],TemplateBox[{b}, Det])=TemplateBox[{g}, Det] 实际上是右侧最大公约数:

此外,史密斯分解的前三行给出裴蜀等式

因此,广义右侧最大公约数如下:

属性和关系  (9)

分解 中的矩阵 是一个 方阵,其中 的行数:

分解 中的矩阵 是一个 方阵,其中 的列数:

分解 中的矩阵 是一个与 维数相同的对角矩形矩阵:

矩阵 都有大小为 1 的行列式:

不管输入是整数还是有理数, 矩阵总是有整数值:

如果输入矩阵是复数,该结论可推广到结果的实部和虚部:

中的每个对角线项都除以其后续项:

如果输入矩阵是整数,则 矩阵具有整数值:

如果输入矩阵有非整数有理元素,那么 矩阵也有:

对于实值输入, 矩阵的对角线项是非负的:

对于复值输入,项的实部和虚部都是非负的:

史密斯分解通常可以通过埃尔米特分解的两次迭代得到:

第一次使用 HermiteDecomposition 给出了 和一个上三角矩阵:

HermiteDecomposition 应用于三角形矩阵的 Transpose 会给出

Wolfram Research (2015),SmithDecomposition,Wolfram 语言函数,https://reference.wolfram.com/language/ref/SmithDecomposition.html.

文本

Wolfram Research (2015),SmithDecomposition,Wolfram 语言函数,https://reference.wolfram.com/language/ref/SmithDecomposition.html.

CMS

Wolfram 语言. 2015. "SmithDecomposition." Wolfram 语言与系统参考资料中心. Wolfram Research. https://reference.wolfram.com/language/ref/SmithDecomposition.html.

APA

Wolfram 语言. (2015). SmithDecomposition. Wolfram 语言与系统参考资料中心. 追溯自 https://reference.wolfram.com/language/ref/SmithDecomposition.html 年

BibTeX

@misc{reference.wolfram_2024_smithdecomposition, author="Wolfram Research", title="{SmithDecomposition}", year="2015", howpublished="\url{https://reference.wolfram.com/language/ref/SmithDecomposition.html}", note=[Accessed: 17-November-2024 ]}

BibLaTeX

@online{reference.wolfram_2024_smithdecomposition, organization={Wolfram Research}, title={SmithDecomposition}, year={2015}, url={https://reference.wolfram.com/language/ref/SmithDecomposition.html}, note=[Accessed: 17-November-2024 ]}