Files
matlab/mc_pd_eval.m
jens 05650c370d - added
git-svn-id: http://moon:8086/svn/matlab/trunk@129 801c6759-fa7c-4059-a304-17956f83a07c
2020-08-09 11:07:18 +00:00

57 lines
1.4 KiB
Objective-C

## 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