添字が不正というエラー

Question

len=512;
order=200;
sfreq=16;
fshift=10;
fwidth=len/sfreq;
[x,y]=loadwave('/Users/inoguchikaede/Sound_data/mhosda01.wav'); 
n_sample=y(8);
disp(n_sample); 
for i=1:n_sample
    ts(i)=i/sfreq;
end
for i=1:len
    win(i)=0.54-0.46*cos(2*%pi*i/len);
end
n_frame=0;
for st=1:sfreq*fshift:n_sample-len
    pwr=0;
    for i=1:len
        x1(i)=x(i+st)*win(i);
        pwr=pwr+x1(i)*x1(i);
    end
    n_frame=n_frame+1;
    tp(n_frame)=fwidth/2+(n_frame-1)*fshift;
    logpwr(n_frame)=10*log10(pwr);
end
ここからエラーが出る--------------------------
thlogpwr=-25;
for n=1:n_frame // n_frameは分析フレーム数
    if logpwr(n) < thlogpwr then
        f0(n)=0;
    else
        st=1+(n-1)*sfreq*fshift; // 現フレームの開始サンプル番号
        j=1;
        for l=0:order
            
           r(j)=0;
           for h=1:len-l
               r(j)=r(j)+x1(h)*x1(h+l);
           end
           j=j+1;
        //
            //
        end
        max=0;
        for j=80:200
           
            argmax=max(r);
            
        end
        period = argmax/16;
        f0(n) = 1000/period;
    end
end

ケプストラム法で音声を解析するのに自己相関関数を求めてピークを検出し,基本周波数を求めるプログラムを書いています.
長くなってしまっているのですが,ここからエラーが出ると書いたところが自己相関関数を定義してそこからピークを検出するプログラムです.
ここから下の部分において

添字が不正です.

というエラーが出ます.具体的にどこが間違っているか分かりません.
エラーを省略しないでくださいとのコメントをいただいたのですが,特にどこが間違っているかエラーには書かれていませんでした.
個人的にmaxのところが違うのではないかと考えてr(j)にしてみたりしたのですが,同じエラーが出ました.
非常に初歩的な質問になるのですが,今のmax関数の定義で言うとr(80)からr(200)の最大値がargmaxに入るという考えてあっていますか?
よろしくお願いいたします.

追記
この質問ですが質問の意図を伝え忘れていました.
max=0がおかしいと言うことですが,最初にmaxに初期値の値を入れておくことはできませんか?
あと今のプログラムだとargmaxに最大値の値が入っていると思うのですが,私が入れたいのは,最大値の値ではなく,そのときのr(j)のjのでした.つまり,r(80)からr(200)の中で最大値の位置です.
自分の理解不足で違う意図の質問をしてしまっていました.申し訳ございません.

max関数は第二引数で最大値のインデックスを返すと言うことなので,プログラムを修正してみたのですが,max=0でエラーが出ます.

         max=0;
       for j=80:200
         
            [z,argmax]=max(r(80:200));
            
        end

エラー
実行されたファイルの   103 行目 
/Users/

Wrong number of output arguments.

私の問題と意図としてはr(j)のj=80から200の最大値の位置をargmaxに入れたいです.

Comments

[yuki23]

エラーメッセージは勝手に省略せず、出力されたメッセージをそのまま質問に記載してください。
「この下のどこか」では何もわかりません

[kaeruuuun]

コメントありがとうございます.
ただ使っているのがscilabだからなのか分かりませんが,エラーメッセージからは場所が特定できなかったからこの部分だけを省略したのです.
ここからエラーが出ると書いた部分より上の部分はそれだけでエラーが出なかったので,それより下だということは確認しています.

--> exec('/Users/', -1)

61648.
組込み exec( /Users/ 行 -1 )

添字が不正です.

[yuki23]

> 今のmax関数の定義で言うとr(80)からr(200)の最大値がargmaxに入るという考えてあっていますか?
自分で書いたコードではないのでしょうか？

[kaeruuuun]

そうです.単純なことをしているはずなのにただ,この計算がうまく行ってないような気がしています.

[kaeruuuun]

