Question

应用matlab进行线性规划遇到的问题

clc,clear
a=0;
hold on
while a<0.05
c=[-0.05,-0.27,-0.19,-0.185,-0.185];
a=[zeros(4:1),diag([0.025,0.015,0.055,0.026])];
b=a*ones(4:1);
aeq=[1,1.01,1.02,1.045,1.065];
beq=1;
LB=zeros(5,1);
[x,Q]=linprog(c,a,b,aeq,beq,LB)；
Q=-Q;
plot(a,Q,'*k');
a=a+0.001;
end
xlabel('a'),ylabel('Q')

运行后显示红字为
[x,Q]=linprog(c,a,b,aeq,beq,LB)；
|
Error: Missing variable or function

Answer 1

存在问题

1、报错的那行，行尾的分号误写为全角符号（；），应改为半角分号（;）。

2、a既作为循环控制变量使用，同时又作为优化问题的系统矩阵使用，造成混乱。

3、有两处zeros的参数本应为(4,1)，被误写为(4:1)。

 

修改建议

a=0;
hold on
while a<0.05
    c=[-0.05,-0.27,-0.19,-0.185,-0.185];
    A=[zeros(4,1),diag([0.025,0.015,0.055,0.026])];
    b=a*ones(4,1);
    aeq=[1,1.01,1.02,1.045,1.065];
    beq=1;
    LB=zeros(5,1);
    [x,Q]=linprog(c,A,b,aeq,beq,LB);
    Q=-Q;
    plot(a,Q,'*k');
    a=a+0.001;
end
xlabel('a'),ylabel('Q')

 

运行结果

Answer 2

help linprog 看是否有该函数