fun.m

来自「本源码用于解决金融领域的指数复制的问题」· M 代码 · 共 14 行

function f=fun(y,x,Q,R,I)
n=size(y,1);
I=rep(I(2:end),[1,n]);
I=I+1;
z=y';
A=Q*z;
B=R*z;
C=A./B-I;
C=C.^2;
f=sum(C);
f=f';
f=sqrt(f);
f=f.*10000;