%對(duì)序列進(jìn)行擴(kuò)頻
%x只能是隨機(jī)點(diǎn)，不能是電瓶形式
function    [,z11,z22]=kuopin(y1,y2,connections,Ts,N)   
%==========================================================================
%產(chǎn)生M序列
%==========================================================================
%m序列的圖形
T=0.001/N;
tdi=0:2*Ts:N*T-Ts;
sqe=m_sequence(connections);

ppp=150;
for i=1:length(sqe),               %計(jì)算碼元得值
   if(sqe(i)==1),                  %如果信息為1
      for j=1:ppp                 %該碼元對(duì)應(yīng)得點(diǎn)取值為1
         y5((i-1)*ppp+j)=1;
      end
   else
      for j=1:ppp                %反之，信息元為0，則變成負(fù)電瓶
        y5((i-1)*ppp+j)=-1;
      end;
   end;
end
figure(10)
n=length(connections);
M=length(y5);
t=linspace(0,2^n-1,M);
plot(t,y5);
axis([0 (2^length(connections)-1) -2 2]);
title('m序列圖形');
xlabel('電瓶數(shù)'); 

%=========================================================================
%實(shí)現(xiàn)擴(kuò)頻
%=========================================================================
%實(shí)現(xiàn)擴(kuò)頻
p1=length(y5);        %150*255
for n=1:length(y1),
    if(y1(n)==1)
        for i=1:p1
            z11((n-1)*p1+i)=y5(i);
        end
    else
          for i=1:p1
              z11((n-1)*p1+i)=-y5(i);  
          end;
    end;
end;

for n=1:length(y2),
    if(y2(n)==1)
        for i=1:p1
            z22((n-1)*p1+i)=y5(i);
        end
    else
          for i=1:p1
              z22((n-1)*p1+i)=-y5(i);  
          end;
    end;
end;

%=======================================================================
%畫圖
%=======================================================================
figure(2); 
subplot(312);
plot(tdi,z11);
axis([0 N*T -1.5 1.5]);
title('I通道擴(kuò)頻后圖形z11');
xlabel('時(shí)間軸                                               t'); 



figure(3)
subplot(312);
plot(tdi,z22);
axis([0 N*T -1.5 1.5]);
title('Q通道擴(kuò)頻后圖形z22');
xlabel('時(shí)間軸                                               t'); 
end