1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
function varargout = calModelPara(varargin)
y = varargin(1)+varargin(2)
N = size(y,'r')
ek = varargin(2)
np = varargin(3)
//pause
mse = sum(ek.^2)/N;//disp(mse)
fpe = mse * (1 + np/N)/(1 - np/N);//disp(fpe)
nrmse = 1 - sqrt(sum(ek^2))/sqrt(sum((y - mean(y))^2));//disp(nrmse)
AIC = N * log(mse) + 2 * np + N * size(y,'c') * (log(2 * %pi) + 1);//disp(AIC)
AICc = AIC * 2 * np * (np + 1)/(N - np - 1);//disp(AICc)
nAIC = log(mse) + 2 * np/N;//disp(nAIC)
BIC = N * log(mse) + N * size(y,'c') * (log(2 * %pi) + 1) + np * log(N);//disp(BIC)
//pause
varargout(1) = struct('MSE',mse,'FPE',fpe,'FitPer',nrmse*100,'AIC',AIC,'AICc',AICc,'nAIC',nAIC,'BIC',BIC)
endfunction
|
