2PSK数字信号的调制与解调-分享版
通信系统实训报告2psk的调制与解调
目录一.摘要和关键词 ..... 错误!未定义书签。
二.小组成员与分工 ... 错误!未定义书签。
三.设计的主要原理 ... 错误!未定义书签。
四．设计的系统仿真 .. 错误!未定义书签。
五．仿真系统的结论 .. 错误!未定义书签。
六.总结和体会： ..... 错误!未定义书签。
七.致谢 ............. 错误!未定义书签。
八.参考文献 ......... 错误!未定义书签。
2PSK的调制与解调一.摘要和关键词2PSK中文是：二进制相移键控，其有两种调制方法，模拟调制法和键控法，解调是用相干解调法。
我们这次做的是2PSK的调制与解调，在实现的过程中，使用了MATLAB的M文件的程序和SIMULINK 实现。
关键词：2PSK 调制解调 MATLAB二.小组成员与分工小组成员分工：确定题目：，查找资料：全部，设计程序： Simulink模拟图：；PPT，演讲：，演示：旁观：三.设计的主要原理二进制相移键控中，通常用相位0和π来分别表示“0”或“1”。
2PSK已调信号的时域表达式为e(t)=s(t)cosωt 。
因此，在某一个码元持续时间内观察时,有0，或π。
当码元宽度为载波周期的整数倍时，2PSK信号的典型波形如下图，2PSK信号的模拟调制法框图；如下图是产生2PSK信号的键控法框图，就模拟调制法而言，与产生2ASK信号的方法比较，只是对s(t)要求不同，因此2PSK信号可以看作是双极性基带信号作用下的DSB调幅信号。
而就键控法来说，用数字基带信号s(t)控制开关电路，选择不同相位的载波输出，这时s(t)为单极性NRZ 或双极性NRZ脉冲序列信号均可。
2PSK信号属于DSB信号,它的解调,不再能采用包络检测的方法,只能进行相干解调。
在这次通信系统仿真实训中，我们使用了MATLAB中的M文件实现，也使用了SIMULINK模块实现了2PSK的调制与解调。
而我负责的是SIMULINK的解调部分，Simulink是MATLAB中的一种可视化仿真工具，是一种基于MATLAB的框图设计环境，是实现动态系统建模、仿真和分析的一个软件包，被广泛应用于线性系统、非线性系统、数字控制及数字信号处理的建模和仿真中。
2PSK调制解调
• 这里，g(t)是脉宽为sT的单个矩形脉冲，
• 上式中而an的统计特性为：
• 即发送二进制符号“0”时（na取＋1），)(2tePSK取0相位；发送 二进制符号“1”时（na取－1），取π相位。这种以载波的不同 相位直接去表示相应二进制数字信号的调制方式，称为二进制绝 对相移方式。调制方法有模拟调制和键控法，解调方法通常采用 的是相干解调法。下面是2PSK的调制解调原理框图。
2PSK调制方法：
• 2PSK信号的解调一般采用相干解调，以下为原理图：
仿真结果：
2PSK调制解调
2PSK调制解调的基本原理
• 相移键控是利用载波的不同相位来传递数字信息，而振幅和 制 “0”和“1”。因此，2PSK 信号的时域表达式为：
• 其中，n表示第n个符号的绝对相位：
• 因此，式子可以改写为：
• 由于表示信号的两种码元的波形相同，极性相反，故2PSK信号一 般可以表述为一个双极性（bipolarity）全占空（100% duty ratio） 矩形脉冲序列与一个正弦载波的相乘，即：
2psk调制解调的原理
2psk调制解调的原理2PSK调制（2-Phase Shift Keying）是一种基本的数字调制方式。
它通过改变载波的相位来传输数字信号，每个数字比特对应两个不同的相位。
以下将详细解析2PSK调制的原理。
2PSK调制主要涉及到两个过程：调制和解调。
调制过程：1. 文字编码：将要传输的信息进行数字编码，例如使用二进制编码方式，将每个数字比特用0和1代表。
2. 符号分配：每个数字比特对应一个相位，通常选择相位0和相位π来表示0和1。
3. 载波生成：产生一个恒定频率和幅度的正弦波，这个波被称为载波信号。
4. 相位调制：根据编码的数字比特，将相应的相位信息融入到载波信号中。
比如，相位0可以对应载波信号的相位不变，而相位π可以对应载波信号的相位反转。
5. 调制信号生成：得到相位调制后的信号，该信号即为调制信号。
解调过程：1. 接收信号采样：接收到经过信道传输的调制信号，并对信号进行采样。
2. 相位判决：根据接收到的信号的相位信息，进行相位判决以确定每个数字比特的数值。
例如，如果接收到的信号相位为0，则判定为0；如果接收到的信号相位为π，则判定为1。
3. 数字解码：将解调的数字比特翻译回原始的信息字符。
2PSK调制的优点：1. 简单性：2PSK调制的实现比较简单，仅需要改变相位即可。
2. 抗噪声性能：2PSK调制的抗噪声性能较好，因为每个数字比特对应的相位差异明显，相位误差引起的误码率较低。
2PSK调制的局限：1. 带宽效率：2PSK调制一次只能传输一个比特，相比其他复杂调制方式，其带宽利用率较低。
2. 灵活性：2PSK调制只能传输二进制信号，不能传输多元信号。
总结：2PSK调制通过改变载波的相位来传输数字信号。
在调制过程中，信号经过文字编码、符号分配、载波生成和相位调制等步骤。
在解调过程中，信号经过接收信号采样、相位判决和数字解码等步骤。
2PSK调制简单易实现，抗噪声性能好，但带宽利用率相对较低，适用于二进制信号的传输。
【2017年整理】2ASK和2PSK调制与解调实验
实验二2ASK和2PSK调制与解调实验（一）实验目的1、掌握振幅键控(ASK)调制与解调的原理，并会用仿真软件绘制仿真的原理图，得出正确的波形图。
2、掌握相移键控(PSK)调制与解调的原理，并会用仿真软件绘制仿真的原理图，得出正确的波形图。
（二）实验设备计算机、SystemView软件（三）实验内容1、振幅键控(ASK)调制与解调：掌握振幅键控(ASK)调制与解调的原理，并用仿真软件绘制仿真的原理图，得出正确的波形图。
2、相移键控(PSK)调制与解调：掌握相移键控(PSK)调制与解调的原理，并用仿真软件绘制仿真的原理图，得出正确的波形图。
（四）实验原理1、2ASK调制部分：二进制幅度键控的调制器可用一个相乘器来实现。
对于2ASK 信号，相乘器则可以用一个开关电路来代替。
调制信号为1时，开关电路导通，为0时切断。
2ASK信号表达式：S（t）=a(n)Acos（ωct）式中：A-载波幅度，ωc -载波频率，a(n)-二进制数字信号2、2PSK二进制相移键控（2PSK ）就是根据数字基带信号的两个电平，使载波相位在连个不同的数值之间不通的数值之间切换的一种相位调制方法。
通常，两个载波相位相差π个弧度。
PSK 信号可以写成如下形式：Spsk （t ）=a(n)Acos （ωct ）1.调制部分：在2PSK 中，通常用相位0°或180°来分别表示1或-1.这里用调相法来生成2PSK ：将数字信号与载波直接相乘。
这也是DSB 信号产生的方法。
S2PSK （t ）=cos(ω0t+φ)， φ=0或πS2PSK （t ）= ACOS(ω0) a(n)=1-ACOS(ω0) a(n)= -12.解调部分2PSK 必须采用相干解调，同步载波是个关键问题。
相干接收2PSK 系统组成如图所示：对2PSK 信号相干接收的前提是首先进行载波提取，可采用平方环或科斯塔斯环来实现。
为分析方便起见，在本实验中可直接在接收端设信道输出图3-2-1 2PSK 系统组成置一个与发送端严格同步的本地载波源。
2PSK数字信号的调制与解调-分享版
信息对抗大作业一、实验目的。
使用MATLAB构成一个加性高斯白噪声情况下的2psk调制解系统，仿真分析使用信道编码纠错和不使用信道编码时，不同信道噪声比情况下的系统误码率。
二、实验原理。
数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。
为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。
这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。
数字调制技术的两种方法：①利用模拟调制的方法去实现数字式调制，即把数字调制看成是模拟调制的一个特例，把数字基带信号当做模拟信号的特殊情况处理；②利用数字信号的离散取值特点通过开关键控载波，从而实现数字调制。
这种方法通常称为键控法，比如对载波的相位进行键控，便可获得相移键控（PSK）基本的调制方式。
图1 相应的信号波形的示例1 0 1数字调相：如果两个频率相同的载波同时开始振荡，这两个频率同时达到正最大值，同时达到零值，同时达到负最大值，它们应处于"同相"状态；如果其中一个开始得迟了一点，就可能不相同了。
如果一个达到正最大值时，另一个达到负最大值，则称为"反相"。
一般把信号振荡一次（一周）作为360度。
如果一个波比另一个波相差半个周期，我们说两个波的相位差180度，也就是反相。
当传输数字信号时，"1"码控制发0度相位，"0"码控制发180度相位。
载波的初始相位就有了移动，也就带上了信息。
相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。
在2PSK中，通常用初始相位0和π分别表示二进制“1”和“0”。
因此，2PSK信号的时域表达式为(t)=Acos t+)其中，表示第n个符号的绝对相位：=因此，上式可以改写为图2 2PSK信号波形解调原理2PSK信号的解调方法是相干解调法。
实验六PSK调制与解调
实验六2PSK调制与解调一、实验目的1、理解二进制移相键控（Phase Shift Keying，PSK）调制和解调的基本原理；2、了解2PSK调制和解调的实现方法。
二、实验原理一个正弦载波。
如果它被一个双极性比特流按照图6-1所示的方案调制，它的极性将在每一次比特流极性改变时跟着改变。
图6-1对正弦波来说，极性的翻转就等价于反相。
因此，乘法器的输出就是BPSK（2PSK）信号。
二进制移相键控的解调可分两个步骤来考虑。
1、限带信号波形的恢复，使其转化到基带信号；2、从基带的限带波形里重建二进制消息比特流。
在本实验中，实现第一步依靠的是一个“窃取”的本地同步载波。
第二步的抽样判决由定标模块实现，最后还应线性解码，重建原始单极性基带信号。
解调原理如图6-2所示。
图6-2三、实验设备1、主机TIMS-301F2、TIMS基本插入模块（1）TIMS-148音频振荡器（Audio Oscillator）（2）TIMS-150乘法器（Multiplier）或TIMS-425正交模块（Quadrature Utilities），此模块集成了2个乘法器和1个加法器（3）TIMS-151移相器（Phase Shifer）（4）TIMS-153序列产生器（Sequence Generator）（5）TIMS-154可调低通滤波（Tuneable LPF）（6）TIMS-402定标模块(decision-maker module)（7）TIMS-406线性编码器（Line Code Encoder）（8）TIMS-407线性译码器（Line Code Decoder）3、计算机4、Pico虚拟仪器四、实验步骤1、将Tims系统中音频振荡（Audio Oscillator）、移相器（Phase Shifter）、序列码产生器（Sequence Generator）、线性编码器（Line-code Encode）、乘法器（Multiplier）按图6-3连接。
2PSK信号的调制与解调
通信原理大作业2PSK 信号的调制与解调3090401014 葛一飞 通信091 一．2PSK 信号的调制原理：当相移常数Kp=π时，当基带数字信号采用幅度为1宽度为TS 的矩形脉冲的双极性非归零码表示时，时域表示式为：受键控的载波相位按基带脉冲而改变的数字调制方式。
()()⎩⎨⎧-===-=P 1,0,0P ,1,1概率为概率为n n s n a a nT t g a t m ()()()t m K t cos t S p c PSK +=ω2()()⎩⎨⎧-==+=P 10P 1以概率以概率,a ,t cos ,a ,K t cos n c n p c ωω()()()⎩⎨⎧-==-=P 1,0,cos P ,1,cos 2以概率以概率n c n c PSK a t a t t S ωωt nT t g a t S c n s n PSK ωcos )()(2⎥⎦⎤⎢⎣⎡-=∑⎩⎨⎧→→”“”“如，相位相位010π二．2PSK 信号相干解调由于绝对移相方式是以某一相位作为基准的，因此解调时在接收端也必须有同样一个固定基准相位作为参考。
即采用相干解调2.波形显示调制：冲 据输出 S 2解调：总结：通过本次大作业，我们在MATLAB平台上对数字信号的传输系统进行了一次仿真，有效的完善了学习过程中实践不足的问题，同时进一步巩固了原先的基础知识。
通过仿真，基本掌握了MATLAB的基本功能和使用方法，对数字基带传输系统有了一定的了解，加深了对2PSK信号的调制原理的认识，理解了如何对他们进行调制，通过使用MATLAB仿真，对个调制和解调电路中各元件的特性有了较为全面的理解。
可以说这次大作业使我收益颇丰，对通信原理也有了新的认识。
附录:Matlab程序代码：clcclose allclear allcodn=60;fc=6e+3;fs=fc*6;bode=1000;code=round(rand(1,codn));code_len=round(1/bode/(1/fs));for i=1:codnx0((i-1)*code_len+1:code_len*i)=code(i);endx=2*x0-1;car=cos(2*pi*fc/fs*(0:length(x0)-1));y=x.*car;figuresubplot(211)plot(x)axis([0 length(x0) -1.5 1.5]) grid onzoom ontitle('原始基带信号')subplot(212)plot(y)zoom ongrid ontitle('2PSK的频谱')ay=abs(fft(y));f=0:fs/length(y):fs/2;ay=ay(1:length(f));figureplot(f,ay)zoom ongrid ontitle('2PSK信号')z=y.*car;fl=fir1(64,fc/fs*2); ¨z1=2*filter(fl,1,z);figuresubplot(211)plot(z)grid ontitle('混频器输出信号') subplot(212)plot(z1)grid ontitle('解调得到的基带信号') zoom onfiguresubplot(211)plot(x)axis([0 length(x0) -1.5 1.5]) grid ontitle('原始基带信号')subplot(212)plot(z1)axis([0 length(x0) -1.5 1.5]) grid ontitle('解调得到的基带信号') zoom on。
2psk解调原理
2psk解调原理
2PSK（相移键控）解调原理如下：
2PSK信号采用相干解调方法解调，解调的关键在于恢复出一个与发送端载波同频同相的本地参考载波。
2PSK信号经过信道传输之后，再和载波相乘，然后经过低通滤波后抽样判决恢复出原始基带码元信号。
在相干解调中，如何得到与接收的2PSK信号同频同相得相干载波是一个关键问题，这涉及到同步原理。
需要注意的是，在实际应用中，由于绝对移相方式存在相位模糊现象（或反向工作现象或倒π现象），因此在相干解调中，通常采用差分移相键控方式。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
信息对抗大作业一、实验目的。
使用 MATLAB构成一个加性高斯白噪声情况下的2psk 调制解系统，仿真分析使用信道编码纠错和不使用信道编码时，不同信道噪声比情况下的系统误码率。
二、实验原理。
数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。
为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。
这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。
数字调制技术的两种方法：①利用模拟调制的方法去实现数字式调制，即把数字调制看成是模拟调制的一个特例，把数字基带信号当做模拟信号的特殊情况处理；②利用数字信号的离散取值特点通过开关键控载波，从而实现数字调制。
这种方法通常称为键控法，比如对载波的相位进行键控，便可获得相移键控（PSK）基本的调制方式。
图 1相应的信号波形的示例101数字调相：如果两个频率相同的载波同时开始振荡，这两个频率同时达到正最大值，同时达到零值，同时达到负最大值，它们应处于" 同相 " 状态；如果其中一个开始得迟了一点，就可能不相同了。
如果一个达到正最大值时，另一个达到负最大值，则称为" 反相 " 。
一般把信号振荡一次（一周）作为360 度。
如果一个波比另一个波相差半个周期，我们说两个波的相位差180 度，也就是反相。
当传输数字信号时， "1" 码控制发 0 度相位， "0" 码控制发 180 度相位。
载波的初始相位就有了移动，也就带上了信息。
相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。
在2PSK 中，通常用初始相位0 和π分别表示二进制“1”和“ 0”。
因此， 2PSK信号的时域表达式为(t)=Acos t+)其中，表示第 n 个符号的绝对相位：=因此，上式可以改写为图 22PSK信号波形解调原理2PSK信号的解调方法是相干解调法。
由于PSK信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息来解调信号。
下图2-3 中给出了一种2PSK信号相干接收设备的原理框图。
图中经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。
判决器是按极性来判决的。
即正抽样值判为1，负抽样值判为0.2PSK信号相干解调各点时间波形如图3所示.当恢复的相干载波产生180°倒相时, 解调出的数字基带信号将与发送的数字基带信号正好是相反, 解调器输出数字基带信号全部出错.图 32PSK信号相干解调各点时间波形这种现象通常称为 " 倒π " 现象 . 由于在 2PSK 信号的载波恢复过程中存在着 180°的相位模糊 , 所以 2PSK信号的相干解调存在随机的 " 倒π " 现象 , 从而使得 2PSK方式在实际中很少采用。
2PSK的基本原理：相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。
在 2PSK中，通常用初始相位为 0 和π表示二进制的“ 1”和“ 0”。
因此 2PSK的信号的时域表达式为：ｅ 2psk (t)=Acos(ω c t+φ n)其中，φ n 表示第0n 个符号的绝对相位：发送“ 0”时φn=π发送“ 1”时因此，上式可改写为Acosω c t概率为Pｅ2psk (t)=- Acosω c t概率为1-P图 4 2PSK信号的时间波形T sA-A100由于表示信号的两种码元的波形相同，记性相反，鼓 2PSK信号一般可以表述为一个双极性全占空矩形脉冲序列与一个正弦载波相乘，即ｅ 2psk(t)=s(t)cosω c t其中s(t)=∑a g(t-nTs )n这里， g(t)是脉宽为Ts 的单个矩形脉冲，而an 得统计特性为1概率为 Pa n=-1概率为 1-P即发送二进制符号“ 0”时（ an 取 +1），ｅ 2psk(t)取 0 相位；发送二进制符号“1”时（ an 取-1 ），ｅ 2psk(t)取π 相位。
2PSK相干解调系统但是由于2PSK信号的载波回复过程中存在着180°的相位模糊，即恢复的本地载波与所需相干载波可能相同，也可能相反，这种相位关系的不确定性将会造成解调出的数字基带信号与发送的基带信号正好相反，即“1”变成“ 0”吗“ 0”变成“ 1”，判决器输出数字信号全部出错。
这种现象称为2PSK方式的“倒π ”现象或“反相工作”。
10011atT sbtctdt10011et图 5 2PSK信号的解调器波形图2PSK信号在一个码元的持续时间Ts 内可以表示为u1T(t)发送“ 1”时S T (t)=u oT(t)=- u1T(t)发送“ 0”时期中Acosω c t0< t < Tsu1T(t)=0其他设发送端发出的信号如上式所示，则接收端带通滤波器输出波形y(t) 为[a+n c (t)]cosω c t-n s(t)sinω c t发送“ 1”时y(t)=[-a+n c(t)]cosω c t-n s(t)sinω c t发送“ 0”时y(t)经过想干解调（相乘—低通）后，送入抽样判决器的输入波形为a+n c(t)发送“ 1”时x(t)=-a+n c(t)发送“ 0”时由于 nc(t)是均值为0，方差为σ 2 的高斯噪声，所以x(t)的一维概率密度函数为1-（x-a）2f1(x)=exp发送“ 1”时2πσn2σ2n21-（x+a）2f(x)=exp发送“ 0”时2πσn22σn由最佳判决门限分析可知，在发送“ 1”和“ 0”概率相等时，即 P(1)=P(0) 时，最佳门限 b*=0. 此时，发“ 1”而错判为“ 0”的概率为P(0/1)=P(x≦ 0)= ∫0-∞ f 1(x)dx=1/2erfc(r )22式中： r=a /2 σn同理，发“ 0 而错判为“ 1”的概率为P(1/0)=P(x＞0)=∫ -∞f0(x)dx=1/2erfc(r )2PSK信号的调制器键控法原理方框图如图：双极性s(t )不归零e2 PSK (t)码型变换乘法器cos c t图 6 2PSK信号的调制器原理方框图说明： 2psk 调制器可以采用相乘器，也可以采用相位选择器就模拟调制法而言，与产生2ASK 信号的方法比较，只是对s(t)要求不同，因此2PSK 信号可以看作是双极性基带信号作用下的DSB调幅信号。
而就键控法来说，用数字基带信号s(t)控制开关电路，选择不同相位的载波输出，这时s(t)为单极性NRZ或双极性NRZ脉冲序列信号均可。
2PSK 信号属于DSB信号，它的解调，不再能采用包络检测的方法，只能进行相干解调。
2PSK信号的解调通常采用相干解调法原理框图如图：带通e2PSK (t )滤波器a相乘器c低通d抽样e滤波器判决器输出cos ctb定时脉冲图 7 2PSK信号的相干解调原理方框图图中，假设相干载波的基准相位与2PSK信号的基准一致（通常默认为0 相位）。
说明：由于PSK信号的功率谱中五载波分量，所以必须采用相干解调的方式。
在相干解调中，如何得到同频同相的本地载波是个关键问题。
只有对PSK信号进行非线性变换，才能产生载波分量。
2PSK信号经过带通滤波器得到有用信号，经相乘器与本地载波相乘再经过低通滤波器得到低频信号v(t)，再经抽样判决得到基带信号。
2PSK相干解调系统性能原理框图如图：带通低通抽样输出发送端信道相乘器滤波器滤波器判决器P e s T (t )y i (t)y(t)x(t )定时n i (t ) 2 cos c t脉冲图 8 2PSK相干解调系统性能原理方框图由最佳判决门限分析可知，在发送“1”符号和发送“ 0”符号概率相等时，最佳判决门限b* = 0。
此时，发“1”而错判为“0”的概率为P ( 0 / 1 )P ( x0 )0f 1 ( x ) dx同理，发送“ 0”而错判为“ 1”的概率为1erf c r 2P ( 1 / 0 )P ( x0 ) f 0 ( x ) dx故 2PSK信号相干解调时系统的总误码率为1erf c r 2P e P(1)P (0 / 1) P( 0) P(0 /1) 1 erfc r2在大信噪比条件下，上式可近似为P e 1e r2r三、实验过程。
根据代码即可实现全部过程，图形如下：代码如下：clc;clear all;close all;max=15;s=randint(1,max);%长度为 max 的随机二进制序列?Sinput=[];for n=1:length(s);if s(n)==0;A=zeros(1,2000);else s(n)==1;A=ones(1,2000);endSinput=[Sinput A];endfigure(1);subplot(211);plot(Sinput);grid onaxis([0 2000*length(s) -2 2]);title('输入信号波形');Sbianma=encode(s,7,4,'hamming');%汉明码编码后序列? a1=[];b1=[];f=1000;t=0:2*pi/1999:2*pi;for n=1:length(Sbianma);if Sbianma(n)==0;B=zeros(1,2000);%每个值 2000 个点 ?????else Sbianma(n)==1;B=ones(1,2000);enda1=[a1 B];%s(t),码元宽度 2000????c=cos(2*pi*f*t);%载波信号 ?b1=[b1 c];%与 s(t)等长的载波信号，变为矩阵形式?endfigure(2);subplot(211)plot(a1);grid on;axis([0 2000*length(Sbianma) -2 2]);title('编码后二进制信号序列 '); a2=[];b2=[];for n=1:length(Sbianma);if Sbianma(n)==0;C=ones(1,2000);%每个值 2000 点 ?????????d=cos(2*pi*f*t);%载波信号 ?????else Sbianma(n)==1;C=ones(1,2000);d=cos(2*pi*f*t+pi);%载波信号 ?????enda2=[a2 C];%s(t)，码元宽度 2000?b2=[b2 d];%与 s(t)等长的载波信号 ?endtiaoz=a2.*b2;%e(t)调制 ?figure(3);subplot(211);plot(tiaoz);grid on;axis([0 2000*length(Sbianma) -2 2]);title('2psk已调制信号 ');figure(2);subplot(212);plot(abs(fft(a1)));axis([0 2000*length(Sbianma) 0 400]);title('编码后二进制信号序列频谱');figure(3);subplot(212);plot(abs(fft(tiaoz)));axis([0 2000*length(Sbianma) 0 400]);title('2psk信号频谱 ')%-----------------带有高斯白噪声的信道 ----------------------?tz=awgn(tiaoz,10);%信号 tiaoz加入白噪声，信噪比为10?figure(4);subplot(211);plot(tz);grid onaxis([0 2000*length(Sbianma) -2 2]);title('通过高斯白噪声后的信号 ');figure(4);subplot(212);plot(abs(fft(tz)));axis([0 2000*length(Sbianma) 0 800]);title('加入白噪声的 2psk 信号频谱 ');%-------------------同步解调-----------------------------?jiet=2*b1.*tz;%同步解调 ?figure(5);subplot(211);plot(jiet);grid onaxis([0 2000*length(Sbianma) -2 2]);title('相乘后的信号波形 ')figure(5);subplot(212);plot(abs(fft(jiet)));axis([0 2000*length(Sbianma) 0 800]);title('相乘后的信号频率 ');%----------------------低通滤波器 ---------------------------fp=500;fs=700;rp=3;rs=20;fn=11025;ws=fs/(fn/2);wp=fp/(fn/2);%计算归一化角频率?[n,wn]=buttord(wp,ws,rp,rs);%计算阶数和截止频率?[b,a]=butter(n,wn);%计算 H（ z ） ?figure(6);freqz(b,a,1000,11025);subplot(211);axis([0 40000 -100 3])title('lpf频谱图 ');jt=filter(b,a,jiet);figure(7);subplot(211);plot(jt);grid onaxis([0 2000*length(Sbianma) -2 2 ]);title('经低通滤波器后的信号波形');figure(7);subplot(212);plot(abs(fft(jt)));axis([0 2000*length(Sbianma) 0 800]); title('经低通滤波器后的信号频率');%-----------------------抽样判决--------------------------? for m=1:2000*length(Sbianma);if jt(m)<0;jt(m)=1;else jt(m)>0;jt(m)=0;endendfigure(8);subplot(211);plot(jt)grid onaxis([0 2000*length(Sbianma) -2 2]);title('经抽样判决后信号jt(t)波形')figure(8);subplot(212);plot(abs(fft(jt)));axis([0 2000*length(Sbianma) 0 800]);title('经抽样判决后的信号频谱');grid on;n=500:2000:2000*length(Sbianma);a5=[];a5=[a5 jt(n)];s1=decode(a5,7,4,'hamming');a6=[];for n=1:length(s1);if s1(n)==0;G=zeros(1,2000);else s1(n)==1;G=ones(1,2000);enda6=[a6 G];endfigure(1);subplot(212);plot(a6);grid onaxis([0 2000*length(s) -2 2]);title('汉明码译码后的波形')grid on%------------------2psk误码率仿真-------------------------? snrdB_min=-10;snrdB_max=10;snrdB=snrdB_min:1:snrdB_max;Nsymbols=200;snr=10.^(snrdB/10);h=waitbar(0,'SNR?Iteration');len_snr=length(snrdB);for j=1:len_snrwaitbar(j/len_snr);sigma=sqrt(1/(2*snr(j)));error_count=0;for k=1:Nsymbolsd=round(rand(1));%随即数据 ?x_d=2*d-1;%0，1分别转化为-1，1?????????n_d=sigma*randn(1);%加噪y_d=x_d+n_d;%加噪后接收?if y_d>0d_est=1;elsed_est=0;endif(d_est~=d)error_count=error_count+1;endenderrors(j)=error_count;endber_sim=errors/Nsymbols;ber_theor=(erfc(sqrt(snr))).*(1-0.5*erfc(sqrt(snr)));figure(9);semilogy(snrdB,ber_theor,'-',snrdB,ber_sim,'*');axis([snrdB_min snrdB_max 0.0001 1]);xlabel('信噪比 ');ylabel('误码率 ');title('2psk信噪比误码率关系图');legend('理论值 ','实际值')。