kuopin.m

就是qpsk用matlab实现调制解调以及扩频接扩

%对序列进行扩频
%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