Experimental Design in Matlab

Can anyone help me with run a complete CCD or BBD experimental design in matlab? Actually I use RStudio for run my DoE but i would like to learn how to perform these in Matlab.

I found this page https://it.mathworks.com/help/stats/design-of-experiments-1.html that give me some information but this is not enough...

I would like to know also the script for give the experimental responses, extimate coefficient, perform the ANOVA and response surface and so on..

Thank you so much!

Answer 1

[Sam Chak]

Hi @Valeria Gallo

You can find some examples in the following links:

Response Surface Designs

Central composite design (CCD)

Box-Behnken design (BBD)

Example 1: A 2-factor central composite design.

dCC = ccdesign(2, 'type', 'circumscribed')

dCC = 16×2

-1.0000 -1.0000 -1.0000 1.0000 1.0000 -1.0000 1.0000 1.0000 -1.4142 0 1.4142 0 0 -1.4142 0 1.4142 0 0 0 0

plot(dCC(:, 1), dCC(:, 2), 'ro', 'MarkerFaceColor', 'b')
X = [1 -1 -1 -1; 1 1 1 -1];
Y = [-1 -1 1 -1; 1 -1 1 1];
line(X, Y, 'Color', 'b')
axis square equal

Example 2: A 3-factor Box-Behnken design.

dBB = bbdesign(3)

dBB = 15×3

-1 -1 0 -1 1 0 1 -1 0 1 1 0 -1 0 -1 -1 0 1 1 0 -1 1 0 1 0 -1 -1 0 -1 1

plot3(dBB(:, 1), dBB(:, 2), dBB(:, 3), 'ro', 'MarkerFaceColor', 'b')
X = [1 -1 -1 -1 1 -1 -1 -1 1 1 -1 -1; ...
     1 1 1 -1 1 1 1 -1 1 1 -1 -1];
Y = [-1 -1 1 -1 -1 -1 1 -1 1 -1 1 -1; ...
     1 -1 1 1 1 -1 1 1 1 -1 1 -1];
Z = [1 1 1 1 -1 -1 -1 -1 -1 -1 -1 -1; ...
     1 1 1 1 -1 -1 -1 -1 1 1 1 1];
line(X, Y, Z, 'Color', 'b')
axis square equal

Answer 2

[Valeria Gallo]

Hi Sam Chak, thank you for your answer. I have already seen the pages, but I need to find the way to run the subsequent steps: coding and uncoding the variables, find the significant factors, find the optimal conditions, find the predicted value of the response under optimal conditions and so on

Answer 3

[Moein Siadaty]

Dear Valeria Gallo

Do you find the answer?

actually I have the same problem and want to dimensionally optimize a heat exchanger as a final part of my matlab code using DOE.