example5_7.m

来自「MATLAB7辅助控制系统设计与仿真 第5章」· M 代码 · 共 14 行

function [eval] = corana(sol)
numv = size(sol,2);
x=sol(1:numv);
d0=[1 1000 10 100 1 10 100 1000 1 10];
d=d0(1:numv);
c=0.15;
s=.2*ones(1,numv);
t=0.05*ones(1,numv);
bk = s.*(round(x./s));
dev= (abs(bk-x)<t) & (bk~=0);
z=c*((bk+sign(bk).*t).^2).*d;
y=x.^2.*d;
eval = sum((dev.*z) + ((~dev).*y));