矩阵论，盖尔圆，MATLAB代码

矩阵论，盖尔圆，MATLAB代码

clc;
clear;
close all;

%%
% A=[1i -1 1 0;1 20 -1 -1;-6 -6 9 3;0 1 0 -5];
% A=[-9 0 -1 1.2;1.1 -3 -1 2;0 -1.4 5 0;0.3 -0.7 0 2];
% A=[-1i 0 -1 1;-1i 7i 1 0;2 0 -20 -1;9 0 -6 -9];
A=[20 8 2;1 0 -2;3 -1 10];
% D=diag([1 2 1]);
% A=D*A/D;

%%
[r,c]=size(A);
R=zeros(r,1);
for k=1:r
   for s=1:c
      if s~=k
          R(k)=R(k)+abs(A(k,s));
      end
   end
end
x=zeros(r,1);
y=zeros(r,1);
for k=1:r
    x(k)=real(A(k,k));
    y(k)=imag(A(k,k));
end

%%
h=figure(1);
hold on
for k=1:r
    para = [x(k)-R(k), y(k)-R(k), 2*R(k), 2*R(k)];
    rectangle('Position', para, 'Curvature', [1 1]);
    plot(x(k),y(k),'r*');
end
axis equal
g=gca;
plot(g.XLim,[0 0],'k-')
plot([0 0],g.YLim,'k-')