你好,你这个代码其中问题很多。注意到你的heavy_wayfact中的数据初始时全部为零。+ p1 c& g+ {- u6 W3 [7 a/ n5 @+ P
- int s = duniform(1,7);' X a, E+ S4 \
- Table L_way = Table("heavy_way");" {) o# ^" [1 c& P2 C3 W: ]* W
- Table L_factway = Table("heavy_wayfact");+ G9 Y# G. f/ V( F( b
- / o. H/ K% V! c7 y. w( W
- int x = 0;
! g: m; ^/ _4 p - int index =0;
$ B* _$ D2 u2 r8 R9 G+ n - double m = 0;! ^+ F" z% t. Y6 }5 r) C
- for(int i=1;i<L_way.numRows;i++)" B/ f6 M( P9 p- i! n) \& d! T* c. U
- {
$ p0 m( H' S" x& L+ u2 Q - double flowrate = L_way[i][s];
6 ]4 R3 l& |, _ - double lfact = L_factway[i][s]; // 每次你这样都是为0,因为你没改变过 r0 k0 v+ ^* N, J( b
- if(flowrate!=0) // 找到非0值
! ~: D) G8 _% g" ^& d4 \ - {% n, S& y8 u! l; T2 z5 ?! c& S
- index++;! M" k1 Z8 e# A. ?4 p! j9 ^
- if(index==1) // 当找到第一个flowrate值不为0# }4 @. Q5 f! C# `8 P+ Y
- {
/ m* ?/ e" B) G - x=i;: |+ s- y3 x+ F. }+ n. K+ {, [
- m=lfact/flowrate; // m还是0,因为lfact始终为0
9 _' q. F, x& W3 l/ t: V - }
- l" g, h7 x- L( M& c% j - if(m>lfact/flowrate) // 所以这个条件始终是没用的2 f% M8 E& H7 R& q5 S
- {
- b6 h; \# G2 g& z0 b - x=i;
. K" W) y3 ^; ^* O - m=lfact/flowrate;
8 t7 B9 p) l+ k9 v5 F9 i - }2 [! Q1 a# [! {1 z0 n" K
- }
2 [- `: L; Z! e, A# ]' d) X! S - }
复制代码 |