Hi Federica!
In the given code the highlighted part code is used create a vector 't' of 'N+1' equally spaced points around the circumference of a circle. Here's what each line does:
- N = 256;: This line sets the value of N to 256. This determines the number of points we want to create on the circle.
- t = (0:N)*2*pi/N;: This line creates the vector t by performing the following steps:
- (0:N) creates a row vector from 0 to N, where each element represents the index of a point on the circle.
- 2*pi/N scales the vector by multiplying each element by 2*pi/N. This step ensures that the points span the full circumference of the circle, from 0 to 2*pi.
So, the resulting vector 't' contains 'N+1' equally spaced points around the circle, covering the range from 0 to 2*pi. These points will be later used to map each point on the circumference of the circle when plotting the Gershgorin Discs.
To plot the row circles and column circles seperately you can introduce two separate loops: one for rows and one for columns.
a sample code would look like this
% Plot row circles
for i = 1:size(A,1)
h = real(A(i,i));
k = imag(A(i,i));
r = 0;
for j = 1:size(A,1)
if i ~= j
r = r + norm(A(i,j));
end
end
N = 256;
t = (0:N)*2*pi/N;
plot(r*cos(t) + h, r*sin(t) + k, '-');
hold on;
plot(h, k, '+');
end
Similarly for the column circles Change the
for i = 1:size(A,2)
...
for j = i:size(A,2)
...
end
end
Hope this helps,
Thank you!
