用MATLAB实现的BPSK调制解调源程序

function  output_frame  =  demodulation(input_modu,  index)%  demodulation  for  IEEE802.11a%  Input:        input_modu,  complex  values  representing  constellation  points%                      index%  Output:      output_frame,  output  bit  stream  (data  unit  is  one  bit)%  In  this  version,  increase  the  quatilization  levels  into  8.%  note:  Matlab  index  starts  from  1Q_length=length(input_modu);QAM_input_I  =  real(input_modu);QAM_input_Q  =  imag(input_modu);output_frame  =  zeros(1,length(input_modu)*index);switch  indexcase  1,        BPSK_Demodu_I  =  [0  1];            %f(m)=(m+1)/2  +  1,  so  I=-1  --->  1,  I=1  --->  2        idx  =  find(QAM_input_I>1);        QAM_input_I(idx)  =  1;        idx  =  find(QAM_input_I<-1);        QAM_input_I(idx)  =  -1;        output_frame  =  BPSK_Demodu_I(round((QAM_input_I+1)/2)  +  1);