24小时热门版块排行榜

返回列表

alex110

铁杆木虫 (正式写手)

应助: 0 (幼儿园)
金币: 6483
散金: 705
帖子: 308
在线: 868.9小时
虫号: 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

» 猜你喜欢

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

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

alex110

铁杆木虫 (正式写手)

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

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

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

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 复试调剂 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考研] 0703化学调剂，六级已过，有科研经历 +11	曦熙兮 2026-03-15	11/550	2026-03-19 17:13 by yywzz
[考研] 材料与化工求调剂 +7	为学666 2026-03-16	7/350	2026-03-19 14:48 by 尽舜尧1
[考研] 求调剂 +3	Mqqqqqq 2026-03-19	3/150	2026-03-19 14:11 by peike
[考研] 328求调剂，英语六级551，有科研经历 +4	生物工程调剂 2026-03-16	12/600	2026-03-19 11:10 by 生物工程调剂
[考研] 求材料调剂 +10	隔壁陈先生 2026-03-12	10/500	2026-03-19 09:15 by 无际的草原
[考研] 0703化学 305求调剂 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 0703化学调剂，求各位老师收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 085700资源与环境308求调剂 +3	墨墨漠 2026-03-18	3/150	2026-03-18 22:35 by bingxueer79
[考研] 354求调剂 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 301求调剂 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 材料工程专硕274一志愿211求调剂 +6	薛云鹏 2026-03-15	6/300	2026-03-17 11:05 by 学员h26Tkc
[考研] 考研调剂 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 一志愿，福州大学材料专硕339分求调剂 +3	木子momo青争 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 283求调剂 +3	听风就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 热情沙漠
[考研] 11408 一志愿西电，277分求调剂 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 东南大学364求调剂 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 321求调剂 +5	大米饭！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 招收0805（材料）调剂 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、