gen_test_sig.m

function [tx] = gen_test_sig(bits, noise_power, channel_bauds, freqs)
configuration;
num_sigs = length(channel_bauds);

for baud=channel_bauds
    if mod(F_S, baud) ~= 0
        disp('Sample Rate must be a multiple of all baud rates!');
        fprintf('%f != k%f\n', F_S, baud);
        tx=0;
        return;
    end
end

samps_per_sym =  F_S./channel_bauds;

%% Run the Transmitter
tx = [];
for idx=1:num_sigs
    % Generate
    sig = gen_sig(bits, ceil(samps_per_sym(idx)), MODULATION);

    if ceil(samps_per_sym(idx)) ~= samps_per_sym(idx)
        sig = resample(sig, 1/ceil(samps_per_sym(idx)), 1/samps_per_sym(idx));
    end

    % Frequency Shift
    t = 0:1/F_S:(length(sig) - 1)/F_S;
    sig = sig.*exp(j*2*pi*freqs(idx)*t);

    % Add to WB signal
    tx = add_and_zero_pad(tx, sig);
end

% Compute Noise variance based on desired noise power
% pretty sure this is wrong
noise_variance = 10.^(0.1*noise_power)/(F_S*2);

% add noise
tx = awgnChannel(tx, noise_variance, F_S, 0, 0, 0, 0);

end

function [samples] = gen_sig(bits, samples_per_sym, modulation)
    configuration;
    SAMPLES_PER_SYMBOL = samples_per_sym;
    recompute_configuration;
    samples = MyTransmitter(bits);
end

function [result] = add_and_zero_pad(left, right)
    if length(left) == length(right)
        result = left + right;
        return
    end
    if length(left) > length(right)
        tmp = left;
        left = right;
        right = tmp;
    end
    left = [left zeros(1, length(right) - length(left))];
    result = left + right;
end