## Copyright (C) 2020 Jens
## 
## 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 3 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, see
## <https://www.gnu.org/licenses/>.

## -*- texinfo -*- 
## @deftypefn {} {@var{retval} =} mc_pd_eval (@var{input1}, @var{input2})
##
## @seealso{}
## @end deftypefn

## Author: Jens <jens@orion>
## Created: 2020-08-07

function retval = mc_pd_eval (delay)

Nc = 4
Nspc = 24;

ka = 0.9;
xmc_i = [ka*ones(1, Nspc/2) -ka*ones(1, Nspc/2)];
xmc_q = [-ka*ones(1, Nspc/4) ka*ones(1, Nspc/2) -ka*ones(1, Nspc/4)];
x = repmat(xmc_i, 1, Nc);

lo_i = repmat(xmc_i, 1, Nc);
lo_q = repmat(xmc_q, 1, Nc);

if delay > 0
	hdly = [zeros(1, delay) 1];
	xdly = filter(hdly, 1, x);
else
	xdly = x
end
	
y_i = lo_i.*xdly;
y_q = lo_q.*xdly;

close all;
subplot(2, 1, 1)
plot(1:Nc*Nspc, lo_i, 1:Nc*Nspc, lo_q); legend('Lo_I','Lo_Q'); grid

subplot(2, 1, 2)
plot(1:Nc*Nspc, y_i, 1:Nc*Nspc, y_q); legend('I','Q'); grid

endfunction