你好,你这个代码其中问题很多。注意到你的heavy_wayfact中的数据初始时全部为零。/ F. j4 s8 p5 G$ a9 g$ v1 e6 B
- int s = duniform(1,7);: F0 p- a+ e9 @1 X. |( i* |5 c. g$ a8 ~
- Table L_way = Table("heavy_way");
5 `- O6 x/ h" a! ]# K - Table L_factway = Table("heavy_wayfact");
* B8 i' |9 T" ]6 l, w; G
$ e0 y+ r5 l3 B( ?. S- int x = 0;
# o* d" k8 n! l( {$ {4 F - int index =0;7 d7 y( R# S5 O# K ^) A2 P$ P; X
- double m = 0;9 v+ q. U1 Z$ f7 h8 Q) [$ k
- for(int i=1;i<L_way.numRows;i++)
8 ~" E2 P9 x# ^9 \& @& n7 u: s - {
* l- `5 R( _8 d1 O2 ? - double flowrate = L_way[i][s];% y0 V3 F) f0 g
- double lfact = L_factway[i][s]; // 每次你这样都是为0,因为你没改变过
. P0 ~5 `( b* L1 \: X8 d! F - if(flowrate!=0) // 找到非0值8 J2 M2 @: I2 |9 M7 [
- {5 u( m8 X1 J) M& |' w7 s+ K
- index++;: a( m4 l( a* b, j
- if(index==1) // 当找到第一个flowrate值不为0
9 Q! R1 c2 Z( a5 N4 n - {
3 ]' B4 O2 K5 y - x=i;( Z, a) A3 p3 F5 m2 P
- m=lfact/flowrate; // m还是0,因为lfact始终为0
( d# n6 o% q) r* X9 i - }
0 y4 ^4 s9 i' q+ A' n - if(m>lfact/flowrate) // 所以这个条件始终是没用的
2 x& @ ~( ]! u4 @6 z w" \ - {% k: q/ k6 A- N; G1 I. [3 W
- x=i;0 B( d' t/ t" ]% ]
- m=lfact/flowrate;# m! i# \3 y c. b; F: s3 I
- }. j4 I$ E7 H" h0 [3 w( Q
- }" _9 R2 E6 [ B6 E1 B
- }
复制代码 |