第一次记录QPSK，BSPK，MPSK，QAM—MATLAB实现

最近有偶然的机会学习了一次QPSK防止以后忘记又得找资料，这里就详细的记录一下

        基于 QPSK 的通信系统如图 1 所示，QPSK 调制是目前最常用的一种卫星数字和数 字集群信号调制方式，它具有较高的频谱利用率、较强的抗干扰性、在电路上实现也较为简单，在后文仅仅使用MATLAB进行模拟。

图 1 基于QOSK调制的通信系统模型

        其相位图如图 2 所示，二进制数 0 和 1 分别表示两个相位，为了提高传输速率，通 常可以采用多项调制的方法，即将待发的数字信号按两比特一组的方式组合，两位二进 制数的组合方式又四种—(00，01，10，11)。每个组合是一个双比特码，通常可以用四 个不同的相位值表示这四组双比特码。在传输过程中，相位改变一次，传输两个二进制 数。这种调相方法成为四项调相或四项调制，广泛广泛应用 于卫星链路、数字集群等 通信业务。

图2相位图

如图 3 所示，QPSK 信号可以采用正交调制器来实现 

图3 QPSK信号生成原理图

 根据图1和图3的原理图，就可以开始着手写MATLAB代码了。 我将其分为以下几个方面

  1. 模拟源信号
  2. 源信号->双极性信号
  3. 双极性信号转为QPSK信号
  4. QPSK信号经过模拟信道传输，引起失真。
  5. 采用相干解调法分离QPSK信号。
  6. 低通滤波过滤噪音
  7. 抽样判决获得信号
  8. 最后将信号极性反转获得源信号

QPSK详细代码如下，在代码中已经做好注释了。

