function [OBJ] = Main_Combined_optimization(Opt_Param)
%% code Description: Main optimization Function
%%Run Model with optimization Parameters
global hot;
Data = xlsread('hotspot.csv');
hot.RecordingTime = Data(:,1);
hot.Loading = Data(:,2);
hot.HotToTopOilGradient = Data(:,3);
t = hot.RecordingTime;
K = hot.Loading;
Hot_Output = hot.HotToTopOilGradient;
%% call the optimizer
[HotSpotTemperture_Optimization] = Combined_Optimization(K,Opt_Param,t);
error_HT = Hot_Output – HotSpotTemperture_Optimization;
OBJ_HT = sum( trapz(t, error_HT.^2));
OBJ = OBJ_HT;
%% Calculate the RMSE
RMSE = sqrt( sum((Hot_Output(:)- HotSpotTemperture_Optimization(:)).^2) / numel(Hot_Output) );
end
Best Answer