你好,你这个代码其中问题很多。注意到你的heavy_wayfact中的数据初始时全部为零。
! H! N6 x8 ?6 C- int s = duniform(1,7);5 ?: l* X, X* y/ ~' C0 S; G
- Table L_way = Table("heavy_way");
# ^1 N' [8 V; R3 D4 W) i$ y% |( i - Table L_factway = Table("heavy_wayfact");
4 `2 q3 g3 \! d) H - & p$ M. |; K- u* A+ P h6 S
- int x = 0;
7 r* N5 V7 A* E! {- a8 q - int index =0;
2 e9 k: `8 `% R - double m = 0;! [3 }! s# L6 M! }( W4 g, B1 g
- for(int i=1;i<L_way.numRows;i++)/ I S1 ?. M. m- c, q
- {
8 {- x% v9 \* q, J - double flowrate = L_way[i][s];# H' i& ]$ e3 p
- double lfact = L_factway[i][s]; // 每次你这样都是为0,因为你没改变过
8 R# S4 i6 ^1 s/ @9 y4 G# ` - if(flowrate!=0) // 找到非0值7 L* t4 Q/ X7 }0 u" F
- { O0 G6 K; m7 s1 j
- index++;: Y/ x; f1 }! d' T
- if(index==1) // 当找到第一个flowrate值不为0: | K4 z; T/ S8 C+ K4 p- {
- {
) x! i1 `, m8 O4 h9 R - x=i;
4 s1 J8 E3 Q9 W6 W - m=lfact/flowrate; // m还是0,因为lfact始终为0
k0 [# r- H* D' L$ D1 R, S" x - }$ t9 y i8 ]; E6 p' q
- if(m>lfact/flowrate) // 所以这个条件始终是没用的
3 ]6 f- }) W5 j5 f( O9 G$ x1 i) P - {
, H5 m" M; q$ |' a1 _ - x=i;) C% Y8 \$ ^ o1 O
- m=lfact/flowrate;9 Z7 P6 Y) o4 D: `2 Y8 a
- }% s/ a; G/ z5 p! }$ G) C' y
- }' `3 Y0 h; x8 ]
- }
复制代码 |