资源描述：
《童玲论文草稿（精荐）》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、「-循环矩阵求逆的几种算法童玲（陕西理工学院数学与计算机科学学院数学与应用数学专业093班，陕西汉屮723100）指导教师：陈露［摘要］厂-循环矩阵逆矩阵是一类特殊的矩阵，及其算法是很重要的一类方法•本文归纳了广-循环矩阵求逆的几种算法，并通过算例对几种算法进行了比较分析.［关键词］r•循环矩阵；逆矩阵；插值法；欧几里德算法；最大公因式.1引言在信号处理、数字图像处理、编码理论、自回归滤波器设计、计算机时序分析以及其他应用学科，常常会遇到各种形式的厂循环矩阵⑴，其中不少问题牵涉到求逆的计算•因此，对它的
2、研究就引起了许多学者的高度重视•近年来，对于「循环短阵求逆的算法的研究呈现出上升的态势，同时也取得一定的成果⑵.木文在查阅文献资料的基础上,对r-循环矩阵求逆算法研究成果进行归纳总结，分析了厂-循环矩阵求逆的初等算法，介绍了厂-循环短阵求逆的新算法，最后利用最大公因式算法介绍r-循环矩阵求逆的快速算法•通过具体例了对儿种算法进行了比较分析.2预备知识定义1⑴若矩阵A具有形状/,0()4…°口-2an-'务）a…仏2A=叫一2叫1G()…S—3•••<5•••ra2•••ra3••••••…叫•••
3、a())(2.1)则称A为厂-循环矩阵.可以看出A完全由厂及®（i=0,/?-1）决定，故可将A写作:AA匚.（佈4，・・.,心_2卫灯）w（“y表示“记为”）.特別地：当厂=1时，就是通常所说的循环矩阵；当r=-l吋，则为反循环矩阵.定义2工称n阶尸■循环矩阵01•■•000000••••••0100为基本r■循环矩阵，简记为J=Cr(0,1,0・・・0).显然JgCMr,且有J1=C,(O...,o/bO,ooo,0),J°=EnJn=rEn.其屮E“为〃阶单位矩阵.上述两定义屮，若r=l时，则/•
4、可以省略不写.3厂-循环矩阵的逆矩阵的初等算法/<0%A=■■■+<畑严1+6Z
5、(0E心0丿+…+E叫］0■■0rE2+…+an-(3.1)3.1定理与算理如果A=a,(i,j=0,1,...一1),则aj_i，i5jran+H^>j・观察(2.1)式对角线及对角线位置上的元素，即任何r-循环矩阵都可以写成:显然有0丿'0eCMr.f0JACr(0,l,0,...,0)=厂=“J1=cr(o,...,o,i,o,...,o)=&Jn+k=rJ(k为非负整数),丿的特征多项式为g(x)=xn-r.由
6、(3.1)式与定义2可知：引理1⑵若AwC"x"，则A-CrgCMrw-1的充分必要条件为A=/(丿)=工％人，其屮/=0〃一1f(x)=^jaix,../=0显然线性无关，因此A的表达式唯一，即任何厂一循环矩阵都可以由丿的多项式唯一表示.引理2⑶设A二C「(d()，e，...,d”_JwCM,.,/二C,0,...,0,T,0,・・・,0)wCM「则有A=/G()■■■+/■■■J+/a、•■■J~+...+z、an-'■■■Jn-[7、B=BAeCA/r.定理1⑷设矩阵A=Cr(^0,6z1,...,^_1)eCMr且A非奇异，则A-1=B=Cr(/?0,^CMr,其中bo，b,・.・bn7是线性方程组"doJd]Qo••••••...ra2...ra}••••••raxra2■■■"X。'•■■<^-l>O的唯一解.证明因为A二C『(Qo4,…4_1)丘CM,,所以A二。0丿°+如丿+...+。心厂一，欲证尸一循环矩阵A的逆矩阵也是厂一循环矩阵,只要找到B=
8、b(/°+bJ+...+b_/"T,注意严=rJk(R为非负数).AB=(%)丿°+d/+…厂一)(»丿°+仇丿+…+b心尸I)=(raibn-l+ra2bn-2+…+®丿
9、+d(0()"°+(ra2bn^+ra3bn^2+…++a}bQ“+•••+g」h、+aQhn_2+...+色_3肉+a』(J广2+So乞一i+a(A—2+・・・+d-2勺+心・仇)厂"仝Co丿。+C]
10、bn-2+也・1=】C]=a}bQ+(70b,+…+ra3bn_2+ra2bn^=0C“_2=+an-3bl+・・・+-2+叫一也刊=0c”_
11、=+an_2b}+...+a{bn_2+a^)n_x=0ran-100an-2^an-5-3an-2ra2、rax“3■■ra2••••0■00■叫1•■06%>O它与线性方程组an-2存3an-2ra25/、兀0■■T0■■■兀_2■0/K->Ora{ra2G()等价.因为A非奇异，所以
12、a7
