這段程序出錯(cuò)，不清楚怎么改

mouseyuanhao

下面是我編的M文件的程序，主要是求一維數(shù)組s的最大值和索引位置，然后加不同寬度的窗，求出窗里信號(hào)能量后與總能量求比值。我編的程序中s是小波包分析的高頻信號(hào)，為了簡化程序段，我用5000個(gè)隨機(jī)數(shù)來代替。編譯的時(shí)候for循環(huán)好像根本就沒有執(zhí)行，不知道是為什么。請前輩們幫忙看看，不勝感激！
9 V& J) y- B/ r. c5 |/ ls=zeros(1,5000);
; X( Q+ t5 P7 Is=rand(1,5000);
6 _- ^# l1 V) w7 o0 g+ L3 y$ Ln=length(s);
+ {" w; I: F- F0 umax=0.0;
5 `0 y. o8 x* k% a; J( d5 J' qindex=0;
for i=1:n
if  max<s(i)
     max=s(i);
' d( i! n: U) K( E! H     index=i;
  else
, m# Q# j: |& |  Q1 G' e/ Q  end
4 L' L& _6 S7 Y# _$ |' ~2 Send   
- q6 B6 z9 n/ u- z2 O( _small=100;
wenergys=sum(s(index-small:index+small).^2);
middle=150;
) c1 r0 U7 @# ]wenergym=sum(s(index-middle:index+middle).^2);
large=200;
wenergyl=sum(s(index-large:index+large).^2);
$ G; l$ E0 J6 y2 W: i, t! xenergy=sum(s.^2);
- l- q+ O# s8 |8 {! Vratios=wenergys/energy;
% T! Y2 |4 G$ |+ j  z+ o5 O) Eratiom=wenergym/energy;
* r0 ?: B. N* L4 g. Y% S2 Mratiol=wenergyl/energy;
m=max*100000;
; l9 Z( a& n2 a' uws=wenergys*100000;
0 H: k" P: a9 Dwm=wenergym*100000;
' Q  T* [0 L; ^( ]wl=wenergyl*100000;
/ u1 u3 x! T2 d+ P" a/ ne=energy*100000;
a=[m,index,e,ws,wm,wl,ratios,ratiom,ratiol];

小土匪

求一維數(shù)組s的最大值和索引位置 max 和index可以求出來呀
# h: \" o& W6 W3 j* ~也就是執(zhí)行了for循環(huán)了
可能計(jì)算機(jī)計(jì)算速度快  我們感覺不到等待的時(shí)間   
/ T" r0 I- A' c8 m+ d如果是5000個(gè)數(shù)字迭代的話肯能時(shí)間很長

mouseyuanhao

其實(shí)我的這個(gè)程序是插在labview里的腳本。我一開始用的就是[max,index]=max(s);但是我是把程序作為腳本插入在labview程序里，labview提示出錯(cuò)：subscript indices must either be real positive integers or logicals. 然后就不給我執(zhí)行了，改成for循環(huán)也一樣。5000個(gè)數(shù)字迭代時(shí)間并不長，我在MATLAB里面試過，瞬間出結(jié)果，但是在labview里編譯就是通不過。我現(xiàn)在懷疑是不是labview和MATLAB的兼容做的不夠好。

小土匪

labview還沒接觸過
4 ?8 F8 P% K- d* B是不是取隨機(jī)數(shù)的時(shí)候  最大值的標(biāo)號(hào)正好小于200了  這樣的話wenergyl=sum(s(index-large:index+large).^2)這步是不是會(huì)出錯(cuò)？
$ S6 N1 f: A' l% }5 I% x, n0 s
提示也正好是indices must either be real positive integers or logicals
你再修改下隨即數(shù)？

小土匪

Subscript indices must either be real positive integers or logicals
R$ P# O9 A3 s
中文解釋：下標(biāo)索引必須是正整數(shù)類型或者邏輯類型
W7 b# p# ]3 I' k8 m
1 g8 j: o8 l2 h) |- j6 A5 f. t: f出錯(cuò)原因：在訪問矩陣（包括向量、二維矩陣、多維數(shù)組，下同）的過程中，下標(biāo)索引要么從 0 開始，要么出現(xiàn)了負(fù)數(shù)。注：matlab 的語法規(guī)定矩陣的索引從 1 開始


這是網(wǎng)上搜到的資料   借鑒下