%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Labskelett för laboration 1b %%%
%%% TSFS05 - Fordonssystem       %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%
%  Denna fil innehåller ett labbskelett som är lämpligt att använda i
%  laboration 1b. Genom att skriva kommandona här istället för vid
%  prompten blir beräkningarna lätta att reproducera och lätta att ändra
%  vid fel. För att hjälpa er att komma igång har vi redan fyllt i ett
%  antal av de kommandon som behövs för uppgift 1. 
%
%  Kontrollera att ni använder rätt enheter! SI-enheter är att föredra i
%  beräkningar, medan andra enheter kan vara att föredra i plottar. Till
%  exempel är det ofta lättare att tyda figurer med vevaxelvinkel angiven
%  i grader än i radianer.
%
%  Glöm inte att kontrollera att de resultat ni får är rimliga! Något som 
%  är ganska lätt att kontrollera är till exempel cylindervolymen.
%  Överensstämmer den slagvolym ni får med den slagvolym som anges i
%  laborationskompendiet och är kompressionsförhållandet rätt? 
%
%  (Filen uppdaterad 2005-09-19)
%
%



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%         
%%%          UPPGIFT 1          %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


load <filnamn>  % Ladda in lämplig mätdatafil

pcyl = ;        % Cylindertryck i [Pa]   
M    = ;        % Moment [Nm]
N    = ;        % Varvtal i [rpm]
ca   = ;        % Vevaxelvinkel i [rad]


cykelNr=1:size(<filnamn>.N,2);   % Välj vilka cykler som ska plottas

figure(1); clf; 
plot(ca*180/pi,pcyl(:,cykelNr)*10^-5);
title('Cylindertryck som funktion av vevaxelvinkel');
xlabel('Vevaxelvinkel [grader]');
ylabel('Tryck [bar]');


% För att lagra figuren som en praktisk eps kan du med fördel använda kommandot:
% print -f1 -depsc2 filnamn.eps
%    -f1    : Figur nummer 1
%    -depsc2: Format eps 2 med färg (c)
% För att lagra filen som jpeg kan du använda:
% print -djpeg filnamn.jpg
% Är du osäker så prova skriv: help print
% I hjälpen finns även andra tips om filformat mm.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%          UPPGIFT 2          %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  

l     = 0.158; 	% Vevstakens längd [m]
a     = 0.044;	% Halva slaglängden [m]
B     = 0.068;	% Cylinderdiameter [m]
ncyl  = 5;      % Antalet motorcylindrar
VD    = 1.6e-3; % Total volym [m³] 

load <filnamn>  % Ladda in lämplig datafil


figure(4); clf; 
plot(...)
title('Momentant moment som funktion av vevaxelvinkel');
xlabel('Vevaxelvinkel []');
ylabel('M [Nm]');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%          UPPGIFT 3          %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



Wcalc = ; % Beräknat arbete
Mcalc = ; % Beräknat moment
Mmeas = ; % Uppmätt medelmoment

disp(['W - beräknat: ', num2str(Wcalc)]) 
disp(['M - beräknat: ', num2str(Mcalc)]) 
disp(['M - uppmätt:  ', num2str(Mmeas)]);


figure(5); clf; 
plot(V*1e3,pcyl*1e-5)
title('pV-diagram (alla cykler)');
xlabel('Volym [dm^3]');
ylabel('Tryck [bar]');


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%          UPPGIFT 4          %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%          UPPGIFT 5          %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%clear all

%%% Motordata för TURBO-motorn

% Tips: Använd gärna funktionen sfc_plot


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%          UPPGIFT 6          %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Tips: Använd gärna funktionen nvol_plot

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%          UPPGIFT 7          %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

