wolffd@0: % filter = SetGain(filter, desired, f, fs)
wolffd@0: % Set the gain of a filter (1x5 vector) to any desired gain
wolffd@0: % at any desired frequency (f).
wolffd@0: 
wolffd@0: % (c) 1998 Interval Research Corporation  
wolffd@0: 
wolffd@0: function filter = SetGain(filter, desired, f, fs)
wolffd@0: oldGain = 10^(FreqResp(filter, f, fs)/20);
wolffd@0: filter(1:3) = filter(1:3)*desired/oldGain;