mumudengyuan
2008-11-11, 21:05
function f = objfun(belta)
f=0;
X=[1.5000 3.5000 1.3000 1.7000
2.5000 4.1000 3.7000 2.6000];
K=[0.0437 0.0359 0.0718 0.0710];
A=[11.1030 17.6180
17.6180 28.4130];
N=length(K);
for i=1:N
temp=1+(K(i)-X(:,i)'*belta/(belta'*A*belta))^2;
f = f+1/temp;
end
f=-f;
[belta,fval,exitflag]=fminsearch(@objfun,[0.2;0.3])
高手帮忙看下。这个规划的函数值怎么也不可能是[0,1]之外的值，但是我运行后出现的是个负值。
困惑....:ft:
f=0;
X=[1.5000 3.5000 1.3000 1.7000
2.5000 4.1000 3.7000 2.6000];
K=[0.0437 0.0359 0.0718 0.0710];
A=[11.1030 17.6180
17.6180 28.4130];
N=length(K);
for i=1:N
temp=1+(K(i)-X(:,i)'*belta/(belta'*A*belta))^2;
f = f+1/temp;
end
f=-f;
[belta,fval,exitflag]=fminsearch(@objfun,[0.2;0.3])
高手帮忙看下。这个规划的函数值怎么也不可能是[0,1]之外的值，但是我运行后出现的是个负值。
困惑....:ft: