プログラムの処理速度を早くする

Question

[x fs]=audioread('-------.wav');
n_length=length(x);
[g fs]=audioread('-------.wav');


g_length=length(g);
y=zeros(1,n_length+g_length-1);
if n_length>g_length
for n = 1:n_length+g_length-1
    if n > n_length
        left = n - n_length + 1;
    else
        left = 1;
    end
    if n < g_length
        right = n;
    else
        right = g_length;
    end
    for k = left:right
        y(n) = y(n) + g(k) * x(n-k+1);
    end
    
end
if n_length<g_length
        for n = 1:n_length+g_length-1
    if n < n_length
         left = n;
    else
        left = g_length;
    end
       
    if n < g_length
         right = n - n_length + 1;
    else
        right = 1;
    end
    for k = right:left
        y(n) = y(n) + g(k) * x(n-k+1);
    end
        end   
end
end

初心者です.
convを使わずに畳み込み積分をしたいです.

前回の質問でこのような回答をいただいたのですが,このプログラムはn_length>g_lengthの時しかうまく動作できません.if文で最初にn_lengthとg_lengthの大小関係を比べてn_length>g_lengthの時と分岐させることもできるのですがすごく複雑になってしまい,処理にすごく時間がかかってしまいます.短く同じ処理をできるプログラムに修正できないでしょうか.初心者なので詳しく説明していだだけると幸いです.

Comments

[K_3578]

その、前の質問で回答された方に聞かれればよろしいのでは？この質問の必要性が分かりません

Answer 1

for i = 1:n_length
  for j = 1:g_length
    f(i + j - 1) = f(i + j - 1) + x(i) * g(j);
  end
end

畳み込みがしたいだけならこれだけです。(MATLABを触ったことがないので正しいコードではないかもしれない)
速度が遅いのはコードが複雑だからではなく単純にこのアルゴリズムが遅いからです。
高速なアルゴリズムとしてはFFTを使ったものが有名でググればいくらでもやり方が出てきますが、それなりに複雑です。

Comments

[kaeruuuun]

できました.ご説明ありがとうございます.