1% Render a triangle roll with default parameters and print diagnostics.
2% Set P.strike_count = 1 for a single strike.
3% Run with: numbl run triangle_demo.m
5P = triangle_defaults();
6tic;
7[y, basis, dinfo] = triangle_render(P);
8elapsed = toc;
10fprintf('modes kept : %d (%d in-plane, %d out-of-plane)\n', basis.n, dinfo.modesIn, dinfo.modesOut);
11fprintf('mode freq : ');
12fprintf('%.0f ', basis.f);
13fprintf('Hz\n');
14fprintf('strikes : %d at %.1f Hz (mean v %.2f m/s, mean contact %.0f us, peak force %.0f N)\n', ...
15 dinfo.strikeCount, P.strike_rate_hz, dinfo.meanV, dinfo.meanTauUs, dinfo.peakForceN);
16fprintf('peak press.: %.4f Pa (%.1f dB SPL) at %.0f cm\n', ...
17 dinfo.peakPa, dinfo.peakSplDb, dinfo.listenerDistM * 100);
18fprintf('length : %.2f s @ %d Hz\n', length(y) / P.sample_rate, P.sample_rate);
19fprintf('render time: %.2f s\n', elapsed);
20if dinfo.listenerInside
21 fprintf('WARNING : listener is at or inside the rod\n');
22end