function [plan_up,plan_down]=setupPSDplan(plan, psd_up, psd_down, method, nosignal,capval);
%% ===========================================================================
%setupPSDplan - Generate a PSD description from a simple list
%
% Parameter:    plan		PSD mask description
% Parameter:    psd_up		PSD-mask for upstream
% Parameter:    psd_down	PSD-mask for downstream
% Parameter:    method  	Interpolation method ('Linear' or 'Log-Linear')
% Parameter:    nosignal  	Signal level when not sending 
% Parameter:    capval  	Maximum signal level
% Returns:      plan_up 	PSD plan for upstream
% Returns:      plan_down 	PSD plan for downstream
%
% Example(s):
%   [tmp_upstream,tmp_downstream]= setupPSDplan({...
%      fmin 'u' .138e6 'd' 2.0e6 'u' 4.3e6, 'd' 7.1e6 ...
%      'u' 10.1e6 'd' 14.35e6 'u' fmax }, ...
%      upstream_val, downstream_val, 'Log-Linear', -120);
%
% Algorithmic details:
%
% Reference:
%% ===========================================================================

%% ===========================================================================
% Copyright (C):                                        
%       1999-2000 by Telia Research AB, Lulea, Sweden;                
%       2000-2003 by Forschungszentrum Telekommunikation Wien, Austria;
%                                                         All rights reserved.
% Project       : FTW's xDSLsimu
% Author(s)     : Tomas Nordstrom (Tomas.B.Nordstrom@@Telia.se)
%               : Daniel Bengtsson (Daniel.J.Bengtsson@Telia.se)
%		          : Frank Sjöberg
%
% CVS:       $Id: setupPSDplan.m,v 3.3 2003/11/03 12:33:22 tono Exp $
%% ===========================================================================
% Change History
%      1999-01-14 (ToNo) Created
%      1999-07-27 (FS)   Modified from setupPSDplan to handle non-flat PSD
%      2000-04-06 (UvAn) Changed default value for "nosignal" to -140
%      2002-10-11 (ToNo) Made internal function external
%                        Fixed off-by-one errors in start and stop index     
%      2003-01-03 (ToNo) Octave to Matlab compatibility port
%      2003-11-03 (ToNo) Lists now use cell arrays for both Octave and Matlab
%% ===========================================================================
global ex;

if nargin<3,
    error('Not enough input arguments.');
end

if nargin<4
    % Default interpolation method is linear
   method = 'Linear';	
end;

if nargin<5,
    % Default value; (before we used nosignal=ex.param.backgroundNoise);
    nosignal=-140;      
end

if nargin<6,
    capval=0;
end

plen=length(plan);
if plen<3,
    error('setupPSDplan must have length > 3');
end;

if ischar(plan{1}) | ischar(plan{plen}),
    error('Start and stop setupPSDplan with a frequency.');
end;

plan_up=[];
plan_down=[];
for i=1:2:plen-2,
    start=plan{i};
    direction=plan{i+1};
    stop=plan{i+2};
% Problem with ischar (actually isalpha) in octave?    
%    if ischar(start) | ischar(stop),
%        error('Frequencies should be around a direction.');
%    end;
    if stop<start
        error('Frequencies must be increasing.');
    end;
    switch direction
        case 'u'
           tmp_psd = setupPSD_def2PSD(psd_up, method, [start+1 stop-1], nosignal,capval);
           if isempty(plan_up)
               plan_up = [start nosignal tmp_psd stop nosignal];
           else
               plan_up = [plan_up start nosignal tmp_psd stop nosignal];
           end;
        case 'd'
           tmp_psd = setupPSD_def2PSD(psd_down, method, [start+1 stop-1], nosignal,capval);
           if isempty(plan_down)
               plan_down = [start nosignal tmp_psd stop nosignal];
           else
               plan_down = [plan_down start nosignal tmp_psd stop nosignal];
           end;
        case 'n'
           if isempty(plan_up)
               plan_up=[start+1 nosignal stop-1 nosignal];
           else
               plan_up=[plan_up start+1 nosignal stop-1 nosignal];
           end;
           if isempty(plan_down)
               plan_down=[start+1 nosignal stop-1 nosignal];
           else
               plan_down=[plan_down start+1 nosignal stop-1 nosignal];
           end;
        otherwise
            estr=sprintf('Unknown direction %s', direction);
            error(estr);       
    end;
end;