EMG - script

Script EMG
clear; clear all; clc; close all;

%% Opdracht 1 - Kalibratie van de krachtopnemer
% Data laden
%1.7806 kg
load('MyData-01-12Mar2019at14h07m30s.mat');
meting06 = DataMinute001;
kalibratie1 = double (meting06);
%4.8956 kg
load('MyData-02-12Mar2019at14h10m48s.mat');
meting02 = DataMinute001;
kalibratie2 = double (meting02);
%9.8366 kg
load('MyData-03-12Mar2019at14h12m33s.mat');
meting03 = DataMinute001;
kalibratie3 = double (meting03);
%14.8636 kg
load('MyData-04-12Mar2019at14h15m32s.mat');
meting04 = DataMinute001;
kalibratie4 = double (meting04);
%19.6916kg
load('MyData-05-12Mar2019at14h17m39s.mat');
meting05 = DataMinute001;
kalibratie5 = double (meting05);

% Kalibratielijn bepalen

m = [1.7806, 4.8956, 9.8366, 14.8636, 19.6916]; g = 9.81;
F_opgelegd = m * g;

gem1 =
mean(kalibratie1(:, 1)); %om 1 getal te krijgen
gem2 =
mean(kalibratie2(:, 1));
gem3 =
mean(kalibratie3(:, 1));
gem4 =
mean(kalibratie4(:, 1));
gem5 =
mean(kalibratie5(:, 1));
V = [gem1, gem2, gem3, gem4, gem5];

coef = polyfit(F_opgelegd, V, 1);
a = coef(1); b = coef(2);

y = a * F_opgelegd + b;

% Plotten kalibratielijn

figure(1)
plot( F_opgelegd, V, 'r', F_opgelegd, y, 'o')
title ('Data en de kalibratielijn')
xlabel ('Kracht [N]')
ylabel ('Spanning [microV]')
legend ('Kalibratielijn', 'Data punten')

% Er is sprake van een offset, omdat b niet nul is.

% Omschrijven
F_opgelegd = (y - b) / a;