Экстраполяция к теоретической зависимости

% очищаем все, что можно
clear all;
clc;
clf;
 
% импортируем данные из файла
Q = importdata('Prakt323.csv');
% вытягиваем оттуда только массив данных
S = Q.data;
% сохраняем данные по своим масивам
U = S(:,1)/1000; % весь первый столбец - напряжение, В
I = S(:,2)/1000; % весь второй столбец - ток, А
dI = S(:,3); % весь третий столбец 
 
% массив положительных значений
x = U(U>=0.0); 
y = I(find(U>=0.0));
% рисуем положительную ветку ВАХ
plot(x,y,'b') 
 
% начинаем усиленно интерполировать
ftype = fittype('a*(exp(b*x)-1)'); % создаем свою ф-цию для интерполяции
fresult1 = fit(x,y,ftype) % интерполируем по оной
% получаем массив результатов интерполяции
fy1 = fresult1.a*(exp(fresult1.b*x)-1); % получаем массив результатов интерполяции
 
In = fresult1.a; %ток насыщения
 
fresult2 = fit(x,y,'exp2') % интерполируем по сумме экспонент
% получаем массив результатов интерполяции
fy2 = fresult2.a*exp(fresult2.b*x)+ fresult2.c*exp(fresult2.d*x);
 
% рисуем все на одном графике
hold on
plot(x,fy1,'*r') % интерполяция своей ф-цией красным
plot(x,fy2,'.g') % интерполяция суммой эксп. - зеленым
grid on
 
disp(['Ток насыщения: ',num2str(In),'А'])