Hi,
This code seems to work for making long PFBs with the green blocks. It
could probably be improved, but at least it's functional. Basically it
assigns the 'buf' variable that contains the PFB coefficients to a
unique name in the base Matlab workspace. I couldn't manage to get it
to work assigning it in any other namespace. The unique name I chose
is just related to the current time, so is meaningless. A more
meaningful unique name could probably be constructed, but this was
easy. Let me know if anyone has a better fix for this problem.
Glenn
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% Center for Astronomy Signal Processing and Electronics Research %
% http://seti.ssl.berkeley.edu/casper/ %
% Copyright (C) 2007 Terry Filiba, Aaron Parsons %
% %
% This program is free software; you can redistribute it and/or modify %
% it under the terms of the GNU General Public License as published by %
% the Free Software Foundation; either version 2 of the License, or %
% (at your option) any later version. %
% %
% This program is distributed in the hope that it will be useful, %
% but WITHOUT ANY WARRANTY; without even the implied warranty of %
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the %
% GNU General Public License for more details. %
% %
% You should have received a copy of the GNU General Public License along %
% with this program; if not, write to the Free Software Foundation, Inc., %
% 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function pfb_coeff_gen_init(blk, varargin)
% Initialize and configure the Polyphase Filter Bank coefficient generator.
%
% pfb_coeff_gen_init(blk, varargin)
%
% blk = The block to configure.
% varargin = {'varname', 'value', ...} pairs
%
% Valid varnames for this block are:
% PFBSize = Size of the FFT (2^FFTSize points).
% CoeffBitWidth = Bit width of coefficients.
% TotalTaps = Total number of taps in the PFB
% CoeffDistMem = Implement coefficients in distributed memory
% WindowType = The type of windowing function to use.
% bram_latency = The latency of BRAM in the system.
% n_inputs = Number of parallel input streams
% nput = Which input this is (of the n_inputs parallel).
% fwidth = The scaling of the bin width (1 is normal).
% Declare any default values for arguments you might like.
defaults = {};
if same_state(blk, 'defaults', defaults, varargin{:}), return, end
check_mask_type(blk, 'pfb_coeff_gen');
munge_block(blk, varargin{:});
disp('Running pfb coeff init8')
PFBSize = get_var('PFBSize', 'defaults', defaults, varargin{:});
CoeffBitWidth = get_var('CoeffBitWidth', 'defaults', defaults, varargin{:});
TotalTaps = get_var('TotalTaps', 'defaults', defaults, varargin{:});
CoeffDistMem = get_var('CoeffDistMem', 'defaults', defaults, varargin{:});
WindowType = get_var('WindowType', 'defaults', defaults, varargin{:});
bram_latency = get_var('bram_latency', 'defaults', defaults, varargin{:});
n_inputs = get_var('n_inputs', 'defaults', defaults, varargin{:});
nput = get_var('nput', 'defaults', defaults, varargin{:});
fwidth = get_var('fwidth', 'defaults', defaults, varargin{:});
% Set coefficient vector
alltaps = TotalTaps*2^PFBSize;
windowval = transpose(window(WindowType, alltaps));
total_coeffs = windowval .*
sinc(fwidth*([0:alltaps-1]/(2^PFBSize)-TotalTaps/2));
buf = zeros(1,alltaps/2^n_inputs);
for i=1:alltaps/2^n_inputs,
buf(i)=total_coeffs((i-1)*2^n_inputs + nput + 1);
end
%%
this_buf_name = ['pfbfirbuf' num2str(round(rem(now,1)*1e8))];
assignin('base',this_buf_name,buf)
%%
delete_lines(blk);
% Add Ports
reuse_block(blk, 'din', 'built-in/inport', 'Position', [235 28 265 42], 'Port',
'1');
reuse_block(blk, 'sync', 'built-in/inport', 'Position', [15 93 45 107], 'Port',
'2');
reuse_block(blk, 'dout', 'built-in/outport', 'Position', [360 28 390 42],
'Port', '1');
reuse_block(blk, 'sync_out', 'built-in/outport', 'Position', [130 28 160 42],
'Port', '2');
reuse_block(blk, 'coeff', 'built-in/outport', 'Position', [450 343 480 357],
'Port', '3');
% Add Static Blocks
reuse_block(blk, 'Delay', 'xbsIndex_r3/Delay', ...
'latency', 'bram_latency+1', 'Position', [65 12 110 58]);
reuse_block(blk, 'Counter', 'xbsIndex_r3/Counter', ...
'cnt_type', 'Free Running', 'n_bits', num2str(PFBSize-n_inputs),
'arith_type', 'Unsigned', ...
'rst', 'on', 'explicit_period', 'on', 'Position', [65 75 115 125]);
reuse_block(blk, 'Delay1', 'xbsIndex_r3/Delay', ...
'latency', num2str(bram_latency+1), 'Position', [290 12 335 58]);
reuse_block(blk, 'Concat', 'xbsIndex_r3/Concat', ...
'num_inputs', num2str(TotalTaps), 'Position', [310 97 365 643]);
reuse_block(blk, 'Register', 'xbsIndex_r3/Register', ...
'Position', [385 325 430 375]);
add_line(blk, 'din/1', 'Delay1/1');
add_line(blk, 'Delay1/1', 'dout/1');
add_line(blk, 'sync/1', 'Counter/1');
add_line(blk, 'sync/1', 'Delay/1');
add_line(blk, 'Delay/1', 'sync_out/1');
add_line(blk, 'Concat/1', 'Register/1');
add_line(blk, 'Register/1', 'coeff/1');
% Add Dynamic Blocks
for a=1:TotalTaps,
blkname = ['ROM', num2str(a)];
v = [this_buf_name, '(', num2str((a-1)*2^(PFBSize-n_inputs) + 1), ':',
num2str(a*2^(PFBSize-n_inputs)), ')'];
%v = ['buf(', num2str(a-1), '*2^(PFBSize-n_inputs) + 1:', num2str(a),
'*2^(PFBSize-n_inputs))'];
%v = mat2str(buf((a-1)*2^(PFBSize-n_inputs)+1 : a*2^(PFBSize-n_inputs)));
reuse_block(blk, blkname, 'xbsIndex_r3/ROM', ...
'depth', num2str(2^(PFBSize-n_inputs)), 'initVector', v, 'arith_type',
'Signed (2''s comp)', ...
'n_bits', num2str(CoeffBitWidth), 'bin_pt', num2str(CoeffBitWidth-1),
...
'latency', num2str(bram_latency), 'use_rpm','on', 'Position', [150
65*(a-1)+74 200 65*(a-1)+126]);
add_line(blk, 'Counter/1', [blkname, '/1']);
reintname = ['Reinterpret', num2str(a)];
reuse_block(blk, reintname, 'xbsIndex_r3/Reinterpret', ...
'Position', [220 65*(a-1)+84 260 65*(a-1)+116]);
add_line(blk, [blkname, '/1'], [reintname, '/1']);
add_line(blk, [reintname, '/1'], ['Concat/', num2str(a)]);
end
% Set coefficient ROMs to use distribute memory (or not).
for a=1:TotalTaps,
blkname = ['ROM', num2str(a)];
if CoeffDistMem,
set_param([blk,'/',blkname], 'distributed_mem', 'on');
else,
set_param([blk,'/',blkname], 'distributed_mem', 'off');
end
end
clean_blocks(blk);
fmtstr = sprintf('PFBSize=%d, n_inputs=%d, taps=%d', PFBSize, n_inputs,
TotalTaps);
set_param(blk, 'AttributesFormatString', fmtstr);
save_state(blk, 'defaults', defaults, varargin{:});
