你好,你这个代码其中问题很多。注意到你的heavy_wayfact中的数据初始时全部为零。! Z+ d1 w, b4 H7 K3 S
- int s = duniform(1,7);
, i' e4 m) z% v8 `. z1 ^ U) m! M - Table L_way = Table("heavy_way");
: q" ` t3 N1 d - Table L_factway = Table("heavy_wayfact");
4 ?5 h7 Q( H6 l3 U: U& F - & Y$ H3 A9 I$ v9 l
- int x = 0;
2 j' B. ?+ i) @0 i T - int index =0;2 ? Y6 Z$ }- G9 W/ l
- double m = 0;- n' @; J4 @; ~4 F
- for(int i=1;i<L_way.numRows;i++)# h+ P" Z' l% _* n4 _
- {
. a- p2 d) x1 z5 a/ ~1 z - double flowrate = L_way[i][s];
1 Z2 d& `! j7 T - double lfact = L_factway[i][s]; // 每次你这样都是为0,因为你没改变过
4 W- Z; v% r$ a& P0 u: ^% l - if(flowrate!=0) // 找到非0值8 J( C: j/ y, X5 L- D3 k" I
- {; i" P0 [$ i V# ]4 `0 c
- index++;7 }4 W0 J* I1 p' N! c9 [
- if(index==1) // 当找到第一个flowrate值不为0% o0 {9 C% ]) d$ z8 a
- {% f3 _9 ]4 y9 b2 `3 Q# P( b9 ^
- x=i;
$ n' F1 l& S* h% y" E* u+ y' j8 ]6 X - m=lfact/flowrate; // m还是0,因为lfact始终为0
5 S! f: h0 _/ X- Z% i - }
7 l. K( s. \9 b0 I- J - if(m>lfact/flowrate) // 所以这个条件始终是没用的9 Y0 S6 {1 y' @ l4 N* }. ?
- {
7 |/ ^8 \ W+ G - x=i; U9 u. `8 i p- T: I8 y9 J7 Q
- m=lfact/flowrate;: H1 F( H" g; m0 h
- }
8 I, t. K. ?! |9 ` - }: }5 D$ G }% n9 G, u, z) {8 Y3 `4 D
- }
复制代码 |