代码问题

lupkpu9

代码如下：
clear all;
8 z, K; k9 `/ D2 Lclc;
randn('state',0);
rand('state',0);
dt = 0.005;
J = 1000000;
) e4 \* d* ^7 a' H7 a# j' W2 A& O5 `for h = 1:100
    x = normrnd(43.68,39.85,[1,J]);%设定分布：正态分布,Avoided fertilizers application
1 {& E( o: Z; A4 R- s* y. J$ `$ y    x1 = find(x>0);
( T- H! L7 w0 m* ]6 t3 D6 s$ u    x2 = x(x1);
3 N, f, Z, j: H% \    X(h, = x2;
: v& B6 o( l, k3 k2 @    A(h, = mean(X(h,);
$ c- W  r8 L$ d' x6 N    B(h,:) = std(X(h,:));
: v" R" A/ p& \2 |# z    y1 = X(:);
' @7 Q: p" x2 w1 u2 L* P0 g    Y1 = sort(y1);
. K$ U+ E; Q* U1 \$ _        W = length(x1);
! a- R. _! @: {) u: A/ {( L    if rem(W,2)==0
        W1 = W;
2 x( x- a* X2 ?9 r7 M) d# {% |    else
0 A) c, O/ K1 P; i! ^& m        W1 = W+1;
    end
" P8 H" F: |# w; o# }; o. |    for r = 1:W1 * 5/100
3 P% m% W. ~0 j        W2(:,r) = X(W1 * 95/100 + r) - X(r);
  S) ^( y6 g% e. H' f9 V" L# h' E    end
2 e) |# j3 W, ^- {6 v    n1 = find(W2 == min(W2));
5 t( J( c6 e0 i6 J/ j- ]0 ~    C(h,:) = Y1(n1);
    D(h,:) = Y1(J * 95/100 + n1);
0 p  v9 r: [9 b    if h > 1
        A1 = 2 * std(A)/sqrt(h);
: F& W  I0 s1 o) V. N  J        B1 = 2 * std(B)/sqrt(h);
. R6 t7 ?6 o! c" [        C1 = 2 * std(C)/sqrt(h);
, A' S+ c- g/ z  _' F% t1 d8 }        D1 = 2 * std(D)/sqrt(h);
3 K& p* V' b1 P" B6 c; Q8 S% x        if (A1<=dt&B1<=dt&C1<=dt&D1<=dt) == 1
: P, E1 v/ n- Z6 B            break
        end
& N' r9 T, B* ]* I1 d# [    end
8 {- b& Z+ o, J1 Y$ ~/ I/ ?end
1 @! `8 Z, i. l7 U

  u% O7 K1 y- \, w4 z4 f7 |

7 }  A1 m+ J; Z2 N# L
上述代码运行后出现“下标索引必须为正整数类型或逻辑类型。”  为什么
7 u' ]' p8 L2 l/ b6 v1 Q3 r& s" [: o

llbnmo

帮你顶一下

kekek

W2(:,r) = X(W1 * 95/100 + r) - X(r);
这句的问题，你的代码前面算出来的W1，并不能够被100整除，导致 W1 * 95/100 就不是一个整数，打个比方，点名时候可以有3个人可以有30个人，但不可能有0.3个人，你这个索引值 W1 * 95/100 + r 不是整数，就相当于点名发现全班有0.3个人一样。

zaiyiaaaa

lupkpu9

kekek 发表于 2020-7-15 14:06
; |8 E/ ?0 C& j; GW2(:,r) = X(W1 * 95/100 + r) - X(r);
这句的问题，你的代码前面算出来的W1，并不能够被100整除，导致 W1 ...

谢谢，我去试一下
6 g# K4 X1 t# n+ H% J) a& \* c2 r" V