При запуске программы в рабочее пространство не записывается ни одна переменная

name_root_dir='D:\Folder2'
[FName, PName] = uigetfile('*.dat', 'Open File',name_root_dir,'MultiSelect', 'off' );
clc
FullName =strcat(PName,FName)
 
%данные из бинарного файла
 
[fid,m]=fopen(FullName,'rb');
[oscil,count]=fread(fid,inf,'bit16');
oscil=oscil./819.2;
fclose('all');
disp('размер массива данных')
%разделение массива
sizeA=size(oscil);
oscil_u_i=reshape(oscil,2,sizeA(1)/2);
oscil_i=oscil_u_i;
oscil_i(1,:)=[];
oscil_i=interp1([1:sizeA(1)/2],oscil_i,[1:0.5:sizeA(1)/2]);
oscil_i=oscil_i';
oscil_u=oscil_u_i;
oscil_u(2,:)=[];
oscil_u=interp1([1:sizeA(1)/2],oscil_u,[1:0.5:sizeA(1)/2]);
oscil_u=oscil_u';
sizeA=size(oscil_i);
c=0;
k=1;
i=1;
j=0;
disp('Нормирование данных')
oscil_mean_i=mean(oscil_i);
oscil_i=oscil_i-oscil_mean_i;
oscil_mean_u=mean(oscil_u);
oscil_u=oscil_u-oscil_mean_u;
oscil_max=max(oscil_i);
for v=1:sizeA(1)
if abs(oscil_i(v))>0.6*oscil_max
if oscil_i(v)>0
oscil_i(v)=0.6*oscil_max;
else
oscil_i(v)=-0.6*oscil_max;
end
end
end
 
disp
 
dd=cumsum(oscil_i);
i=1;
ff=0;
k1=0;
k2=0;
m=100;
d_z=zeros(3800,1);
dd=[dd;d_z];
oscil_i=[oscil_i;d_z];
oscil_u=[oscil_u;d_z];
while ff<1
d=mean(abs(dd(i:i+m)));
try 
d_r=[d_r,d];
catch
d_r=d;
end
if (d>20 && k1==0)
k1=i;
m=1800;
end
if (d<5 && k1>0)
k2=i;
ff=2;
end
i=i+m;
k=i;
end
 
oscil_i=oscil_i(k1-450:k2+450)*r_osc;
oscil_u=oscil_u(k1-450:k2+450)*0.25;
A=zeros(1900,1);
oscil_i=cat(1,oscil_i,A);
oscil_i=filtfilt(bL,aL,oscil_i);
oscil_u=cat(1,oscil_u,A);
oscil_u=filtfilt(bL,aL,oscil_u);
oscil_max=size(oscil_i);
kk=oscil_max(1);
disp('')
disp('')
oscil_int=cumsum(oscil_i);
oscil_int=oscil_int-mean(oscil_int);
AB=size(oscil_int);
a=1;
b=AB(1);
for i=1:AB(1)
if abs(oscil_int(i))<50
oscil_int(i)=0;
end
if abs(oscil_u(i))<0.05
oscil_u(i)=0;
end
end
oscil_int=oscil_int;
disp('')
%начальные значения R
Rmin=0;
Rmax=0.5;
R=(Rmax-Rmin)/2;
 
%главное окно GUI
hwin=figure('NumberTitle','off',...
'Name','Resist', ...
'Tag','mainWindow',...
'MenuBar','figure',...
'HandleVisibility','callback',...
'IntegerHandle','off',...
'Renderer','painters',...
'Units', 'pixels',...
'HandleVisibility','on',...
'ResizeFcn',@figResize);
 
