www.gusucode.com > signal 工具箱matlab源码程序 > signal/@dfilt/@parallel/dgdfgen.m
function msg = dgdfgen(Hd,hTar,doMapCoeffsToPorts,pos)
%DGDFGEN
% Author(s): V. Pellissier
% Copyright 2004-2006 The MathWorks, Inc.
narginchk(4,4);
sys = hTar.system;
t=regexp(sys, '/'); % to get the system and the block
if length(t)==1 % eg: sys =Filter/stage
% Find out if a block with the same name already exist
h= find_system(sys(1:t-1),'LookUnderMasks','all','Name', sys(t+1:end));
% checks to see whether the FILTER blk needs to be redrawn
[redraw,nstages_equal,h] = chkIfRedrawNecessary(Hd,h,sys,'dfilt.parallel');
elseif length(t)>1 % eg: sys =Filter/stage1/stage1
% if n Stages are present inside FILTER BLOCK
h= find_system(sys(1:t(end)-1),'Searchdepth',1,'LookUnderMasks','all','Name', sys(t(end)+1:end));
% checks to see whether the FILTER blk needs to be redrawn
[redraw,nstages_equal,h] = chkIfRedrawNecessary(Hd,h,sys,'dfilt.parallel');
end
p = positions;
if isempty(h) || ... % If the destination is not found add the block
redraw || ... % If the above algorithm says to redraw, add a new block
strcmp(hTar.OverwriteBlock, 'off') % When we are not overwriting blocks
% we have a unique subsystem name
%Add SubSystem
hsubsys = add_block('built-in/subsystem',hTar.system, 'Tag', 'FilterWizardSubSystem');
% Restore position of the block
set_param(hsubsys,'Position', pos);
% Determine horizontal and vertical offset for input:
[x_offset, y_offset] = lcloffset(1);
offset = [x_offset, y_offset, x_offset, y_offset];
p = positions;
% Input:
blk = 'Input';
[hndl,hndl1]= hTar.inport(blk, Hd);
pos = p.input + offset;
set_param([sys '/' blk], 'Position',pos);
end
blk='Input';
last_conn{1,1} = [blk '/1'];
last_conn{1,2} = '';
set_param(sys, 'UserDataPersistent', 'on');
% store filter and mapcoeffstoports information
CurrentFilter.filter = class(Hd);
CurrentFilter.nstages = nstages(Hd);
CurrentFilter.mapcoeffstoports = hTar.MapCoeffsToPorts;
set_param(sys,'UserData',CurrentFilter);
% Sections
nsections = length(Hd.Stage);
for k=1:nsections,
% Force sys to be the current system
idx = findstr(sys, '/');
set_param(0,'CurrentSystem',sys(1:idx(end)-1));
% Add new subsystem for filter realization:
section_name = ['Stage', sprintf('%d',k)];
hTarSection = copy(hTar);
hTarSection.destination = 'current';
% Get section coefficient names and variables
if doMapCoeffsToPorts
hTarSection.CoeffNames = hTar.CoeffNames.(sprintf('%s%d','Stage',k));
coeffvar = hTar.privCoefficients.(sprintf('%s%d','Stage',k));
setprivcoefficients(hTarSection,coeffvar);
end
idx = findstr(sys, '/');
if length(idx)==1,
blockpath = hTar.blockname;
else
blockpath = sys(idx(end)+1:end);
end
hTarSection.blockname = [blockpath '/' section_name];
pos = createmodel(hTarSection);
subsys = hTarSection.system;
% Obtain state information for each stage
hTarSection = parse_filterstates(Hd.Stage(k),hTarSection);
% Realize this section
msg='';
if isa(Hd.Stage(k),'dfilt.multistage') || ...
isa(Hd.Stage(k),'dfilt.coupledallpass') || ...
isa(Hd.Stage(k),'mfilt.abstractiirmultirate'),
msg = dgdfgen(Hd.Stage(k),hTarSection,doMapCoeffsToPorts,pos);
% Add an extra input port for Farrow filters
nports = privnports(Hd.Stage(k));
joffset = length(find_system(sys,'SearchDepth',1,'LookUnderMasks','all', ...
'Regexp','on','Name','^FracDelay'));
if nports>1,
for j=1:nports-1,
blk = ['FracDelay',num2str(joffset+j)];
[hndl,hndl1]= hTar.inport(blk, Hd.Stage(k));
pos = p.input + offset+(joffset+j)*[0 30 0 30];
set_param([hTar.system '/' blk], 'Position',pos);
interstg(hTar, {[blk '/1']},{[hTarSection.blockname '/' num2str(j+1)]}, 12);
end
end
else
DGDF = dgdfgen(Hd.Stage(k),hTarSection,doMapCoeffsToPorts);
DG = expandToDG(DGDF,doMapCoeffsToPorts);
check_if_optimizezeros_possible(Hd.Stage(k), hTarSection);
% Optimisations for this section
if isfield(get(Hd.Stage(k)), 'Arithmetic'),
optimize(DG,...
strcmpi(hTarSection.optimizeones,'on'),...
strcmpi(hTarSection.optimizenegones,'on'),...
strcmpi(hTarSection.optimizezeros,'on'),...
strcmpi(hTarSection.optimizedelaychains,'on'),...
strcmpi(hTarSection.mapcoeffstoports,'on'),...
Hd.Stage(k).Arithmetic);
else
optimize(DG,...
strcmpi(hTarSection.optimizeones,'on'),...
strcmpi(hTarSection.optimizenegones,'on'),...
strcmpi(hTarSection.optimizezeros,'on'),...
strcmpi(hTarSection.optimizedelaychains,'on'), ...
strcmpi(hTarSection.mapcoeffstoports,'on'),...
'double');
end
% Garbage Collection (clean up)
DG = gc(DG);
% generate mdl model
dg2mdl(DG,hTarSection,pos);
% Add an extra input port for Farrow filters
if isa(Hd.Stage(k),'dfilt.abstractfarrowfd'),
koffset = length(find_system(sys,'SearchDepth',1,'LookUnderMasks','all', ...
'Regexp','on','Name','^FracDelay'));
blk = ['FracDelay',num2str(koffset+1)];
[hndl,hndl1]= hTar.inport(blk, Hd.Stage(k));
pos = p.input + (koffset+1)*(offset+[0 30 0 30]);
set_param([hTar.system '/' blk], 'Position',pos);
interstg(hTar, {[blk '/1']},{[hTarSection.blockname '/2']}, 12);
end
end
% Determine horizontal and vertical offset for rendering filter stage:
[x_offset, y_offset] = lcloffset(k);
offset = [x_offset, y_offset, x_offset, y_offset];
pos = p.section + offset;
set_param(subsys,'Position', pos);
if ~nstages_equal
new_conn{1,1} = [section_name '/1'];
% Last stage: no summer - use output of Section block for inter-stage connections
new_conn{1,2} = new_conn{1,1};
new_conn{2,2}='';
if k<nsections,
% Use a SUMMER :
if k==1,
sum_str = '|++';
orient = 'right';
else
sum_str = '++|';
orient = 'up';
end
blk = ['Sum' sprintf('%d', k)];
sum(hTar, blk, sum_str, 'double');
pos = p.sum + offset;
set_param([sys '/' blk], 'Position', pos, ...
'Orientation', orient, 'ShowName','off');
% %ic{k} = blk;
% end
ic{k} = blk;
% Internal connection:
add_line(sys, [section_name '/1'], [ic{k} '/1'],'autorouting','on'); % Connect Section subsystem to summer
new_conn{1,2}=[blk '/1']; % output of the summer
new_conn{2,2}=[blk '/2']; % 2nd input of the summer
elseif k==1,
ic{1} = section_name;
end
% Connect input to this stage:
interstg(hTar, last_conn,new_conn, [12 21]);
last_conn{1,2} = new_conn{2,2};
end
end
if ~nstages_equal
% Determine horizontal and vertical offset for output:
[x_offset, y_offset] = lcloffset(1);
offset = [x_offset, y_offset, x_offset, y_offset];
p = positions;
blk = 'Output';
outport(hTar, blk);
set_param([sys '/' blk], 'Position', p.output + offset);
% Connect last stage to output:
add_line(sys, [ic{1} '/1'], [blk '/1'],'autorouting','on');
end
% --------------------------------------------------------------
% Utility functions
% --------------------------------------------------------------
function [x_offset, y_offset] = lcloffset(stage)
x_offset = 0;
y_offset = 80*stage-80;
% --------------------------------------------------------------
function p = positions
p.input = [90, 52, 120, 67];
p.section = [225, 40, 300, 80];
p.sum = [340, 47, 365, 72];
p.output = [440, 52, 470, 67];
% [EOF]