clc
clear all
 
% graphics_toolkit("gnuplot")
N=8 
num_ch=2 
num_fig=0;
cl=['b','r','g','k','c','m','y'];
%
% name_root_dir='D:\drive-download-20180510T144415Z-001\2 канала 1-ток 2- напряжение'
%[FName, PName] = uigetfile('*.dat', 'Open File',name_root_dir,'MultiSelect', 'on' );
[FName, PName] = uigetfile('*.dat', 'Open File');
FullName =strcat(PName,FName);
 
[fid,m]=fopen(FullName,'rb');
[oscil,count]=fread(fid,inf,'int16');
fclose('all');
 
oscil=oscil./819.2;
sizeA=size(oscil);
oscil=reshape(oscil,num_ch,sizeA(1)/num_ch);
L=N*1000+750+2000;
u2=oscil(1,1:L);
i2=oscil(2,1:L);
t=[0:L-251];
t=t*0.00002;%час
U2=u2;
clear oscil;
figure('numbertitle', 'off','name','u2,i2')
plot(u2,'b')
axis ([0 L -2 2])
hold on
plot(i2,'r')
axis ([0 L -2 2])
title('u2,i2')
grid on
hold off
 
% 
for i=250:500:L-525
    ii2=cumsum(i2(i:i+525)-sum(i2(i:i+525))/525);
    try
        I2=[I2,ii2(1:500)];
    catch
    I2=ii2(1:500);
    end
end
a=zeros(1,250);
U2=[a,u2];
U2=-U2(1:L-250);%
a=zeros(1,500);
I2=[a,I2];
I2=I2*0.925;%
%
%clear i2;
%clear u2;
%
figure('numbertitle', 'off','name','U2,I2')
plot(t,U2*50,'b')
hold on
plot(t,I2,'r')
axis ([0 max(t) -100 100])
hold off
grid on
title('U2,I2 від t')
for i=1:L-250
if i<=1000
ief=0;
uef=0;
else
ief=sqrt(sum(I2(i-1000:i).^2)/1000);
uef=sqrt(sum(U2(i-1000:i).^2)/1000);
end
 
try
Ief=[Ief,ief];
Uef=[Uef,uef];
catch
Ief=ief;
Uef=uef;
end
end
disp('max I2')
max(abs(I2))
disp('max U2')
max(abs(U2))
%
figure('numbertitle', 'off','name','t,Ief,Uef')
plot(t,Ief,'r')
hold on 
plot(t,Uef,'b')
axis ([0 max(t) 0 45])
grid on
title('Ief,Uef')
hold off
 
%
P2=abs(U2.*I2);
figure('numbertitle', 'off','name','P2')
plot(t,P2,'g')
axis ([0 max(t) 0 120])
grid on
title('Потужність')
%
ed=ones(1,500);
for i=500:500:L-500
e2=cumsum(P2(i:i+500));
%e2=ed.*ee2;
try
E2=[E2,e2(1:500)];
catch
E2=e2(1:500);
end
end 
a=zeros(1,500);
E2=[E2,a];
figure('numbertitle', 'off','name','E2')
plot(t,E2,'r')
axis ([0 max(t) 0 2500])
grid on
title('Енергія за кожен період') 
hold off
figure('numbertitle', 'off','name','u2=f(i2)')
plot(abs(I2/75),abs(U2/2))
grid on
title('ВАХ') 
figure('numbertitle', 'off','name','R(t)')
R=Uef./Ief;
plot(t,R)
grid on
title('R')
figure('numbertitle', 'off','name','du/di')
A=size(I2);
for j=1:A(2)
if abs(I2(j))<0.4
r(j)=0;
else
r(j)=U2(j)/I2(j);
end
end
plot(t,r)
grid on
title('r')