2楼: Originally posted by hzlhm at 2021-12-20 13:14:32
我们的眼睛视力不好，看不清图片的内容，请将代码直接复制到贴内。

clear;clc;

format shortG

options = optimset('Display','off');

H = [1.32706 -0.26629 -0.10485 2.37964
    -2.14291 1.11803 0.89604 0.92137
    -0.27878 -0.86637 -0.43404 0.18012
    -0.35438 0.09057 0.22529 -0.09379
    0.42988 -1.50226 2.11735 -0.66289
    0.34362 0.63772 -0.49262 -0.80495
    -0.00402 -0.09919 -1.29207 0.04081
    -0.60345 -0.65027 -0.68906 -0.44298
    -0.04899 -0.38018 -0.83804 -0.83195
    1.33196 1.91824 0.61199 -0.68537]

q = zeros(10,4);
w = zeros(10,1);

for j=1:10
   
    c=[-H(j,4),0,0,0];
   
    a=[H(1,4),-H(1,1:3);
        H(2,4),-H(2,1:3);
        H(3,4),-H(3,1:3);
        H(4,4),-H(4,1:3);
        H(5,4),-H(5,1:3);
        H(6,4),-H(6,1:3);
        H(7,4),-H(7,1:3);
        H(8,4),-H(8,1:3);
        H(9,4),-H(9,1:3);
        H(10,4),-H(10,1:3);
        H(j,4),0,0,0];
   
    b=[0;0;0;0;0;0;0;0;0;0;1];
    aeq=[0,H(j,1:3)];
    beq=[1];
    lb=[0,0,0,0];
    ub=[1,1,1,1];
   
    [xx,fval,exitflag,output] = linprog(c,a,b,aeq,beq,lb,ub,[],options);

    q(j,:) = xx';
    w(j,:) = exitflag';
   
end

q
w

5楼: Originally posted by 小灰灰SHG at 2021-12-21 20:34:04
您好，源代码已发布在评论，感谢您的回帖，麻烦您帮忙啦⊙▽⊙
...

6楼: Originally posted by zyh700600 at 2021-12-23 11:27:45
clear;clc;

format shortG

options = optimset('Display','off');

H =

q = zeros(10,4);
w = zeros(10,1);

for j=1:10
   
    c=;
   
    a=;
   
    b=;
    aeq=;
    be ...

7楼: Originally posted by hzlhm at 2021-12-24 12:33:28
根据你的代码运行，得到是空解。所以软件会报错。
有具体的方程和条件吗？...

9楼: Originally posted by 小灰灰SHG at 2021-12-25 16:50:40
程序逻辑上是通的，是因为数据问题没有解吗？
...