h = get(hwin,'OuterPosition');
%кнопка выхода
uicontrol('Style','PushButton',...
'Tag','btnExit',...
'Parent',hwin,...
'Position',[h(3)-100 20 80 30],...
'String','Exit',...
'CallBack', 'close');
h = get(hwin,'OuterPosition');
%кнопка расчета сопротивления
uicontrol('Style','PushButton',...
'Tag','btnResist',...
'Parent',hwin,...
'Position',[h(3)-200 20 80 30],...
'String','R',...
'CallBack', @btnResist_value);
uicontrol('Style','Slider',...
'Tag','btnBegin',...
'Parent',hwin,...
'Position',[h(3)-10 60 20 h(4)-10],...
'Min',1,'Max',20000,'Value',1,...
'SliderStep',[0.0001 0.1],...
'CallBack', @btnBegin_value);
uicontrol('Style','Slider',...
'Tag','btnEnd',...
'Parent',hwin,...
'Position',[h(3)-40 60 20 h(4)-200],...
'Min',1,'Max',20000,'Value',20000,...
'SliderStep',[0.0001 0.1],...
'CallBack', @btnEnd_value);
%слайдер радиуса
uicontrol('Style','Slider',...
'Tag','btnCircle',...
'Parent',hwin,...
'Position',[40 5 floor(h(3)-260) 20],...
'Min',Rmin,'Max',Rmax,'Value',R,...
'SliderStep',[0.0001 0.01],...
'CallBack', @btnCircle_value);
 
%
hA = axes('Parent',hwin,...
'Tag','Axes_UI',...
'Units', 'pixels',...
'OuterPosition', [0 20 floor(h(3)-260) floor(h(4)-150)],...
'XGrid','on',...
'YGrid','on');
 
set(get(hA,'XLabel'),'String','i, kA',...
'FontName','times',...
'FontAngle','italic',...
'FontSize',14)
set(get(hA,'YLabel'),'String','u, V',...
'FontName','times',...
'FontAngle','italic',...
'FontSize',14) 
dat.a=a;
dat.b=b;
dat.u=oscil_u;
dat.i=oscil_int;
dat.R=R;
dat.Rmin=Rmin;
dat.Rmax=Rmax;
figure(hwin)
axes(hA)
i=dat.i;
u=dat.u;
axes(hA)
plot(hA,i(a:b)./1000,u(a:b))
set(hA,'XGrid','on','YGrid','on')
child_handles = allchild(hwin);
guidata(hwin,dat);
guihandles(hwin);
 
 
function figResize(src,evt)
get(gca,'Children');
allchild(gca);
handles = guihandles(src);
dat = guidata(src);
try
R=dat.R;
Rmin=dat.Rmin;
Rmax=dat.Rmax;
catch
Rmin=0;
Rmax=2;
R=(Rmax-Rmin)/2;
end
h = get(handles.mainWindow,'OuterPosition');
hF = get(0, 'CurrentFigure');
hA = get(hF, 'CurrentAxes');
set(hA,'Units', 'pixels')
set(handles.btnExit,'Position',[h(3)-99 20 80 30]);
set(handles.btnResist,'Position',[h(3)-200 20 80 30]);
set(handles.btnBegin,'Position',[h(3)-100 60 20 h(4)-180]);
set(handles.btnEnd,'Position',[h(3)-40 60 20 h(4)-180]);
set(handles.Text_1,'Position',[h(3)-120 h(4)-110 50 20]);
set(handles.Text_2,'Position',[h(3)-60 h(4)-110 50 20]);
set(handles.Text_k,'Position',[h(3)-200 h(4)-110 50 20]);
set(handles.Text_R,'Position',[floor(R/(Rmax-Rmin)*(floor(h(3)-260))) 30 60 20]);
set(handles.btnCircle,'Position',[20 5 floor(h(3)-260) 20]);
set(hA,'Position',[20 90 floor(h(3)-260) floor(h(4)-150)]);
 
function [a]=btnBegin_value(src,~)
clc
handles = guihandles(src);
dat = guidata(src);
hA = gca;
a=floor(get(src,'Value'));
dat.a=a;
b=dat.b;
guidata(src,dat);
set(handles.Text_1,'String',num2str(a));
set(handles.btnEnd,'Min',a+1);
i=dat.i;
u=dat.u;
axes(hA) 
plot(hA,i(a:b)./1000,u(a:b))
set(gca,'XGrid','on','YGrid','on')
set(get(hA,'XLabel'),'String','i, kA',...
'FontName','times',...
'FontAngle','italic',...
'FontSize',14)
set(get(hA,'YLabel'),'String','u, V',...
'FontName','times',...
'FontAngle','italic',...
'FontSize',14) 
 
