各位大神，求教。目标函数中包含6个变量：x,y,z,w,g,h。于此同时，这6个变量都是关于时间变化的，即T=1：24. 相当于，确定6*24组变量的值，来使最终的目标函数f最小。

难点在于：目标函数f(t),和约束条件中，不仅包含x(t),还包含x(t-1)。

具体描述如下：

function f=objfun(x,y,z,w,g,h)

for t=1:24

if (x(t)~=0&&x(t-1)~=0)||(x(t)==0&&x(t-1)==0)

output=0

elseif (x(t)~=0&&x(t-1)==0)

a=2

else

b=3

end

if (y(t)~=0&&y(t-1)~=0)||(y(t)==0&&y(t-1)==0)

output=0

elseif (y(t)~=0&&y(t-1)==0)

c=2

else

d=3

end

if (z(t)~=0&&x(t-1)~=0)||(z(t)==0&&x(t-1)==0)

output=0

elseif (z(t)~=0&&x(t-1)==0)

e=2

else

f=3

end

f=2.03*(1/9.7)*x(t)/((0.0753*(x(t)/65)^3)-0.3095*(x(t)/65)^2+0.4174*(x(t)/65)+0.1068)...%微型燃气轮机燃料成本

+(2.03*(1/9.7)*y(t)/(-0.0023*y(t)+0.6735))...

+(z(t)^2+2*z(t)+3)...

+0.165179*x(t)...

+0.103616*y(t)...

+2.5*w(t)...

+0.68*g(t)...%

+5*(a+c+e)...

+4*(b+d+f) %调用上边的函数。

end

end

function [c,ceq]=confun(x,y,z,w,g,h)

c=[x(t)-x(t-1)-3;x(t-1)-x(t)-3;y(t)-y(t-1)-3;y(t-1)-y(t)-3;z(t)-z(t-1)-3;z(t-1)-z(t)-3]

ceq=[x(t)+y(t)+z(t)+w(t)+g(t)+h(t)-19]

end

clear

x0=[0,0,0,0,0,0,0,0,0,5,30,30,30,30,30,30,0,5,15,15,17,15,5,4];%24

y0=[0,0,0,0,0,0,0,30,30,30,30,30,30,30,30,30,30,30,30,30,30,30,30,30]

z0=[0,0,0,0,0,0,0,30,30,30,30,30,30,30,30,30,30,30,30,30,30,30,30,30]

w0=[27,23,28,30,23,12,20,0,11,5,-25,-27,-21,-30,-27,-25,15,30,30,30,30,30,30]

g0=[0,0,0,0,0,0,0,0,0,0,0,9,12,13,15,9,8,4,4,3,3,0,0,0]

h0=[0,0,0,0,0,0,0,0,0,0,0,-10,-5,-8,-5,0,10,10,0,-8,-10,-2,0,0]

%上下限约束的值不一定是这些，暂时这样

xM=[12,13,13,5,12,13,13,5,15,12,13,13,5,15,12,13,13,5,12,13,13,5,5,5]

xm=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]

yM=[12,13,13,5,12,13,13,5,15,12,13,13,5,15,12,13,13,5,12,13,13,5,5,5]

ym=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]

zM=[12,13,13,5,12,13,13,5,15,12,13,13,5,15,12,13,13,5,12,13,13,5,5,5]

zm=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]

wM=[12,13,13,5,12,13,13,5,15,12,13,13,5,15,12,13,13,5,12,13,13,5,5,5]

wm=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]

gM=[12,13,13,5,12,13,13,5,15,12,13,13,5,15,12,13,13,5,12,13,13,5,5,5]

gm=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]

hM=[12,13,13,5,12,13,13,5,15,12,13,13,5,15,12,13,13,5,12,13,13,5,5,5]

hm=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]

options=optimset('largescale','off','display','iter');

[x,y,z,w,g,h,fval,exitflag,output]=fmincon(@objfun,x0,y0,z0,w0,g0,h0,[],[],[],[],xM,xm,yM,ym,zM,zm,wM,wm,gM,gm,hM,hm @confun,options)

完全运行不出来，错误：

??? 'confun' is not a valid base class.

求各位大神解答，应该错误挺多的。多谢，多谢！