reduce_diag.m

MATLAB科学计算与工程分析源代码源程序4

function [P,D]=reduce_diag(A)
% 该函数用来将一个矩阵A对角化
% 输出变量为矩阵P，满足inv(P)*A*P=diag(lambda_1,...,l.lambda_n)
% 编写日期：2007-5-12

if ~isdiag(A)       % 判断矩阵A是否可化为对角矩阵
   error('该矩阵不能对角化！');
else    
    disp('注意：将下面的矩阵P的任意列乘以任意非零数所得矩阵仍满足inv(P)*P*A=D');
    [P,D]=eig(A);      
end