lab script

clc;
clear all;
close all;
%% Channels
Time = 1;
Accelerometer = 2;
LED = 3;
IR = 4;
Voltage = 5;
Current = 6;
Impedance = 7;
Freq = [1 5 10 15 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 50 100 200 300 400 500 600 700 800 900 1000 2000 3000 4000 5000];
Freq2= [1 5 10 15 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 50 100];
%% Set the Dir for Part B step 3
folder = 'Data/Part B Step 3';
filetype='*.xls';
f=fullfile(folder,filetype);
d=dir(f);

%% Set the Dir for Part B step 5
folder = 'Data/Part B Step 5';
filetype='*.xls';
f=fullfile(folder,filetype);
d=dir(f);
%% Part B laod files based on the input freq
for k=1:numel(d);
  data{k}=xlsread(fullfile(folder,d(k).name));
end

%% Offset the voltage data

for k=1:numel(d);
  data{k}(5:end)-0.749;
end

%% Getting the Impedance Matrix
for k=1:numel(d);
    for i = 1:200;
        data{k}(i,Impedance) = data{k}(i,Voltage)/data{k}(i,Current);
    end
end
%% Getting the average impedacne
for j = 1:k;
    Average_Impedance(j)= abs(mean(data{j}(Impedance:end)));
end
%% Plot the impedance over frequency for Part B step 3
loglog(Freq, Average_Impedance);
xlabel('Freqency in log scale');
ylabel('Impedance in lof scale');
title('Impedance & Frequency');

%% Plot the impedance over frequency for Part B step 5
loglog(Freq2, Average_Impedance);
xlabel('Freqency in log scale');
ylabel('Impedance in lof scale');
title('Impedance & Frequency with Extra Weights');