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

mouseyuanhao

下面是我編的M文件的程序，主要是求一維數組s的最大值和索引位置，然后加不同寬度的窗，求出窗里信號能量后與總能量求比值。我編的程序中s是小波包分析的高頻信號，為了簡化程序段，我用5000個隨機數來代替。編譯的時候for循環好像根本就沒有執行，不知道是為什么。請前輩們幫忙看看，不勝感激！
, Y5 B' j# L. M' |3 Ds=zeros(1,5000);
" S; Y3 }1 i& n$ P3 os=rand(1,5000);
" S3 n) G  i8 [% k% k* gn=length(s);
7 z0 S* A- n: ]) @8 w# _, Cmax=0.0;
index=0;
6 W) l0 C! A7 m6 _& \% F/ r! ufor i=1:n
$ g! P$ K* m2 k6 Q if  max<s(i)
. B/ p; U1 u4 l: W$ p3 G/ i     max=s(i);
7 M9 s; ]- X2 A. n# \+ s- E+ p     index=i;
' k# y8 z% D# K: r" m9 D  else
  end
end   
small=100;
wenergys=sum(s(index-small:index+small).^2);
middle=150;
wenergym=sum(s(index-middle:index+middle).^2);
' p) u2 c' ?+ b( I+ M1 p3 Glarge=200;
3 ?4 L( ?. ]2 d% J% g1 cwenergyl=sum(s(index-large:index+large).^2);
energy=sum(s.^2);
ratios=wenergys/energy;
ratiom=wenergym/energy;
ratiol=wenergyl/energy;
m=max*100000;
ws=wenergys*100000;
+ A: |+ q+ H& j1 c1 Rwm=wenergym*100000;
5 R% c' M6 c7 ~+ Bwl=wenergyl*100000;
e=energy*100000;
a=[m,index,e,ws,wm,wl,ratios,ratiom,ratiol];

小土匪

求一維數組s的最大值和索引位置 max 和index可以求出來呀
' N- J5 m; f/ ]也就是執行了for循環了
4 T  `! x' m: C# a- ^. c4 ?可能計算機計算速度快  我們感覺不到等待的時間   
如果是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還沒接觸過
是不是取隨機數的時候  最大值的標號正好小于200了  這樣的話wenergyl=sum(s(index-large:index+large).^2)這步是不是會出錯？
+ ^. V* T+ D& y( V5 {
提示也正好是indices must either be real positive integers or logicals
你再修改下隨即數？

小土匪

Subscript indices must either be real positive integers or logicals
. B, N5 m& P# r+ f' ER$ P# O9 A3 s
# x: O. j: t7 B# f0 F/ t" e中文解釋：下標索引必須是正整數類型或者邏輯類型
W7 b# p# ]3 I' k8 m
出錯原因：在訪問矩陣（包括向量、二維矩陣、多維數組，下同）的過程中，下標索引要么從 0 開始，要么出現了負數。注：matlab 的語法規定矩陣的索引從 1 開始
9 j/ `: _5 v* w& F. w6 O$ ^

這是網上搜到的資料   借鑒下