代码问题

lupkpu9

代码如下：
clear all;
* C7 s# r# \2 H. Sclc;
7 o( N% J5 O* T+ |; Srandn('state',0);
rand('state',0);
dt = 0.005;
' }7 {3 z2 c3 a$ B! Z# i' ~J = 1000000;
for h = 1:100
    x = normrnd(43.68,39.85,[1,J]);%设定分布：正态分布,Avoided fertilizers application
    x1 = find(x>0);
    x2 = x(x1);
    X(h, = x2;
! F8 {+ G3 c: |8 R    A(h, = mean(X(h,);
! X: R, K) z- H# k7 w! P    B(h,:) = std(X(h,:));
+ r7 i% g3 {5 \6 W    y1 = X(:);
    Y1 = sort(y1);
3 W: S: n; N. C5 t        W = length(x1);
    if rem(W,2)==0
7 J6 n2 W, h! o# g+ x# H- }, ~        W1 = W;
    else
8 B. ]' H$ E0 `6 X        W1 = W+1;
    end
    for r = 1:W1 * 5/100
        W2(:,r) = X(W1 * 95/100 + r) - X(r);
% D, v9 I+ q/ [2 ^9 [* z    end
! f5 y$ O1 y8 |" p" z' T' z    n1 = find(W2 == min(W2));
    C(h,:) = Y1(n1);
    D(h,:) = Y1(J * 95/100 + n1);
" T) p# g/ [- h& ?    if h > 1
) F: T8 s0 k4 N" B- L& _        A1 = 2 * std(A)/sqrt(h);
9 X3 o7 H; {- ^6 t7 b- Y        B1 = 2 * std(B)/sqrt(h);
8 D" {; e  o9 ~, \; @) b        C1 = 2 * std(C)/sqrt(h);
; e) b0 z7 I& z3 ]- |3 u        D1 = 2 * std(D)/sqrt(h);
        if (A1<=dt&B1<=dt&C1<=dt&D1<=dt) == 1
( _  M! l/ i* [) |! L            break
7 \3 u/ C& M: d3 ]- l        end
    end
end


8 x& J( f- y" U0 U- Y: ?


上述代码运行后出现“下标索引必须为正整数类型或逻辑类型。”  为什么
2 ^( Y8 K: N6 w5 C4 u9 o. @: G

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
W2(:,r) = X(W1 * 95/100 + r) - X(r);
这句的问题，你的代码前面算出来的W1，并不能够被100整除，导致 W1 ...

9 O  w4 y" x- Y7 {3 c2 o+ T谢谢，我去试一下
3 _# A+ e+ g' y