通信原理 简易QPSK AB方法的调制解调及星座图的matlab实现

已于 2022-06-18 21:38:28 修改
阅读量6.1k 收藏 105
点赞数 8
文章标签： matlab
于 2022-03-02 20:34:37 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_45471796/article/details/123238496
版权

QPSK调制解调方框图如下所示

A方式 框图如下所示：

A方式代码如下所示

clc
clear

Nm = 100;%100个码元
snrdB = 20;%信噪比20dB
Ts = 2;%原始码元间隔
N_sample = 1000;
dt = Ts/N_sample;%采样间隔

Msg = randi(4,1,Nm);
Msg_q = Msg-1;%转换为四进制,为产生的信号
sig_Msg_q = [];
for ik = 1:Nm
    sig_Msg_q = [sig_Msg_q,Msg_q(ik)*ones(1,2*N_sample)];
end
t =0:dt/2:Nm*Ts-dt/2;
figure(1);
plot(t,sig_Msg_q);axis([0,10,-1,4]);%四进制消息码电平
Msg_bin = [];
for ik = 1:Nm
    Msg_bin = [Msg_bin,fix(Msg_q(ik)/2)];
    Msg_bin = [Msg_bin,Msg_q(ik)-2*fix(Msg_q(ik)/2)];
end
Msg_bin_a = [];
Msg_bin_b = [];
for ik = 1:Nm
    Msg_bin_a = [Msg_bin_a,Msg_bin(2*ik-1)];
    Msg_bin_b = [Msg_bin_b,Msg_bin(2*ik)];
end%分成a,b两列
sig_Msg_bin = [];
for ik = 1:Nm
    sig_Msg_bin = [sig_Msg_bin,Msg_bin(2*ik-1)*ones(1,N_sample)];
    sig_Msg_bin = [sig_Msg_bin,Msg_bin(2*ik)*ones(1,N_sample)];
end
figure(2);
plot(t,sig_Msg_bin);axis([0,10,-1,2]);

wavNum = 20;%wavNum单码元三角函数波周期数
Ph = 0:2*wavNum*pi/N_sample:2*pi*wavNum*(Ts-1/N_sample);%Ph为相位
gt11 = cos(Ph.*ones(1, Ts*N_sample)+pi/2);
gt10 = cos(Ph.*ones(1, Ts*N_sample));
gt01 = cos(Ph.*ones(1, Ts*N_sample)+pi);
gt00 = cos(Ph.*ones(1, Ts*N_sample)+pi*3/2);
ph_out_sig= [3/2*pi,pi,0,pi/2];
sc_sig_out = [];
for ik = 1:50
    sc_sig_out =[sc_sig_out,cos(ph_out_sig(Msg(ik)))+i*sin(ph_out_sig(Msg(ik)))];
end
sc_out = sc_sig_out;

sig_out = [];
for ik = 1:Nm
    if Msg_bin_a(ik)==1
        if Msg_bin_b(ik)==1
            sig_out = [sig_out,gt11];
        else
            sig_out = [sig_out,gt10];
        end
    else
        if Msg_bin_b(ik)==1  
            sig_out = [sig_out,gt01];
        else
            sig_out = [sig_out,gt00];
        end
    end
end

figure(3);plot(t,sig_out);axis([0,1,-1.5,1.5]);%总输出波形

sig_channel = awgn(sig_out,snrdB);
figure(4);plot(t,sig_channel);axis([0,1,-1.5,1.5]);%添加噪声之后的信号，噪声单位SNRdB

wp=[wavNum-5 wavNum+5]*2/N_sample;
ws=[wavNum-10 wavNum+10]*2/N_sample;
rp=3;
rs=30;
[N,wn]=cheb2ord(wp,ws,rp,rs);
[b,a]=cheby2(N,rp,wn,'bandpass');%切比雪夫带通滤波
sig_bandpass = filter(b,a,sig_channel);
figure(5);plot(t,sig_bandpass);axis([0,1,-1.5,1.5]);%经过带通滤波器之后的信号

Ph1 = 0:2*wavNum*pi/N_sample:2*pi*wavNum*(Nm*Ts-1/N_sample);
sig_local1 = cos(Ph1.*ones(1, Nm*Ts*N_sample));

