你好,你这个代码其中问题很多。注意到你的heavy_wayfact中的数据初始时全部为零。
6 ~+ N5 A1 G" m( q0 S8 D- int s = duniform(1,7);
1 y1 N& ]8 q) c+ ?+ y* L2 S7 Z - Table L_way = Table("heavy_way");2 W. Y$ z: V$ P! c6 v u
- Table L_factway = Table("heavy_wayfact");
, r U/ |% N8 G0 }$ N0 G" u% S
2 z% r- |5 ]/ s6 x- int x = 0;
& v* T" W+ G' Y; g, ~ - int index =0;
% i m$ h# V* V - double m = 0;! h2 m0 o3 ~; c6 s: Q' y+ B% e
- for(int i=1;i<L_way.numRows;i++)
f1 E. K- v' Q4 S" r" W5 u, } - {9 W2 w2 F7 T- h# h: w
- double flowrate = L_way[i][s];
5 O0 i6 `% K& e' p* A - double lfact = L_factway[i][s]; // 每次你这样都是为0,因为你没改变过; D6 l! V8 G" S% S1 z% h0 Q
- if(flowrate!=0) // 找到非0值& B, j7 @5 u" i A6 t+ g; u5 Z
- {
- `" A9 ^5 \3 f: _ - index++;# C4 _7 t7 \/ l# d) R) @ m
- if(index==1) // 当找到第一个flowrate值不为0% l7 P* q% U' z# [, C6 _1 n+ s
- {% Q# _& I% L8 S( ]
- x=i;
3 _$ L2 v' o& T6 H - m=lfact/flowrate; // m还是0,因为lfact始终为0- U- ^2 Q' t& b3 R& T2 n
- }
/ N1 b7 C4 L" H6 U2 Q# D - if(m>lfact/flowrate) // 所以这个条件始终是没用的
$ P. Y0 x1 v/ W# m% A, Q - {/ ~+ r- {' R* {# S" T% ]! Z& s
- x=i;5 ~* [( O6 z- S
- m=lfact/flowrate;1 w5 e; I% }0 _' z# Z8 D; l b
- }
& z$ b* w2 A5 \0 f: l - }
$ X3 {7 x9 K; G. V( F - }
复制代码 |