这段程序出错，不清楚怎么改

mouseyuanhao

下面是我编的M文件的程序，主要是求一维数组s的最大值和索引位置，然后加不同宽度的窗，求出窗里信号能量后与总能量求比值。我编的程序中s是小波包分析的高频信号，为了简化程序段，我用5000个随机数来代替。编译的时候for循环好像根本就没有执行，不知道是为什么。请前辈们帮忙看看，不胜感激！
s=zeros(1,5000);
s=rand(1,5000);
! [* P  X2 Z: E, E5 _* o8 wn=length(s);
7 V! [  v3 X) nmax=0.0;
index=0;
- j  _8 k3 C  wfor i=1:n
if  max<s(i)
( b4 r! h* z* l9 ~     max=s(i);
, I1 m! Q: S) ^9 K     index=i;
# p! A6 o* n# j5 }$ {4 H  else
  end
end   
small=100;
wenergys=sum(s(index-small:index+small).^2);
0 o. I1 [$ h, F2 Amiddle=150;
5 s) o3 u& m. J# Z$ r: gwenergym=sum(s(index-middle:index+middle).^2);
large=200;
( h- g6 T  {8 A) Kwenergyl=sum(s(index-large:index+large).^2);
energy=sum(s.^2);
$ k* _4 _+ o3 Yratios=wenergys/energy;
ratiom=wenergym/energy;
ratiol=wenergyl/energy;
8 f7 P6 `! w  J  {4 Jm=max*100000;
0 B+ ^- d( y' c: a0 v0 aws=wenergys*100000;
, q  i& @; d) Ywm=wenergym*100000;
' Y0 N3 j7 T$ N" ^6 e* v4 k$ _wl=wenergyl*100000;
- x: |1 e; k9 k9 W' M$ f1 fe=energy*100000;
a=[m,index,e,ws,wm,wl,ratios,ratiom,ratiol];

小土匪

求一维数组s的最大值和索引位置 max 和index可以求出来呀
, }/ K; r' A' a  s* T2 Y也就是执行了for循环了
2 c# {2 r( x; m可能计算机计算速度快  我们感觉不到等待的时间   
% R  Q1 ^6 @+ v4 z% S. a如果是5000个数字迭代的话肯能时间很长

mouseyuanhao

其实我的这个程序是插在labview里的脚本。我一开始用的就是[max,index]=max(s);但是我是把程序作为脚本插入在labview程序里，labview提示出错：subscript indices must either be real positive integers or logicals. 然后就不给我执行了，改成for循环也一样。5000个数字迭代时间并不长，我在MATLAB里面试过，瞬间出结果，但是在labview里编译就是通不过。我现在怀疑是不是labview和MATLAB的兼容做的不够好。

小土匪

labview还没接触过
! I- }: _1 |' P& `: r! j是不是取随机数的时候  最大值的标号正好小于200了  这样的话wenergyl=sum(s(index-large:index+large).^2)这步是不是会出错？
! R( u1 ~: }* F  O; q
提示也正好是indices must either be real positive integers or logicals
( c5 h/ x5 _3 P9 C$ z' @) O$ x你再修改下随即数？

小土匪

Subscript indices must either be real positive integers or logicals
R$ P# O9 A3 s
* ?0 b/ x; ^3 L' j- c; w# F中文解释：下标索引必须是正整数类型或者逻辑类型
W7 b# p# ]3 I' k8 m
) s, t1 j) }# c" Y! n5 ]出错原因：在访问矩阵（包括向量、二维矩阵、多维数组，下同）的过程中，下标索引要么从 0 开始，要么出现了负数。注：matlab 的语法规定矩阵的索引从 1 开始
2 T; F1 j; ^8 m* O

这是网上搜到的资料   借鉴下