关于fmincon函数，求大神帮忙看看

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
出现问题如下：
>> rgsj_xym_lwltff
未定义函数或变量 'w1'。出错 rgsj_xym_lwltff>stBplm (line 23)
Z=simuff(XX,w1,b1,'logsig',w2,b2,'purelin');出错 rgsj_xym_lwltff>st (line 19)
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.