QPSK调制系统matlab仿真函数

function [y,signal] = qpsk_function(data_length,fs,fc,EbN0,interpolation_rate,Sample_position)
%qpsk 调制系统
%y误码率
%signal 发送信号

%data_length 数据长度
%fs  数据频率
%fc  载波频率
%interpolation_rate 内插系数
%Sample_position 采样点
w_carrier = fc/fs;
SNR = EbN0 +10*log10(4/interpolation_rate); 

%% 发送端
data_send = randi([0,1],data_length,1);  %0/1序列
data_bio = 1-2*data_send;                %对极序列0->1  1->-1
data_i = data_bio(1:2:end);%奇数位
data_q = data_bio(2:2:end);%偶数位
data_upsample_i = upsample(data_i,interpolation_rate);  %内插（每个数据间插入interpolation_rate个0）
data_upsample_q = upsample(data_q,interpolation_rate);  %内插

alpha = 0.35;
span = 6;
hrc1 = rcosdesign(0.35,span,interpolation_rate,'sqrt'); 
hrc = ones(1,interpolation_rate);                    %方形滤波器

%脉冲成型
baseband_i = conv(data_upsample_i,hrc)';                   %成型
baseband_q = conv(data_upsample_q,hrc)';                   %成型

%% 调制
carrier_i = cos(2*pi*w_carrier*[0:length(baseband_i)-1]);   %载波
carrier_q = -sin(2*pi*w_carrier*[0:length(baseband_q)-1]);   %载波
send_qpsk_i = baseband_i.*carrier_i;
send_qpsk_q = baseband_q.*carrier_q;
send_qpsk = send_qpsk_i +send_qpsk_q;  %BPSK调制
 carrier = carrier_i + j*carrier_q;                      %合成复载波


 %% 信道
 recv = awgn(send_qpsk,SNR,'measured');                %加性高斯白
    
 %% 接收

 recv_carrier_i =  carrier_i.*recv;
 recv_carrier_q =  carrier_q.*recv;

 recv_bio_i = conv(recv_carrier_i,hrc);              %匹配滤波
 recv_bio_q = conv(recv_carrier_q,hrc);              %匹配滤波器
 recv_bio_i = recv_bio_i(length(hrc):end-length(hrc)+1);
 recv_bio_q = recv_bio_q(length(hrc):end-length(hrc)+1);

  recv_data_i = downsample(recv_bio_i,interpolation_rate,Sample_position);%下采样——恢复（数据，内插插入点，最佳采样点）
  recv_data_q = downsample(recv_bio_q,interpolation_rate,Sample_position);


 %画星座图 
   scatterplot(recv_data_i+1j*recv_data_q);

  
%判决 计算误码
demod_data = demodulatorQPSK (recv_data_i, recv_data_q,data_length);%demodulatorQPSK（I,Q,符号数） qpsk判决函数
errors = sum(data_send ~= demod_data);%~=统计不相等元素个数
BER =errors /data_length;


 y=BER;
 signal = send_qpsk;

function [ demodulatedBitArrayQPSK ] = demodulatorQPSK (x, y, BIT_NUMBER)
% x 奇数位  y 偶数位     BIT_NUMBER——符号数
demodulatedBitArrayQPSK = zeros(BIT_NUMBER, 1);
%-1->1 1->0
for i = 1:BIT_NUMBER/2

    if x(i) < 0
        
        if y(i) < 0
            
            demodulatedBitArrayQPSK((2*i) - 1) = 1;
            demodulatedBitArrayQPSK(2*i) = 1;
            
        else
            
            demodulatedBitArrayQPSK((2*i) - 1) = 1;
            demodulatedBitArrayQPSK(2*i) = 0;
            
        end
       
    else
        
        if y(i) < 0
            
            demodulatedBitArrayQPSK((2*i) - 1) = 0;
            demodulatedBitArrayQPSK(2*i) = 1;
            
        else
            
            demodulatedBitArrayQPSK((2*i) - 1) = 0;
            demodulatedBitArrayQPSK(2*i) = 0;
            
        end
        
    end
    
end

end