function [RX_I,RX_Q,data]=generate_random_data(N) 
%**************************************************************************
% generate N random datas
%**************************************************************************
rand('state',sum(100*clock));
data=rand(1,N);
for i=1:N
    if (data(i)<=0.5)
        data(i)=0;
    else
        data(i)=1;
    end;
end;
data(N+1:N+8)=zeros(1,8);                  %生成八個零在隨機數的后面，
                                           %是為了仿真“點亮”的過程
%**************************************************************************
% compute the data transmitted in the channel
%**************************************************************************
G1=171;                                                 %G1,G2均為八進制
G2=133;                                                 %并且約束長度為七

G1_vector=fliplr(dec2binvec(oct2dec(G1)));
G2_vector=fliplr(dec2binvec(oct2dec(G2)));

shift_d=zeros(1,6);
N=length(data);

for i=1:N
    t1=0;
    t2=0;
    for j=1:6
        t1=xor(t1,and(shift_d(j),G1_vector(j+1)));
        t2=xor(t2,and(shift_d(j),G2_vector(j+1)));
        
        RX_I(i)=xor(data(i),t1);
        RX_Q(i)=xor(data(i),t2);
    end;
       
    T_vector=dec2binvec((2*binvec2dec(shift_d)),6);     %狀態值右移的過程
    shift_d=T_vector(1:6);
    shift_d(1)=data(i);
end;