Ось фигуры на графике

sum(abs(diff(complex(x,y))))

m = Q(1);                    % текущее максимальное значение 
for i=1:length(Q)            % цикл от 1 до конца вектора с 
                             % шагом 1 (по умолчанию) 
    if m < Q(i)                   % если a(i) > m, 
         m = Q(i);           % то m = a(i) 
    end 
end                          % конец цикла for

clear all; close all; clc;
load('matlab.mat')
x = Q(1:2:length(Q)-1);
y = Q(2:2:length(Q));
 
ax(1) = min(x);
ax(2) = max(x);
ay(1) = min(y);
ay(2) = max(y);
 
L = sum(abs(diff(complex(x,y))))
a = sum(abs(diff(complex(ax,ay))))
 
b = 2*(L/pi - a/2)

clear all; close all; clc;
load('matlab.mat')
load('matlab1.mat', 'D')
 
ax(1) = min(Q);
ax(2) = max(Q);
ay(1) = min(D);
ay(2) = max(D);
 
L = sum(abs(diff(complex(Q,D))))
a = sum(abs(diff(complex(ax,ay))))
 
b = 2*(L/pi - a/2)
plot(Q,D,ax,ay)

clear all; close all; clc;
load('matlab.mat')
load('matlab1.mat', 'D')
 
j = 1;
for i = 1:length(Q)
    if or(Q(i) <= -5.2, Q(i) >= 1.8)
        ind(j) = i;
        j = j+1;
    end
end % удалил странные значения посредине. искажали решение
 
Q = Q(ind);
D = D(ind);
Q(length(Q)+1) = Q(1);
D(length(D)+1) = D(1);
 
ax(1) = min(Q);
ax(2) = max(Q);
ay(1) = min(D);
ay(2) = max(D);
 
L = sum(abs(diff(complex(Q,D))))
a = sum(abs(diff(complex(ax,ay)))); % большая ось
a = a/2 % большая полуось
 
b = 1:200;
bb = fzero(@(b)(4*(pi*a*b + (a-b)^2)/(a+b)-L),[0,200]); % малая полуось
% была неправильная формула для длины
b = bb*2 % малая ось
 
plot(Q,D,ax,ay)
grid on