function zd
global Mu Sigama %定義全域變數
Mu=[10 2.5]; %輸入平均值
Sigama=[2 0.375]; %輸入標準差
X0=[10; 2.5]; %輸入初始迭代點
A=[];b=[];Aeq=[];beq=[];
lb=[];ub=[];
options=optimset('Algorithm','active-set');
[X,fval,exitflag,output]=fmincon(@bata2,X0,A,b,Aeq,beq,lb,ub,@st,options) %呼叫優化工具箱求解
bata=sqrt(fval)
Pf=cdf('norm',-bata,0,1) %計算失效概率
end
function CC=bata2(X) %目標函式子函式
global Mu Sigama %定義全域變數
CC=((X(1)-Mu(1))/Sigama(1))^2+((X(2)-Mu(2))/Sigama(2))^2; %X均為正態分布
end
function [c,ceq]=st(X) %約束條件子函式
c=[]; %非線性不等式約束
ceq=stBplm(X); %非線性等式約束
end
function Z = stBplm(X)
W1=[0.1242 -1.9466; 67.3014 22.5133; -90.9424 -20.8143];
B1=[-1.5036; 4.0199; 3.0905];
W2=[-552.2048 28.3283 -130.4954];
B2=[-10.8555];
XX=[X(1);X(2)];
Z=simuff(XX,w1,b1,'logsig',w2,b2,'purelin');
end
function[w1,b1,w2,b2]=trainBplm(X);
X1=[10 16 4 10 10];
X2=[2.5 2.5 2.5 3.625 1.375];
p=[X1;X2];
t=[13.6728;10.8005;16.5451;16.8897;-10.3136];
s1=3;
tp=[1 1000 0.0000001];
[w1,b1,w2,b2] = initff(p,s1,'logsig',t,'purelin');
[w1,b1,w2,b2] = trainlm(w1,b1,'logsig',w2, b2,'purelin',p,t,tp);
end
》》出現錯誤如下:
ceq=stBplm(X); %非線性等式約束
出錯 fmincon (line 622)
[ctmp,ceqtmp] = feval(confcn{3},X,varargin{:});
出錯 rgsj_xym_lwltff (line 9)
[X,fval,exitflag,output]=fmincon(@bata2,X0,A,b,Aeq,beq,lb,ub,@st,options) %呼叫優化工具箱求解
原因:
Failure in initial nonlinear constraint function evaluation. FMINCON cannot continue.
轉載請註明出處,本文鏈接:https://www.uj5u.com/qita/89319.html
標籤:其他技術專區
上一篇:410c的串口通訊
下一篇:填坑Ⅱ | 簡單的資料結構
