(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211118814.6 (22)申请日 2022.09.13 (71)申请人 陕西黄河集团有限公司 地址 710043 陕西省西安市新城区幸福北 路21号 (72)发明人 张晓伟　晁智郁　杨咚咚　 (74)专利代理 机构 西安亚信智佳知识产权代理 事务所(普通 合伙) 61241 专利代理师 段国刚 (51)Int.Cl. G06F 17/16(2006.01) (54)发明名称 一种基于并行计算雅克比旋转的矩阵特征 分解方法 (57)摘要 本公开是关于一种基于并行计算雅克比旋 转的矩阵特征分解方法， 该方法包括： 按照给定 的处理顺序， 选择待处理的第一矩阵中的对应元 素， 得到雅克比旋转因子矩阵； 根据所述雅克比 旋转因子矩 阵中每行和每列只有两个非零元素 的特点， 进行矩阵相乘运算， 分别得到第二矩阵 和第三矩阵； 迭代更新所述第一矩阵中的对应元 素， 重复上述所有步骤， 直至处理完所述第一矩 阵中的所有元素， 该方法能够 有效降低计算过程 中重构误差幅值和正交矩 阵与单位矩 阵之间的 绝对误差幅值。 权利要求书4页 说明书8页 附图4页 CN 115455354 A 2022.12.09 CN 115455354 A 1.一种基于并行计算 雅克比旋转的矩阵特 征分解方法， 其特 征在于， 包括以下步骤： 按照给定的处理顺序， 选择待处理的第一矩阵中的对应元素， 得到雅克比旋转因子矩 阵； 根据所述雅克比旋转因子矩阵中每行和每列只有两个非零元素的特点， 进行矩阵相乘 运算， 分别得到第二矩阵和第三矩阵； 迭代更新所述第一矩阵中的对应元素， 重复上述所有步骤， 直至处理完所述第一矩阵 中的所有元 素。 2.根据权利要求1所述基于并行计算雅克比旋转的矩阵特征分解方法， 其特征在于， 按 照给定的处理顺序， 选择待处理的第一矩阵中的对应元素， 得到雅克比旋转因子矩阵的步 骤包括： 按照给定的处理顺序， 选 择待处理的第一矩阵A0中的对应元素， 计算得到雅克比旋转因 子矩阵 其中， c， s和 ‑s分别为雅克比旋转因子 。 3.根据权利要求1所述基于并行计算雅克比旋转的矩阵特征分解方法， 其特征在于， 根 据所述雅克比旋转因子矩阵中每行和每列只有两个非零元素 的特点， 进行矩阵相乘运算， 得到第二矩阵的步骤 包括： 将所述雅克比 旋转因子矩阵 进行转置 ， 得到转置矩阵 并将所述转置矩阵 和所述第一矩阵A0中对 应的行向量进行相乘运 算， 得到矩阵B0； 将所述矩阵B0对应的列向量和所述雅克比旋转因子矩阵 进行相 乘运算， 得到第二矩阵Ai； 其中， T表示矩阵的转置操作。权　利　要　求　书 1/4 页 2 CN 115455354 A 24.根据权利要求3所述基于并行计算雅克比旋转的矩阵特征分解方法， 其特征在于， 将 所述转置矩阵 和所述第一矩阵A0中对应的行向量进行相乘运算， 得 到矩阵B0的公式包括： 其中， c， s和 ‑s分别为雅克比旋转因子； ai,j表示所述第一矩阵中的元 素。 5.根据权利要求 4所述基于并行计算 雅克比旋转的矩阵特 征分解方法， 其特 征在于， 将所述矩阵B0对应的列向量和所述雅克比旋转因子矩阵 进行相 乘运算， 得到第二矩阵Ai的公式包括： 其中， 所述 雅克比旋转因子矩阵D中非零元 素的位置为每 个i＝j的位置 。 6.根据权利要求3所述基于并行计算雅克比旋转的矩阵特征分解方法， 其特征在于， 根 据所述雅克比旋转因子矩阵中每行和每列只有两个非零元素 的特点， 进行矩阵相乘运算， 得到第三矩阵的步骤 包括： 引入单位矩阵I0， 将所述转置矩阵 和所述单位矩阵I0中对应的行 向量相乘， 得到矩阵C0；权　利　要　求　书 2/4 页 3 CN 115455354 A 3