China-pub.CoM 第I2章MATLAB程序设计 171 下载 if rem(i,m)==0 A,i,nd end end 下面的程序是用来完成非移位的QR方法（见8.2节）30次迭代，每隔15次输出结果： A0=[-9-3-16;13716;3310]; m=15;n=30; format long; A=A0; qrmethod; A= 9.9899746707437722.62301237506363-15.53274662438004 0.00708686385759-5.985685125529255,77401643542405 00.00741470005235 3.99571045478546 i= 15 nd 0.01025677416162 A= 10.0000047162466022.6274399374496715.51339551121122 -0.00000333488655-6.00001449452640-5.77348898879412 0 0.000016936546124.00000977827978 i= 30 nd 1.726175143943722e-05 (e)在7.5节中定义了命令planerot。这个算法是使用该命令来返回一个矩阵，该矩阵的 0元素都在任何作为参数输入的m×n矩阵的主对角线的下方。 function B Givens(A) 号如果将A矩阵乘函数B返回的矩阵，就能将大小为m×的A矩阵变成上三角阵。也就是说，根据Q=和· R=B+A,其中满足Q+R=A,该函数可以对A进行QR分解。 [m,n]size(A); =eye(m); for j=1:n for i j:m for k (i+1):m G =eye(m); Plan planerot([A(j,j)A(k,j)]'); %找出2×2的矩阵 G(j,j)=P1an(1,1);G(k,j)=P1an(2,1); G(j,k)=P1an(1,2);G(k,k)=P1an(2,2);下面的程序是用来完成非移位的 Q R方法(见8 . 2节) 3 0次迭代，每隔1 5次输出结果： (e) 在7 . 5节中定义了命令p l a n e r o t。这个算法是使用该命令来返回一个矩阵，该矩阵的 0元素都在任何作为参数输入的 m×n矩阵的主对角线的下方。 function B = Givens(A) % % 如果将 A矩阵乘函数 B返回的矩阵，就能将大小为 m×n的A矩阵变成上三角阵。也就是说，根据 Q = B ´ 和 R = B*A，其中满足Q*R = A，该函数可以对A进行Q R分解。 第1 2章 M AT L A B程序设计 1 7 1 下载 ％找出2×2的矩阵