本发明属于信号处理领域,尤其涉及种基于jacobi迭代算法的高精度矩阵特征值分解实现方法。
背景技术:
在信号处理中,矩阵的特征值分解evd是一个应用广泛的矩阵运算。如数据压缩、噪声去除、数值分析,包括近几年兴起的机器学习、深度学习其基本核心操作也包括矩阵特征值分解。实现矩阵特征值分解的常用方法有gauss变换、householder变换、jacobi迭代等,其中,jacobi迭代是精度较高的方法,并且很适合在fpga中实现。因此一种基于jacobi迭代算法的高精度矩阵特征值分解实现技术在实际工程中具有很高的应用价值。
经典的jacobi迭代算法计算共轭矩阵a∈cn×n的特征值分解如图1所示,这种经典的迭代算法虽然有较快收敛速度,但是该算法需要在矩阵a的众多元素中选取aij,使得aij为非对角元素中绝对值最大的一个,再进行后面的计算操作。这样每一步都要寻找绝对值最大的非对角元,比较费时也不适合在fpga实现,因此经典的jacobi迭代算法在实际工程中并不实用。
目前实际工程中多数采用如图2所示的循环jacobi迭代算法,通过逐行扫描遍历法选取aij,这样避免了寻找最大绝对值的非对角元的复杂繁琐步骤。这样选取aij的方式,在aij数值比较大时,fpga中使用cordic算法计算φ、
技术实现要素:
发明的目的在于解决在循环jacobi迭代算法进行矩阵特征值分解过程中,因为逐行扫描遍历法选取aij在aij比较小甚至接近于0时,导致fpga中使用cordic算法计算φ、
一种基于jacobi迭代算法的高精度矩阵特征值分解实现方法,包括如下步骤:
s1、设数据矩阵a∈cn×n为共轭矩阵,并且设定最大遍历次数为t、最小清扫门限a、扩位门限b和算术左移位数m,其中,最小清扫门限a应小于要求计算结果的精度一个数量级,扩位门限b与算术左移位数m与fpga实现里数据位宽size有关,满足b×2m<2size-4保证计算结果不溢出,n为不为零的自然数,共轭矩阵a中的元素为aij,i=1,2,3,...,n,j=1,2,3,...,n,1≤t≤t且t为自然数;
s2、初始化遍历次数计数器,令t=0,
初始化特征向量初始矩阵,令v=e,其中,e为单位阵;
s3、在s1所述共轭矩阵a中选取aij,初始化清扫元aij行列下标,令i=1,j=2;
s4、判断aij是否满足跳过清扫条件|real(aij)|<a&|imag(aij)|<a,若满足则转入s10,如不满足则转入s4;
s5、判断aij是否满足扩展条件|real(aij)|<b&|imag(aij)|<b,若满足则转入s6,若不满足则转入s7;
s6、进行位扩展,即计算a′ij=aij×2m,转入s7;
s7、令a′ij=aij,进入s8;
s8、计算
s9、计算a=qhaq和v=qhv,其中,q∈cn×n为复数域内的平面旋转矩阵,
s10、判断j=n是否成立,是则进入s11,否则j=j+1后跳转到s4;
s11、判断i=n-1是否成立,是则进入s12,否则i=i+1,j=i+1后跳转到s4;
s12、判断t=t是否成立,是则进入s13,否则t=t+1后跳转到s3;
s13、输出迭代计算结果a和v,其中,a的对角元数位s1输入数据矩阵a的特征值,v为对应的特征向量矩阵。
进一步地,s1所述遍历次数t越大迭代次数越多,计算越精确,但计算时间越长,为了取得速度与精度的平衡,当n≤8时t=3,当n>8时t=6。
进一步地,s1所述a=e×10-1。
本发明的有益效果是:
在不明显增加算法复杂度、fpga实现难度与硬件资源消耗、计算消耗时间的情况下,提高基于循环jacobi迭代算法的fpga实现矩阵特征值分解的计算精度和计算速度,在实际工程中具有重要价值。
附图说明
图1为经典jacobi迭代算法流程。
图2为循环jacobi迭代算法流程。
图3为本发明算法流程。
具体实施方式
下面将结合实施例和附图,对本发明方法进行进一步说明。
本发明应用于估计信号与噪声对应的盖氏圆盘半径,提高圆盘半径计算精度,和计算速度。
实施例1、
接收阵列为8阵元组成的均匀线阵。
如图3所示,考虑n=1个载波频率为
在估计性能包括计算精度、计算速度和资源消耗,具体用下面指标评价:
1.计算精度:
(1).圆盘半径计算精度:
(2).圆盘平均计算精度:
2.计算速度:
(1).计算消耗的时钟数nclk,越小表示计算消耗时间越少,计算速度越快。
3.资源消耗:
(1).寄存器消耗数量nreg,越小对应寄存器资源消耗越少。
(2).逻辑门消耗数量nlut,越小对应逻辑门资源消耗越少。
应用特征值分解估计信号与噪声对应的盖氏圆盘半径包括,a.仿真接收信号数据建模、b.应用本发明进行特征值分解、c.计算圆盘半径,具体为以下步骤:
a.仿真接收信号数据建模。
a1.由下式产生阵元数n=8的阵列接收信号向量x(k)=[x1(k)x2(k)…x8(k)]h进入步骤a2。
x(k)=a(γ)s(k)+n(k),k=1,2,…,l
式中,n(k)为8×1为均值为零、方差σ2=1的高斯白噪声矢量;远场接收信号s(k)=ass(k),其中其幅度as=10snr/20;a(γ)=[1e-jφ…e-j(n-1)φ]t,
a2.由
a3.根据
b.应用本发明对块矩阵r′进行特征值分解,计算r′的特征值矩阵d与特征向量矩阵v。
b1.进行初始化,具体方法为:
b11.设数据矩阵a=r′为共轭矩阵,并且设定遍历次数t=5、最小清扫门限a=10-8、扩位门限b=10-5和算术左移位数m=8,进入步骤b12。
其中t越大迭代次数越多计算越精确但计算时间越长,根据矩阵维数n进行选取,当n≤8时t=4,n>8时t=6可以取得速度与精度的平衡;最小清扫门限a与计算精度要求e有关,a=e×10-1,比如要求计算精度为10-5则a≈10-6。扩位门限b与算术左移位数m与fpga实现里数据位宽size有关,满足b×2m<2size-4保证计算结果不溢出。
b12.初始化遍历次数计数器和特征向量初始矩阵,t=0、v=e,其中,e为单位阵,进入步骤b13。
b13.初始化清扫元aij的行列下标,i=1、j=2,进入步骤b2。
b2.进行jaocbi旋转,具体方法如下:
b21.在矩阵a中选取aij,进入步骤b22。
b22.判断是否满足跳过清扫条件|real(aij)|<a&|imag(aij)|<a,是则跳转到步骤b3,否则进入步骤b23。
b23.判断是否满足扩展条件|real(aij)|<b&|imag(aij)|<b,是则进入步骤b24,否则进入步骤b25。
b24.进行位扩展,即a′ij=aij×2m,进入步骤b26。
b25.不进行位扩展,即a′ij=aij,进入步骤b26。
b26.计算相角和模值,即
b27.计算平面旋转角,即
b28.进行jacobi旋转,即计算a=qhaq和v=qhv,其中q∈cn×n为复数域内的平面旋转矩阵。
即q的对角元素中除了qii=ejφcosθ、qjj=e-jφcosθ之外其他均为1,非对角元素中除了qij=-ejφsinθ、qji=e-jφsinθ其他元素均为0,进入步骤b3。
b3.对迭代过程进行判断。
b31.判断j=n是否成立,是则进入步骤b32,否则j=j+1后跳转到步骤b21。
b32.判断i=n-1是否成立,是则进入步骤b33,否则i=i+1,j=i+1后跳转到步骤b21。
b33.判断t=t是否成立,是则进入步骤b4,否则t=t+1后跳转到步骤b13。
b4.输出迭代计算结果a和v,其中a的对角元就是数据矩阵a的特征值,v为对应的特征向量矩阵,进入步骤c。
c.对数据协方差矩阵r进行酉变换,计算信号与噪声对应的圆盘半径。
c1.由下式构造酉变换矩阵t∈cn×n,进入步骤c2。
其中,v∈c(n-1)×(n-1)为前面计算块矩阵r′的特征向量,满足vvh=e,e为单位阵。
c2.进行酉变换得到圆盘半径,即计算下式,进入步骤c3。
式中,λi,i=1,2,…,n-1为块矩阵r′的特征值。
c3.由ri=|ρi|,i=1,2,…,n-1计算圆盘半径ri,进入步骤c4。
c4.计算圆盘半径计算精度:
c5.统计计算消耗的时钟数nclk、寄存器消耗数量nreg和逻辑门消耗数量nlut,算法结束。
仿真结果为:
计算精度:
计算速度:nclk=11710
资源消耗:nreg=29104、nlut=30254
此时,估计信号所对应的圆盘半径计算精度为ε1≈10-9;估计噪声所对应的圆盘半径计算精度为εi≈10-4,i=2,3,…,7;圆盘平均计算精度
实施例2、
经典方法循环jacobi算法应用于估计信号与噪声对应的盖氏圆盘半径的估计性能,作为实施例1的对比例。
实施例2的方法如附图2所示,其余仿真条件与实施例1的相同,进行信号与噪声对应的盖氏圆盘半径的估计。
实施例2的评价标准与实施例1的一致。
仿真结果为:
计算精度:
计算速度:n′clk=17960
资源消耗:n′reg=29101、n′lut=29998
本发明此时,估计信号所对应的圆盘半径计算精度为ε1≈10-8;估计噪声所对应的圆盘半径计算精度为εi≈10-1,i=2,3,…,7;圆盘平均计算精度
综上所述,对比实施例1与实施例2,本发明相对于经典方法在增加(nreg-n′reg)/n′reg×%≈0.01%寄存器资源消耗,(nlut-n′lut)/n′lut×%≈0.85%查找表资源消耗的情况下,平均计算精度从10-1提高到了10-4数量级,提高了3个数量级,同时计算速度提高了|nclk-n′clk|/n′clk×%≈34.8%。
所以,本发明在基本不增加资源消耗的情况下,不仅可以提高计算精度,还可以提高计算速度,在实际工程中具有重要价值。
资讯网魔兽图可以制作海报的网站胡兵演宋尧的电视剧网站排名优化公司哪如何制作社交网站年生的女孩子起名做直饮水给公司起名字企业文化网站设计周公能解梦大全电影客栈最新茶叶店起名重生文推荐女装品牌名字起名周易道德经黄帝内经算命一般姻缘怎么说美容养生排行榜876av电影天堂酒店起名字大全免费吗公司取名起名大全凶吉温州网站建设好的公司硬盘分区合并企业网企业网站建设公司网页seo优化技巧周易起名软件手机版穿越之废材七小姐阳字起名大全线下营销推广方式企业英文网站建设seo技术培训南阳睢县海港城少年生前被连续抽血16次?多部门介入两大学生合买彩票中奖一人不认账让美丽中国“从细节出发”淀粉肠小王子日销售额涨超10倍高中生被打伤下体休学 邯郸通报单亲妈妈陷入热恋 14岁儿子报警何赛飞追着代拍打雅江山火三名扑火人员牺牲系谣言张家界的山上“长”满了韩国人?男孩8年未见母亲被告知被遗忘中国拥有亿元资产的家庭达13.3万户19岁小伙救下5人后溺亡 多方发声315晚会后胖东来又人满为患了张立群任西安交通大学校长“重生之我在北大当嫡校长”男子被猫抓伤后确诊“猫抓病”测试车高速逃费 小米:已补缴周杰伦一审败诉网易网友洛杉矶偶遇贾玲今日春分倪萍分享减重40斤方法七年后宇文玥被薅头发捞上岸许家印被限制高消费萧美琴窜访捷克 外交部回应联合利华开始重组专访95后高颜值猪保姆胖东来员工每周单休无小长假男子被流浪猫绊倒 投喂者赔24万小米汽车超级工厂正式揭幕黑马情侣提车了西双版纳热带植物园回应蜉蝣大爆发当地回应沈阳致3死车祸车主疑毒驾恒大被罚41.75亿到底怎么缴妈妈回应孩子在校撞护栏坠楼外国人感慨凌晨的中国很安全杨倩无缘巴黎奥运校方回应护栏损坏小学生课间坠楼房客欠租失踪 房东直发愁专家建议不必谈骨泥色变王树国卸任西安交大校长 师生送别手机成瘾是影响睡眠质量重要因素国产伟哥去年销售近13亿阿根廷将发行1万与2万面值的纸币兔狲“狲大娘”因病死亡遭遇山火的松茸之乡“开封王婆”爆火:促成四五十对奥巴马现身唐宁街 黑色着装引猜测考生莫言也上北大硕士复试名单了德国打算提及普京时仅用姓名天水麻辣烫把捣辣椒大爷累坏了