Next  |  Prev  |  Up  |  Top  |  Index  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search

Matlab for Passive Reflectance Synthesis Method 2

... as in Method 1 for poles ...

% Construct a resonator as a sum of arbitrary modes 
% with unit residues, adding a near-zero at dc.

B = zeros(1,2*nsec+1);
impulse = [1,zeros(1,2*nsec)];
for i=1:nsec, 
    % filter() does polynomial division here:
    B = B + filter(A,denoms(i,:),impulse);
end;

% add a near-zero at dc 
B = filter([1 -0.995],1,B);

... as in Method 1 for display ...


Next  |  Prev  |  Up  |  Top  |  Index  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search

[How to cite this work]  [Order a printed hardcopy]  [Comment on this page via email]

``Physical Audio Signal Processing'', by Julius O. Smith III, W3K Publishing, 2010, ISBN 978-0-9745607-2-4
Copyright © 2023-08-20 by Julius O. Smith III
Center for Computer Research in Music and Acoustics (CCRMA),   Stanford University
CCRMA