%% Discrete modelling
Ae=(1-(K*ts)/Cth);Be=(ts*COP*HP)/Cth;D=ts*K/Cth;Tin=zeros(1,i);u=zeros(1,i);for j=1:i % f=@(u) sum(u(j)*gamma(j));
f=u(j); %%%Starting point%%%
u0=0; if j>1 Tin(j+1)=Ae*Tin(j)+Be*u(j)+D*Tout(j); else Tin(j)=0; end %%%Constraints%%%
%%%Equality%%%%
Aeq=Be; beq=(Tin(j+1)-Ae*Tin(j)-D*Tout(j)); %%%Inequality%%%
A=[(Tin(j)*Be)/(Tin(j+1)-Ae*Tin(j)-D*Tout(j)) ;(-(Tin(j)*Be)/(Tin(j+1)-Ae*Tin(j)-D*Tout(j)))]; b=[Tmax ; -Tmin]; options = optimoptions(@linprog,'Display', 'off'); u(j)=linprog(f,u0,A,b,Aeq,beq,umin,umax,options);end
I am getting the error: The number of rows in A must be the same as the number of elements of b.
Does anyone know why?
Thanks!
Best Answer