代码问题

lupkpu9

代码如下：
clear all;
8 W* u: H2 d: V  Gclc;
, G* Y6 Y" @. m' g) M" zrandn('state',0);
9 R% Q5 q. {( ^6 Nrand('state',0);
dt = 0.005;
1 ~2 x3 n1 H% t6 P1 c; HJ = 1000000;
, {: C* M+ f. P5 xfor h = 1:100
    x = normrnd(43.68,39.85,[1,J]);%设定分布：正态分布,Avoided fertilizers application
8 v! O7 O3 S% K7 n5 O6 k" c    x1 = find(x>0);
  T' G2 ?" h1 I$ p8 m& ~    x2 = x(x1);
- C" ~8 P' v# {: \' h$ P4 p5 w$ |    X(h, = x2;
    A(h, = mean(X(h,);
    B(h,:) = std(X(h,:));
0 H3 Y; g9 Y5 I$ H( z5 l+ f    y1 = X(:);
    Y1 = sort(y1);
& u; y- n& t2 {        W = length(x1);
    if rem(W,2)==0
        W1 = W;
7 i2 t  t4 x0 h# S    else
        W1 = W+1;
/ @+ [, w9 ]- z5 }" Y& w( ~. u9 |    end
9 h, H1 i9 w& H  r" s; j/ E    for r = 1:W1 * 5/100
8 V5 @# Z. z0 }6 q  k9 z* a        W2(:,r) = X(W1 * 95/100 + r) - X(r);
# r- Z3 f  @( I6 H* Q6 c    end
    n1 = find(W2 == min(W2));
3 O- D5 o8 Y6 r4 Q    C(h,:) = Y1(n1);
    D(h,:) = Y1(J * 95/100 + n1);
4 k* p; X, I, k/ Y    if h > 1
        A1 = 2 * std(A)/sqrt(h);
        B1 = 2 * std(B)/sqrt(h);
0 I) C7 @8 F6 p  C" W% j        C1 = 2 * std(C)/sqrt(h);
/ t, L5 [& u3 [  \5 g        D1 = 2 * std(D)/sqrt(h);
        if (A1<=dt&B1<=dt&C1<=dt&D1<=dt) == 1
9 x# v6 @/ Y, {' @            break
& ~) D" W+ E: U9 A( x, |! `        end
) @: M4 c9 [  x2 D* {( c    end
end


1 A* U3 t" J3 @


/ U, r4 a" C. Q1 U3 |( }上述代码运行后出现“下标索引必须为正整数类型或逻辑类型。”  为什么
; W$ A* d6 W: |8 C! J0 R

lupkpu9

kekek 发表于 2020-7-15 14:06
W2(:,r) = X(W1 * 95/100 + r) - X(r);
) I) Q0 B0 f8 s5 n, @5 W这句的问题，你的代码前面算出来的W1，并不能够被100整除，导致 W1 ...

谢谢，我去试一下
: L8 k6 K4 }) P7 Z5 e, u

zaiyiaaaa

kekek

W2(:,r) = X(W1 * 95/100 + r) - X(r);
6 g2 }9 Q1 P+ C3 M' Z这句的问题，你的代码前面算出来的W1，并不能够被100整除，导致 W1 * 95/100 就不是一个整数，打个比方，点名时候可以有3个人可以有30个人，但不可能有0.3个人，你这个索引值 W1 * 95/100 + r 不是整数，就相当于点名发现全班有0.3个人一样。

llbnmo

帮你顶一下