代码问题

lupkpu9

代码如下：
clear all;
8 Q6 Q- o; L2 p" f. W2 {clc;
randn('state',0);
rand('state',0);
dt = 0.005;
J = 1000000;
9 d6 N' M  S" M! e+ E( r4 J* Ifor h = 1:100
    x = normrnd(43.68,39.85,[1,J]);%设定分布：正态分布,Avoided fertilizers application
    x1 = find(x>0);
: [3 F% C8 K3 o    x2 = x(x1);
% {. O+ i. y8 W) l: g    X(h, = x2;
    A(h, = mean(X(h,);
    B(h,:) = std(X(h,:));
: |" x6 X% u* w! e  O5 r    y1 = X(:);
    Y1 = sort(y1);
9 _  q5 P+ f: X$ c; V        W = length(x1);
    if rem(W,2)==0
        W1 = W;
1 ?" p3 @7 `6 g% @! h7 h/ d    else
        W1 = W+1;
    end
    for r = 1:W1 * 5/100
        W2(:,r) = X(W1 * 95/100 + r) - X(r);
; E  e6 w: e8 S1 u    end
    n1 = find(W2 == min(W2));
3 _9 z, w0 t3 u* l4 Y; W    C(h,:) = Y1(n1);
    D(h,:) = Y1(J * 95/100 + n1);
6 G/ i/ w  ?7 i' v    if h > 1
( C+ C" q( Q2 Q" ]        A1 = 2 * std(A)/sqrt(h);
8 q( F. o3 s) g! u# o+ [, H4 g        B1 = 2 * std(B)/sqrt(h);
1 h2 v9 B6 O& F$ ]# j; d        C1 = 2 * std(C)/sqrt(h);
  j  \# a1 P2 T3 d0 s( u        D1 = 2 * std(D)/sqrt(h);
. E6 a7 l' k1 N* h9 t9 T        if (A1<=dt&B1<=dt&C1<=dt&D1<=dt) == 1
6 U+ a8 L3 C7 q0 s$ B            break
        end
    end
end
1 D/ b+ ?& `" t3 u+ v

1 M! g  e) s7 N4 Z' S

0 i0 E$ |. w' b7 K
' J0 O1 P9 Z4 `# A! m+ B上述代码运行后出现“下标索引必须为正整数类型或逻辑类型。”  为什么
# L% x3 s5 E6 D/ v# z( l( J' o% [1 S

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

谢谢，我去试一下