你好,你这个代码其中问题很多。注意到你的heavy_wayfact中的数据初始时全部为零。/ h8 `; L1 I3 t& ^+ ~/ N8 u
- int s = duniform(1,7);
) n% S; }& b& ?/ [ - Table L_way = Table("heavy_way");
; F# y9 X) p6 G t7 G5 ? - Table L_factway = Table("heavy_wayfact");
& v# U, O6 s/ z
2 @! o4 V0 \, |3 u4 M- o- int x = 0;4 ~2 K$ r/ s$ t4 o/ g
- int index =0;1 @5 W& G% k3 I6 Y) |
- double m = 0;( v, d. J. h: u: L( i. {3 y3 z
- for(int i=1;i<L_way.numRows;i++)
: R- h8 V! `8 n; h - {7 R: h+ d# m3 |3 _0 R4 z. Q( o! E
- double flowrate = L_way[i][s];
4 [$ B) n$ p- w, A - double lfact = L_factway[i][s]; // 每次你这样都是为0,因为你没改变过
1 g2 t" k9 e9 z4 z# O. A3 p' ? - if(flowrate!=0) // 找到非0值' q) D' [, \! m/ [7 n
- {
- O5 x& I0 _" k0 | R* c+ ?( N) J - index++;
+ _, @. U' h2 B% f - if(index==1) // 当找到第一个flowrate值不为0
) X+ u# |3 ~9 V# g' X5 L - {! e5 k+ H ~: y7 a, U' |: V
- x=i;
. \! R' g2 m; }% z' V' {' ]% X - m=lfact/flowrate; // m还是0,因为lfact始终为0: v2 B/ Y: E" a) D/ e
- }# j( c( T8 i( R+ K+ u
- if(m>lfact/flowrate) // 所以这个条件始终是没用的6 b. L. J1 T7 }8 r4 y
- {
2 h# g0 O/ V7 K3 W - x=i;
* ~" M) k0 ?. h, v - m=lfact/flowrate;
, X6 ]; D0 _9 m" ?/ T: O - }1 t7 Y2 y& b" q- |) y" Q
- }/ u$ l0 G2 Y! f: [. ~
- }
复制代码 |