内容简介作者简介 - 科学与工程计算系

11.3.2 由 M 文件创建外部应用程序的入门算例
【例 11.3.2-1】建立一个脱离 MATLAB 环境,可独立运行的外部程序。该程序的功能是:
对于给定矩阵 A,如果存在 S 使得 S -1 AS=Λ,则要求出一个 S,否则给出信息说明所给的矩
阵 A 不能对角化。
(1)
[exm2.m]
function exm2
A=[4,0,0;0,3,1;0,1,3];
S=exm2_f(A)
[exm2_f.m]
function S=exm2_f(A)
[m,n]=size(A);
if m~=n
error('输入矩阵应是方阵!');
end;
e=eig(A);
%
same=0;
for i=1:m-1
for j=(i+1):m
if e(j)==e(i)
same=1;
end
end
end
%
if any(any((A'-A)))&(same==1)
error('矩阵无法对角化!');
end
[v,d]=eig(A);
S=v;
(2)
exm2
S =
0 0 1.0000
-0.7071 0.7071 0
0.7071 0.7071 0
(3)
mcc -m exm2 exm2_f
(4)
5