%% 采用代码实现的4PSK通信系统仿真
% 假设在T=1,加入高斯噪声
clc;
clear all;
close all;
%% 发端
% 1.调制,生成二进制信号
bit_in = randi([0 1],1000,1); 
% 2.变为双极性码
data = -2*bit_in+1;        
% 3.串并转换模块：奇数位为I，偶数为为Q
data_I  = data(1:2:1000);   % 间断获取 I
data_Q  = data(2:2:1000);   % 间断获取 Q
data_I1=repmat(data_I',20,1);
data_Q1=repmat(data_Q',20,1);

% 按列优先将data_I1中的数据存入data_I2
for i=1:1e4
    data_I2(i)=data_I1(i);
    data_Q2(i)=data_Q1(i);
end

% 4.产生升余弦的基带信号
f=0:0.1:1;
xrc=0.5+0.5*cos(pi*f); 

data_I2_rc=conv(data_I2,xrc)/5.5;
data_Q2_rc=conv(data_Q2,xrc)/5.5;

figure
subplot(2,2,1)
stem(bit_in(1:20),'black','LineWidth',2);
axis([0,20,0,1]);
title("发送的消息序列");
subplot(2,2,2)
plot(f,xrc,'black','LineWidth',2);
title("升余弦信号");
subplot(2,2,3)
plot(data_I2_rc(1:20),'black','LineWidth',2);
title("升余弦I信号");
subplot(2,2,4)
plot(data_Q2_rc(1:20),'black','LineWidth',2);
title("升余弦Q信号");

% 5. 正交调制（调相法：将基带数字信号（双极性）与载波信号直接相乘的方法）
f1=1;                        % 载波频率
t1=0:0.1:1e3+0.9;

I_rc=data_I2_rc.*cos(2*pi*f1*t1);
Q_rc=data_Q2_rc.*sin(2*pi*f1*t1);
x=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);

figure(1)
subplot(2,1,1);
plot(t1,x,'black','LineWidth',2); xlabel('t'); ylabel('幅度'); 
grid on; 
axis([0 1/f1*10 -1.2 1.2]);  % 输出2个周期的信号
title('QPSK信号'); 

%% 6.仿真信道噪声
n0=rand(size(t1))/2;
y=x+n0;

subplot(2,1,2);
plot(t1,y,'black','LineWidth',2); xlabel('t'); ylabel('幅度'); 
grid on; 
axis([0 1/f1*10 -2 2]);  % 输出2个周期的信号
title('带噪声的QPSK信号'); 


%% 仿真接收端
% 7.正交解调：只能采用相干解调
I_demo=y.*cos(2*pi*f1*t1);
Q_demo=y.*sin(2*pi*f1*t1);
% 8.低通滤波
I_recover=conv(I_demo,xrc);    
Q_recover=conv(Q_demo,xrc);
I=I_recover(11:10010);
Q=Q_recover(11:10010);

% 9.抽样判决
data_recover=[];
for i=1:20:10000
    data_recover=[data_recover I(i:1:i+19) Q(i:1:i+19)];
end
bit_recover=[];
for i=1:20:20000
    if sum(data_recover(i:i+19))>0
        data_recover_a(i:i+19)=1;
        bit_recover=[bit_recover 1];
    else
        data_recover_a(i:i+19)=-1;
        bit_recover=[bit_recover -1];
    end
end

% 10.变为单极性码
bit_recovered=(1-bit_recover)/2; 
figure(2)
subplot(2,1,1)
stem(bit_in(1:20),'black','LineWidth',2);
axis([0,20,0,1]);
title("发送的消息序列");

subplot(2,1,2)
stem(bit_recovered(1:20),'black','LineWidth',2);
axis([0,20,0,1]);
title("接收的消息序列");

figure(3)
subplot(2,1,1)
stem(bit_in(1:20),'black','LineWidth',2);
axis([0,20,0,1]);
title("发送的消息序列");

subplot(2,1,2)
stem(data(1:20),'black','LineWidth',2);
axis([0,20,-1,1]);
title("双极性码");
set(gcf,'color','w')

BSPK的代码也贴在这里

clc
clear
close all;

%%初始化参数设置
data_len = 100000;                       % 原始数据长度
SNR_dB = 0:10;                           % 信噪比 dB形式
SNR = 10.^(SNR_dB/10);                   % Eb/N0
Eb = 1; % 每比特能量
N0 = Eb./SNR ; %噪声功率
error2 = zeros(1,length(SNR_dB));          % 码元错误个数
simu_ber_BPSK = zeros(1,length(SNR_dB));         % 仿真误误码率
theory_ber_BPSK = zeros(1,length(SNR_dB));   % BPSK理论误码率
demod2_signal= zeros(1,data_len);         % 解调信号

%%基带信号产生
data_source = round(rand(1,data_len));  % 二进制随机序列

%%BPSK基带调制   
send_signal2 = (data_source - 1/2)*2; % 双极性不归零序列 

%%高斯信道无编码
for z = 1:length(SNR_dB)
     noise2 = sqrt(N0(z)/2) * randn(1,data_len); %高斯白噪声
     receive_signal2 = send_signal2 + noise2;
     demod_signal2 = zeros(1,data_len);
        for w = 1:data_len
                if (receive_signal2(w) > 0)
                demod_signal2(w) = 1;              % 接收信号大于0  则判1
                else
                demod_signal2(w) = 0;              % 接收信号小于0  则判0
                end
        end
        %统计错误码元个数
       for w = 1:data_len
           if(demod_signal2(w) ~=data_source(w) )
                  error2(z) = error2(z) + 1;    % 错误比特个数
           end
       end
           %计算误码率
        simu_ber_BPSK(z) = error2(z) / data_len;         % 仿真误比特率
        theory_ber_BPSK(z) = qfunc(sqrt(2*SNR(z)));   % 理论误比特率
end

    

%%二进制序列、基带信号图像
figure(1);
stem(data_source);
title("二进制随机序列");
axis([0,50,0,1]);
figure(2);
stem(send_signal2);
title("BPSK基带调制--发送信号");
axis([0,50,-1.5,1.5]);

figure(4);
stem(noise2);
title("高斯白噪声");
axis([0,50,-0.5,0.5]);

figure(5)
stem(receive_signal2);
title("接收信号");
axis([0,50,-1.5,1.5]);

figure(7)
stem(demod_signal2);
title("解调信号");
axis([0,50,0,1]);

figure(8);
 semilogy(SNR_dB,simu_ber_BPSK,'M-X',SNR_dB,theory_ber_BPSK,'k-s');     

grid on;                                      
axis([0 10 10^-5 10^-1])                      
xlabel('Eb/N0 (dB)');                     
ylabel('BER');                                  
 legend('BPSK仿真误码率','BPSK理论误码率');  

%%画星座图
scatterplot(send_signal2);
title('发送信号星座图');
scatterplot(receive_signal2);
title('接收信号星座图');
scatterplot(demod_signal2);
title('解码信号星座图');

MPSK代码

clc;
clear all;
close all;
%% 调用库函数实现MPSK的通信系统仿真

M=4;
cycl=80;       % 运行次数
SNR=0:1:30;    % 信噪比
s=randi([0 M-1],1,1000);  % 输入信息一行1000列;

BER1=zeros(cycl,length(SNR));

for  n=1:cycl
    for k=1:length(SNR)
        x=pskmod(s,M,pi/4);               % M进制PSK
        y=awgn(x,SNR(k),'measured');      % 在传输序列中加入噪声
        r=pskdemod(y,M,pi/4);             % 解调
        r1=reshape(r',1,[]);
        [num,rat]=biterr(r1,s,log2(M));           % 误码率计算
        BER1(n,k)=rat;
    end
end

figure(1)
subplot(2,1,1)
stem(s(1:20),'black','LineWidth',2);
axis([0,20,0,M]);
title("发送的消息序列");

subplot(2,1,2)
stem(r1(1:20),'black','LineWidth',2);
axis([0,20,0,M]);
title("接收的消息序列");


%% 8PSK
M=8;
cycl=80;       % 运行次数
SNR=0:1:30;    % 信噪比
s=randi([0 M-1],1,1000);  % 输入信息一行1000列;

BER2=zeros(cycl,length(SNR));

for  n=1:cycl
    for k=1:length(SNR)
        x=pskmod(s,M,pi/4);               % M进制PSK
        y=awgn(x,SNR(k),'measured');      % 在传输序列中加入噪声
        r=pskdemod(y,M,pi/4);             % 解调
        r1=reshape(r',1,[]);
        [num,rat]=biterr(r1,s,log2(M));           % 误码率计算
        BER2(n,k)=rat;
    end
end

figure(2)
subplot(2,1,1)
stem(s(1:20),'black','LineWidth',2);
axis([0,20,0,M]);
title("原始消息序列");

subplot(2,1,2)
stem(r1(1:20),'black','LineWidth',2);
axis([0,20,0,M]);
title("传递消息序列");

% 统计平均误码率
figure(3)

BER1=mean(BER1);
subplot(2,1,1)
semilogy(SNR,BER1,'k-o','LineWidth',2);
xlabel('SNR/dB'); ylabel('BER');
title("4PSK");
grid on

BER2=mean(BER2);
subplot(2,1,2)
semilogy(SNR,BER2,'k-o','LineWidth',2);
xlabel('SNR/dB'); ylabel('BER');
title("8PSK");
grid on

QAM代码也在这里

clc;
clear all;
close all;
%% 基于16QAM的通信系统仿真
% 发端
nbit=10000; 
M=16;                                              % M表示QAM调制的阶数
k=log2(M);
graycode=[0 1 3 2 4 5 7 6 12 13 15 14 8 9 11 10];  % 格雷映射编码规则
EsN0=5:20;                         % 信噪比范围
snr=10.^(EsN0/10);                 % 将db转换为线性值
% 产生16进制的消息符号
s=randi([0,1],1,nbit);         
s_reshape=reshape(s,k,nbit/k)';    % 对数据流进行分组，对于16QAM，则每4位一组
msg=bi2de(s_reshape,'left-msb');   % 转化成10进制，作为qammod的输入
% 进行格雷映射
msg1=graycode(msg+1);        
% 调制
r=qammod(msg1,M);         % 调用matlab中的qammod函数，16QAM调制方式的调用(输入0到15的数，M表示QAM调制的阶数)得到调制后符号
spow1=norm(r).^2/nbit;    % 取a+bj的模.^2得到功率除整个符号得到每个符号的平均功率
for i=1:length(EsN0)
    % 信道
    sigma=sqrt(spow1/(2*snr(i)));                          % 16QAM根据符号功率求出噪声的功率
    x=r+sigma*(randn(1,length(r))+1i*randn(1,length(r)));  % 16QAM混入高斯加性白噪声
    % 16QAM的解调
    y1=qamdemod(x,M);             
    % 格雷逆映射
    y2=graycode(y1+1);        % 返回译码出来的信息，十进制
    test=de2bi(y2,k,'left-msb');
    y3=reshape(test',1,nbit);
    [err1,ber1(i)]=biterr(s,y3); 
end
%% 绘图
figure(1)
subplot(2,1,1)
stem(s(1:20),'black','LineWidth',2);
axis([0,20,0,1.2]);
title("发送的消息序列");

subplot(2,1,2)
stem(y3(1:20),'black','LineWidth',2);
axis([0,20,0,1.2]);
title("接收的消息序列");

scatterplot(r);           % 调用matlab中的scatterplot函数,画星座点图
scatterplot(x);           % 调用matlab中的scatterplot函数,画rx星座点图

% 16QAM调制信号在AWGN信道的性能
figure( )
semilogy(EsN0,ber1,'black','LineWidth',2);                            % ber ser比特仿真值 ser1理论误码率 ber1理论误比特率
title('16QAM调制信号在AWGN信道的性能分析');grid;
xlabel('Es/N0(dB)');                      
ylabel('误比特率');

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：/a/270020.html

如若内容造成侵权/违法违规/事实不符，请联系我们进行投诉反馈qq邮箱809451989@qq.com，一经查实，立即删除！

相关文章

Vue3设计目标和优化

Vue3设计目标和优化

vue3.0的设计目标是什么？做了哪些优化？ Vue3的设计目标：更小、更快、更友好、优化方案 一、设计目标 Vue3之前我们会面临的许多问题： 随着功能的增长，复杂组件的代码变得越来越难以维护起来缺少一种比较【干净】的…
阅读更多...
Redis基础-Redis概念及常见命令

Redis基础-Redis概念及常见命令

1.nosql数据库 NoSQL数据库是一种提供了非关系型数据存储的数据库系统，与传统的关系型数据库（如SQL数据库）不同。NoSQL数据库的特点是灵活性高，能够处理结构化、半结构化或非结构化数据。它们通常用于大数据和实时Web应用。NoSQL数…
阅读更多...
平衡二叉树的构建（递归

平衡二叉树的构建（递归

目录 1.概念：2.特点：3.构建方法：4.代码：小结： 1.概念： 平衡二叉树（Balanced Binary Tree），也称为AVL树，是一种二叉树，它满足每个节点的左子树和右…
阅读更多...
如何分析 JVM 内存瓶颈浅谈

如何分析 JVM 内存瓶颈浅谈

背景： 当操作系统内存出现瓶颈时，我们便会重点排查那个应用占用内存过大。对于更深一步分析内存的使用，就进一步去了解内存结构，应用程序使用情况，以及内存如何分配、如何回收，这样你才能更好地确定内存的…
阅读更多...
图像九宫格切分1x3、3x3 Python

图像九宫格切分1x3、3x3 Python

文章目录 1、需求2、实现2-1 贴图、切分2-2 GUI 3、运行效果4、代码 1、需求 把一个图像切分成 1x3 或者 3x3切分出来的图像比例希望都是 1:1 正方形如果图像尺寸满足 切分条件，自动填充一些“白边”然后继续切分如果填充了白边的话，希望能够调整原图像…
阅读更多...
挖到宝了，大数据分析工具做分析真的太快了

挖到宝了，大数据分析工具做分析真的太快了

随着企业越做越大，累积数据的速度越来越快，但分析的效率却不升反降，不利于数字化运营决策。但大数据分析工具的出现让这一现象成为过去，无他，就是大数据分析工具做分析的真的太快了，可在任意终端上随时按需…
阅读更多...
网络编程--socket编程

网络编程--socket编程

这里写目录标题 套接字概念通信原理总结 预备知识网络字节序简介字节转换函数 IP地址转换函数为什么单独列出函数原型sockaddr结构体 一级目录二级目录二级目录二级目录 一级目录二级目录二级目录二级目录 套接字 概念 Socket本身有插座的意思，但他是进程之间网络通…
阅读更多...
《LIO-SAM阅读笔记》1.IMU预积分模块

《LIO-SAM阅读笔记》1.IMU预积分模块

前言： LIO-SAM是一个多传感器融合的紧耦合SLAM框架，融合的传感器类型有雷达、IMU和GPS，其中雷达和IMU在LIO-SAM框架中必须使用的。LIO-SAM的优化策略采用了GTSAM库，GTSAM库采用了因子图的优化方法，其提供了一些列C的外…
阅读更多...
K8S理论

K8S理论

kubernetes：8个字母省略，就是k8s 自动部署自动扩展和管理容器化部署的应用程序的一个开源系统 k8s是负责自动化运维管理多个容器化程序的集群，是一个功能强大的容器编排工具 分布式和集群化的方式进行容器化管理 版本有1.15 .1.18 .1.20 …
阅读更多...
力扣-收集足够苹果的最小花园周长[思维+组合数]

力扣-收集足够苹果的最小花园周长[思维+组合数]

题目链接 题意： 给你一个用无限二维网格表示的花园，每一个 整数坐标处都有一棵苹果树。整数坐标 (i, j) 处的苹果树有 |i| |j| 个苹果。 你将会买下正中心坐标是 (0, 0) 的一块 正方形土地 ，且每条边都与两条坐标轴之一平行。 给你一个整…
阅读更多...
【flink番外篇】7、flink的State（Keyed State和operator state）介绍及示例 - 完整版

【flink番外篇】7、flink的State（Keyed State和operator state）介绍及示例 - 完整版

Flink 系列文章 一、Flink 专栏 Flink 专栏系统介绍某一知识点，并辅以具体的示例进行说明。 1、Flink 部署系列 本部分介绍Flink的部署、配置相关基础内容。 2、Flink基础系列 本部分介绍Flink 的基础部分，比如术语、架构、编程模型、编程指南、基本的…
阅读更多...
dpdk原理概述及核心源码剖析

dpdk原理概述及核心源码剖析

dpdk原理 1、操作系统、计算机网络诞生已经几十年了，部分功能不再能满足现在的业务需求。如果对操作系统做更改，成本非常高，所以部分问题是在应用层想办法解决的，比如前面介绍的协程、quic等，都是在应用层重新开发的框…
阅读更多...
docker 私有仓库

docker 私有仓库

Docker 私有仓库 一、私有仓库搭建 # 1、拉取私有仓库镜像 docker pull registry # 2、启动私有仓库容器 docker run -id --nameregistry -p 5000:5000 registry # 3、打开浏览器 输入地址http://私有仓库服务器ip:5000/v2/_catalog，看到{"repositories&quo…
阅读更多...
Linux操作系统——进程(三) 进程优先级

Linux操作系统——进程(三) 进程优先级

进程优先级 首先呢，我们知道一个进程呢（或者也可以叫做一个任务），它呢有时候要在CPU的运行队列中排队，要么有时候阻塞的时候呢又要在设备的等待队列中排队，其实我们排队的本质就是：确认优先级。…
阅读更多...
【数据结构】什么是二叉树?

【数据结构】什么是二叉树?

🦄个人主页:修修修也 🎏所属专栏:数据结构 ⚙️操作环境:Visual Studio 2022 目录 📌二叉树的定义 📌二叉树的特点 📌特殊二叉树 📌二叉树的性质 📌二叉树的存储结构 📌二叉树…
阅读更多...
大语言模型说明书

大语言模型说明书

在浩瀚的信息宇宙中，大语言模型如同一颗璀璨的星星正在熠熠生辉。21世纪以来，人工智能可谓是飞速发展，从简单的神经网络到大语言模型、生成式AI，这并非仅仅是一种技术的进步，更是人类智慧的飞跃。大语言模型不仅仅是语…
阅读更多...
opencv入门到精通——形态学转换

opencv入门到精通——形态学转换

目录 目标 理论 1. 侵蚀 2. 扩张 3. 开运算 4. 闭运算 5. 形态学梯度 6. 顶帽 7. 黑帽 结构元素 目标 在这一章当中， 我们将学习不同的形态学操作，例如侵蚀，膨胀，开运算，闭运算等。我们将看到不同的功能&…
阅读更多...
期末复习【计算机图像处理】

期末复习【计算机图像处理】

期末复习【计算机图像处理】 前言版权推荐期末复习期末考点内容期末考试题型一、填空 2分*10 概念二、简答 5分*2三、计算 10分*6四、绘图 10分*1 测验二测验三最后 前言 2023-12-25 15:09:07 昨天没有睡好 0点~3点看B站 Google模拟面试 3点~5点没睡着 5~6点睡了一会 6~12点终…
阅读更多...
isp代理/双isp代理/数据中心代理的区别？如何选择？

isp代理/双isp代理/数据中心代理的区别？如何选择？

本文我们来详细科普一下几种不同的代理类型：isp代理/双isp代理/数据中心代理，了解他们的区别，选择更适合自己的代理类型。 在讲述这几种代理类型之前，我们先复习一下代理大类有哪几种。 一、机房代理和非机房代理 在做代理ip选…
阅读更多...
《试题与研究》期刊发表投稿方式

《试题与研究》期刊发表投稿方式

《试题与研究》杂志是面向全国公开发行的国家CN级权威教育期刊。创刊以来一直以服务教育服务学生为办刊宗旨，以优秀的内容质量和编校质量深受广大读者好评，其权威性、导向性、针对性、实用性在全国教育期刊中独树一帜。为推动教育科研事业的发展&#xf…
阅读更多...
最新文章

Copyright @ 2022~2023