2

sph2cartによって返されるデータをベクトル要素を持つ行列として整理し、この行列の各要素を操作したいと考えています(ベクトル-ベクトルまたはベクトル-スカラー計算)。これを達成する例を次に示します。

lightV = zeros(1, 1, 3);
lightV(1,1,1) = 0.5;
lightV(1,1,2) = 0.4;
lightV(1,1,3) = 0.7;
[Az El] = meshgrid(0:60:360, 0:15:90);
[x y z] = sph2cart(Az*pi/180, El*pi/180, 1);
refV = zeros(size(Az,1), size(Az,2), 3);
radius = zeros(size(Az,1), size(Az,2));
for i = 1:size(Az,1)
    for j = 1:size(Az,2)
       refV(i,j,1) = -x(i,j);
       refV(i,j,2) = -y(i,j);
       refV(i,j,3) = z(i,j);
       radius(i,j) = dot(refV(i,j,:), lightV(1,1,:));
    end
 end

ただし、これはやや冗長に見えますが、どうすればもっと簡潔にできますか?

4

1 に答える 1

0

コードを次のようにベクトル化できます。

S.lightV = [0.5, 0.4, 0.7];
[Az, El] = meshgrid(0:60:360, 0:15:90);
[S.x, S.y, S.z] = sph2cart(Az * pi/180, El * pi/180, 1);
S.refV = cat(3, -x, -y, z);
S.radius = sum(bsxfun(@times, S.refV, reshape(S.lightV, 1, 1, [])), 3);

cat3 番目の次元に沿って連結するための の使用と、ネストされた for ループ内の内積を置き換えるためのbsxfunとの組み合わせに注意してください。sumまた、すべてを 1 つの struct にバインドしましたS。

于 2013-05-30T09:01:54.267 に答える