//
//
// L03 lab01
//
clear; //czysc tablice zmiennych
clc; //czysc konsole
close(winsid()); //zamknij aktywne okna graficzne
//
// Deklaracje
//
wzo=9.81; //[s] sygnal kalibrujacy
A=1.2; //amplituda sygnalu
t0=0; //[s] czas poczatkowy
krp=0.01; // [s] krok pomiaru/symulacji
tk=20; //czas koncowy [s]
t=(t0:krp:tk); //[s] sciezka czasu
disp('t0=', t0); //
syggen=A*sin(t); //sygnal sinusoidalny o amplitudzie A
//
//
scf(0); //ustaw aktywne okno graficzne
clf(0); //czysc aktywne okno graficzne
plot2d([t],[syggen],1);
plot2d([t],[syggen*2/3],2);
xtitle('','t[s]','A[m/s^2]');
legends(['sin(t)', '2/3*sin(t)'],[1,2],opt='ur');
xs2pdf(0,'okno0.pdf'); //ss do pdf
//
//Wczytanie danych pomiarowych
//
kalibracja=read('kalibracja.txt',-1,2); //wczytaj dane
scf(1); //ustaw aktywne okno graficzne
clf(1); //czysc aktywne okno graficzne
plot2d([kalibracja(:,1)],[kalibracja(:,2)],3);
xtitle('','t[s]','A[m/s^2]');
legends(['kalibracja'],[3],opt='ur');
xs2pdf(1,'okno1.pdf'); //ss do pdf
wspkX=wzo/max(kalibracja(:,2));
 
scf(2); //ustaw aktywne okno graficzne
clf(2); //czysc aktywne okno graficzne
plot2d([kalibracja(:,1)],[kalibracja(:,2)*wspkX],5);
xtitle('','t[s]','A[m/s^2]');
legends(['kalibracja'],[5],opt='ur');
xs2pdf(2,'okno2.pdf'); //ss do pdf
disp(max(kalibracja(:,2)*wspkX));
 
//pomiar
 
//wczytanie pomiaru do macierzy pomiar
pomiar=read('pomiar.txt',-1,2);
//kalibracja danych pomiarowych
pomiar=pomiar*wspkX;
//prezentacja i analiza wynikow pomiarowych
scf(3); //ustaw aktywne okno graficzne
clf(3); //czysc aktywne okno graficzne
plot2d([pomiar(:,1)],[pomiar(:,2)*wspkX],6);
xtitle('','t[s]','A[m/s^2]');
legends(['pomiar'],[6],opt='ur');
xs2pdf(3,'okno3.pdf'); //ss do pdf

//pomiar2

//wczytanie pomiaru do macierzy pomiar
pomiar2=read('pomiar2.txt',-1,2);
//kalibracja danych pomiarowych
pomiar2=pomiar2*wspkX;
//prezentacja i analiza wynikow pomiarowych
scf(4); //ustaw aktywne okno graficzne
clf(4); //czysc aktywne okno graficzne
plot2d([pomiar2(:,1)],[pomiar2(:,2)*wspkX],7);
xtitle('','t[s]','A[m/s^2]');
legends(['pomiar2'],[7],opt='ur');
xs2pdf(4,'okno4.pdf'); //ss do pdf