clear all;  % Clear workspace
close all;  % Close plot windows
addpath('../utility');% Add path to utilities

dataBRSE=importdata('../../SuspensionRequirement/BRSE/DisplacementNoiseRequirementCopy.dat');
dataBRSE=[dataBRSE(:,1),dataBRSE(:,2),dataBRSE(:,3),dataBRSE(:,4),dataBRSE(:,5),dataBRSE(:,6),...
	      dataBRSE(:,7),dataBRSE(:,8),dataBRSE(:,9)];
dataDRSE=importdata('../../SuspensionRequirement/DRSE/DisplacementNoiseRequirementCopy.dat');
dataDRSE=[dataDRSE(:,1),dataDRSE(:,2),dataDRSE(:,3),dataDRSE(:,4),dataDRSE(:,5),dataDRSE(:,6),...
	      dataDRSE(:,7),dataDRSE(:,8),dataDRSE(:,9)];
dataTF_LBF2LTM=importdata('sumconplot/160620_typeBp_forceTF_LBF2LTM.dat');
dataTF_LBF2LTM=[dataTF_LBF2LTM(:,1),dataTF_LBF2LTM(:,2),dataTF_LBF2LTM(:,3)];
dataTF_VBF2VTM=importdata('sumconplot/160620_typeBp_forceTF_VBF2VTM.dat');
dataTF_VBF2VTM=[dataTF_VBF2VTM(:,1),dataTF_VBF2VTM(:,2),dataTF_VBF2VTM(:,3)];

dataTF_LBF2LTM_int=interp1(dataTF_LBF2LTM(:,1),dataTF_LBF2LTM(:,2),dataBRSE(:,1));
dataTF_VBF2VTM_int=interp1(dataTF_VBF2VTM(:,1),dataTF_VBF2VTM(:,2),dataBRSE(:,1));


% ALL
fig=figure;
loglog(dataBRSE(:,1),dataBRSE(:,8),'-r',...
	   dataDRSE(:,1),dataDRSE(:,8),'-b',...
	   dataBRSE(:,1),dataTF_LBF2LTM_int(:,1),'-k',...
	   dataBRSE(:,1),dataTF_VBF2VTM_int(:,1)/100.0,'--k',...
	   'LineWidth',1)
grid on
title('BF actuator requirement (memo)','FontSize',12,'FontWeight','bold','FontName','Times New Roman',...
    'interpreter','none')
ylabel('Magnitude [m/N] or [m/rtHz]','FontSize',12,'FontWeight','bold','FontName','Times New Roman')
xlabel('Frequency [Hz]','FontSize',12,'FontWeight','bold','FontName','Times New Roman')
%ylim([1e-20 1e-5])
xlim([1e-1 1e3])
set(gca,'FontSize',12,'FontName','Times New Roman')
legend('PRM(BRSE)','PRM(DRSE)','force TF LBF2LTM','force TF VBF2VTM')
set(fig,'Color','white')
%export_fig('figure/BFactuatorRequirement.pdf')

% BRSE
fig=figure;
loglog(dataBRSE(:,1),dataBRSE(:,8)./dataTF_LBF2LTM_int(:,1),'-r',...
	   dataDRSE(:,1),dataDRSE(:,8)./dataTF_LBF2LTM_int(:,1),'-b',...
	   dataBRSE(:,1),dataBRSE(:,8)./dataTF_VBF2VTM_int(:,1)*100.0,'--m',...
	   dataDRSE(:,1),dataDRSE(:,8)./dataTF_VBF2VTM_int(:,1)*100.0,'--c',...
	   'LineWidth',1)
grid on
title('BF actuator requirement','FontSize',12,'FontWeight','bold','FontName','Times New Roman',...
    'interpreter','none')
ylabel('Magnitude [N/rtHz]','FontSize',12,'FontWeight','bold','FontName','Times New Roman')
xlabel('Frequency [Hz]','FontSize',12,'FontWeight','bold','FontName','Times New Roman')
%ylim([1e-20 1e-5])
xlim([5e-1 1e3])
set(gca,'FontSize',12,'FontName','Times New Roman')
legend('PRM(BRSE) in L DoF','PRM(DRSE) in L DoF','PRM(BRSE) in V DoF','PRM(DRSE) in V DoF','Location','northwest')
set(fig,'Color','white')
%export_fig('figure/BFactuatorRequirementBRSE.pdf')