医药健康 基础科学 工程技术 资源环境 前沿交叉 政策管理
 
论文作者:Steven Delvaux 期刊:《线性代数及应用》 发布时间:2008-2-2 9:50:7
ELSEVIER
柯西型矩阵的快速QR分解算法
 
对于形如A=UTV(其中U和V为归一化矩阵,T为上三角矩阵)的位移结构矩阵,设计其正交分解的有效算法,目前存在很多困难。然而,这类问题的求解对于最小二乘和基于秩的QR分解等问题是非常关键的。
 
该问题的求解可转化为计算一个相关的柯西型矩阵C(其中,数据点分布在复平面内的实轴或单位圆上)的正交分解。基于柯西型矩阵的Cholesky分解和QR分解等技巧,目前已设计出一些相对有效的算法,但存在计算精度差和正交性退化等问题。
 
设H为与柯西型矩阵C相关的拟可分解矩阵,对于C的QR分解C=QR,Q的每一行为r阶矩阵H对应的特征向量,通过求解形如H—aIn的拟可分离线性方程可相应地得到R的每一列。因此,C的QR分解可通过求解H的逆特征值问题进一步得到。为求解逆特征值问题,可对C做出合理假设,利用H的一些简化形式,设计出基于发生器的算法,但缺乏数值稳定性。还可以通过求解H的QR分解迭代序列,设计基于QL的算法求解逆特征值问题,该算法具有较高的数值精度。从而,在此算法基础上可设计出柯西型矩阵C的快速QR分解算法。
 
相关论文1月15日发表在爱思唯尔期刊《线性代数及应用》(Linear Algebra and its Applications)上。(科学新闻杂志 常红旭/编译)
 
(《线性代数及应用》(Linear Algebra and its Applications),Volume 428, Issues 2-3, 15 January 2008, Pages 697-711,Luca Gemignani,Steven Delvaux)
 
更多阅读(英文)
 
发E-mail给: 
    
| 打印 | 评论 | 论坛 | 博客 |
相关论文 一周论文排行

小字号

中字号

大字号

Baidu
map