max(r(j))とmax(r)は同じですよね？これを治しても同じエラーが出るので

[kaeruuuun]

rまでは値が入っているのですがargmaxには値がどうしても入りません.

[jbpb0]

質問内容とは違いますが、
n_sample=y(8);
って、合ってますか？

https://help.scilab.org/docs/6.1.0/en_US/loadwave.html
を見ると、8番目は「bytes per sample (per channel)」ってなってますが、n_sampleに入れる数値の単位は「バイト」でいいのですか？

Scilab使ってないので、外してたらごめんなさい

【追記】 やっぱり外してるみたいなので、スルーしてください

[kaeruuuun]

>jbpb0さん
ありがとうございます.
承知いたしました.
実はこのプログラムはここからエラーが出ると書いたところからmax=0の上までで自己相関係数を求めています.ただ私が実現したいのは自己相関係数ではなく自己相関関数です.自己相関関数は自己相関係数を逆離散フーリエ変換したら出てくると言う認識だったので,r(j)の式の下にy(j)=ifft(r(j))としたのですが,上手く実現できませんでした.もしわかるようでしたら教えていただきたいです.

[jbpb0]

> 自己相関関数は自己相関係数を逆離散フーリエ変換したら出てくる

違いますよ
自己相関関数は、タイムラグ(遅延)を変えながら自己相関係数を計算するようなものです
(いろんなタイムラグに対して、それぞれ値があるので、関数と呼ばれる)
https://oceanone.hatenablog.com/entry/2020/04/13/235830

ただし、正規化するかしないかとか、細部には差異はありますが
https://tmytokai.github.io/open-ed/activity/acr/text01/page01.html

[kaeruuuun]

ありがとうございます.
>自己相関関数は、タイムラグ(遅延)を変えながら自己相関係数を計算するようなものです
と言うことは自己相関関数を用いて実現はできると言うことですよね？

[kaeruuuun]

自己相関係数が自己相関関数の離散的なものというイメージだったので逆離散フーリエ変換を取れば良いのではないかと考えていました.

[jbpb0]

パワースペクトルをフーリエ逆変換すると自己相関関数になります
https://oresoken.blog.ss-blog.jp/2015-03-17
そのことと混同してませんか？

あと、Scilabではxcorr()で自己相関関数が計算できるのでは？
https://help.scilab.org/docs/6.1.0/ja_JP/xcorr.html

[kaeruuuun]

混同していたかもしれません.
>Scilabではxcorr()で自己相関関数が計算できるのでは？
r(j)=xcorr(x1(h));
このように定義し直してみたのですが,計算量が多いのはあるのですがかれこれ10分以上クルクルしている状態で,間違っているような気がしています.

[jbpb0]

この質問のもともとの内容とかなりずれてきてるので、別の質問にしてください

【追記】 https://help.scilab.org/docs/6.1.0/ja_JP/xcorr.html
の「例」の計算式
[c, ind] = xcorr(y, "biased");
と、

> r(j)=xcorr(x1(h));

を見比べてください
全然違いますよね
「例」の計算式には、(j), (h)みたいなもの付いてないですよね

一度、「例」の計算をそのまま自分でもやってみて、その元データ「y」の形や、計算結果「c」と「ind」の形を確認して、どのように計算させたらいいのかを理解してください

[kaeruuuun]

承知しました.
別の質問にします.

[kaeruuuun]

新しく質問にしました.良ければ見ていただきたいです.

[jbpb0]

> 違う意図の質問をしてしまっていました

ということなら、別の質問にした方がいいと思います
「解決済」だと、あんまり見てもらえませんし、質問文の最後の追記まで読まないと質問の意図が分からないとなると、なおさら読まれません

[kaeruuuun]

わかりました.別の質問にします.教えていただきありがとうございます.

Answer 1

        max=0;
        for j=80:200

            argmax=max(r);

        end

ここが間違っています。maxに0を代入しているので関数ではなくなってしまっています。

そもそもベクトルの最大値を求めるのにfor文を使う必要はありません。

argmax = max(r(80:200));

Comments

[kaeruuuun]

ありがとうございます.
実行はできました.