bugfix modal strength

9f36ac5a · Marco Berzborn · 0e753a93 · 9f36ac5a
Commit 9f36ac5a authored Dec 12, 2016 by Marco Berzborn
Hide whitespace changes
Inline Side-by-side

Showing with 1 addition and 1 deletion

applications/SphericalHarmonics/ita_sph_modal_strength.m applications/SphericalHarmonics/ita_sph_modal_strength.m +1 -1

No files found.
--- a/applications/SphericalHarmonics/ita_sph_modal_strength.m
+++ b/applications/SphericalHarmonics/ita_sph_modal_strength.m
@@ -90,7 +90,7 @@ for idxRad = 1:numel(uniqueRad)
    switch sArgs.transducer
        case {'microphone','mic'}
            if strcmp(type,'rigid') && isempty(sArgs.scatterer)
-                bn = bn*diag(4*pi./(kVec.*uniqueRad(idxRad)));
+                bn = bn.*(4*pi*(1i.^(n-1))*(1./(kVec.*uniqueRad(idxRad)).^2));
            else
                bn = diag(4*pi*1i.^n)*bn;
            end