Commit c53c9324 authored by Aytekin Gel's avatar Aytekin Gel
Browse files

Script file to run PSUADE in batch mode to construct surrogate model 

This script is provided as input to PSUADE to first construct a quadratic polynomial based response surface and then perform 120 fold cross-validation error assessment. Finally generate several response surface plots.
parent bc84bcaf

Case1_ParticleSettling/DeterministicCalibration/PSUADE/1_ResponseSurface/Quadratic/PSbuildRS_2.in

0 → 100644
+33 −0
Original line number Diff line number Diff line 
load ../../psData_OLH_n120_i6_o1_y2 # load data file

ana_expert      # psuade> analysis expert mode on

printlevel 3    # psuade> printlevel set from 0 to 3

rscheck         # Check for best data-fitted surrogate model

2               # Select Linear regression as response surface method

1               # Quantity of interest no

y               # Perform cross-cross_validation

120             # Enter the number of groups to validate : (2 - 120)

n               # Random selection of leave-out groups ? (y or n)

sys mv RSFA_CV_err.m R2_QUAD_120F_RSFA_CV_err.m # Rename the CV error file before next one

sys mv RSFA_training_err.m R2_QUAD_RSFA_training_err.m # Rename the training error file before next one

rs2

256             # Grid resolution ? (32 - 256)

2               # Available response surface tools:  select Quadratic 

6               # Enter the input for x axis (1 - 7) :

5               # Enter the input for y axis (1 - 7), not 1 : 2

y               # Set other nominal values automatically ? (y or n)

1               # Enter output number (1 - 1) :

n               #  Ymin and Ymax found = 2.037723e-01 3.421112e-01. Set lower threshold ? (y or n) :

n               # Set upper threshold ? (y or n) : n

sys mv matlabrs2.m R2_QUAD_rs2_x6x5.m    # Rename the response surface constructed

rs2

256             # Grid resolution ? (32 - 256)

2               # Available response surface tools:  select LINEAR

6               # Enter the input for x axis (1 - 7) :

3               # Enter the input for y axis (1 - 7), not 1 : 2

y               # Set other nominal values automatically ? (y or n)

1               # Enter output number (1 - 3) :

n               #  Ymin and Ymax found = 2.037723e-01 3.421112e-01. Set lower threshold ? (y or n) :

n               # Set upper threshold ? (y or n) : n

sys mv matlabrs2.m R2_QUAD_rs2_x6x3.m    # Rename the response surface constructed

sys ls -alrt R2_QUAD*  # Check most one created after mv

q