24小时热门版块排行榜

返回列表

alex110

铁杆木虫 (正式写手)

应助: 0 (幼儿园)
金币: 6533
散金: 705
帖子: 307
在线: 868.7小时
虫号: 1160169
注册: 2010-12-01
专业: 岩土力学与岩土工程

[求助] lingo计算求解，一直不出现结果，请高手帮忙看看如何解决得到结果，程序如下。

MODEL:
TITLE DEMOPWL_V4;

SETS:
FEELING/ 1..166 /: X;                                                                                        ! Different types of feeling;
PWLPOINTS/ 1..7 /;                                                                                                                               ! PWL points;
PWLPARS/ ROW, COL /;                                                                                                    ! PWL parameters;
PWLFUNS( PWLPARS, PWLPOINTS ): veryUnCom, unCom, normCom, com, veryCom;                      ! Types of PWL functions;
WEIGHTS( FEELING, PWLPOINTS ): W;                                                                                              ! Weights of PWL functions;
COST( FEELING ): SA, SB, SC, SD, SE, SAC, SBC, SCC, SDC, SEC;
ENDSETS

DATA:
X = 0.0119, 0.0114, 0.0134, 0.0147, 0.0147,                                                                         ! very Uncom;
      0.0078, 0.0087, 0.0071, 0.0088, 0.0096, 0.0088, 0.0083, 0.0089,
      0.0095, 0.0076, 0.0094, 0.0081, 0.0073, 0.0067, 0.0073, 0.0061,
      0.0088, 0.0053, 0.0083, 0.0065, 0.0088, 0.0068, 0.0059, 0.0083, 0.0166, 0.0074, 0.0081, 0.0136                   ! uncom;
      0.0045, 0.0043, 0.0061, 0.0065, 0.0048, 0.0041, 0.0067, 0.0051,
      0.0039, 0.0046, 0.0052, 0.0048, 0.0052, 0.0053, 0.0038, 0.0051,
      0.0065, 0.0061, 0.0047, 0.0057, 0.0056, 0.0048, 0.0051, 0.0067,
      0.0046, 0.0052, 0.0067, 0.0058, 0.0037, 0.0065, 0.0052, 0.0076,
      0.0078, 0.0054, 0.0044, 0.0081, 0.0073, 0.0067, 0.0048, 0.0053,
      0.0078, 0.0073, 0.0061, 0.0087, 0.0047, 0.0051, 0.0068, 0.0039,
      0.0047, 0.0051, 0.0066, 0.0074, 0.0058, 0.0049, 0.0038, 0.0041,
      0.0052, 0.0065, 0.0088, 0.0058, 0.0049, 0.0097, 0.0065, 0.0078,
      0.0047, 0.0074, 0.0073, 0.0098, 0.0074                               ! normal com;
      0.0021, 0.0009, 0.0004, 0.0049, 0.0037, 0.0027, 0.0037, 0.0018,
      0.0025, 0.0033, 0.0027, 0.0026, 0.0029, 0.0034, 0.0038, 0.0029,
      0.0039, 0.0033, 0.0019, 0.0018, 0.0029, 0.0014, 0.0034, 0.0027,
      0.0011, 0.0021, 0.0016, 0.0019, 0.0018, 0.0014, 0.0019, 0.0018,
      0.0039, 0.0061, 0.0027, 0.0036, 0.0019, 0.0046, 0.0048, 0.0051,
      0.0061, 0.0069, 0.0061, 0.0058, 0.0037, 0.0074, 0.0034, 0.0076,
      0.0054, 0.0109, 0.0055                                                                   ! com;
      0.0016, 0.0007, 0.0008, 0.0006, 0.0014, 0.0018, 0.0008, 0.0015,
      0.0008, 0.0007, 0.0007, 0.0089, 0.0032;                                                                                           ! very com;
ENDDATA

DATA:
veryUnCom = 0.0001,  , 0.0026, 0.0056, 0.0079,       , 0.02,
                     ,  ,  ,  ,  ,  ,  ;
unCom = 0.0001,  , 0.0026, 0.0056, 0.0079,       , 0.02,
               ,  ,  ,  ,  ,  ,  ;
normCom = 0.0001,  , 0.0026, 0.0056, 0.0079,  , 0.02,
               ,  ,  ,  ,  ,  ,  ;
com = 0.0001,  , 0.0026, 0.0056, 0.0079,  , 0.02,
         ,  ,  ,  ,  ,  ,  ;
veryCom = 0.0001,  , 0.0026, 0.0056, 0.0079,  , 0.02,
               ,  ,  ,  ,  ,  ,  ;
ENDDATA

DATA:
EPS = 1E-8;
ENDDATA

