dolden.txt

一种在扩频领域应用十分广泛的编码序列,GOLDEN码序列的MATLAB仿真实现!

%function goldpn=goldgen(GoldGenPoly) %GoldGenPoly是由两个m序列的生成多项式组成的矩阵
GoldGenPoly=[1 0 0 1 ;0 0 1 1 ];
%pn_len = 2^(size(GoldGenPoly,2))-1; %size(GoldGenPoly,2)为阶数
%ini1_state = [zeros(1,3),1]; % m1序列的初始状态
%ini2_state = [1,zeros(1,3)]; % m2序列的初始状态
ini1_state = [0 0 1 1];
ini2_state = [0 1 0 1];
%--------------------------------------------------------------------------------------
% -- to generate the 1st and 2nd m sequences. then the gold sequence is generated
%for i=1:pn_len
for i=1:15
m1_out(i) = ini1_state(length(ini1_state)); 
m2_out(i) = ini2_state(length(ini2_state));
%for i01=1size(GoldGenPoly,2)-1)
s1_step = sum(GoldGenPoly(1,1:4).*ini1_state);
s2_step = sum(GoldGenPoly(2,1:4).*ini2_state);

for j=4:-1:2
ini1_state(j) = ini1_state(j-1);
ini2_state(j) = ini2_state(j-1); 
end

ini1_state(1) = mod(s1_step,2);
ini2_state(1) = mod(s2_step,2);

goldpn(i) = xor(m1_out(i),m2_out(i));
end
m1_out;
m2_out;
m1_out=m1_out.*2-1;
m2_out=m2_out.*2-1;
% m1_out=[m1_out m1_out];
% m2_out=[m2_out m2_out];
% 
% for indx=1:15
% tmp=0;
% for i=1:15
% tmp=tmp+m1_out(i)*m2_out(i+indx);
% end
% tcf(indx)=tmp;
% end
goldpn = goldpn.*2-1;