bpsk.m

采用MATLAB仿真BPSK-OFDM的BER特性。从理论和实际两方面把输出波形作为对比。

%********************************************************
% Assignment 6
% bpskmod.m
% Function to perform BPSK modulation
% Programmed by HZL
%  2007/12/14
%****************** variables *************************
% out :output data
% num_data : Number of input data
% num_samp : number of sampling
% mes_orig  :data sequence
% *******************BPSK Modulation ************************

 %**********************原始数据产生**************************************
 num_data =100;
 num_samp=100;
 msg_orig=randint(1,num_data);  %生成数据序列
 data_1=msg_orig*2-1;        %单/双极性变换
 %************************************************************
 
 %**********************抽样*********************************
 for ii=1: num_data
     for iii=num_samp*(ii-1)+1:num_samp*ii;
     data_2(iii)=data_1(ii);
     end
 end
 %************************************************************
 
%****************modulation***********************************
for jj=1:num_data*num_samp
	carr(jj)=cos((jj-1)*pi/50);  %载波产生
end
  data_3=data_2.*carr;            % modulation
  out=data_3;                      %输出信号
%**********************************************************
data4=awgn(data_2,10,'measured');           %经过AWGN信道
%******************AWGN************************************

%******************解调************************************
data5=data_3.*carr;
%**********************************************************

%*********************判决*************************************
out=zeros(1,num_data);
for i=1:num_data;
    data_pan(i)=sum(data4((i-1)*num_samp+1:i*num_samp));
   if data_pan(i)>0;
       out(i)=1;
   else out(i)=-1;
   end
end
%****************************************************** 
%***********************绘图******************************
for i=1:100*num_samp 
   data2_g(i)=data_2(i);          %抽样数据
   data3_g(i)=data_3(i);          %调制数据
end 
figure(1);                   %调制波形
x=0+0.01:0.01:100;
subplot(2,1,1);
plot(x,data2_g);
axis([0 10 -1.4 1.4]);  
title('线性调制');
subplot(2,1,2);
plot(x,data3_g);
axis([0 10 -1.4 1.4]);
title('BPSK调制信号');
  
figure(2);                           %画经过高斯信道后的信号波形和频谱图
plot(x,data4);
axis([0 10 -1.4 1.4]);
title('经过高斯信道后的信号波形');

figure(3);                           %解调后的接受到的信号波形和频谱
plot(x,data5);
axis([0 10 -1.4 1.4]);
title('解调后的信号波形');
%******************************************************************

%******************** end of file *************************************