function [b]=btnEnd_value(src,evt)
handles = guihandles(src);
dat=guidata(src);
hA = gca;
set(src,'Min',dat.a);
b=floor(get(src,'Value'));
a=dat.a;
dat.b=b;
guidata(src,dat);
set(handles.Text_2,'String',num2str(b))
set(handles.btnBegin,'Max',b-1);
i=dat.i;
u=dat.u;
axes(hA)
plot(hA,i(a:b)./1000,u(a:b));
set(gca,'XGrid','on','YGrid','on')
set(get(hA,'XLabel'),'String','i, kA',...
'FontName','times',...
'FontAngle','italic',...
'FontSize',14)
set(get(hA,'YLabel'),'String','u, V',...
'FontName','times',...
'FontAngle','italic',...
'FontSize',14) 
 
function btnCircle_value(src,evt)
handles = guihandles(src);
dat=guidata(src);
h = get(handles.mainWindow,'OuterPosition');
try
p=dat.p;
delete(p)
catch
end
hF = get(0, 'CurrentFigure');
hA = get(hF, 'CurrentAxes');
hA = gca;
R=get(src,'Value');
Rmin=dat.Rmin;
Rmax=dat.Rmax;
set(handles.Text_R,'Position', [floor(R/(Rmax-Rmin)*(floor(h(3)-250))) 30 60 20]);
set(handles.Text_R,'String',num2str(R));
axes(hA)
hold on
[X Y]=calculateEllipse(0,0,R,R,0,360);
p=plot(X,Y,'r');
dat.p=p;
dat.R=R;
guidata(src,dat);
hold off
set(gca,'XGrid','on','YGrid','on')
set(get(hA,'XLabel'),'String','i, kA',...
'FontName','times',...
'FontAngle','italic',...
'FontSize',14)
set(get(hA,'YLabel'),'String','u, V',...
'FontName','times',...
'FontAngle','italic',...
'FontSize',14) 
 
function [X Y] = calculateEllipse(x, y, a, b, angle, steps)
%# This functions returns points to draw an ellipse
%#
%# @param x X coordinate
%# @param y Y coordinate
%# @param a Semimajor axis
%# @param b Semiminor axis
%# @param angle Angle of the ellipse (in degrees)
%#
 
error(nargchk(5, 6, nargin));
if nargin<6, steps = 36; end
 
beta = -angle * (pi / 180);
sinbeta = sin(beta);
cosbeta = cos(beta);
 
alpha = linspace(0, 360, steps)' .* (pi / 180);
sinalpha = sin(alpha);
cosalpha = cos(alpha);
 
X = x + (a * cosalpha * cosbeta - b * sinalpha * sinbeta);
Y = y + (a * cosalpha * sinbeta + b * sinalpha * cosbeta);
 
if nargout==1, X = [X Y]; end
 
function [b]=btnResist_value(src,evt)
 
handles = guihandles(src);
dat=guidata(src);
h = get(handles.mainWindow,'OuterPosition');
begin_R=dat.a
end_R=dat.b;
R=dat.R;
i=dat.i./1000;
u=dat.u;
 
figure
plot(i(begin_R:end_R),u(begin_R:end_R))
line(i(begin_R:end_R),u(begin_R:end_R),'Color',[1 0 0],'LineWidth',2)
hold on
for k=begin_R+1:end_R
 
if abs(i(k))>0
A=(u(k)-u(k-1))/(i(k)-i(k-1));
xdata=i(k:k+10,1);
ydata=u(k:k+10,1);
I=mean(xdata);
 
