代码问题

lupkpu9

代码如下：
% t5 A& P+ V2 V3 ^, Bclear all;
) w  h9 H+ \4 }7 uclc;
randn('state',0);
rand('state',0);
) m7 w& F2 U6 V: Y4 _6 E' I6 s2 Wdt = 0.005;
J = 1000000;
for h = 1:100
7 M' K# S9 F7 p8 S    x = normrnd(43.68,39.85,[1,J]);%设定分布：正态分布,Avoided fertilizers application
    x1 = find(x>0);
! k  U( u+ r$ x  o9 ^3 t    x2 = x(x1);
6 _% k. F5 r! U; j4 ~    X(h, = x2;
    A(h, = mean(X(h,);
    B(h,:) = std(X(h,:));
    y1 = X(:);
    Y1 = sort(y1);
( A8 a5 C* s5 C* @9 E5 B        W = length(x1);
* o; J# J0 V6 x    if rem(W,2)==0
        W1 = W;
    else
( R( @+ R) K6 l5 v  ]. u+ J        W1 = W+1;
, F6 }. x1 R! g0 l  u, c9 V# A    end
    for r = 1:W1 * 5/100
# D! Y; e+ M3 j1 T2 N4 z        W2(:,r) = X(W1 * 95/100 + r) - X(r);
. F8 S# `; W" [* v6 S    end
3 f, w! X+ R; C  f. Y! g$ P# E    n1 = find(W2 == min(W2));
    C(h,:) = Y1(n1);
3 u& n4 d2 ~& J" |    D(h,:) = Y1(J * 95/100 + n1);
1 E8 Y" c) k7 s  f    if h > 1
2 h# u4 ~. x) ]  Y/ ?0 f$ P2 ?5 l        A1 = 2 * std(A)/sqrt(h);
2 q& a  \8 T; N' R1 K. b/ [        B1 = 2 * std(B)/sqrt(h);
2 B5 d, S( k' O& p! [        C1 = 2 * std(C)/sqrt(h);
  f6 Y1 P2 Z* b/ c0 z2 E9 c        D1 = 2 * std(D)/sqrt(h);
        if (A1<=dt&B1<=dt&C1<=dt&D1<=dt) == 1
) j7 I8 \/ @- `' X3 x            break
3 p' N- a+ W1 e! J7 h! ~        end
. O1 X3 g) O2 I1 v9 H- ?# x# c) k% P    end
end
; J0 O% l7 p9 {. g9 w


! T0 f4 ?7 ]& z- W

上述代码运行后出现“下标索引必须为正整数类型或逻辑类型。”  为什么

llbnmo

帮你顶一下

kekek

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

zaiyiaaaa

lupkpu9

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

0 Y+ |( Q( N; J+ R4 _谢谢，我去试一下
2 n" n0 `- a# x8 s* t/ a+ A