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

[SPClean] = fft(PClean_f);

% Output

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

[SPdisturbed] = fft(Pdisturbed_f);

reading=20

for q = 1:readings

filename = sprintf('Tfunct%d_%d.mat',sn,q);
TF = importdata(filename);
TFtrans=transpose(TF(1:end,1));
T(q,1:end)=TFtrans;
end

% TF=T';

T_f_trans=mean(T); %Taking mean of T i.e. all TF from 1 TF1 to TF5
TFinal=transpose(T_f_trans);
A=TF.*SPdisturbed;

 B=ifft(A);