代码问题

lupkpu9

代码如下：
2 B$ p! ]$ j* V1 d. d5 |& w/ h/ `clear all;
' c) @( l, ]$ l5 Q( \- Aclc;
randn('state',0);
rand('state',0);
dt = 0.005;
* Y% k9 u# w1 \7 T$ IJ = 1000000;
for h = 1:100
    x = normrnd(43.68,39.85,[1,J]);%设定分布：正态分布,Avoided fertilizers application
0 R1 B: {/ n  _3 @    x1 = find(x>0);
3 @' F% k- Y6 |' n8 m3 l' O  R9 Z    x2 = x(x1);
( Y8 b. ~9 i! O6 K    X(h, = x2;
    A(h, = mean(X(h,);
    B(h,:) = std(X(h,:));
; ]4 W* l% ?$ K6 K7 E    y1 = X(:);
8 h& c( E; A' p- j' Q7 K) A% K    Y1 = sort(y1);
        W = length(x1);
' M( E; U& n, }    if rem(W,2)==0
        W1 = W;
    else
        W1 = W+1;
    end
+ L) `/ V. c& b    for r = 1:W1 * 5/100
        W2(:,r) = X(W1 * 95/100 + r) - X(r);
    end
    n1 = find(W2 == min(W2));
    C(h,:) = Y1(n1);
    D(h,:) = Y1(J * 95/100 + n1);
    if h > 1
# {- z. S6 e$ _$ B" o3 H* n        A1 = 2 * std(A)/sqrt(h);
        B1 = 2 * std(B)/sqrt(h);
        C1 = 2 * std(C)/sqrt(h);
        D1 = 2 * std(D)/sqrt(h);
        if (A1<=dt&B1<=dt&C1<=dt&D1<=dt) == 1
            break
, n# k9 W( ]/ N3 D' l3 K4 G        end
( S8 m$ b+ R* R* ~1 U3 P    end
end




7 y5 {, S; U6 z: b
上述代码运行后出现“下标索引必须为正整数类型或逻辑类型。”  为什么

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 ...

谢谢，我去试一下
7 ^/ K! L- [/ y4 w! k