a11_3b.m

来自「Matlab numerical methods,examples of mat」· M 代码 · 共 121 行

echo on;
%---------------------------------------------------------------------------
%A11_3B   MATLAB script file for implementing Algorithm 11.3
%
% NUMERICAL METHODS: MATLAB Programs, (c) John H. Mathews 1995
% To accompany the text:
% NUMERICAL METHODS for Mathematics, Science and Engineering, 2nd Ed, 1992
% Prentice Hall, Englewood Cliffs, New Jersey, 07632, U.S.A.
% Prentice Hall, Inc.; USA, Canada, Mexico ISBN 0-13-624990-6
% Prentice Hall, International Editions:   ISBN 0-13-625047-5
% This free software is compliments of the author.
% E-mail address:      in%"[e-mail]"
%
% Algorithm 11.3 (Jacobi Iteration for Eigenvalues and Eigenvectors).
% Section	11.3, Jacobi's Method, Page 571
%---------------------------------------------------------------------------

clc; clear all; format long; delete output;

% - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%
%        JACOBI METHOD FOR FINDING EIGEN-PAIRS
%
% Assume that A is an n by n symmetric real matrix, thus
%
% it has a full set of eigenvectors  V , V  ,..., V .
%                                     1   2        n
%
% The cyclic Jacobi`s method of iteration is used
%
% to find all of the eigen-pairs.
%
% Remark. jacobi2.m is used for Algorithm 11.3b

pause % Press any key to continue.

clc;

% - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%
%             JACOBI`S  METHOD FOR EIGENVECTORS
%
%      Assume that A is an n by n real symmetric matrix.
% Then A has a full set of eigenvectors:
% V , V ,..., V .  Jacobi`s method of iteration is used
%  1   2       n
% to find all the eigenvalues and eigenvectors of A. Let
% A = A , and construct a sequence of orthogonal
%      1
%                                    T
% matrices { R  }  such that  D  =  R  A  R .
%             j                j     j  j  j
% { D  } converges to the diagonal matrix D ,
%    j
% of eigenvalues, and the sequence {V  = R R ...R }
%                                    j    1 2    j
%
% converges to the matrix of eigenvectors.
%
pause % Press any key to continue.

clc;

% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%
%      We use the cyclic Jacobi strategy for selecting the off
% diagonal element to annihilate in the construction process.
% Sweep through the matrix and annihilate elements in the
% strict order
%
% a  ,a  ,...,a  ; a  ,a  ,...,a  ; ...;a
%  12  13      1n   23  24      2n       n-1,n
%
%      For each sweep throughout the matrix, the computed
% value t is the [R.M.S.] average of the diagonal elements
% of A.  The user must supply the error tolerance for
% annihilating the off diagonal elements
%
%        |d   |  >  t*epsilon       for all  q > p.
%          p,q

pause % Press any key to continue.

clc;

% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%
% Example 11.7, page 569.  Finding all eigenvalues of the matrix A.
%
% Enter the matrix in  A
%
% Enter the tolerance in   epsilon
%
% Enter the maximum cycles in  max1

A = [ 8   -1    3   -1;
     -1    6    2    0;
      3    2    9    1;
     -1    0    1    7];

epsilon = 1e-16;

[V,D] = jacobi2(A,epsilon,1);

pause % Press any key to continue.

clc;
%............................................
% Begin section to print the results.
% Diary commands are included which write all
% the results to the Matlab textfile   output
%............................................
Mx1 = 'Implementation of the cyclic Jacobi method.';
Mx2 = 'The matrix  A  is:';
Mx3 = 'The eigenvalues of  A  are:';
Mx4 = 'The eigenvectors of  A  are:';
clc,echo off,diary output,...
disp(''),disp(Mx1),disp(''),disp(Mx2),disp(A),
disp(Mx3),disp(diag(D)),disp(Mx4),disp(V),...
diary off,echo on