最大配列サイズの基本設定を超えてしまう

Question

%%テストデータの読み込み
[x2,fs]=audioread('sound.wav');

d=zeros(length(x2),1); %d:エコー信号
x=zeros(length(x2),1) %x:参照信号
e=zeros(length(x2),1); %誤差信号
%wave信号を参照信号とエコー信号に分ける
x(:)=x2(:,1);
d(:)=x2(:,2);
myu=0.5; %ステップサイズパラメータ
filter_length=400;
w=zeros(filter_length,1);%フィルタ係数初期化
N=length(x2)-filter_length;

for i=1:N;
    xx=x(i:i+filter_length-1);
    y(i)=w.'.*xx;
    e(i)=d(i)-y(i);
    w=w+(myu.*((e(i).*x)./(x.'.*x)));
end
figure(1)
plot(flipud(w))

音源x2(2,80000)の中身のwavファイルでサンプリング周波数は8000Hzに設定しています.
学習同定法を用いてインパルス応答を推定するプログラムをリンク内容を参考にして作成しています.
x2の長さの80000までwを計算し,最後にwをプロットすることで音源のインパルス応答が取得できます.
for文のところを作成してみたのですが,
左辺と右辺の要素数が異なるため、代入を実行できません。

左辺と右辺の要素数が異なるため、代入を実行できません。

エラー:  (行 16)
    y(i)=w.'.*xx;
 
>> 

このようなエラーが出ます.

>> size(y)

ans =

       80000         400

>> size(w.'.*x)

ans =

       80000         400

>> 

と確認しましたが,sizeは合っているように見えます.
w,xxは400×1の行列です.

>> size(y(1))

ans =

     1     1

y=w.'.*x;
e=d-y;
w=w+(myu.*((e.*x)/(x.'.*x)));

またfor文をこのように直すと以下のようなエラーが出ます.

要求された 80000x80000 (47.7GB) 配列は、最大配列サイズの基本設定 (8.0GB) を超
えています。これにより、MATLAB は反応しなくなる可能性があります。

エラー:  (行 17)
w=w+(myu.*((e.*x)/(x.'.*x)));

Comments

[jbpb0]

質問のコードでは、yはforループ内の
y(i)=w.'.*xx;
でいきなり出てきてますけど、

> size(y)
ans = 80000 400

ということは、実際はそれよりも前に
y=zeros(length(x2),filter_length)
みたいにして、形状を定義してるのでしょうか？

> w,xxは400×1の行列です.

の場合、
size(w.'.*xx)
ans = 400 400
なので、
y(i)=w.'.*xx;
とするには、y(1), y(2)...それぞれが400x400のサイズで、それがN個あるのがyとなり、yは80000x400x400という、かなりサイズがでかい3次元配列になるはずで、

> size(y)
ans = 80000 400

とは矛盾します

y(i, :)=(w.*xx)';
ならば、式の右辺のサイズは
size((w.*xx)')
ans = 1 400
と1x400になるので、yはNx400の2次元配列になりますけど

[jbpb0]

あと
> w=w+(myu.*((e(i).*x)./(x.'.*x)));
の「(x.'.*x)」は、80000x80000というサイズの配列なので、サイズがでかすぎて領域確保できないと思う

そこを「(x.*x)」に変えたら80000x1のサイズになるけど、
size(myu.*((e(i).*x)./(x.*x)))
ans = 80000 1
と
size(w)
ans = 400 1
が違うので、結局
w=w+(...
の計算はできない

Answer 1

以下でどうですか？

for i=1:N;
xx=x(i:i+filter_length-1);
y(i)=w'xx;
e(i)=d(i)-y(i);
w=w+ (myu. (e(i)*xx)/(xx'*xx) );
end