sig_local2 = -1*sin(Ph1.*ones(1, Nm*Ts*N_sample));

sig_mul_a = sig_local1.*sig_bandpass;
sig_mul_b = sig_local2.*sig_bandpass;
figure(6);
subplot(2,1,1);plot(t,sig_mul_a);axis([0,0.1,-1.5,1.5]);
subplot(2,1,2);plot(t,sig_mul_b);axis([0,0.1,-1.5,1.5]);

wp=(0.5)*2/N_sample; 
ws=(1.5)*2/N_sample; 
Rp=3;Rs=30; 
[N,wn]=cheb2ord(wp,ws,rp,rs);
[b,a]=cheby2(N,rp,wn,'low');
sig_lp_a = filter(b,a,sig_mul_a);
sig_lp_b = filter(b,a,sig_mul_b);
figure(7);
subplot(2,1,1);plot(t,sig_lp_a);axis([0,1,-1.5,1.5]);
subplot(2,1,2);plot(t,sig_lp_b);axis([0,1,-1.5,1.5]);%切比雪夫低通滤波

t1 = Ts*N_sample/2:Ts*N_sample:2*Nm*N_sample-Ts*N_sample/2;
Msg_out_a = [];
Msg_out_b = [];
lvl_rec_a = [];
lvl_rec_b = [];
for ik = 1:length(t1)
    Msg_out_a = [Msg_out_a,sign(sig_lp_a(t1(ik)))];
    Msg_out_b = [Msg_out_b,sign(sig_lp_b(t1(ik)))];
    lvl_rec_a = [lvl_rec_a,sig_lp_a(t1(ik))];
    lvl_rec_b = [lvl_rec_b,sig_lp_b(t1(ik))];
end
sc_rec = lvl_rec_a+lvl_rec_b*i;
Msg_rec_bin = [];
for ik = 1:Nm
    Msg_rec_bin = [Msg_rec_bin,Msg_out_a(ik)];
    Msg_rec_bin = [Msg_rec_bin,Msg_out_b(ik)];
end
Msg_rec_q = (Msg_out_a+1)/2*2+(Msg_out_b+1)/2;

sig_Msg_rec_bin = [];
sig_Msg_rec_q = [];
for ik = 1:Nm
    sig_Msg_rec_bin = [sig_Msg_rec_bin,Msg_rec_bin(2*ik-1)*ones(1,N_sample)];
    sig_Msg_rec_bin = [sig_Msg_rec_bin,Msg_rec_bin(2*ik)*ones(1,N_sample)];
    sig_Msg_rec_q = [sig_Msg_rec_q,Msg_rec_q(ik)*ones(1,2*N_sample)];
end
figure(8);
subplot(211);plot(t,sig_Msg_rec_bin);axis([0,10,-2,2]);
subplot(212);plot(t,sig_Msg_rec_q);axis([0,10,-1,4]);
       
scatterplot(sc_out);
scatterplot(sc_rec);

count = 0;
for i = 1:Nm
    if Msg_q(i) == Msg_rec_q(i)
        count = count + 1;
    end
end
Pe = (Nm-count)/Nm;

运行结果如下所示

 

 

B方式框图如下所示：

B方式代码如下所示：

clc
clear

Nm = 1000;
snrdB = -10;%信噪比20dB
Ts = 2;%原始码元间隔
N_sample = 1000;
dt = Ts/N_sample;%采样间隔
Msg = randi(4,1,Nm);
Msg_q = Msg-1;%转换为四进制,为产生的信号
sig_Msg_q = [];
for ik = 1:Nm
    sig_Msg_q = [sig_Msg_q,Msg_q(ik)*ones(1,2*N_sample)];
end
t =0:dt/2:Nm*Ts-dt/2;
figure(1);
plot(t,sig_Msg_q);axis([0,10,-1,4]);%四进制消息码电平
Msg_bin = [];
for ik = 1:Nm
    Msg_bin = [Msg_bin,fix(Msg_q(ik)/2)];
    Msg_bin = [Msg_bin,Msg_q(ik)-2*fix(Msg_q(ik)/2)];
