通信原理与matlab仿真v3 第十一章 QPSK调制解调器仿真(12)

算法工匠

2月前浏览1006

关注

继续讲解！本文给出优化后的完整程序。之前的文章详细讲解了学习QPSK调制解调器的各个过程，希望大家能够真正的消化和理解，这样才能学习更难的知识！本文内容很长，算是这个系列文章的总结，希望读者能够耐心的看下去！只要看完，肯定有很大收获！当然这也是含有重磅代码的文章！希望读者们能够珍惜！！！程序开始慢慢实用化，学习代码渐进的过程非常难得！

 

归档版本的程序！

赶快上硬货！

一睹为快吧！

 

%%%%%%     程序说明

% 中频信号仿真，分析误码性能。

% 修订程序完成！

% 测试不同的信噪比下的误码性能！

% 注意以下内容：

% 本程序没有将信噪比和比特信噪比进行等效换算！！！

% 本程序采用相干解调，没有深入载波恢复的内容！

% 没有进行组帧内容的仿真！

%%%                       软件环境 

% 软件版本：R2015b

clear all,close all;

%-------------            程序主体            ---------------------%

%%%  初始化参数   %%%

T = 1;                        % 基带信号宽度

fc = 10;                      % 载波频率

ml = 2;                       % QPSK

nb = 900;                  % 传输的比特数

fs = 200;                     % 采样频率 Hz

SNR = 100;                 % 信噪比

%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%             调制部分              %%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 信源 基带信号的产生

data = randi([0,1],1,nb);        

% 产生随机信号

% 串并转换

idata = data(1:ml:(nb-1));  

qdata = data(2:ml:nb);

%%% 星座图映射

idatamap = zeros(1,nb/2);        

% 存放映射之后的信号

qdatamap = zeros(1,nb/2);

% QPSK星座图 正方形 即为pi/4QPSK A方式

% 实际中不用菱形的星座图

for i=1:nb/2

    if idata(i)==1&&qdata(i)==1  

        idatamap(i) = 1;

        qdatamap(i) = 1;

    elseif idata(i)==0&&qdata(i)==1

        idatamap(i)=-1;

        qdatamap(i)=1;

    elseif idata(i)==0&&qdata(i)==0

        idatamap(i)=-1;

        qdatamap(i)=-1;

    elseif idata(i)==1&&qdata(i)==0

        idatamap(i)=1;

        qdatamap(i)=-1;

    end

end

% 内插

Izero = zeros(1,nb/2);        

Qzero = zeros(1,nb/2);

sample_zero = 10;

for i = 1:sample_zero*nb/2

    if rem(i,sample_zero)==1

        Izero(i) = idatamap(fix((i-1)/sample_zero)+1);

        Qzero(i) = qdatamap(fix((i-1)/sample_zero)+1);

    else

        Izero(i) = 0; % 插0

        Qzero(i) = 0;

    end

end

%%% 脉冲成型滤波器

% 设计参数

rollfactor = 1;

psf = rcosdesign(rollfactor, 5, 10,'sqrt');

% 之前用过rcosflt函数，该函数在2015年后被软件剔除！

来源：通信工程师专辑