I try to use reliability analysis by using PCK model. 200 experimental data is availabe and there is no explicit function. I write the code as follows:
FILELOCATION = fullfile(…
uq_rootPath, ‘Examples’, ‘SimpleDataSets’, ‘Truss_Matlab_FEM’);
load(fullfile(FILELOCATION,‘data.mat’), ‘X’, ‘Y’);
InputOpts.Marginals(i).Name = sprintf(‘E%d’,i);
InputOpts.Marginals(i).Type = ‘Gaussian’;
InputOpts.Marginals(i).Moments = [0.25 1];
APCKMCSOpts.Type = ‘Reliability’;
APCKMCSOpts.Method = ‘AKMCS’;
AKOpts.AKMCS.MaxAddedED = 200;
AKOpts.AKMCS.IExpDesign.X = X;
AKOpts.AKMCS.IExpDesign.Y = Y;
APCKMCSOpts.AKMCS.MetaModel = ‘PCK’;
APCKMCSOpts.AKMCS.PCK.Kriging.Corr.Family = ‘Gaussian’;
APCKMCSOpts.AKMCS.IExpDesign.N = 5;
APCKMCSOpts.Simulation.MaxSampleSize = 1e6;
myAPCKMCSAnalysis = uq_createAnalysis(APCKMCSOpts);
I get this error:
No model defined!
Error in uq_evalLimitState (line 67)
M_X = uq_evalModel(FullModel,X);
Could you help me to find my mistake?
In fact, I decide to predict the response of generated random variables and then estimate the failure probability. I didn’t find any example which has such process. If there is any example in UQLab examples, please let me know. It could help me to get idea for improving aformentioned code.