代码问题

lupkpu9

代码如下：
clear all;
& T) t0 K/ T! ~  v* v, m5 m, yclc;
" u* m5 P8 e8 E: srandn('state',0);
5 m; o7 V! b5 {- D8 S8 B2 \+ Lrand('state',0);
dt = 0.005;
J = 1000000;
- k5 u2 U! ^% _' \# R3 O7 Y5 ~for h = 1:100
    x = normrnd(43.68,39.85,[1,J]);%设定分布：正态分布,Avoided fertilizers application
/ X- x5 T% k& K, l    x1 = find(x>0);
    x2 = x(x1);
$ ?. z) q8 W6 v7 j& c+ i) F4 z    X(h, = x2;
' O5 \8 h. ?8 T/ f* t* z" h    A(h, = mean(X(h,);
    B(h,:) = std(X(h,:));
( D2 S* R( u5 S' ?# T    y1 = X(:);
    Y1 = sort(y1);
8 O3 e* y/ Q" D7 N% A        W = length(x1);
    if rem(W,2)==0
6 x- C( j! t( i5 X0 ^        W1 = W;
    else
        W1 = W+1;
' i% B# }7 y( `% o    end
# R9 |: d3 j6 G& x    for r = 1:W1 * 5/100
; C5 `1 [6 h" ~* \: E! c) ]  L$ y        W2(:,r) = X(W1 * 95/100 + r) - X(r);
% ~9 H" f! Z- P( x* @    end
' m/ O  G" q0 P9 N' |" i6 H; `* t    n1 = find(W2 == min(W2));
6 w. P& T1 _7 }    C(h,:) = Y1(n1);
    D(h,:) = Y1(J * 95/100 + n1);
0 |+ [6 w$ |. V1 {$ i. w0 a$ w    if h > 1
- W6 F0 Z& N) P4 c. _) n! d        A1 = 2 * std(A)/sqrt(h);
        B1 = 2 * std(B)/sqrt(h);
# W  t1 z+ M2 g7 R        C1 = 2 * std(C)/sqrt(h);
* K6 A- q/ o5 M. G* S$ @  t        D1 = 2 * std(D)/sqrt(h);
        if (A1<=dt&B1<=dt&C1<=dt&D1<=dt) == 1
            break
        end
    end
8 I' e3 E( F' I- x* Oend
# o+ C3 B; {( X& C, Z
! H* Z- h, x. `

# ?1 e7 O4 o" Z* J! ~5 j
: i* O! S5 ^+ t4 p: M
4 ~; y* }4 p5 H# }2 h3 t% @3 [上述代码运行后出现“下标索引必须为正整数类型或逻辑类型。”  为什么

llbnmo

帮你顶一下

kekek

W2(:,r) = X(W1 * 95/100 + r) - X(r);
8 u9 M5 @. T  m这句的问题，你的代码前面算出来的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);
8 \* V3 c. C7 T" x3 \这句的问题，你的代码前面算出来的W1，并不能够被100整除，导致 W1 ...

3 F! Y( k' K! J9 f/ q2 U" V谢谢，我去试一下