Вот код двух последовательностей (Голда и Франка), нужно их сложить любым способом, и получить новую последовательность

goldseq1 = comm.GoldSequence('FirstPolynomial', 'x^6+x+1', 'SecondPolynomial', 'x^6+x^5+x^2+x+1', ...
'FirstInitialConditions', [0 0 0 0 0 1], ...
'SecondInitialConditions', [0 0 0 0 0 1], ...
'Index', 0,'SamplesPerFrame', 64); % последовательность Голда 1
goldseq2 = comm.GoldSequence('FirstPolynomial', 'x^6+x+1', 'SecondPolynomial', 'x^6+x^5+x^2+x+1', ...
'FirstInitialConditions', [0 0 0 0 0 1], ...
'SecondInitialConditions', [0 0 0 0 0 1], ...
'Index', 1,'SamplesPerFrame', 64); % последовательность Голда 1
x = goldseq1 ();
y = goldseq2 ();
a = xcov(goldseq2());
figure ()
 plot (abs (a)/max(a))
 % Исходные данные
N = 64; % Длина многофазной последовательности
mfprs_type = 'Frank';
% mfprs_type = 'Chu';
% mfprs_type = 'Milewski';
 
% Frank sequence
if (strcmp(mfprs_type, 'Frank'))
    Q = sqrt(N);
    r = euler (Q); 
    % length(r) - количество последовательностей Фрэнка
    fprintf( 'Всего последовательностей в ансамбле Фрэнка = %d\r', length(r) );
    fprintf('Длина последовательности Фрэнка = %d\r', N);
    k = 0:1:Q-1;
    j = 0:1:Q-1;
    kj = zeros(Q,Q);
    for ii = 1:1:Q
        for iii = 1:Q
            kj(ii, iii) = k(ii)*j(iii);
        end
    end
    seqs_amount = length(r);
%     seqs_amount = 1;
    source_multiphase_seqs_bit = zeros(seqs_amount, N);
 
    for kk = 1:1:seqs_amount
        alpha = exp(2*pi*r(kk)*1i/Q);
        for ii = 1:1:N
            source_multiphase_seqs_bit(kk, ii) = alpha^kj(  floor( (ii-1)/Q )+1, rem( (ii-1), Q )+1  );
        end
    end
end
   fprintf( 'Формирующее число = %d\rДлина ПСП = %d\rВсего ПСП = %d\r', Q, N, r );
  
kk = 1; % номер ПСП в ансамбле
z = xcov (source_multiphase_seqs_bit(kk,:));
figure ()
plot (abs (z)/ max(z))