x0=[1, 1,];
f= @(x,xdata)x(1).*xdata+x(2);
[x,resnorm] = lsqcurvefit(f,x0,xdata,ydata);
U=x(1)*I+x(2);
dYdX=x(1);
%dYdX=A;
a=R;
%if i(k-2)<i(k+2)
%dYdX=x(1); 
phi= 2*atan((a^2 - U^2 + (a^4*dYdX^2 + a^4 - 2*a^2*dYdX^2*I^2 + 4*a^2*dYdX*I*U - 2*a^2*U^2 + dYdX^2*I^4 + dYdX^2*I^2*U^2 - 2*dYdX*I^3*U - 2*dYdX*I*U^3 + I^2*U^2 + U^4)^(1/2) + dYdX*I*U)/(dYdX*a^2 - dYdX*I^2 + U*I));
b=((U - dYdX*I)*(-(- a^2 + I^2 + U^2)/(a^2*dYdX^2 + a^2 - dYdX^2*I^2 + 2*dYdX*I*U - U^2))^(1/2));
try
phi_m=[phi_m,phi];
b_m=[b_m,b];
dYdX_m=[dYdX_m,dYdX];
catch 
phi_m=phi;
b_m=b;
dYdX_m=dYdX;
end
set(handles.Text_k,'String',num2str(x(1)));
else
A=-(u(k)-u(k-1))/(i(k)-i(k-1));
xdata=i(k:k+10,1);
ydata=u(k:k+10,1);
I=mean(xdata);
 
x0=[1, 1,];
f= @(x,xdata)x(1).*xdata+x(2);
[x,resnorm] = lsqcurvefit(f,x0,xdata,ydata);
U=x(1)*I+x(2);
dYdX=x(1);
%dYdX=A;
a=R; 
phi= -2*atan((U^2 - a^2 + (a^4*dYdX^2 + a^4 - 2*a^2*dYdX^2*I^2 - 4*a^2*dYdX*I*U - 2*a^2*U^2 + dYdX^2*I^4 + dYdX^2*I^2*U^2 + 2*dYdX*I^3*U + 2*dYdX*I*U^3 + I^2*U^2 + U^4)^(1/2) + dYdX*I*U)/(- dYdX*a^2 + dYdX*I^2 + U*I))
phi=2*atan((a^2 - U^2 + (a^4*dYdX^2 + a^4 - 2*a^2*dYdX^2*I^2 - 4*a^2*dYdX*I*U - 2*a^2*U^2 + dYdX^2*I^4 + dYdX^2*I^2*U^2 + 2*dYdX*I^3*U + 2*dYdX*I*U^3 + I^2*U^2 + U^4)^(1/2) - dYdX*I*U)/(- dYdX*a^2 + dYdX*I^2 + U*I)) 
b=((U + dYdX*I)*((- a^2 + I^2 + U^2)/(-a^2*dYdX^2 - a^2 + dYdX^2*I^2 + 2*dYdX*I*U + U^2))^(1/2)) ; 
 
end
try
phi_m=[phi_m,phi];
b_m=[b_m,b];
dYdX_m=[dYdX_m,dYdX];
catch 
phi_m=phi;
b_m=b;
dYdX_m=dYdX;
end
[X Y] = calculateEllipse(0,0, R, b, phi*180/pi,360);
plot(X,Y,'b')
 
 
 
end
 
R_button=tan(pi/2-phi_m)/1000;
R_button=R_button';
sizeR=size(R_button);
for k=1:sizeR(1)
if abs(R_button(k))>0.5
R_button(k)=0;
end
end
t=(0:0.00001:(sizeR(1)-1)*0.00001)';
u=u(1:sizeR(1));
i=i(1:sizeR(1));
E=cumsum(u.*i);
save ('phi','R_button','u','i','E','t','sizeR')
disp('OK')
[X Y] = calculateEllipse(0,0, R, b, phi*180/pi,360);
plot(X,Y,'b') 
X0 = 1:length(u);
% интерполировать в X, Y 
X = linspace(min(X), max(X), 100);
Y = interp1(X0, Y0, X);
% найдём производную dy/dx = (y(x+dx) - y(x-dx)) / (2*dx)
dX = X(2:end-1);
dY = (Y(3:end) - Y(1:end-2)) ./ (X(3:end) - X(1:end-2));