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

mouseyuanhao

下面是我編的M文件的程序，主要是求一維數組s的最大值和索引位置，然后加不同寬度的窗，求出窗里信號能量后與總能量求比值。我編的程序中s是小波包分析的高頻信號，為了簡化程序段，我用5000個隨機數來代替。編譯的時候for循環好像根本就沒有執行，不知道是為什么。請前輩們幫忙看看，不勝感激！
+ C! X8 @$ s0 `s=zeros(1,5000);
s=rand(1,5000);
n=length(s);
max=0.0;
) X2 ?+ I, N7 \2 f+ k$ M: ^3 tindex=0;
for i=1:n
: g4 u6 T0 y: S9 A; y. s if  max<s(i)
     max=s(i);
     index=i;
1 q% p- `$ z/ t& X' M6 E' [  else
6 `2 ]  R6 m0 @  end
/ C2 m. W4 B8 A4 D3 c5 Wend   
small=100;
  }9 n* |" n0 t+ {4 K6 nwenergys=sum(s(index-small:index+small).^2);
8 m7 Y' c  Y8 F* F, C- Imiddle=150;
wenergym=sum(s(index-middle:index+middle).^2);
' v2 |5 P; b6 E% C3 M0 nlarge=200;
wenergyl=sum(s(index-large:index+large).^2);
energy=sum(s.^2);
% G9 z, M2 v5 V5 J( Z7 l+ t- |ratios=wenergys/energy;
ratiom=wenergym/energy;
ratiol=wenergyl/energy;
m=max*100000;
4 l8 n/ p# @3 F7 E! tws=wenergys*100000;
6 ~1 M. h8 g* c5 [# v9 Bwm=wenergym*100000;
% m! e# G* O$ h2 n) L7 Kwl=wenergyl*100000;
$ H8 R! F3 {! o1 q; Le=energy*100000;
a=[m,index,e,ws,wm,wl,ratios,ratiom,ratiol];

小土匪

求一維數組s的最大值和索引位置 max 和index可以求出來呀
也就是執行了for循環了
; `* [' p6 A" l" T/ T可能計算機計算速度快  我們感覺不到等待的時間   
如果是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還沒接觸過
  M, U, ?' G* A  J是不是取隨機數的時候  最大值的標號正好小于200了  這樣的話wenergyl=sum(s(index-large:index+large).^2)這步是不是會出錯？
9 J; ]7 [! p( `9 b) H
提示也正好是indices must either be real positive integers or logicals
  |- t+ l+ b* s8 [4 R7 l- O8 t; o你再修改下隨即數？

小土匪

Subscript indices must either be real positive integers or logicals
R$ P# O9 A3 s
中文解釋：下標索引必須是正整數類型或者邏輯類型
W7 b# p# ]3 I' k8 m
出錯原因：在訪問矩陣（包括向量、二維矩陣、多維數組，下同）的過程中，下標索引要么從 0 開始，要么出現了負數。注：matlab 的語法規定矩陣的索引從 1 開始

+ E! m' M' V7 x# T
這是網上搜到的資料   借鑒下