% Example: Create a space-filling design factors = 'Speed', 'Load', 'Timing'; range = [800 6000; % RPM 20 120; % Load (%) -10 30]; % Timing (deg) des = xydesign(factors, range, 'NumPoints', 50); scatter(des); % Visualize Load measured data from engine tests.
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.
% 3. Build knock model (binary: 0=no knock, 1=knock) knock_model = mbcgp(data, 'Knock', 'Speed','Load','Timing', 'Distribution','binomial'); knock_model = fit(knock_model);