%
% Evaluation of natural frequencies of pressure sensors for calibration
%
% 01.8.2014

clc;
close all hidden;
clear ;

addpath('Y:AmitAmit_readingstest')
%% Readings Laden

Pfad = 'Y:AmitAmit_readingstest';
% Pfad_1 = reference sensor 1cm to the speaker - measuring sensor 10cm in Box + 100cm for speakers

PfadSave='Y:AmitAmit_readingstest';

sn = 1;  % Sensornumber already checked
readings =20;

for ii=1:readings

    assignin('base', 'ActualReading', dlmread([Pfad, 'Sensor1_', num2str(ii),'.lvm'], 't'))

    t = ActualReading(1:end,1);           %Time where 1st column corresponds to time

    freq = 1/(t(2)-t(1));                    %Frequeny

    % Inputvalue

    assignin('base', ['PClean_',num2str(ii)],ActualReading(1:end,17));  % Reference sensor: Always number ////sensor no:16///// !!!!

    PClean_mean=eval(['mean','([PClean_',num2str(ii),'])']);

    assignin('base', ['PClean_',num2str(ii)],eval(['PClean_',num2str(ii)])-PClean_mean);

    % Outputvalue

    assignin('base', ['Pdisturbed_',num2str(ii)],ActualReading(1:end,2)); % Measuring sensor number ///// 1 + sensor number///// !!!!as 1st column is time

    Pdisturbed_mean=eval(['mean','([Pdisturbed_',num2str(ii),'])']);

    assignin('base', ['Pdisturbed_',num2str(ii)],eval(['Pdisturbed_',num2str(ii)])-Pdisturbed_mean);

    % Reference
    assignin('base', ['Imp_',num2str(ii)],['Pdisturbed_',num2str(ii)]);

    % Representation of the pressure measurements
    assignin('base', 'PClean',eval(['PClean_',num2str(ii)]));
    assignin('base', 'Pdisturbed',eval(['Pdisturbed_',num2str(ii)]));

    % Fast Fourier Analyse
    NFenster=2^12;

    %  Input
    assignin('base', 'PClean_f',eval(['PClean_',num2str(ii)]));
    [SPClean]      = fft(PClean_f,NFenster);

    %  Output
    assignin('base', 'Pdisturbed_f',eval(['Pdisturbed_',num2str(ii)]));
    [SPdisturbed]      = fft(Pdisturbed_f,NFenster);

    TFunct= SPClean./SPdisturbed;

    SP=TFunct.*SPdisturbed;

    % SPi=ifft(SP);
    % Z=norm(abs(PClean)-abs(SPi));
    filename = sprintf('Tfunct%d_%d.mat',sn,ii);
    save(filename,'TFunct');

end

close all;