data = mbcdata.import('engine_test.csv'); % Remove outliers data = removeoutliers(data, 'Response', 'BSFC'); % Split into training/validation [train, val] = splitdata(data, 0.8); Use mbcmodels to create response surface models.
quad = mbcquadratic(train, 'Response', 'NOx', 'Predictors', 'Speed','Load','Timing'); quad = fit(quad); summary(quad); % See coefficients Use the calset object to optimize lookup tables for multiple responses (e.g., minimize BSFC, keep NOx < limit). mcc toolbox
% 1. Load data load engine_data.mat % contains Speed, Load, Timing, Torque, Knock % 2. Build torque model torque_model = mbcgp(data, 'Torque', 'Speed','Load','Timing'); torque_model = fit(torque_model); data = mbcdata
(best for non-linear):
