clear all
close all
clc
load bifurkation.txt
data=bifurkation;

stabil_tid=5; %venter til systemet har stabilisert seg etter ventilbruk [sek]
TS=data(:,1); %samler sammen data kolonnene
% S1=data(:,2);
% S2=data(:,3);
PP1=data(:,3);
% P2=data(:,5);
Q=data(:,4);%volumstrom vann
% SP=data(:,7);
z=data(:,7);
tid=cumsum(TS); %samler sammen tidsvektor
[ant_rader,ant_kol]=size(data);

P1_max=[];%tomme vektorer
P1_min=[];
P2_max=[];
P2_min=[];

    X=[]; %tom vektor
    for i=1:length(z)-1  %gir vektor av enere hvor ventilposisjonene i z er like 
    x=isequal(z(i),z(i+1));
    X=[X x];  %lagrer null og en
    end

Ind=find(X==0); %finner indekser hvor ventilpos. er endret 
Ind=[0 Ind ant_rader]; % brukes for aa splitte opp hoveddatafil 

i=1; %initsierer
m=length(Ind); %antall forskjellige ventil-posisjoner

while i<m
ind_tid=find(abs(tid-stabil_tid)==min(abs(tid-stabil_tid))); %finner hvor vi skal begynne aa analyse data
    
DATA=data(Ind(i)+1:Ind(i+1),:);  %lagrer matrise med samme ventilposisjon
Z(i)=round(DATA(1,7)); %prosent aapninng

P1=DATA(:,3);
P2=DATA(:,4);
TS=DATA(:,1);
n=length(TS);
TS=TS(ind_tid:n);
P1=P1(ind_tid:n);
P2=P2(ind_tid:n);
P1=sgolayfilt(P1,5,21); % Savitzky Golay filter
P2=sgolayfilt(P2,5,21); % Savitzky Golay filter
tid=cumsum(TS)/1000;%TS oppgis i msek.

figure(i)
subplot(2,1,1)
plot(tid,P1)
xlabel('time [s]')
ylabel('P1 [barg]')
AXIS([0 120 0.06 0.18])
title(['',num2str(Z(i)),'  % Valve opening'])
subplot(2,1,2)
plot(tid,P2)
xlabel('time [s]')
ylabel('P2 [barg]')
AXIS([0 120 0 0.18])

a='trykk_mot_tid_data3';
filnavn=[a '_' num2str(Z(i)) '.eps'];
g=gcf;
set(g,'PaperUnits','centimeters','PaperType','A4')
opts = struct('FontMode','fixed','FontSize',9,'Width',14,'Color','rgb');
exportfig(g,filnavn,opts)
%previewfig(g,opts)

%[ZZ(i).P1_max ZZ(i).P1_min]=gjsn_bunn_topp(tid,P1,20);
%[ZZ(i).P2_max ZZ(i).P2_min]=gjsn_bunn_topp(tid,P2,20);

%P1_max=[P1_max ZZ(i).P1_max]; 
%P1_min=[P1_min ZZ(i).P1_min];
%P2_max=[P2_max ZZ(i).P2_max];
%P2_min=[P2_min ZZ(i).P2_min];

i=i+1;
end

break
save eksperimentelle_data P1_max P1_min P2_max P2_min Z
%Trykk mot ventil-posisjon:
figure(i)
subplot(2,1,1)
plot(Z,P1_max,Z,P1_min)
ylabel('P_1 [barg]')
legend('Maks trykk','Min. Trykk')
subplot(2,1,2)
plot(Z,P2_max,Z,P2_min)
xlabel('VentilÅpning % ')
ylabel('P_2 [barg]')
% g=gcf;
% set(g,'PaperUnits','centimeters','PaperType','A4')
% opts = struct('FontMode','fixed','FontSize',9,'Width',15,'Color','rgb');
% exportfig(g,'trykk_mot_ventil.eps',opts)
% % previewfig(g,opts)

% For data3.txt hvor et par negative punkter oppstod
%Fjerner disse negative pkt.
% X=find(Q<=0);
% a=Q(1:(X(1)-1));
% b=Q((X(1)+1):(X(2)-1));
% c=Q((X(2)+1):ant_rader);
% QQ=[a ;b; c];
% a=PP1(1:(X(1)-1));
% b=PP1((X(1)+1):(X(2)-1));
% c=PP1((X(2)+1):ant_rader);
% PP1=[a ;b; c];

% figure(i+1)
% plot(PP1,QQ,'*')
% xlabel('Trykk innløp føderør')
% ylabel('Volumstrøm vann innløp føderør')
% pause
% g=gcf;
% pause
% set(g,'PaperUnits','centimeters','PaperType','A4')
% opts = struct('FontMode','fixed','FontSize',9,'Width',10,'Color','rgb');
% exportfig(g,'volumstrom_mot_trykk.eps',opts)
% %previewfig(g,opts)