%SKRYPT start.m DO URUCHAMIANIA
%MODELU "tarcie.mdl"
%SKRYPT UMOŻLIWA ANALIZĘ POPRAWNOŚCI ZAMODELOWANIA
%SIŁ DZIAŁAJĄCYCH NA OBIEKT

%**********************************************************
%DANE DO NORMLNEJ SYMULACJI
P=10			%przyłożona siła
m1=3;			%masa suwaka
m2=2;			%masa walca
mi=0.5; 		%współczynnik tarcia kinetycznego
b=5;			%współczynnik oporów "kwadratowych"
g=10;			%przyśpieszenie ziemskie
R=0.1;		%promień walca
f=0.01		%współczynnik tarcia tocznego
v0=0;			%prędkość początkowa 
x0=0;			%położenie początkowe 

%**************************************************************
% symulacja "normalnego" układu:
[t,x,y]=sim('tarcie',[0 25]);
figure(1)		%nowy wykrs
plot(t,y);		%naniesienie wyników na wykres - 
					%czas(wektor t), położenie(1 kolumna macierzy y), prędkość(2 kol. macierzy y)
grid;				%rysowanie siatki
%axis([0 17.5 -10 10]);%skalowanie osi - axis([Xmin Xmax Ymin Ymax])
title('Zmiana położenia suwaka i walca w czasie');%tytuł wykresu
xlabel('czas [s]'); 			%etyxieta osi poziomej
ylabel('położenie [m]');	%etykieta osi pionowej
legend('położenie suwaka','położenie walca');	%wyświetla legendę charakterystyk na wykresie
													%w kolejności kolumn macierzy y
                                       
%****************************************************************************                                       
%zmiana parametrów - zerowe tłumienia
b=0;
mi=0;
f=0;
[t,x,y]=sim('tarcie',[0 25]);
figure(2)	
plot(t,y);	
grid;	
title('Zmiana położenia suwaka i walca bez tłumienia');
xlabel('czas [s]');
ylabel('położenie [m]');
legend('położenie suwaka','położenie walca');	

%***********************************************************************************
%zmiana przyłozonej siły
b=5;
mi=0.5;
f=0.01;
P=3;
[t,x,y]=sim('tarcie',[0 50]);
figure(3)	
plot(t,y(:,1),'g'); %wyświetla tylko pierwszą kolumnę macierzy y w kolorze zielonym - 'g'reen 
hold on;
plot(t,y(:,2),'b'); %b-niebieski
P=30;
[t,x,y]=sim('tarcie',[0 50]);
hold on;			%zachowując stary wykres umożliwia naniesienie na układ współrzędnych kolejnych wyników
plot(t,y(:,1),'r');  %r-czerwony
hold on;
plot(t,y(:,2),'m');   %m-purpurowy
grid;
title('Porównanie zmian przyłozonej siły P');
xlabel('czas [s]');
ylabel('położenie [m]');
legend('suwak - P=3','walec - P=3','suwak - P=30','walec - P=30');
hold off;
%*********************************************************
% badanie tarcia
P=10;
mi=0.2;
[t,x,y]=sim('tarcie',[0 50]);
figure(4)	
plot(t,y(:,1),'g');% wyświetla tylko pierwszą kolumnę macierzy y w kolorze zielonym - 'g'reen 
hold on;
plot(t,y(:,2),'b');
mi=0.8;
[t,x,y]=sim('tarcie',[0 50]);
hold on;
plot(t,y(:,1),'r');
hold on;
plot(t,y(:,2),'m');
grid;
title('Porównanie zmian współczynnika tarcia');
xlabel('czas [s]');
ylabel('położenie [m]');
legend('suwak - mi=0.2','walec - mi=0.2','suwak - mi=0.8','walec - mi=0.8');
hold off;
%*********************************************************
%zmiana masy obiektów
mi=0.5;
m1=4;
m2=2;
[t,x,y]=sim('tarcie',[0 50]);
figure(5)	
plot(t,y(:,1),'g');
hold on;
plot(t,y(:,2),'b');
m1=40;
m2=20;
[t,x,y]=sim('tarcie',[0 50]);
hold on;
plot(t,y(:,1),'r');
hold on;
plot(t,y(:,2),'m');
grid;
title('Porównanie zmian masy obiektów');
xlabel('czas [s]');
ylabel('położenie [m]');
legend('położenie suwaka dla m1=4 i m2=2','położenie walca dla m1=4 i m2=2','położenie suwaka dla m1=40 i m2=20','położenie walca dla m1=40 i m2=20');
hold off;
%*********************************************************
%zmiana masy obiektów
m1=2;
m2=4;
[t,x,y]=sim('tarcie',[0 50]);
figure(6)	
plot(t,y(:,1),'g');% wyświetla tylko pierwszą kolumnę macierzy y w kolorze zielonym - 'g'reen 
hold on;
plot(t,y(:,2),'b');
m1=20;
m2=40;
[t,x,y]=sim('tarcie',[0 50]);
hold on;%zachowując stary wykres umożliwia naniesienie na układ współrzędnych kolejnych wyników
plot(t,y(:,1),'r');
hold on;
plot(t,y(:,2),'m');
grid;
title('Porównanie zmian masy obiektów');
xlabel('czas [s]');
ylabel('położenie [m]');
legend('położenie suwaka dla m1=2 i m2=4','położenie walca dla m1=2 i m2=4','położenie suwaka dla m1=20 i m2=40','położenie walca dla m1=20 i m2=40');
hold off;
%*********************************************************