這段程序出錯，不清楚怎么改

mouseyuanhao

下面是我編的M文件的程序，主要是求一維數組s的最大值和索引位置，然后加不同寬度的窗，求出窗里信號能量后與總能量求比值。我編的程序中s是小波包分析的高頻信號，為了簡化程序段，我用5000個隨機數來代替。編譯的時候for循環好像根本就沒有執行，不知道是為什么。請前輩們幫忙看看，不勝感激！
s=zeros(1,5000);
/ C( {2 W! ^/ f# @$ p% ws=rand(1,5000);
- [# ^- E, }" Q) Q$ u; m! T, [2 R4 F  Zn=length(s);
max=0.0;
$ q! l' i3 H* aindex=0;
2 j: T% u# J& u  yfor i=1:n
3 W/ M" e, m+ j5 W if  max<s(i)
     max=s(i);
     index=i;
. S: f: O/ A& H& b  else
  end
end   
9 S$ Z' g9 J1 E" m1 gsmall=100;
( v  R" ^$ u% K# N. A8 Qwenergys=sum(s(index-small:index+small).^2);
) ?  V2 ]1 J( M. e; r/ `; Imiddle=150;
wenergym=sum(s(index-middle:index+middle).^2);
1 |; ?- k" H( w  m$ y' Q4 |; j* I- dlarge=200;
wenergyl=sum(s(index-large:index+large).^2);
energy=sum(s.^2);
# ^. H; o! B& D( x' nratios=wenergys/energy;
8 I. [, R7 o+ ?- N6 v; f! Q6 @ratiom=wenergym/energy;
ratiol=wenergyl/energy;
" I9 x+ c5 u8 h$ j8 Mm=max*100000;
ws=wenergys*100000;
wm=wenergym*100000;
9 V- }. @3 I8 p0 W- ~) s5 Cwl=wenergyl*100000;
( q# {0 j5 l0 x+ je=energy*100000;
a=[m,index,e,ws,wm,wl,ratios,ratiom,ratiol];

小土匪

求一維數組s的最大值和索引位置 max 和index可以求出來呀
& l0 H4 E+ p. r/ E也就是執行了for循環了
& |' Z7 Y  @* Y; B  [可能計算機計算速度快  我們感覺不到等待的時間   
如果是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還沒接觸過
% }% N) ~; [) O( d$ E" R是不是取隨機數的時候  最大值的標號正好小于200了  這樣的話wenergyl=sum(s(index-large:index+large).^2)這步是不是會出錯？

2 P% m! M3 L. \7 y提示也正好是indices must either be real positive integers or logicals
你再修改下隨即數？

小土匪

Subscript indices must either be real positive integers or logicals
/ V9 F. O( l; ^4 R! xR$ P# O9 A3 s
1 G8 U) h1 N- J2 L) y" h中文解釋：下標索引必須是正整數類型或者邏輯類型
7 W2 P3 E7 E' }; N1 N) a( t; D1 KW7 b# p# ]3 I' k8 m
出錯原因：在訪問矩陣（包括向量、二維矩陣、多維數組，下同）的過程中，下標索引要么從 0 開始，要么出現了負數。注：matlab 的語法規定矩陣的索引從 1 開始
: g$ n  d1 t" I8 B' I6 d

! N' `. k3 F7 b這是網上搜到的資料   借鑒下