%GTopt.m
global Z

format short g
format compact
%format long g

%load results/GTopt/optGT;
load optGT;

%Pådrag 
    %u = [Ws_CH4kompressor; Ws_LUFTkompressor]
    %u = [8.2E5; 2.6E8];     % J/s (verdier hentet fra høstprosjektet)

%Forstyrrelser
    %d = [p_CH4; p_LUFT; T_CH4; T_LUFT]
    d = [20E5; 1E5; 283; 283];  %Pa (d(1:2)) og K (d(3:4)) (verdier hentet fra Skogestad(2000))
        
%Initialbetingelser i brennkammeret 
    %X_BK =[n_O2; n_N2; n_CO2; n_H20(g); n_CH4; T_BK]
        
    uX = uXGT;  %Henter initialverdier fra tidligere optimalisering
    uXsc = [1E4; 1E7; 1; 1; 1; 1; 1E1]; %Skaleringsvektor
    uX = uX./uXsc;                      %Skalering

%Options for optimalisering (fmincon)
    options = optimset;
    options.Display = 'iter';
    options.Diagnostics = 'on';
    options.LargeScale = 'off';
    options.GradConstr = 'off';
    options.GradObj = 'off';
    options.MaxFunEvals = 10000;
    options.MaxIter = 2000;

%fmincon
    [uX,fval,exitflag,output,lambda,grad,hessian] = fmincon('GTfmincon',uX,[],[],[],[],[],[],'GTnonlcon',options,d,uXsc)

%d82=d; 
%uX82=uX.*uXsc; 
%save results/GTopt/opt82 d82 uX82 exitflag fval %01 henviser til forstyrrelse og opt.nr. (ingen forstyrrelse og første opt.)
                                                              %11 betyr redusert metantrykk og første opt.
Z