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

Question

このような式をmatlabで書いていて

f=1000;%対象とする周波数
c=340;%音速
dist=0.04;%評価音場の観測点間隔 
rs=1.5; %スピーカアレーの半径
k=2*pi*f/c; %波数
L=64; %スピーカの数
dr=zeros(1,L); %L個のスピーカ駆動信号
phi_s=(0:L-1)/L*2*pi; %スピーカ角度
N=fix((L-1)/2); %展開次数の最大値
phi_t=pi/3; %平面波の到来方向
for s=1:L
    for n=-N:N
        dr(s)=dr(s)+((1i^(n+1)*exp(-1i*n*phi_t))/(besselh(n,2,k*rs))*exp(1i*n*phi_s));
    end
    dr(s)=dr(s)*2/pi;
end

のように書いてみたのですが,

左辺と右辺の要素数が異なるため、代入を実行できません。
dr(s)=dr(s)+((1i*exp(-1i*n*phi_t))/(besselh(n,2,k*rs))*exp(1i*n*phi_s));

のようなエラーが出て困っています.

Answer 1

「for s=1:L」以降の代わりに

s=1;
n=-N;
size(dr(s))
size(((1i^(n+1)*exp(-1i*n*phi_t))/(besselh(n,2,k*rs))*exp(1i*n*phi_s)))

を実行して、サイズを確認してください
二つの「size()」の結果が一致しないとダメです

Comments

[kaeruuuun]

ご回答ありがとうございます.
s=1;
n=-N;
size(dr(s))
size(((1i^(n+1)*exp(-1i*n*phi_t))/(besselh(n,2,k*rs))*exp(1i*n*phi_s)))

ans =

1 1

ans =

1 64

>>
for s=1:L以降ということなのでこのようにしてみましたが,これだとsが1のままだと思います.

[jbpb0]

「for s=1:L」は、
s=1で計算
s=2で計算
...
s=Lで計算
をするという意味です
だから、ループの1回目はs=1で計算されます

[kaeruuuun]

そうでした.そもそもこの1回目を計算するところからでした.
私が求めたいのは1×64の行列なので,dr(s)が1行1列になってしまっているということですね

[kaeruuuun]

左辺を1×64の行列にできますか?

[jbpb0]

size(dr)

[kaeruuuun]

drは1×64の行列です.という結果が表示されました.

[jbpb0]

> drは1×64の行列

size(((1i^(n+1)*exp(-1i*n*phi_t))/(besselh(n,2,k*rs))*exp(1i*n*phi_s)))
と同じですよね

[kaeruuuun]

そうです.

[kaeruuuun]

なぜdr(s)は1×1なのかがわかりません.

[jbpb0]

s=1のとき、dr(s)=dr(1)

[kaeruuuun]

であれば,s=1の時は右辺も1×1にならないといけないのですか?

[kaeruuuun]

このままだと行列のサイズが違うので計算できないですよね？

[kaeruuuun]

dr=zeros(1,L);と定義しているのでdrは1×64の行列になると思います.

[jbpb0]

> dr=zeros(1,L);と定義しているのでdrは1×64の行列になる

そうです

なので、
%for s=1:L
for n=-N:N
dr=dr+((1i^(n+1)*exp(-1i*n*phi_t))/(besselh(n,2,k*rs))*exp(1i*n*phi_s));
end
dr=dr*2/pi;
%end
でいいのでは？？
いちいちdrを分解しなくて、64個を一括で計算すれば

[kaeruuuun]

>そうです
なので、
%for s=1:L
for n=-N:N
dr=dr+((1i^(n+1)*exp(-1i*n*phi_t))/(besselh(n,2,k*rs))*exp(1i*n*phi_s));
end
dr=dr*2/pi;
%end
でいいのでは？？
いちいちdrを分解しなくて

それはそうなのですが,このようなスタイルで与えられていたのでsを使ってできないかと考えていましたがsを使うと無理なのでしょうか.

[jbpb0]

size(((1i^(n+1)*exp(-1i*n*phi_t))/(besselh(n,2,k*rs))*exp(1i*n*phi_s(s))))

[kaeruuuun]

やってみます

[jbpb0]

ただし、matlabでは、なるべくfor...使わないでベクトルや配列を一括計算する方がいいです

参考
https://jp.mathworks.com/help/matlab/matlab_prog/vectorization.html

[kaeruuuun]

ans =

1 1
のような結果でした.

[kaeruuuun]

なのでこれだとdr(s)に代入できます.

[kaeruuuun]

>ただし、matlabでは、なるべくfor...使わないでベクトルや配列を一括計算する方がいいです
計算量の問題ですか?
>size(((1i^(n+1)*exp(-1i*n*phi_t))/(besselh(n,2,k*rs))*exp(1i*n*phi_s(s))))
なぜこのphi_s(s)にしたら1行1列になったのか頭が追いついていません.

[jbpb0]

> 計算量の問題ですか?

上記参考Webページ見てください

> なぜこのphi_s(s)にしたら1行1列になったのか

size(phi_s)
size(phi_s(s))

[kaeruuuun]

一つ一つsizeを用いてその中身が何行何列か確認したら,phi_sでエラーが出ていることがわかりました.
ありがとうございます.