EDA365电子论坛网

标题: 代码问题 [打印本页]

---

作者: lupkpu9    时间: 2020-7-13 11:11
标题: 代码问题
代码如下：
clear all;
. x) j- ~+ ?9 Z7 F) Gclc;
. M. O( F: C) M# m* m( ]randn('state',0);
0 G" p9 h/ ^( v# V6 `2 Q4 Rrand('state',0);
dt = 0.005;
6 J( Z, k7 ~; F  XJ = 1000000;
# t6 S& n) ^0 v6 w7 ^for h = 1:100
    x = normrnd(43.68,39.85,[1,J]);%设定分布：正态分布,Avoided fertilizers application
; K9 U; c! ~0 g  V    x1 = find(x>0);
- K# B  Z. x/ v' ?    x2 = x(x1);
! W, x+ I9 C6 m0 j7 D    X(h, = x2;
    A(h, = mean(X(h,);
; T3 k9 m* n4 T, D    B(h,:) = std(X(h,:));
    y1 = X(:);
; o! E" j6 {; k9 i    Y1 = sort(y1);
        W = length(x1);
    if rem(W,2)==0
        W1 = W;
    else
        W1 = W+1;
    end
    for r = 1:W1 * 5/100
* ]. G3 |" R2 r, u) z        W2(:,r) = X(W1 * 95/100 + r) - X(r);
    end
    n1 = find(W2 == min(W2));
    C(h,:) = Y1(n1);
( W/ ]( ]1 p; O+ Q: o. T    D(h,:) = Y1(J * 95/100 + n1);
    if h > 1
( b4 i8 P  p$ q+ c/ b7 h        A1 = 2 * std(A)/sqrt(h);
' W9 a, M  [" i& Y        B1 = 2 * std(B)/sqrt(h);
        C1 = 2 * std(C)/sqrt(h);
) g  n. I  c6 h1 K( Q' b7 f        D1 = 2 * std(D)/sqrt(h);
        if (A1<=dt&B1<=dt&C1<=dt&D1<=dt) == 1
; @, O( n0 S' @; u, O            break
8 Q$ i0 ^+ R7 ?8 Y! ]: M3 E+ F# b: B        end
& X( y0 |8 L4 z! @) {# w    end
- e0 e2 x$ H+ W0 K# uend


" P; D3 ~7 U) t' {( Y2 L) O- g; @
# M  E5 d' n- p0 p" D' u
" |/ x5 T* F' j
上述代码运行后出现“下标索引必须为正整数类型或逻辑类型。”  为什么

---

作者: llbnmo    时间: 2020-7-13 13:07
帮你顶一下

---

作者: kekek    时间: 2020-7-15 14:06
W2(:,r) = X(W1 * 95/100 + r) - X(r);
& u7 v- t( ]* Z+ k: d3 \这句的问题，你的代码前面算出来的W1，并不能够被100整除，导致 W1 * 95/100 就不是一个整数，打个比方，点名时候可以有3个人可以有30个人，但不可能有0.3个人，你这个索引值 W1 * 95/100 + r 不是整数，就相当于点名发现全班有0.3个人一样。

---

作者: zaiyiaaaa    时间: 2020-7-15 15:53

---

作者: lupkpu9    时间: 2020-7-16 15:05

kekek 发表于 2020-7-15 14:06
$ e; c' o7 u2 ]0 R; E/ h+ UW2(:,r) = X(W1 * 95/100 + r) - X(r);
这句的问题，你的代码前面算出来的W1，并不能够被100整除，导致 W1 ...

谢谢，我去试一下

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2