Поиск прямолинейного участка на графике по экспериментальным данным

% Очистка переменных и рабочего пространства
clc;
close all;
clear ALL;
 
% форматы отображения данных
format shortEng
format compact
% Чтение из файла
proh_har = dlmread('HAT-CH_V5_d.txt','', 2);
% разбив массива по столбцам
Vds = proh_har(:,1);
Igs = proh_har(:,2);
Time_1 = proh_har(:,3);
Vgs = proh_har(:,4);
Ids = proh_har(:,5);
Time_2 = proh_har(:,6);
% абсолютные значения токов
Igs_abs = abs(Igs);
Ids_abs = abs(Ids);
% корень квадратный токов
Igs_abs_sqrt = sqrt(Igs_abs);
Ids_abs_sqrt = sqrt(Ids_abs);
% берем производную от тока и напряжения
Vds_diff = diff(Vds);
Ids_abs_sqrt_diff = diff(Ids_abs_sqrt);
Igs_abs_sqrt_diff = diff(Igs_abs_sqrt);
dydx = Ids_abs_sqrt_diff./Vds_diff;
% определение основных величин
Vds_max = max(Vds);
Vds_min = min(Vds);
Vgs_mean = mean(Vgs);
Igs_max = max(Igs_abs);
Ids_max = max(Ids_abs);
Ion_Ioff = Ids_max/Igs_max;
% расчет подвижности
% параметры транзистора
slope = mean(dydx);
L = 7.5e-5;
W = 2000e-6;
C = 360e-9;
mobiliti = (((slope^2)*L*2)/(W*C));
mobiliti_tab = (((dydx.^2)*L*2)/(W*C));
max_mobiliti_elem = max(mobiliti_tab)
min_mobiliti_elem = min(mobiliti_tab);
mean_mobiliti_elem = mean(mobiliti_tab);
 
    for i = 1:length(Ids_abs)
 
    end
 
% построение характеристик
semilogy(Vds, Ids_abs);
grid on;
title('Transfer characteristics'); 
xlabel('V_G_S'); 
ylabel('I_D_S'); 
legend('I_D_S');