八、矩阵特征值的jacobi方法.doc_第1页
八、矩阵特征值的jacobi方法.doc_第2页
八、矩阵特征值的jacobi方法.doc_第3页
八、矩阵特征值的jacobi方法.doc_第4页
八、矩阵特征值的jacobi方法.doc_第5页
已阅读5页，还剩2页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、用jacobi方法计算对称矩阵A的特征值和对应的特征向量。function k,Bk,V,D,Wc=jacobite(A,jd,max1)n,n=size(A);P0=eye(n);Vk=eye(n);Bk=A;k=1;state=1;while (kjd) state=1; end k,k=k+1;Pk,Vk,Bk=B2,Wc,endif (kmax1) disp(迭代次数k已经达到最大迭代次数ma1，迭代次数k，对称矩阵Bk，以特征向量为列向量的矩阵V，特征值为对角元的对角矩阵D如下：)else disp(迭代次数k，对称矩阵Bk，以特征向量为列向量的矩阵V，特征值为对角元的对角矩阵D如下：)endk=k-1;Bk=B2;V=Vk;D=diag(diag(Bk);Wc;V1,D1=eig(A,nobalance) A=12 -56 3 -1;-56 7 2 0;3 2 5 1;-1 0 1 12; k,Bk,V,D,Wc=jacobite(A,0.0001,100)k = 1Pk = 0.7227 0.6912 0 0 -0.6912 0.7227 0 0 0 0 1.0000 0 0 0 0 1.0000Vk = 0.7227 0.6912 0 0 -0.6912 0.7227 0 0 0 0 1.0000 0 0 0 0 1.0000Bk = 65.5558 0 0.7858 -0.7227 -0.0000 -46.5558 3.5189 -0.6912 0.7858 3.5189 5.0000 1.0000 -0.7227 -0.6912 1.0000 12.0000Wc = 56k = 2Pk = 1.0000 0 0 0 0 0.9977 0.0678 0 0 -0.0678 0.9977 0 0 0 0 1.0000Vk = 0.7227 0.6896 0.0468 0 -0.6912 0.7210 0.0490 0 0 -0.0678 0.9977 0 0 0 0 1.0000Bk = 65.5558 -0.0533 0.7840 -0.7227 -0.0533 -46.7948 0 -0.7574 0.7840 0.0000 5.2391 0.9509 -0.7227 -0.7574 0.9509 12.0000Wc = 3.5189k = 3Pk = 1.0000 0 0 0 0 1.0000 0 0 0 0 0.9906 0.1367 0 0 -0.1367 0.9906Vk = 0.7227 0.6896 0.0464 0.0064 -0.6912 0.7210 0.0485 0.0067 0 -0.0678 0.9883 0.1364 0 0 -0.1367 0.9906Bk = 65.5558 -0.0533 0.8754 -0.6088 -0.0533 -46.7948 0.1035 -0.7502 0.8754 0.1035 5.1079 -0.0000 -0.6088 -0.7502 -0.0000 12.1312Wc = 0.9509k = 4Pk = 0.9999 0 -0.0145 0 0 1.0000 0 0 0.0145 0 0.9999 0 0 0 0 1.0000Vk = 0.7233 0.6896 0.0359 0.0064 -0.6904 0.7210 0.0585 0.0067 0.0143 -0.0678 0.9882 0.1364 -0.0020 0 -0.1367 0.9906Bk = 65.5685 -0.0518 -0.0000 -0.6087 -0.0518 -46.7948 0.1043 -0.7502 -0.0000 0.1043 5.0952 0.0088 -0.6087 -0.7502 0.0088 12.1312Wc = 0.8754k = 5Pk = 1.0000 0 0 0 0 0.9999 0 -0.0127 0 0 1.0000 0 0 0.0127 0 0.9999Vk = 0.7233 0.6896 0.0359 -0.0024 -0.6904 0.7211 0.0585 -0.0025 0.0143 -0.0660 0.9882 0.1372 -0.0020 0.0126 -0.1367 0.9905Bk = 65.5685 -0.0595 -0.0000 -0.6080 -0.0595 -46.8044 0.1044 0.0000 -0.0000 0.1044 5.0952 0.0075 -0.6080 0.0000 0.0075 12.1407Wc = 0.7502k = 6Pk = 0.9999 0 0 0.0114 0 1.0000 0 0 0 0 1.0000 0 -0.0114 0 0 0.9999Vk = 0.7233 0.6896 0.0359 0.0059 -0.6903 0.7211 0.0585 -0.0103 0.0127 -0.0660 0.9882 0.1374 -0.0132 0.0126 -0.1367 0.9905Bk = 65.5754 -0.0595 -0.0001 -0.0000 -0.0595 -46.8044 0.1044 -0.0007 -0.0001 0.1044 5.0952 0.0075 -0.0000 -0.0007 0.0075 12.1338Wc = 0.6080k = 7Pk = 1.0000 0 0 0 0 1.0000 0.0020 0 0 -0.0020 1.0000 0 0 0 0 1.0000Vk = 0.7233 0.6895 0.0373 0.0059 -0.6903 0.7209 0.0600 -0.0103 0.0127 -0.0680 0.9881 0.1374 -0.0132 0.0129 -0.1366 0.9905Bk = 65.5754 -0.0595 -0.0002 -0.0000 -0.0595 -46.8046 -0.0000 -0.0007 -0.0002 -0.0000 5.0954 0.0075 -0.0000 -0.0007 0.0075 12.1338Wc = 0.1044k = 8Pk = 1.0000 0.0005 0 0 -0.0005 1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000Vk = 0.7229 0.6899 0.0373 0.0059 -0.6907 0.7206 0.0600 -0.0103 0.0128 -0.0680 0.9881 0.1374 -0.0133 0.0129 -0.1366 0.9905Bk = 65.5754 0.0000 -0.0002 0.0000 -0.0000 -46.8046 -0.0000 -0.0007 -0.0002 -0.0000 5.0954 0.0075 0.0000 -0.0007 0.0075 12.1338Wc = 0.0595k = 9Pk = 1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000 0.0011 0 0 -0.0011 1.0000Vk = 0.7229 0.6899 0.0373 0.0059 -0.6907 0.7206 0.0600 -0.0103 0.0128 -0.0680 0.9880 0.1384 -0.0133 0.0129 -0.1377 0.9903Bk = 65.5754 0.0000 -0.0002 0.0000 -0.0000 -46.8046 0.0000 -0.0007 -0.0002 0.0000 5.0954 0.0000 0.0000 -0.0007 -0.0000 12.1338Wc = 0.0075k = 10Pk = 1.0000 0 0 0 0 1.0000 0 -0.0000 0 0 1.0000 0 0 0.0000 0 1.0000Vk = 0.7229 0.6899 0.0373 0.0059 -0.6907 0.7206 0.0600 -0.0103 0.0128 -0.0680 0.9880 0.1384 -0.0133 0.0129 -0.1377 0.9903Bk = 65.5754 0.0000 -0.0002 0.0000 0.0000 -46.8046 0.0000 0.0000 -0.0002 0.0000 5.0954 -0.0000 0.0000 -0.0000 -0.0000 12.1338Wc = 6.9206e-004k = 11Pk = 1.0000 0 0.0000 0 0 1.0000 0 0 -0.0000 0 1.0000 0 0 0 0 1.0000Vk = 0.7229 0.6899 0.0373 0.0059 -0.6907 0.7206 0.0600 -0.0103 0.0128 -0.0680 0.9880 0.1384 -0.0133 0.0129 -0.1377 0.9903Bk = 65.5754 -0.0000 0.0000 0.0000 -0.0000 -46.8046 0.0000 0.0000 -0.0000 0.0000 5.0954 -0.0000 0.0000 -0.0000 -0.0000 12.1338Wc = 2.0482e-004k = 12Pk = 1.0000 0 0 0 0 1.0000 0.0000 0 0 -0.0000 1.0000 0 0 0 0 1.0000Vk = 0.7229 0.6899 0.0373 0.0059 -0.6907 0.7206 0.0600 -0.0103 0.0128 -0.0680 0.9880 0.1384 -0.0133 0.0129 -0.1377 0.9903Bk = 65.5754 -0.0000 0.0000 0.0000 -0.0000 -46.8046 -0.0000 0.0000 -0.0000 0.0000 5.0954 -0.0000 0.0000 -0.0000 -0.0000 12.1338Wc = 6.2740e-007迭代次数k，对称矩阵Bk，以特征向量为列向量的矩阵V，特征值为对角元的对角矩阵D如下：V1 = 0.6899 -0.0373 0.0059 -0.7229 0.7206 -0.0600 -0.0103 0.6907 -0.0680 -0.9880 0.1384 -0.0128 0.0129 0.1377 0.9903 0.0133D1 = -46.8046 0 0 0 0 5.0954 0 0 0 0 12.1338 0 0 0 0 65.5754k = 12Bk = 65.5754 -0.0000 0.0000 0.0000 -0.0000 -46.8046 -0.0000 0.0000 -0.0000 0.000
人人文库> 全部分类> 应用文书 > 技术指导

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论