end
Msg_bin_a = [];
Msg_bin_b = [];
for ik = 1:Nm
    Msg_bin_a = [Msg_bin_a,Msg_bin(2*ik-1)];
    Msg_bin_b = [Msg_bin_b,Msg_bin(2*ik)];
end%分成a,b两列
sig_Msg_bin = [];
sig_Msg_bin_a = [];
sig_Msg_bin_b = [];
for ik = 1:Nm
    sig_Msg_bin = [sig_Msg_bin,Msg_bin(2*ik-1)*ones(1,N_sample)];
    sig_Msg_bin = [sig_Msg_bin,Msg_bin(2*ik)*ones(1,N_sample)];
    sig_Msg_bin_a = [sig_Msg_bin_a,Msg_bin_a(ik)*ones(1,2*N_sample)];
    sig_Msg_bin_b = [sig_Msg_bin_b,Msg_bin_b(ik)*ones(1,2*N_sample)];
end
figure(2);
subplot(311);plot(t,sig_Msg_bin);axis([0,10,-1,2]);title('二进制码');
subplot(312);plot(t,sig_Msg_bin_a);axis([0,10,-1,2]);title('二进制A通道码');
subplot(313);plot(t,sig_Msg_bin_b);axis([0,10,-1,2]);title('二进制B通道码');
    
wavNum = 20;%wavNum单码元三角函数波周期数
Ph = 0:2*wavNum*pi/N_sample:2*pi*wavNum*(Ts-1/N_sample);%Ph为相位
gta1 = cos(Ph.*ones(1, Ts*N_sample)); %a通道消息码为1
gta0 = -1*cos(Ph.*ones(1, Ts*N_sample));%a通道消息码为0
gtb1 = -1*sin(Ph.*ones(1, Ts*N_sample)); %b通道消息码为1
gtb0 = sin(Ph.*ones(1, Ts*N_sample));%b通道消息码为0
sc_out = (Msg_bin_a*2-1)+(Msg_bin_b*2-1)*i;

sig_a = [];
sig_b = [];
for ik = 1:Nm
    if Msg_bin_a(ik) == 1
        sig_a = [sig_a,gta1];
    else
        sig_a = [sig_a,gta0];
    end
    if Msg_bin_b(ik) == 1
        sig_b = [sig_b,gtb1];
    else
        sig_b = [sig_b,gtb0];
    end
end
sig_out = sig_a+sig_b;
figure(3);
subplot(311);plot(t,sig_out);axis([0,10,-1,2]);
subplot(312);plot(t,sig_a);axis([0,10,-1,2]);
subplot(313);plot(t,sig_b);axis([0,10,-1,2]);

sig_channel = pesnrf(sig_out,snrdB);

wp=[wavNum-0.5 wavNum+0.5]*2/N_sample;
ws=[wavNum-1.5 wavNum+1.5]*2/N_sample;
rp=3;
rs=30;
[N,wn]=cheb2ord(wp,ws,rp,rs);
[b,a]=cheby2(N,rp,wn,'bandpass');%切比雪夫带通滤波
sig_bandpass = filter(b,a,sig_channel);
figure(4);
subplot(311);plot(t,sig_out);axis([0,10,-1,2]);title('输出');
subplot(312);plot(t,sig_channel);axis([0,10,-1,2]);title('信道输出');
subplot(313);plot(t,sig_bandpass);axis([0,10,-1,2]);title('带通滤波器输出');

Ph1 = 0:2*wavNum*pi/N_sample:2*pi*wavNum*(Nm*Ts-1/N_sample);
sig_local1 = cos(Ph1.*ones(1, Nm*Ts*N_sample));
sig_local2 = -1*sin(Ph1.*ones(1, Nm*Ts*N_sample));

sig_mul_a = sig_local1.*sig_bandpass;
sig_mul_b = sig_local2.*sig_bandpass;
figure(5);
subplot(211);plot(t,sig_mul_a);axis([0,10,-2,2]);title('A通道相乘信号');
subplot(212);plot(t,sig_mul_b);axis([0,10,-2,2]);title('B通道相乘信号');

