%對序列進行擴頻
%x只能是隨機點，不能是電瓶形式
function    [sqe,z11,z22]=kuopin(y1,y2,connections,T,Ts)  
N=0.001/T;
fs=1/Ts;
tdi=0:2*Ts:(N*T)-Ts; 
%==========================================================================
%產生M序列
%==========================================================================
%m序列的圖形
sqe=m_sequence(connections);
ppp=fs/(N*255*1000);
for i=1:length(sqe),               %計算碼元得值
   if(sqe(i)==1),                  %如果信息為1
      for j=1:ppp                %該碼元對應得點取值為1
         y5((i-1)*ppp+j)=1;
      end
   else
      for j=1:ppp,                %反之，信息元為0，則變成負電瓶
         y5((i-1)*ppp+j)=-1;
      end;
   end;
end
figure(10)
n=length(connections);
t=linspace(0,(2^n-1),length(y5));
plot(t,y5);
axis([0 (2^length(connections)-1) -2 2]);
title('m序列圖形');
xlabel('電瓶數'); 

%=========================================================================
%實現擴頻
%=========================================================================
%實現擴頻
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(7); 
subplot(211);
plot(tdi,z11);
axis([0 N*T -1.5 1.5]);
title('I通道擴頻后圖形z11');
xlabel('時間軸                                               t'); 



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