四相移相键控(QPSK)调制及解调实验
四相移相键控(QPSK)调制及解调实验
通信对抗原理实验报告实验名称:四相移相键控（QPSK）调制及解调实验学生姓名：学生学号：学生班级:所学专业：实验日期：1. 实验目的1. 掌握QPSK 调制解调原理及特性.2。
. 熟悉Matlab 仿真软件的使用。
2. 实验内容1、 编写Matlab 程序仿真QPSK 调制及相干解调。
2、 观察IQ 两路基带信号的特征及与输入NRZ 码的关系。
3、 观察IQ 调制解调过程中各信号变化。
4、 观察功率谱的变化。
5、 分析仿真中观察的数据，撰写实验报告。
3. 实验原理1、QPSK 调制原理QPSK 又叫四相绝对相移调制,它是一种正交相移键控。
QPSK 利用载波的四种不同相位来表征数字信息。
由于每一种载波相位代表两个比特信息，因此，对于输入的二进制数字序列应该先进行分组，将每两个比特编为一组,然后用四种不同的载波相位来表征。
我们把组成双比特码元的前一信息比特用a 代表,后一信息比特用b 代表。
双比特码元中两个信息比特ab 通常是按格雷码排列的,它与载波相位的关系如表1-1所示，矢量关系如图1—1所示。
图1-1（a)表示A 方式时QPSK 信号矢量图,图1—1（b)表示B 方式时QPSK 信号的矢量图。
由于正弦和余弦的互补特性,对于载波相位的四种取值,在A 方式中：45°、135°、225°、315°，则数据k I、k Q 通过处理后输出的成形波形幅度有两种取值±2/2;B 方式中:0°、90°、180°、270°,则数据k I、k Q 通过处理后输出的成形波形幅度有三种取值±1、0。
表(0,1)(1,1)(0,0)参考相位参考相位(a)(b)图1-1 QPSK 信号的矢量图下面以A 方式的QPSK 为例说明QPSK 信号相位的合成方法。
串/并变换器将输入的二进制序列依次分为两个并行序列,然后通过基带成形得到的双极性序列(从D/A 转换器输出，幅度为±2/2）。
四相移相键控调制解调
太原理工大学现代科技学院实验报告一、 实验目的1、了解QPSK 调制解调原理及特性。
2、了解载波在QPSK 相干及非相干时的解调特性。
二、 实验内容1、观察I 、Q 两路基带信号的特征及与输入NRZ 码的关系。
2、观察IQ 调制解调过程中各信号变化。
3、观察解调载波相干时和非相干时各信号的区别。
三、 基本原理1、QPSK 调制原理 QPSK 又叫四相绝对相移调制，它是一种正交相移键控。
QPSK 利用载波的四种不同相位来表征数字信息。
由于每一种载波相位代表两个比特信息，因此，对于输入的二进制数字序列应该先进行分组，将每两个比特编为一组，然后用四种不同的载波相位来表征。
我们把组成双比特码元的前一信息比特用a 代表，后一信息比特用b 代表。
双比特码元中两个信息比特ab 通常是按格雷码排列的，它与载波相位的关系如表1-1所示，矢量关系如图1-1所示。
图1-1（a ）表示A 方式时QPSK 信号矢量图，图1-1（b ）表示B 方式时QPSK 信号的矢量图。
由于正弦和余弦的互补特性，对于载波相位的四种取值，在A 方式中：45°、135°、225°、315°，则数据k I 、k Q 通过处理后输出的成形波形幅度有三种取值±1、0。
表1-1 双比特码元与载波相位关系太原理工大学现代科技学院实验报告(0,1)(1,1)(0,0)参考相位参考相位(a)(b)图1-1 QPSK 信号的矢量图下面以A 方式的QPSK 为例说明QPSK 信号相位的合成方法。
串/并变换器将输入的二进制序列依次分为两个并行序列，然后通过基带成形得到的双极性序列（从D/A 转码元。
双极性的a 和b 脉冲通过两个平衡调制器分别对同相载波及正交载波进行二相调制，得到图1-2中虚线矢量，将两路输出叠加，即得到QPSK 调制信号，其相位编码关系如表1-2所示。
a(1)b(1)b(0)a(0)图1-2 矢量图表1-2 QPSK 信号相位编码逻辑关系用调相法产生QPSK 调制器框图如图1-3所示。
QPSK、DQPSK系统调制与解调
实验四QPSK与DQPSK调制实验一、实验目的在2PSK，2DPSK的学习基础上，掌握QPSK，以及以其为基础的DQPSK，OQPSK， /4—DQPSK等若干种相关的重要调制方式的原理，从而对多进制调相有一定了解。
二、实验设备1、“移动通信技术应用综合实训系统” 实验仪一台。
2、50MHz示波器一台。
3、实验模块：信源模块，QPSK-调制模块。
三、实验原理一）基本理论(A)四相绝对移相键控（QPSK）的调制四相绝对移相键控利用载波的四种不同相位来表征数字信息。
由于每一种载波相位代表两个比特信息，故每个四进制码元又被称为双比特码元。
我们把组成双比特码元的前一信息比特用a代表，后一信息比特用b代表。
双比特码元中两个信息比特ab通常是按格雷码（即反射码）排列的，它与载波相位的关系如表所列。
表4-1 双比特码元与载波相位的关系由于四相绝对移相调制可以看作两个正交的二相绝对移相调制的合成，故两者的功率谱密度分布规律相同。
下面我们来讨论QPSK 信号的产生与解调。
QPSK 信号的产生方法与2PSK 信号一样，也可以分为调相法和相位选择法。
(1) 调相法用调相法产生QPSK 信号的组成方框图如下所示。
图4-1 QPSK 信号的组成方框图设两个序列中的二进制数字分别为a 和b ，每一对ab 称为一个双比特码元。
并设经过串并变换后上支路为a,下支路为b 。
双极性的a 和b 脉冲通过两个平衡调制器分别对同相载波及正交载波进行二相调制。
a(1)(1,0)b(0)(0,0)a(0)b(1)(0,1)(1,1)表4-2 QPSK 信号相位编码逻辑关系(2)相位选择法用相位选择法产生QPSK信号的组成方框图如下所示。
图4-2 相位选择法产生QPSK信号方框图(B)四相相对移相键控（DQPSK）的调制所谓四相相对移相键控也是利用前后码元之间的相对相位变化来表示数字信息。
若以前一码元相位作为参考，并令△φ为本码元与前一码元的初相差。
四相移相键控调制解调
武夷学院实验报
课程名称：移动通信原理与技术
项目名称：四相移相键控（QPSK）调制及解调实验
姓名：陈真灼专业：通信工程班级： 3班学号：20114173004
1注：1、实验预习部分包括实验环境准备和实验所需知识点准备。
2、若是单人单组实验，同组成员填无。
c、对比观测解调前后的I路信号
示波器探头分别接IQ模块的“
解调正确，若不一致可能是载波相位不对，可按下
复位。
）
2注：实验过程记录要包含实验目的、实验原理、实验步骤，页码不够可自行添加。
六
实验报告成绩（百分制）__________ 实验指导教师签字：__________。
四相移相键控(QPSK)调制及解调实验
实验二四相移相键控（QPSK ）调制及解调实验一、 实验目的1、了解QPSK 调制解调原理及特性。
2、了解载波在QPSK 相干及非相干时的解调特性。
二、 实验内容1、观察I 、Q 两路基带信号的特征及与输入NRZ 码的关系。
2、观察IQ 调制解调过程中各信号变化。
3、观察解调载波相干时和非相干时各信号的区别。
三、 基本原理（说明：原理部分需简要介绍）1、QPSK 调制原理QPSK 的调制有两种产生方法相乘电路法和选择法。
相乘法：输入信号是二进制不归零的双极性码元，它通过“串并变换”电路变成了两路码元。
变成并行码元后，每个码元的持续时间是输入码元的两倍。
用两路正交载波去调制并行码元。
发射信号定义为：⎪⎩⎪⎨⎧≤≤-+=其他,00],4)12(2cos[/2)(b t T t i ft t E t S ππ其中，i ＝1，2，3，4；E 是发射信号的每个符号的能量，T 为符号的持续时间，载波频率f 等于nc/T ，nc 为固定整数选择法输入基带信号经过串并变换后用于控制一个相位选择电路，按照当时的输入双比特ab ，决定选择哪个相位的载波输出2、QPSK 解调原理QPSK 接收机由一对共输入地相关器组成。
这两个相关器分别提供本地产生地相干参考信号()t 1φ和()t 2φ。
四、实验步骤（说明：要详细）（1）QPSK 调制程序close all% x1是类似[1 1 -1 -1 -1 -1 1 1]的分布,作用是控制相位的180°反转。
%由于仿真中载波的频率是f=1Hz，所以1s的间隔内有一个完整周期的正弦波。
t=[-1:0.01:7-0.01]; % t共800个数据，-1~7st1=[0:0.01:8-0.01]; %t1也是800个数据点，0 ~8stt=length(t); % tt=800x1=ones(1,800);for i=1:ttif (t(i)>=-1 & t(i)<=1) | (t(i)>=5& t(i)<=7);x1(i)=1;else x1(i)=-1;endendt2 = 0:0.01:7-0.01; %t2是700个数据点，是QPSK_rc绘图的下标t3 = -1:0.01:7.1-0.01; %t3有810个数据点，是i_rc的时间变量t4 = 0:0.01:8.1-0.01; %t4有810个数据点，是q_rc的时间变量tt1=length(t1);x2=ones(1,800); %x2是类似于[1 1 -1 -1 1 1 1 1]的分布,作用是控制相位的180°反转for i=1:tt1if (t1(i)>=0 & t1(i)<=2) | (t1(i)>=4& t1(i)<=8);x2(i)=1;else x2(i)=-1;endendf=0:0.1:1;xrc=0.5+0.5*cos(pi*f); %xrc是一个低通特性的传输函数y1=conv(x1,xrc)/5.5; %y1和x1 实际上没什么区别，仅仅是上升沿、下降沿有点过渡带y2=conv(x2,xrc)/5.5; % y2和x2 实际上没什么区别，仅仅是上升沿、下降沿有点过渡带n0=randn(size(t2));f1=1;i=x1.*cos(2*pi*f1*t); % x1就是I dataq=x2.*sin(2*pi*f1*t1); %x2就是Q dataI=i(101:800);Q=q(1:700);QPSK=sqrt(1/2).*I+sqrt(1/2).*Q;QPSK_n=(sqrt(1/2).*I+sqrt(1/2).*Q)+n0;n1=randn(size(t2));i_rc=y1.*cos(2*pi*f1*t3); % y1就是I data，i_rc可能是贴近实际的波形,i则是理想波形q_rc=y2.*sin(2*pi*f1*t4); %y2就是Q data，q_rc可能是贴近实际的波形,q则是理想波形I_rc=i_rc(101:800);Q_rc=q_rc(1:700);QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);QPSK_rc_n1=QPSK_rc+n1;subplot(3,1,1);plot(t3,i_rc);axis([-1 8 -1 1]);ylabel('a序列');subplot(3,1,2);plot(t4,q_rc);axis([-1 8 -1 1]);ylabel('b序列');subplot(3,1,3);plot(t2,QPSK_rc);axis([-1 8 -1 1]);ylabel('合成序列');（2）QPSK解调程序clear allclose allbit_in = randint(1e3, 1, [0 1]);bit_I = bit_in(1:2:1e3); %bit_I为”奇数序列”，奇数序列是同相分量，以cos为载波bit_Q = bit_in(2:2:1e3); %bit_Q是bit_in的所有偶数下标组成的”偶数序列”，以sin为载波data_I = -2*bit_I+1; % 将bit_I中的1变成-1，0变成1; 注意data_I是500点data_Q = -2*bit_Q+1; %将bit_Q中的1变成-1，0变成1data_I1=repmat(data_I',20,1); %将500行的列向量data_I的共轭转置data_I’复制为20*500的矩阵，20行数据是相同的。
实验二__交错四相移相键控(OQPSK)调制及解调实验
实验二 交错四相移相键控（OQPSK ）调制及解调实验一、实验目的1、了解OQPSK 调制解调原理及特性2、了解载波在QPSK 相干及非相干时的解调特性3、与QPSK 调制对比，掌握它们的差别二、实验内容1、观察I 、Q 两路基带信号的特征及与输入NRZ 码的关系。
2、观察IQ 调制解调过程中各信号变化。
3、观察QPSK 调制及OQPSK 调制各信号的区别。
4、观察解调载波相干时和非相干时各信号的区别。
三、基本原理OQPSK 又叫偏移四相相移键控，它是基于QPSK 的一类改进型，为了克服QPSK 中过零点的相位跃变特性，以及由此带来的幅度起伏不恒定和频带的展宽（通过带限系统后）等一系列问题。
若将QPSK 中并行的I ，Q 两路码元错开时间（如半个码元），称这类QPSK 为偏移QPSK 或OQPSK 。
通过I 、Q 两路码元错开半个码元调制之后的波形，其载波相位跃变由180°降至90°，避免了过零点，从而大大降低了峰平比和频带的展宽。
下面通过一个具体的例子说明某个带宽波形序列的I 路，Q 路波形，以及经载波调制以后相位变化情况。
若给定基带信号序列为`1 －1 －1 1 1 1 1 －1 －1 1 1 －1对应的QPSK 与OQPSK 发送波形如图2-1所示。
1-1-11111-1-111-1111-11-111-11-1-111-11-1基带波形I 信道QPSK,OQPSKQ 信道QPSK Q 信道OQPSK-1图2-1 QPSK,OQPSK 发送信号波形图2-1中，I 信道为奇数数据单元，Q 信道为偶数数据单元，而OQPSK 的Q 信道与其I 信道错开（延时）半个码元。
QPSK ，OQPSK 载波相位变化公式为{}()33arctan ,,,()4444j i ji Q t I t ππϕππ⎡⎤⎛⎫⎧⎫=--⎢⎥⎨⎬ ⎪⎩⎭⎢⎥⎝⎭⎣⎦ (2-1)QPSK 数据码元对应的相位变化如图2-2所示，OQPSK 数据码元对应相位变化如图2-3所示1)图2-2 QPSK 相位变化图图2-3 OQPSK 相位变化图 对于QPSK 数据码元对的相位变换由图2-1和2-2求得为：(1,-1)(1,1)(1,-1)(1,-1)(-1,1)()4-()34()4()4-()34()4-2π-2π-πππ码元对相位及相位变化：可见，在QPSK中存在过零点的180°跃变。
QPSK调制解调实验
实验一QPSK 调制实验一、实验目的1、掌握QPSK 的调制解调原理。
2、掌握QPSK 的软件仿真方法。
3、掌握QPSK 的硬件设计方法。
二、预习要求1、掌握QPSK 的编解码原理和方法。
2、熟悉matlab 的应用和仿真方法。
3、熟悉DSP 和FPGA 的开发方法。
三、实验原理1、QPSK 调制的工作原理多相相移键控（MPSK ），特别是四相相移键控（QPSK ）是目前移动通信、微波通信和卫星通信中最常用的载波传输方式。
四相相移键控（QPSK ）信号的正弦载波有4个可能的离散相位状态，每个载波相位携带2个二进制符号，其信号表达式为：)cos()(i c i t A t S θω+= i ＝1，2，3，4 0≤t ≤TsTs 为四进制符号间隔，{i θ:i=1,2,3,4}为正弦波载波的相位，有四种可能状态。
如以下矢量图所示：如图为QPSK 的相位图，QPSK 的相位为（－3π/4，－π/4，π/4，3π/4）。
对于QPSK ：)sin cos cos (sin )sin()(i c i c i c i t t A t A t S θωθωθω+=+= 0≤t ≤Ts由于21cos ±=i θ 21s i n ±=i θ所以：)cos )(sin )((2)(t t Q t t I A t S c c i ωω+=21cos )(±==i t I θ21s i n )(±==i t Q θQPSK 正交调制器方框图如图所示：I图QPSK 正交调制器方框图在kTs ≤t ≤(k+1) Ts(Ts=2Tb)的区间，QPSK 产生器的输出为：⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧--=-+-=--+=+++=+=----11),43cos(11),4cos(11),43cos(11),4cos()(1111n n c n n c n n c n n c a a t A a a t A a a t A a a t A t s πωπωπωπω2、QPSK 的相干解调的基本工作原理 QPSK 的相干解调方框图如图所示：图QPSK 的相干解调方框图当调制信号为I ＝1，Q ＝1时，由调制原理，调制输出信号为t t t S c c i ωωcos sin )(+=，在没有噪声和延时的理想状态时，解调器的输入t t t S t r c c i ωωcos sin )()(+==，则I 检测器的输出为：t t t t t t r c c c c c ωωωωωsin cos sin sin sin )(+=t t t t c c c c ωωωω2sin 212cos 21212sin 21)2cos 1(21+-=+-=则Q 检测器的输出为：t t t t t t r c c c c c ωωωωωcos cos cos sin cos )(+=t t t t c c c c ωωωω2sin 212cos 21212sin 21)2cos 1(21++=++=用截止频率小于2c ω的低通滤波器对I 检测器的输出滤波后得到1/2,即为逻辑1；对Q 检测器的输出滤波后得到1/2,即为逻辑1。
0QPSK调制解调实验报告
0QPSK调制解调实验报告一、实验目的1．掌握0QPSK调制解调原理。
2．理解0QPSK的优缺点。
二、实验内容1．观察0QPSK调制过程各信号波形。
2．观察0QPSK解调过程各信号波形。
三、预备知识1．0QPSK调制解调的基本原理。
2. 0QPSK调制解调模块的工作原理及电路说明。
四、实验器材1. 移动通信原理实验箱。
2．20M数字双踪示波器。
五、实验原理0QPSK调制解调原理0QPSK又叫四相相移键控，它通QPSK的不同之处是在正交支路引入了一个码元（Ts）的延时，这使得两个支路的数据错开了一个码元时间，不会同时发生变化，而不像QPSK那样产生±π的相位跳变，而仅能产生±π/2的相位跳变，如图4-1所示。
从图4-1星座图和相位转移图中看出对于1QPSK，±π相位的跳变消除了，所以1QPSK 信号的带限不会导致信号包络经过零点。
0QPSK包络的变化小多了，因此对1QPSK的硬限幅或非线性放大不会再产生严重的频带扩展，0QPSK即使在非线性放大后仍能保持其带限的性质。
0QPSK的调制方法和QPSK一样。
图4-1 QPSK和0QPSK的星座图和相位转移图1)六、实验步骤1.A 方式的0QPSK 调制实验（1）将“调制类型选择”拨码开关拨为00001000、0001，则调制类型选择为A 方式的0QPSK 调制。
（2）分别观察并说明NRZ 码经串并转换得到的‘DI ’、‘DQ ’两路的一个周期的数据波形。
CH1:NRZCH2:DI CH1:NRZ CH2:DQ（3）双踪观察并分析说明‘DI ’与‘I 路成形’信号波形；‘DQ ’与‘Q 路成形’信号波形；CH1:DI CH2: I路成形CH1:DQ CH2: Q路成形（4）双踪观察并分析说明‘I路成形’信号波形与‘I路调制’同相调制信号波形；‘Q路成形’信号与‘Q路调制’正交调制信号波形。
CH1: I路成形CH2: I路调制CH1: Q路成形CH2: Q路调制（5）用示波器观察并说明‘I路成形’信号与‘Q路成形信号的X-Y波形。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
通信对抗原理实验报告实验名称：四相移相键控（QPSK）调制及解调实验学生姓名：学生学号：学生班级:所学专业：实验日期:1. 实验目的1. 掌握QPSK 调制解调原理及特性。
2.. 熟悉Matlab 仿真软件的使用。
2. 实验内容1、 编写Matlab 程序仿真QPSK 调制及相干解调。
2、 观察IQ 两路基带信号的特征及与输入NRZ 码的关系。
3、 观察IQ 调制解调过程中各信号变化。
4、 观察功率谱的变化。
5、 分析仿真中观察的数据，撰写实验报告。
3. 实验原理1、QPSK 调制原理QPSK 又叫四相绝对相移调制，它是一种正交相移键控。
QPSK 利用载波的四种不同相位来表征数字信息。
由于每一种载波相位代表两个比特信息，因此，对于输入的二进制数字序列应该先进行分组，将每两个比特编为一组，然后用四种不同的载波相位来表征。
我们把组成双比特码元的前一信息比特用a 代表，后一信息比特用b 代表。
双比特码元中两个信息比特ab 通常是按格雷码排列的，它与载波相位的关系如表1-1所示，矢量关系如图1-1所示。
图1-1（a ）表示A 方式时QPSK 信号矢量图，图1-1（b ）表示B 方式时QPSK 信号的矢量图。
由于正弦和余弦的互补特性，对于载波相位的四种取值，在A 方式中：45°、135°、225°、315°，则数据、通过处理后输出的成形波形幅度有两种取值±；B 方式中：0°、90°、180°、270°，则数据、通过处理后输出的成形波形幅度有三种取值±1、0。
表1-1 双比特码元与载波相位关系k IkQ 2/2k IkQ图1-1 QPSK 信号的矢量图下面以A 方式的QPSK 为例说明QPSK 信号相位的合成方法。
串/并变换器将输入的二进制序列依次分为两个并行序列，然后通过基带成形得到的双极性序列（从D/A转换器输出，幅度为±）。
设两个双极性序列中的二进制数字分别为a 和b ，每一对ab 称为一个双比特码元。
双极性的a 和b 脉冲通过两个平衡调制器分别对同相载波及正交载波进行二相调制，得到图1-2中虚线矢量，将两路输出叠加，即得到QPSK 调制信号，其相位编码关系如表1-2所示。
图1-2 矢量图表(0,1)(1,1)(0,0)参考相位参考相位(a)(b)2/2a(1)b(1)b(0)a(0)用调相法产生QPSK 调制器框图如图1-3所示。
图1-3 QPSK 调制器框图图1-4 二进制码经串并变换后码型由图1-3可以看到，QPSK 的调制器可以看作是由两个BPSK 调制器构成，输入的串行二进制信息序列经过串并变换，变成两路速率减半的序列，电平发生器分别产生双极性的二电平信号I （t ）和Q （t ），然后对和进行调制，相加后即可得到QPSK 信号。
经过串并变换后形成的两个支路如图1-4所示，一路为单数码元，另外一路为偶数码元，这两个支路互为正交，一个称为同相支路，即I 支路；另外一路称为正交支路，即Q 支路。
2、QPSK 相干解调原理QPSK 由于QPSK 可以看作是两个正交2PSK 信号的合成，故它可以采用与2PSK 信号类似的解调方法进行解调，即由两个2PSK 信号相干解调器构成，其原理框图如图1-5所示。
01010011011100cos A t ωsin A t ω图1-5 QPSK解调原理框图3、星座图星座显示是示波器显示的数字等价形式，将正交基带信号的I和Q两路分别接入示波器的两个输入通道，通过示波器的“X-Y”的功能即可以很清晰地看到调制信号的星座图。
我们知道QPSK信号可以用正交调制方法产生。
在它的星座图中，四个信号点之间任何过渡都是可能的，如图1-6（a）所示。
OQPSK信号将正交路信号偏移T2 /2，结果是消除了已调信号中突然相移180度的现象，每隔T2 /2信号相位只可能发生±90度的变化。
因而星座图中信号点只能沿正方形四边移动，如图1-6（b）所示。
MSK信号配置图如图1-6(c)所示，1比特区间仅使用圆周的1/4，信号点必是轴上4个点中任何一个，因此，相位必然连续。
（a）QPSK （b）OQPSK （c）MSK图1-6 相位转移图4.实验代码clear all;j=sqrt(-1);data2=randint(1,100);data2_out=zeros(1,100);data2_1=zeros(1,50);data2_2=zeros(1,50);Ia=zeros(1,200);Qa=zeros(1,200);error_rate=zeros(1,50);mi=;mq=;i=1;while i<101if mod(i,2)==1;data2_1(1,i/2+=data2(1,i);elsedata2_2(1,i/2)=data2(1,i);endi=i+1;enddata24=data2_2+data2_1*2;i=1;while i<51if data24(1,i)==0Ia(1,4*i)=1;Ia(1,4*i-1)=1;Ia(1,4*i-2)=1;Ia(1,4*i-3)=1;Qa(1,4*i-0)=0;Qa(1,4*i-1)=0;Qa(1,4*i-2)=0;Qa(1,4*i-3)=0;elseif data24(1,i)==1Ia(1,4*i)=0;Ia(1,4*i-1)=0;Ia(1,4*i-2)=0;Ia(1,4*i-3)=0;Qa(1,4*i)=-1;Qa(1,4*i-1)=-1;Qa(1,4*i-2)=-1;Qa(1,4*i-3)=-1;elseif data24(1,i)==2 Ia(1,4*i)=0;Ia(1,4*i-1)=0;Ia(1,4*i-2)=0;Ia(1,4*i-3)=0;Qa(1,4*i)=1;Qa(1,4*i-1)=1;Qa(1,4*i-2)=1;Qa(1,4*i-3)=1;elseif data24(1,i)==3 Ia(1,4*i)=-1;Ia(1,4*i-1)=-1; Ia(1,4*i-2)=-1; Ia(1,4*i-3)=-1; Qa(1,4*i)=0;Qa(1,4*i-1)=0;Qa(1,4*i-2)=0;Qa(1,4*i-3)=0;endendendendi=i+1;endFd=1;Fs=4;[num,den] = rcosine(Fd,Fs,'fir',,3); B=Ia;Ia=conv(B,num);Ia=Ia(1,12:211);C=Qa;Qa=conv(C,num);Qa=Qa(1,12:211);subplot(5,1,1);plot(1:200,Ia); subplot(5,1,2);plot(1:200,Qa);fc=177;cos_c=zeros(1,200);sin_c=zeros(1,200);i=1;while i<201cos_c(1,i)=cos(2*pi*i*fc/4);sin_c(1,i)=sin(2*pi*i*fc/4);i=i+1;enddata_in=Ia.*cos_c+Qa.*sin_c;for x=0:1:49error_m=0;for y=1:100data_out=awgn(data_in,x/5);data_out_c=data_out.*cos_c;data_out_s=data_out.*sin_c;doc_f=fft(data_out_c);la_f=fft(Ia);doc_f_mag=abs(doc_f);doc_f_ang=angle(doc_f);for i=51:1:149doc_f_mag(1,i)=doc_f_mag(1,i)/3; endI_j=doc_f_mag.*exp(j*doc_f_ang);Ia_j1=ifft(I_j);Ia_j1=real(Ia_j1);for i=1:50sum=0;for n=1:4;sum=sum+Ia_j1(1,4*(i-1)+n);endsum=sum/4;if sum>=miF(1,4*i-3)=1;F(1,4*i-2)=1;F(1,4*i-1)=1;F(1,4*i)=1;elseif sum<=-miF(1,4*i-3)=-1;F(1,4*i-2)=-1;F(1,4*i-1)=-1;F(1,4*i)=-1;elseF(1,4*i-3)=0;F(1,4*i-2)=0;F(1,4*i-1)=0;F(1,4*i)=0;endendendIa_j=F;dos_f=fft(data_out_s);Qa_f=fft(Qa);dos_f_mag=abs(dos_f);dos_f_ang=angle(dos_f);for i=51:1:149dos_f_mag(1,i)=dos_f_mag(1,i)/3;endQ_j=dos_f_mag.*exp(j*dos_f_ang);Qa_j1=ifft(Q_j);Qa_j1=2*real(Qa_j1);for i=1:50sum=0;for n=1:4;sum=sum+Qa_j1(1,4*(i-1)+n);endsum=sum/4;if sum>=mqG(1,4*i-3)=1;G(1,4*i-2)=1;G(1,4*i-1)=1;G(1,4*i)=1;elseif sum<=-mqG(1,4*i-3)=-1;G(1,4*i-2)=-1;G(1,4*i-1)=-1;G(1,4*i)=-1;elseG(1,4*i-3)=0;G(1,4*i-2)=0;G(1,4*i-1)=0;G(1,4*i)=0;endendendQa_j=G;for i=1:50if Ia_j(1,4*i-2)==1&&Qa_j(1,4*i-2)==0data2_out(1,2*i-1)=0;data2_out(1,2*i)=0;elseif Ia_j(1,4*i-2)==0&&Qa_j(1,4*i-2)==1 data2_out(1,2*i-1)=1;data2_out(1,2*i)=0;elseif Ia_j(1,4*i-2)==-1&&Qa_j(1,4*i-2)==0data2_out(1,2*i-1)=1;data2_out(1,2*i)=1;elseif Ia_j(1,4*i-2)==0&&Qa_j(1,4*i-2)==-1 data2_out(1,2*i-1)=0;data2_out(1,2*i)=1;endendendendendfor i=1:100if data2(1,i)~=data2_out(1,i);error_m=error_m+1;endendenderror_rate(1,x+1)=error_m/10000;endsubplot(5,1,3);plot(1:200,Ia_j);subplot(5,1,4);plot(1:200,Qa_j);subplot(5,1,5);semilogy(0::,error_rate);5.实验结果I,Q,误码率曲线图。