CALC:
@FOR(PWLFUNS(I, J) | J #LE# 2: veryUnCom(2, J) = 1 - EPS);
@FOR(PWLFUNS(I, J) | J #GE# 3: veryUnCom(2, J) = EPS);
@FOR(PWLFUNS(I, J) | J #EQ# 3: unCom(2, J) = 1 - EPS);
@FOR(PWLFUNS(I, J) | J #NE# 3: unCom(2, J) = EPS);
@FOR(PWLFUNS(I, J) | J #EQ# 4: normCom(2, J) = 1 - EPS);
@FOR(PWLFUNS(I, J) | J #NE# 4: normCom(2, J) = EPS);
@FOR(PWLFUNS(I, J) | J #EQ# 5: com(2, J) = 1 - EPS);
@FOR(PWLFUNS(I, J) | J #NE# 5: com(2, J) = EPS);
@FOR(PWLFUNS(I, J) | J #LE# 5: veryCom(2, J) = EPS);
@FOR(PWLFUNS(I, J) | J #GE# 6: veryCom(2, J) = 1 - EPS);
ENDCALC

MIN = Z;
Z = -COSTA - COSTB - COSTC - COSTD - COSTE;

@FOR( COST:
SAC = 1 - SA;
SBC = 1 - SB;
SCC = 1 - SC;
SDC = 1 - SD;
SEC = 1 - SE;
@BND( EPS, SA, 1 - EPS );
@BND( EPS, SB, 1 - EPS );
@BND( EPS, SC, 1 - EPS );
@BND( EPS, SD, 1 - EPS );
@BND( EPS, SE, 1 - EPS );
@BND( EPS, SAC, 1 - EPS );
@BND( EPS, SBC, 1 - EPS );
@BND( EPS, SCC, 1 - EPS );
@BND( EPS, SDC, 1 - EPS );
@BND( EPS, SEC, 1 - EPS );
   );

COSTA = @SUM( COST: SA*@LOG( SA ) + SAC*@LOG( SAC ) );
COSTB = @SUM( COST: SB*@LOG( SB ) + SBC*@LOG( SBC ) );
COSTC = @SUM( COST: SC*@LOG( SC ) + SCC*@LOG( SCC ) );
COSTD = @SUM( COST: SD*@LOG( SD ) + SDC*@LOG( SDC ) );
COSTE = @SUM( COST: SE*@LOG( SE ) + SEC*@LOG( SEC ) );
@FREE(COSTA); @FREE(COSTB); @FREE(COSTC); @FREE(COSTD); @FREE(COSTE); @FREE(Z);

@FOR( COST( I ): X( I )  = @SUM( PWLPOINTS( J ): veryUnCom( 1, J )*W( I, J ) ) );
@FOR( COST( I ): SA( I ) = @SUM( PWLPOINTS( J ): veryUnCom( 2, J )*W( I, J ) ) );
@FOR( COST( I ): SB( I ) = @SUM( PWLPOINTS( J ): unCom( 2, J )*W( I, J ) ) );
@FOR( COST( I ): SC( I ) = @SUM( PWLPOINTS( J ): normCom( 2, J )*W( I, J ) ) );
@FOR( COST( I ): SD( I ) = @SUM( PWLPOINTS( J ): com( 2, J )*W( I, J ) ) );
@FOR( COST( I ): SE( I ) = @SUM( PWLPOINTS( J ): veryCom( 2, J )*W( I, J ) ) );
@FOR( COST( I ): @SUM( PWLPOINTS( J ): W( I, J ) ) = 1 );
@FOR( WEIGHTS( I, J ): @SOS2( 'SOSSETS_' + FEELING( I ), W( I, J ) ) );

veryUnCom( 1, 2 ) - unCom( 1, 2 ) = 0;
unCom( 1, 2 ) - normCom( 1, 2 ) = 0;
normCom( 1, 2 ) - com( 1, 2 ) = 0;
com( 1, 2 ) - veryCom( 1, 2 ) = 0;
veryUnCom( 1, 6 ) - unCom( 1, 6 ) = 0;
unCom( 1, 6 ) - normCom( 1, 6 ) = 0;
normCom( 1, 6 ) - com( 1, 6 ) = 0;
com( 1, 6 ) - veryCom( 1, 6 ) = 0;
@BND( 0.0018, veryUnCom( 1, 2 ), 0.0035 );  @BND( 0.0084, veryUnCom( 1, 6 ), 0.0132 );

END

回复此楼

» 本帖@通知

@NicoleLas

» 猜你喜欢

AI论文写作工具：是科研加速器还是学术作弊器？已经有5人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有7人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

1楼2018-12-26 11:42:19

已阅关注TA 给TA发消息送TA红花 TA的回帖

alex110

铁杆木虫 (正式写手)

2楼2018-12-26 14:34:33

已阅关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 alex110 的主题更新

返回列表