wp=(0.25)*2/N_sample; 
ws=(1)*2/N_sample; 
Rp=3;Rs=30; 
[N,wn]=cheb2ord(wp,ws,rp,rs);
[b,a]=cheby2(N,rp,wn,'low');
sig_lp_a = filter(b,a,sig_mul_a);
sig_lp_b = filter(b,a,sig_mul_b);
figure(6);
subplot(211);plot(t,sig_lp_a);axis([0,10,-2,2]);title('A通道低通滤波输出信号');
subplot(212);plot(t,sig_lp_b);axis([0,10,-2,2]);title('B通道低通滤波输出信号');


t1 = Ts*N_sample/2:Ts*N_sample:2*Nm*N_sample-Ts*N_sample/2;
Msg_out_a = [];
Msg_out_b = [];
lvl_rec_a = [];
lvl_rec_b = [];
for ik = 1:length(t1)
    Msg_out_a = [Msg_out_a,sign(sig_lp_a(t1(ik)))];
    Msg_out_b = [Msg_out_b,sign(sig_lp_b(t1(ik)))];
    lvl_rec_a = [lvl_rec_a,sig_lp_a(t1(ik))];
    lvl_rec_b = [lvl_rec_b,sig_lp_b(t1(ik))];
end
sc_rec = lvl_rec_a+lvl_rec_b*i;
Msg_rec_bin = [];
for ik = 1:Nm
    Msg_rec_bin = [Msg_rec_bin,Msg_out_a(ik)];
    Msg_rec_bin = [Msg_rec_bin,Msg_out_b(ik)];
end
Msg_rec_q = (Msg_out_a+1)/2*2+(Msg_out_b+1)/2;

sig_Msg_out_a = [];
sig_Msg_out_b = [];
sig_Msg_rec_bin = [];
sig_Msg_rec_q = [];
for ik = 1:Nm
    sig_Msg_out_a = [sig_Msg_out_a,Msg_out_a(ik)*ones(1,2*N_sample)];
    sig_Msg_out_b = [sig_Msg_out_b,Msg_out_b(ik)*ones(1,2*N_sample)];
    sig_Msg_rec_bin = [sig_Msg_rec_bin,Msg_rec_bin(2*ik-1)*ones(1,N_sample)];
    sig_Msg_rec_bin = [sig_Msg_rec_bin,Msg_rec_bin(2*ik)*ones(1,N_sample)];
    sig_Msg_rec_q = [sig_Msg_rec_q,Msg_rec_q(ik)*ones(1,2*N_sample)];
end
figure(7);
subplot(411);plot(t,sig_Msg_out_a);axis([0,10,-2,2]);title('A通道输出码元');
subplot(412);plot(t,sig_Msg_out_b);axis([0,10,-2,2]);title('B通道输出码元');
subplot(413);plot(t,sig_Msg_rec_bin);axis([0,10,-2,2]);title('接收到的二进制码元');
subplot(414);plot(t,sig_Msg_rec_q);axis([0,10,-1,4]);;title('接收到的四进制码元');
       
scatterplot(sc_out);title('理论星座图');
scatterplot(sc_rec);title('实际星座图');

count = 0;
for i = 1:Nm
    if Msg_q(i) == Msg_rec_q(i)
        count = count + 1;
    end
end
Pe = (Nm-count)/Nm;

运行结果如下所示：

 

 

 

出错时把pesnrf换成awgn即可

