代码问题

lupkpu9

代码如下：
clear all;
& m6 q" H1 p# X6 @! w$ [clc;
+ K0 b" {! |) T7 arandn('state',0);
6 E$ F1 h7 o$ |* T/ v! p! _rand('state',0);
" J; L' X  @5 _( ~4 g2 wdt = 0.005;
0 W3 x6 u1 {/ [6 G) P. BJ = 1000000;
  J" E6 I$ g( |0 g  q& Nfor h = 1:100
1 {( E/ f5 W! i! p    x = normrnd(43.68,39.85,[1,J]);%设定分布：正态分布,Avoided fertilizers application
    x1 = find(x>0);
, }* Y2 ?9 Q; E" M    x2 = x(x1);
    X(h, = x2;
! {/ ?6 P! N+ j) y4 S; b    A(h, = mean(X(h,);
& z3 R0 j/ b; M. K% k    B(h,:) = std(X(h,:));
    y1 = X(:);
    Y1 = sort(y1);
        W = length(x1);
    if rem(W,2)==0
- L* i7 q0 M. m3 U        W1 = W;
    else
7 U# f2 K# M8 j1 X9 L5 G6 [        W1 = W+1;
5 b5 E5 D0 a  C+ l- [    end
$ S6 m. Y/ t6 r! D    for r = 1:W1 * 5/100
        W2(:,r) = X(W1 * 95/100 + r) - X(r);
    end
    n1 = find(W2 == min(W2));
# `" Y$ z+ N6 B+ L: ~& V' t3 i( o    C(h,:) = Y1(n1);
    D(h,:) = Y1(J * 95/100 + n1);
    if h > 1
        A1 = 2 * std(A)/sqrt(h);
+ l, S/ z8 M+ L2 Y* ^        B1 = 2 * std(B)/sqrt(h);
! j; n+ s" A! A. I1 r        C1 = 2 * std(C)/sqrt(h);
+ h) Y9 N9 n% ]4 i        D1 = 2 * std(D)/sqrt(h);
        if (A1<=dt&B1<=dt&C1<=dt&D1<=dt) == 1
            break
: l  Q% T6 p9 U6 y, K        end
6 K, P0 v4 C3 ?( [! {9 X    end
end
! v$ u1 J2 F" T5 k8 y
& |/ |2 a" l# f9 ~9 M5 M8 U, Z$ }

- p2 p0 `. H& C% g
9 q: R, {& s0 a+ Q" Y1 T' i
8 }( x  f$ o3 g& B! i上述代码运行后出现“下标索引必须为正整数类型或逻辑类型。”  为什么

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
' v* q6 g) {, m7 i+ r- qW2(:,r) = X(W1 * 95/100 + r) - X(r);
$ U1 k/ w. M* m/ g+ D, C这句的问题，你的代码前面算出来的W1，并不能够被100整除，导致 W1 ...

* b/ _" j  e+ i) h$ V! L谢谢，我去试一下