% steady state plot 17 dec

clear,clc,clf,close all,

load('xPC_Data_20071217_202522.mat');
xpcdata6.Blow=(-0.0655*xpcdata6.Blow.^2+5.1865*xpcdata6.Blow); %omregning

%% Til plot af steadystate
% figure(5)
% plot(xpcdata6.time,[xpcdata6.Tr0xF8g,xpcdata6.Temp,xpcdata6.Retur,...
%     xpcdata6.Stok,xpcdata6.Blow,xpcdata6.Lambda,...
%     xpcdata6.Fan],...
%     'Linewidth',1.5);
% 
% xlabel('Tid [s]','FontSize',12')
% ylabel('Amplitude','FontSize',12')
% 
% xlim([2600 3049]);
% ylim([0 115]);
% 
% legend({'T_{rg} [{\circ}C]','T_{frem} [{\circ}C]',...
%     'T_{retur} [{\circ}C]','Stoker [%]','Blæser [%]','\lambda [%]',...
%     'Kalorifer [Hz]','Temp Forskel'},...
%     'FontSize',12','location','EO')
% grid on;

%Udregner middelværdierne ved steady state
% AvgStoker = xpcdata6.Stok(27000:29490);
% x1 = mean(AvgStoker);
% disp('Middelværdi stoker (2700-2949 sek) = '), disp(x1);
% 
% AvgTemp = xpcdata6.Temp(27000:29490);
% x1 = mean(AvgTemp);
% disp('Middelværdi fremløbstemperatur (2700-2949 sek) = '), disp(x1);
% 
% AvgRetur = xpcdata6.Retur(27000:29490);
% x1 = mean(AvgRetur);
% disp('Middelværdi returvand (2700-2949 sek) = '), disp(x1);
% 
% AvgRg = xpcdata6.Tr0xF8g(27000:29490);
% x1 = mean(AvgRg);
% disp('Middelværdi røggastemperatur (2700-2949 sek) = '), disp(x1);
% 
% AvgBlow = xpcdata6.Blow(27000:29490);
% x1 = mean(AvgBlow);
% disp('Middelværdi blæser (2700-2949 sek) = '), disp(x1);
% 
% AvgLambda = xpcdata6.Lambda(27000:29490);
% x1 = mean(AvgLambda);
% disp('Middelværdi Lambda (2700-2949 sek) = '), disp(x1);

% 
% line([2700 2700],[0 150], 'Color','k','LineWidth',2)
% line([2949 2949],[0 150], 'Color','k','LineWidth',2)
%Saving plot
% print -depsc2 ../../rapport/Fig/StokerSteadyState.eps

%% Til plot af steps til verification af dynamik - blower step

% figure(6)
% plot(xpcdata6.time,[xpcdata6.Tr0xF8g,xpcdata6.Temp,xpcdata6.Retur,...
%     xpcdata6.Stok,xpcdata6.Blow,xpcdata6.Lambda,...
%     xpcdata6.Fan],...
%     'Linewidth',1.5);
% 
% xlabel('Tid [s]','FontSize',12')
% ylabel('Amplitude','FontSize',12')
% 
% xlim([2900 4350]);
% ylim([0 125]);
% 
% legend({'T_{rg} [{\circ}C]','T_{frem} [{\circ}C]',...
%     'T_{retur} [{\circ}C]','Stoker [%]','Blæser [%]','\lambda [%]',...
%     'Kalorifer [Hz]','Temp Forskel'},...
%     'FontSize',12','location','EO')
% grid on;

% line([10720 10720],[0 150], 'Color','k','LineWidth',1.5)
% line([12720 12720],[0 150], 'Color','k','LineWidth',1.5)
% line([14720 14720],[0 150], 'Color','k','LineWidth',1.5)
% line([16720 16720],[0 150], 'Color','k','LineWidth',1.5)

%Saving plot
% print -depsc2 ../../rapport/Fig/StokerBlowStepMeas.eps

%% Til plot af steps til verification af dynamik - stoker step
% også udsnit af plot til zoom på lambda osc.
% figure(7)
% subplot(2,1,1),plot(xpcdata6.time(29000:80000)-2900,[xpcdata6.Tr0xF8g(29000:80000),xpcdata6.Temp(29000:80000),xpcdata6.Retur(29000:80000),...
%     xpcdata6.Stok(29000:80000),xpcdata6.Blow(29000:80000),xpcdata6.Lambda(29000:80000),...
%     xpcdata6.Fan(29000:80000)],...
%     'Linewidth',1.5);
% 
% xlabel('Tid [s]','FontSize',12')
% ylabel('Amplitude','FontSize',12')
% 
% xlim([0 5100]);
% ylim([0 135]);
% 
% legend({'T_{rg} [{\circ}C]','T_{frem} [{\circ}C]',...
%     'T_{retur} [{\circ}C]','Stoker [%]','Blæser [%]','\lambda [%]',...
%     'Kalorifer [Hz]'},...
%     'FontSize',12','location','EO')
% grid on;
% 
% %Boks omkring lambda oscillering
% line([2480 2480],[1 11], 'Color','k','LineWidth',1.5)
% line([3520 3520],[1 11], 'Color','k','LineWidth',1.5)
% line([2480 3520],[1 1], 'Color','k','LineWidth',1.5)
% line([2480 3520],[11 11], 'Color','k','LineWidth',1.5)
% 
% %%Plot af lambda osc.
% subplot(2,1,2),plot(xpcdata6.time(29000:80000)-2900,xpcdata6.Lambda(29000:80000),...
%     'Linewidth',1.5,'Color',[0.75,0.75,0]);
% 
% xlabel('Tid [s]','FontSize',12')
% ylabel('Amplitude','FontSize',12')
% 
% legend({'\lambda [%]     '},...
%     'FontSize',12','location','EO')
% 
% xlim([2480 3520]);
% ylim([1 11]);
% grid on;
% 
% %Saving lambda osc. plot
% % print -depsc2 ../../rapport/Fig/LambdaOsc.eps
% 
% 
% % line([10720 10720],[0 150], 'Color','k','LineWidth',1.5)
% % line([12720 12720],[0 150], 'Color','k','LineWidth',1.5)
% % line([14720 14720],[0 150], 'Color','k','LineWidth',1.5)
% % line([16720 16720],[0 150], 'Color','k','LineWidth',1.5)
% 
% %Saving plot
% % print -depsc2 ../../rapport/Fig/StokerStokStepMeas.eps
% 
% %% FFT af Lambda osc.
% 
% %omskrivning af eksempel i Matlab hjælp funktionen
% signal = xpcdata6.Lambda(70000:90000);
% L = length(signal);
% Fs= 10;
% NFFT = 2^nextpow2(L); % Next power of 2 from length of y
% signal2=signal-mean(signal);
% Y = fft(signal2,NFFT)/L;
% f = Fs/2*linspace(0,1,NFFT/2);
% 
% % Plot single-sided amplitude spectrum.
% figure(2)
% plot(f,2*abs(Y(1:NFFT/2)),'LineWidth',1.5,'Color','k') 
% title('Enkelt-Sidet amplitude spektrum of \lambda_{rg}','FontSize',12)
% xlabel('Frekvens [Hz]','FontSize',12)
% ylabel('Amplitude |Y(f)|','FontSize',12)
% 
% xlim([0 0.05]);
% ylim([0 2.5]);
% grid on;

%Saving plot
% print -depsc2 ../../rapport/Fig/LambdaOscFFT.eps