量大龙神魔
关注
  • 8
    点赞
  • 踩
  • 105
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
高码率QPSK调制解调方案（2.1 QPSK调制）
sdau2005的博客
02-03 1226
高码率QPSK调制解调方案，第2篇，2.1 QPSK调制，全文共分4部分，分成数篇文章发出。
2PSK,4PSK,8PSK,16PSK四个调制方式的误码率及星座图
平时工作较为繁忙,私信消息一般晚上回复,谢谢大家~~
01-19 6326
MPSK - multiple phase shift keying 多进制数字相位调制 , 又称多相制，是二相制的推广。我们知道PSK调制是将调制信号的信息附加在载波的相位上，例如载波的0相位表示调制信号的“1”电平，180°相位表示调制信号的“0”电平。如果将载波的相位进一步细分，分成m个不同的相位，一般m=2N，N是正整数。常见的有m=4，8，16，32，分别称为4PSK、8PSK等。这样，每个具有特定相位的载波就代表N比特的信息量。它是利用载波的多种不同相位状态来表征数字信息的调制方式。
QPSK星座图原理
热门推荐
王敏
04-19 2万+
文章目录QPSK调制星座图原理python如何实现星座图作图遇到的问题 QPSK调制星座图原理 1、QPSK将输入序列每两个bit分为一组（码元），前一个记做a，后一个记做b。这样就把输入序列分成了2个子序列。也就是说，将原串行序列转化成两组并行序列。对这两组序列中的每个bit进行极性转换，即一种电平转换，将0->-1，1->+1。这样原来的两组2进制01序列就转化成了+1，-1序列。...
【MATLAB源码-第125期】基于matlab的QPSK系统IQ调制方式仿真，输出每一个节点的波形。
2401_83582367的博客
06-14 730
这样，每个符号携带两个比特的信息，相比于传统的二进制相移键控（BPSK），QPSK在相同的带宽条件下能传输更多的数据。1. 偏移QPSK（OQPSK）：在OQPSK中，I和Q两个通道的数据是分别调制的，Q通道的数据相对于I通道有半个符号周期的偏移。理论上，QPSK的误码率性能与BPSK相当，但由于它在相同带宽下传输的数据量是BPSK的两倍，因此在实际应用中更有效率。在这些应用中，QPSK能够在保证较高数据传输速率的同时，有效抵抗信道中的干扰和噪声，提高传输的可靠性。这一步通常使用同相和正交混频器来完成。
基于matlab的qpsk与bpsk信号性能比较仿真,基于matlab的QPSK与BPSK信号性能比较仿真...
weixin_36032446的博客
03-24 3764
-欢迎下载-目 录第一章概述2第二章QPSK通信系统原理与仿真22.1 QPSK系统框图介绍22.2QPSK信号的调制原理32.2.1QPSK信号产生方法32.2.2QPSK星座图42.3QPSK解调原理及误码率分析42.3.1QPSK解调方法42.3.2QPSK系统误码率52.4QPSK信号在AWGN信道下仿真5第三章BPSK通信系统原理与仿真63.1BPSK信号的调制原理63.2BPSK...
程序阅读_全面详解LTE：MATLAB建模仿真与实现_自学笔记（1）调制与编码_程序阅读
王少靖的博客
03-06 3350
程序阅读_全面详解LTE：MATLAB建模仿真与实现_自学笔记（1）调制与编码_程序阅读 在粗浅地掌握了LTE知识后，从今天开始对《全面详解LTE：MATLAB建模仿真与实现》一书的学习。在此做个记录，供今后学习与参考。 开始之前，附上源程序的下载路径（源自《全面详解LTE：MATLAB建模仿真与实现》一书） https://pan.baidu.com/s/11oyJoH8pjUfZvPhD8M3E2A?pwd=1111 （含电子版教材） 写在开始：由于《全面详解LTE：MATLAB建模仿真与实现》一书已经
pi/4dqpsk的matlab及FPGA仿真
Adolfor的博客
06-21 5915
基于八相偏移调制的原理，本着资料稀少的原则，分享一下自己的经验。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % PIQPSKMODEM使用清单 % 卷积编码/编码的Ser/Ber~Snr曲线 % 仿真星座图，pi/4_QPSK信号的产生方法基带及调制； % 仿真差分解调 QPSK信号的过程及结果； % 符号速率 ＝1Mbps； % 基带成形滤波器滚降系数a=0.5； % 采样速率 ＝8Mbps； % 绘制Hilbert滤波器及普通带通滤波器的频率响应图 % 绘制 p
基于matlab的qpsk与bpsk信号性能比较仿真,基于matlab的QPSK与BPSK信号性能比较仿真.doc...
weixin_30788829的博客
03-24 1282
基于matlab的QPSK与BPSK信号性能比较仿真.doc2装订线目 录第一章概述2第二章QPSK通信系统原理与仿真22.1 QPSK系统框图介绍22.2QPSK信号的调制原理32.2.1QPSK信号产生方法32.2.2QPSK星座图42.3QPSK解调原理及误码率分析42.3.1QPSK解调方法42.3.2QPSK系统误码率52.4QPSK信号在AWGN信道下仿真5第三章BPSK通信系统原理与...
MPSK图像传输系统
weixin_45736693的博客
02-21 3116
（1）构建MPSK数字通信系统Clarke基带仿真模型， （2）分别在理想信道、AWGN信道、平坦Rayleigh衰落信道中，完成基于Matlab的MPSK基带调制与解调， （3）完成图像传输系统的分析与计算，绘制理论和仿真误码率曲线， （4）完成不同接收信噪比下图像信号的星座图演示， （5）设计完成基于Matlab的图像传输系统图形用户界面（GUI）。
通信笔记
qq_36222714的博客
07-04 1万+
1. OFDM OFDM 在FDM频分复用方法的基础上，采用多个正交的子载波来传输信号，极大限度地提升了系统的频谱利用率。 它将高速率数据流进行串并转换，分配到N个相互正交的子载波上，转换成并行的低速子数据流，所以调制符号的持续间隔远大于信道的时间扩散，从而有效地减少由于无线信道时间弥散所带来地ISI。 当信号带宽大于信道相干带宽时，就会产生频率选择性衰落，而OFDM可以将一个频率选择性衰落信道...
基于MATLAB的QPSK调制解调全过程仿真,输出基带,眼图,星座图,根升余弦滤波等+代码仿真操作视频
07-02
2.内容：基于MATLAB的QPSK调制解调全过程仿真,输出基带,眼图,星座图,根升余弦滤波等+代码仿真操作视频 3.用处：用于QPSK调制解调算法编程学习 4.指向人群：本硕博等教研学习使用 5.运行注意事项： 使用matlab...
通过matlab调用simulink进行QPSK调制解调仿真,并输出基带数据和星座图+matlab操作视频
06-09
2.内容：通过matlab调用simulink进行QPSK调制解调仿真,并输出基带数据和星座图+matlab操作视频 3.用处：用于QPSK调制解调算法编程学习 4.指向人群：本硕博等教研学习使用 5.运行注意事项： 使用matlab2021a或者...
qpsk调制解调与星座图映射
05-23
理解QPSK调制解调的工作原理和实现方法对于通信工程和信号处理领域的研究者以及从业者至关重要，它涉及到数字信号处理、通信理论和信息论等多个方面的知识。通过MATLAB这样的工具，我们可以对这些概念进行直观的可视...
基于Matlab的QPSK调制解调星座图
06-25
载波频率2.4GHZ 码元速率500MHz 含有加噪处理过程 代码可运行 无问题 学渣草草之作 作业用 有点水
MATLAB基础：函数与函数控制语句
Zjs15630471263的博客
07-26 1027
今天我们继续学习Matlab中函数相关知识。help 命令是最基本的查询方法，可查询所有目录、指定目录、命令、函数。我们直接点击帮助菜单即可查询所需的API函数。lookfor 关键字用于搜索相关的命令和函数。如，我们输入lookfor input，MATLAB就会显示所有关于input的API函数type M 文件名用于显示文件代码。如果想显示以.m扩展名结尾的就要加通配符 which M 文件名用于显示指定文件的路径。返回一个路径demo 用于演示matlab功能。如图，我们可以获取各种演示代码来解
MATLAB怎么实现多条曲线共用一个图例
最新发布
qq_37040743的博客
07-26 491
例子：如果我们有4条曲线，分布是cos类和sin类，我们的图例就想区分是cos类还是sin类。
Radon(拉当) 变换：超详细讲解（附MATLAB，Python 代码）
qq_39538718的博客
07-26 682
Radon 变换是数学上用于函数或图像的一种积分变换，广泛应用于图像处理领域，尤其是在计算机断层成像 (CT) 中。本文档将详细介绍 Radon 变换的数学含义及其在图像处理中的应用。
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值