Re: [Faudiostream-devel] plotting frequency response

2013-09-15 Thread urban
Thanks very much, Next time I get into an argument whether we really need a 5 band EQ I can make my point with the support of frequency plots! I'll add it to my debug script, too. -Urban -- LIMITED TIME SALE - Full Year

Re: [Faudiostream-devel] plotting frequency response

2013-09-14 Thread Julius Smith
Hi Urban, You surmise correctly. Here is my implementation: function [xdb] = db(x,mindb) % [y] = db(x,mindb) % convert x to db, optionally clipped to mindb ax = abs(x); xmax = max(ax(:)); if nargin<2, xmin = xmax*eps; else xmin = 10^(mindb/20); end xdb = 20*log10(max(abs(x),xmin*ones(size(x

Re: [Faudiostream-devel] plotting frequency response

2013-09-14 Thread urban
Hi Julius, there is an undefined function called "db" (marked below). I guessed it could be implemented like: %function db(a,b) % min(20*log10(a),b) ; %endfunction but I dont know octave. Can you provide the function definition? Urban > I always use the "pfo.m" script below ("plot faustout"

[Faudiostream-devel] plotting frequency response

2013-09-12 Thread urban
Hi, I learned that I can view the frequency response of a filter like this: octave:1> plot(20*log10(abs(fft(faustout,1024))(1:512,:))); but how do I get the frequency axis to display Hz(log) ? Can someone provide me with a octave script? Thanks, Urban --

Re: [Faudiostream-devel] plotting frequency response

2013-09-12 Thread Julius Smith
Hi Urban, I always use the "pfo.m" script below ("plot faustout") - NOT pretty, but handy. - Julius % input parameters: if !exist('fs') fs = 1; freqstr = 'Normalized Frequency (cycles/sample)'; else if fs<1000 % (Hack) freqstr = 'Frequency (kHz)'; else freqstr = 'Frequency