警告从;
% gerschgorin.m
一个=魔法(5);
gerschgorin (A)
函数gerschgorin (A)
如果大小(1)~ =大小(A, 2)
错误(“La矩阵deve essere quadrata。”);%的矩阵必须是正方形。
返回;
结束
图1:%计算和代表与矩阵相关的圈子。
图;
为i = 1:尺寸(1)
%圆的中心就是h (h, k)的实部(我)和k的虚部(我)。
h =真正的((我));
k =图像放大((我));
%发现圆的半径。
r = 0;
为j = 1:尺寸(1)
如果我~ = j
r = r + abs ((i, j));
结束
结束
t = 0:0.01:2 *π;
情节(r * cost + h r * sin (t) + k,“- - -”);%使用极坐标。
持有在;
c =情节(h, k,“波”);
结束
%绘制的实际特征值矩阵。
ev = eig(一个);
牧师=情节(真正的(ev),图像放大(ev),“罗”);
参照线(0);
yline (0);
轴平等的;
传奇([c]牧师,“Centri”,“Autovalori”);
网格在;
包含(“杰克逊Reale”);
ylabel (“杰克逊Immaginaria”);
标题(“矩阵A的圈子”);
%图2:计算和代表与矩阵相关的圈子^ T。
图;
为i = 1:尺寸(1)
h =真正的((我));
k =图像放大((我));
r = 0;
为j = 1:尺寸(1)
如果我~ = j
r = r + abs ((j, i));
结束
结束
t = 0:0.01:2 *π;
情节(r * cost + h r * sin (t) + k,“- - -”);
持有在;
c =情节(h, k,“波”);
结束
%绘制的实际特征值矩阵。
ev = eig(一个);
牧师=情节(真正的(ev),图像放大(ev),“罗”);
参照线(0);
yline (0);
轴平等的;
传奇([c]牧师,“Centri”,“Autovalori”);
网格在;
包含(“杰克逊Reale”);
ylabel (“杰克逊Immaginaria”);
标题(“矩阵^ T圈”);
%图3:十字路口联盟的圆和一个^ T。
图;
为i = 1:尺寸(1)
h =真正的((我));
k =图像放大((我));
r1 = 0;
r2 = 0;
x = [];
y = [];
为j = 1:尺寸(1)
如果我~ = j
r1 = r1 + abs ((i, j));
r2 = r2 + abs ((j, i));
x1 = r1 * cost + h;
日元= r1 * sin (t) + k;
x2 = r2 * cost + h;
y2 = r2 * sin (t) + k;
x = [x, x1);
y = [y, y₁];
结束
结束
polyout1 = polyshape (x, y);
polyout2 = polyshape (x2, y2);
十字路口=相交(polyout1 polyout2);
情节(十字路口,“EdgeColor”,“红色”);
持有在;
c =情节(h, k,“波”);
结束
%绘制的实际特征值矩阵。
ev = eig(一个);
牧师=情节(真正的(ev),图像放大(ev),“罗”);
参照线(0);
yline (0);
轴平等的;
网格在;
包含(“杰克逊Reale”);
ylabel (“杰克逊Immaginaria”);
标题(十字路口的联盟的圆和一个^ T ');
传奇([c]牧师,“Centri”,“Autovalori”);
结束