求助：关于linmap模型的程序编写问题!（程序已写，可以运行，可是结果都是0，为什么）

jqylb

以下程序应该怎么编写？？？我研究了好多天还是一头雾水，求大神帮帮忙啊！！！

其中，


设决策人最偏好的方案用指标空间的理想点表示，则指标空间任一样品点 到理想点的加权欧几里德距离平方Si为


式中Wj(j=1，2，L，p)为第j个指标的权重平方。

Q为方案的有序对集

具体数据如下：

   3.2 4.0 3.2 （方案1）
   4.0 4.8 3.8 （方案2）
X= 2.8 3.8 4.0 （方案3）
   3.6 4.2 4.2 （方案4）

方案的优劣顺序O={方案2，方案4，方案1，方案3}

方案的有序对集Q={（2,4），（2,1），（2,3），（4,1），（4,3），（1,3）}

要求的是权重Wj

jqylb

求权重W

jqylb

图片看不到？？？

jqylb

我自己写了一下，程序可以运行了，可是计算出来的结果都是0，这是为什么啊？？？？


model:
sets:
fa/s1 s2 s3 s4/;
zb/zb1 zb2 zb3/:w,v;
bzh(fa,zb):z;
q(fa,fa)/s2,s4 s2,s1 s2,s3 s4,s1 s4,s3 s1,s3/:r;
endsets
data:
z=-0.39,-0.46,-1.39
1.16,1.39,0.00
-1.16,-0.93,0.46
0.39,0.00,0.93;
enddata
min=@sum(q(k,l):r(k,l));
@for(q(k,l):@sum(zb(j):w(j)*(z(l,j)^2-z(k,j)^2))-2*@sum(zb(j):v(j)*(z(l,j)-z(k,j)))+r(k,l)>0);
@for(q(k,l):@sum(zb(j):w(j))*@sum(zb(j):z(l,j)^2-z(k,j)^2)-2*@sum(zb(j):v(j))*@sum(zb(j):z(l,j)-z(k,j))=h);
@for(q(k,l):r(k,l)>0);
@for(zb:w>0);
h>0;
@for(zb:@free(v););
end




Global optimal solution found.
Objective value: 0.000000
Infeasibilities: 0.000000
Total solver iterations: 0

Variable Value Reduced Cost
H 0.000000 0.000000
W( ZB1) 0.000000 0.000000
W( ZB2) 0.000000 0.000000
W( ZB3) 0.000000 0.000000
V( ZB1) 0.000000 0.000000
V( ZB2) 0.000000 0.000000
V( ZB3) 0.000000 0.000000
Z( S1, ZB1) -0.3900000 0.000000
Z( S1, ZB2) -0.4600000 0.000000
Z( S1, ZB3) -1.390000 0.000000
Z( S2, ZB1) 1.160000 0.000000
Z( S2, ZB2) 1.390000 0.000000
Z( S2, ZB3) 0.000000 0.000000
Z( S3, ZB1) -1.160000 0.000000
Z( S3, ZB2) -0.9300000 0.000000
Z( S3, ZB3) 0.4600000 0.000000
Z( S4, ZB1) 0.3900000 0.000000
Z( S4, ZB2) 0.000000 0.000000
Z( S4, ZB3) 0.9300000 0.000000
R( S2, S4) 0.000000 1.000000
R( S2, S1) 0.000000 1.000000
R( S2, S3) 0.000000 1.000000
R( S4, S1) 0.000000 1.000000
R( S4, S3) 0.000000 1.000000
R( S1, S3) 0.000000 1.000000