var = symvar(f);
varc = symvar(varvec);
s1 = length(var);
s2 = length(varc);
m =floor((s1-1)/3+1);
varv = zeros(1,m);
if s1 ~= s2
for i=0: ((s1-1)/3)
k = symvar(varc,var(3*i+1));
index = (k-1)/3;
varv(i+1) = varval(index+1);
end
fv = subs(f,var,varv);
else
fv = subs(f,varvec,varval);
end
Inconsistency between sizes of second and third arguments.
出错 sym/subs>mupadsubs (line 157)
[X2,Y2,symX,symY] = normalize(X,Y); %#ok
出错 sym/subs (line 145)
G = mupadsubs(F,X,Y);
出错 Funval (line 17)
fv = subs(f,varvec,varval);
出错 minNT (line 18)
v=Funval(gradf,var,x0);
出错 minNF (line 27)
[x2,minf]=minNT(SumF,transpose(x1),var); %用牛顿法求解无约束规划
0.0001
minf=1.0004