想用Fmincon函数求最小值,但是不会弄非线性约束函数,
x1=[4,0,0,0,0];
x1_step1.xoffset = [8;0.41;0.58;4;1.63];
x1_step1.gain = [0.5;10;10;0.5;2];
x1_step1.ymin = -1;
% Layer 1
b1 = [-0.043565515598342514;-0.56099578077009082;0.59013546057547384;-0.017459943136561403;-1.0124263684328507;-0.20878731073027418;0.031244278125342556;-0.3101619375736791];
IW1_1 = [-0.16719132346005225 -0.026797544292805463 -0.034384828437724134 0.022541034715052075 0.018886125813680096;-0.063362066358251978 -0.25930242971889805 0.29688282151075523 0.35410679513986487 0.14418009876171486;9.4989521505021819e-05 0.0099221821076234881 -0.77543394074759731 -0.038686727118812617 0.017095274010232038;0.12181365143990305 -0.35570716555781268 -0.19905444993113527 -0.24824990651980408 0.074840080298249817;-0.0074142602772374368 -0.046797482159538099 0.61170233925146211 -0.54744875708124652 0.17371231944270074;-0.044509415064994493 0.29184314667334654 -0.12555067739168302 -0.051228805557555615 -0.47004367575757849;-0.085451487877427776 0.088274551700570122 -0.10838745085041833 0.16233000995335692 0.44859889392386987;0.39508031931455723 -0.17487120854028659 0.10912099796712307 0.1770023953000347 0.098048086256165776];
% Layer 2
b2 = -0.096297419579917456;
LW2_1 = [0.15733080823956919 0.5582584860812797 -0.71403216826159432 -0.36506387740905455 -0.79069223936174493 -0.36151802575898395 0.34596579661852539 0.30739705348959051];
% Output 1
y1_step1.ymin = -1;
y1_step1.gain = 0.00140254904551141;
y1_step1.xoffset = 2031.0241;
% ===== SIMULATION ========
% Dimensions
Q = size(x1,2); % samples
% Input 1
xp1 = mapminmax_apply(x1,x1_step1);
% Layer 1
a1 = tansig_apply(repmat(b1,1,Q) + IW1_1*xp1);
% Layer 2
a2 = repmat(b2,1,Q) + LW2_1*a1;
% Output 1
y1 = mapminmax_reverse(a2,y1_step1);
ZN=[0.99];
CX(x1)=ZN - y1;
end
% ===== MODULE FUNCTIONS ========
% Map Minimum and Maximum Input Processing Function
function y = mapminmax_apply(x,settings)
y = bsxfun(@minus,x,settings.xoffset);
y = bsxfun(@times,y,settings.gain);
y = bsxfun(@plus,y,settings.ymin);
end
% Sigmoid Symmetric Transfer Function
function a = tansig_apply(n,~)
a = 2 ./ (1 + exp(-2*n)) - 1;
end
% Map Minimum and Maximum Output Reverse-Processing Function
function x = mapminmax_reverse(y,settings)
x = bsxfun(@minus,y,settings.ymin);
x = bsxfun(@rdivide,x,settings.gain);
x = bsxfun(@plus,x,settings.xoffset);
end
想用0.99-原函数,但是一直行不通