二进制移相键控(2PSK)仿真实验报告
一、实验任务
使用SystemView软件对二进制移相键控(2PSK)的调制与解调进行仿真。
二、实验要求
-
获取系统各点时域波形,波形、坐标、标题等要清楚;滤波器的单位冲击相应和幅频特性曲线;
-
获取主要信号的频谱;
-
获取眼图;
-
提取相干载波;
-
完成实验报告。
三、实验说明
-
传码率:20kbit/s;
-
调制可采用模拟调制法或数字键控法;
-
解调为相干解调;
-
要求报告内容完整,结构清晰,撰写规范,波形充分、正确、清晰,便于分析。仿真模型及其说明、图符参数表、仿真波形、波形分析等内容完整;
-
报告要对数据、波形进行解释、分析、说明,获取正确信息,得出合理、正确结论。
-
本题适用于本校非文科类所有专业。
四、实验原理
2PSK原理:
用两个初相相隔为 π \pi π 的载波来传递二进制信息,初相为 0 0 0的载波表示二进制的0,初相为** π \pi π的载波表示二进制的1**,两个载波只有相位不同,频率和振幅都一致;
可用时域表达式: e 2 p s k ( t ) = A c o s ( ω c t + θ n ) e_{2psk}(t)=Acos(\omega_ct+ \theta_n) e2psk(t)=Acos(ωct+θn)来表示,其中第 n n n个二进制信息的绝对相位 θ n = { 0 发送 0 时 π 发送 1 时 \theta_{n}=\left\{ \begin{aligned} 0 \ \ 发送0时\\ \pi \ \ 发送1时\end{aligned} \right. θn={0 发送0时π 发送1时;
即可化简为该式: e 2 p s k ( t ) = { A c o s ( ω c t ) 发送 0 时 − A c o s ( ω c t ) 发送 1 时 e_{2psk}(t)=\left\{ \begin{aligned} Acos(\omega_ct) \ \ 发送0时\\-Acos(\omega_ct) \ \ 发送1时\end{aligned} \right. e2psk(t)={Acos(ωct) 发送0时−Acos(ωct) 发送1时;
由于表示信号的两种码元的波形相同,极性相反,故2PSK信号一般可以表述为一个双极性全占空矩形脉冲序列与一个 余弦载波的相乘,即 e 2 p s k ( t ) = s ( t ) c o s ( ω c t ) e_{2psk}(t)=s(t)cos(\omega_ct) e2psk(t)=s(t)cos(ωct);其中 s ( t ) s(t) s(t)是双极性非归零矩形脉冲信号;
调制原理:
调制分为模拟调制和数字调制(本实验采用模拟调制);
原理框图:
相干解调原理:
例如:基带信号 A A A和载波信号 c o s ( ω t + θ ) cos(\omega t+\theta) cos(ωt+θ)调制得到调制信号 A c o s ( ω t + θ ) Acos(\omega t +\theta) Acos(ωt+θ),乘法器将调制信号和一同频同相的相干载波信号相乘,得到 A c o s 2 ( ω t + θ ) = A 2 + A 2 c o s ( 2 ω t + 2 θ ) Acos^2(\omega t+\theta)={A\over 2}+{A\over 2}cos(2\omega t+2 \theta) Acos2(ωt+θ)=2A+2Acos(2ωt+2θ);通过低通滤波器将高频信号滤除,得到 A 2 {A\over 2} 2A,即可得到原始信号 A A A;
载波频率应该远高于调制信号频率,这样才能比较容易的解调出来;
为了防止噪声干扰,且电压并不是理想的阶跃信号,所以添加抽样判决器,其基本功能是:信号经过抽样判决器,即可确定接收码元是“1”还是“0”。假设抽样判决门限为b,当信号抽样值大于b时,判为“1”码;信号抽样值小于b时,判为“0”码。
解调过程波形图:
波形图标号 | 信号 | 表达式 |
---|---|---|
a | 调制信号 | v a = { A s i n ( ω c t ) 发送 0 时 − A s i n ( ω c t ) 发送 1 时 v_a=\left\{ \begin{aligned} Asin(\omega_ct) \ \ 发送0时\\-Asin(\omega_ct) \ \ 发送1时\end{aligned} \right. va={Asin(ωct) 发送0时−Asin(ωct) 发送1时 |
b | 相干载波 | v b = s i n ( ω c t ) v_b=sin(\omega_ct) vb=sin(ωct) |
c | 含高频分量的解调信号 | v c = v a v b = ± A s i n 2 ( ω c t ) = ± A 2 ∓ A 2 c o s ( 2 ω c t ) v_c=v_av_b=\pm Asin^2(\omega_c t)=\pm {A\over 2}\mp {A\over 2}cos(2\omega_c t) vc=vavb=±Asin2(ωct)=±2A∓2Acos(2ωct) |
d | 不含高频分量的解调信号 | v c = { A 2 发送 0 时 − A 2 发送 1 时 v_c=\left\{ \begin{aligned} {A\over 2} \ \ 发送0时\\-{A\over 2} \ \ 发送1时\end{aligned} \right. vc=⎩ ⎨ ⎧2A 发送0时−2A 发送1时 |
e | 最终解调信号 | v d = { A 2 发送 0 时 − A 2 发送 1 时 v_d=\left\{ \begin{aligned} {A\over 2} \ \ 发送0时\\-{A\over 2} \ \ 发送1时\end{aligned} \right. vd=⎩ ⎨ ⎧2A 发送0时−2A 发送1时 |
科斯塔斯环法原理:
测试点 | 电压值 | 含义 |
---|---|---|
a | v a = e 2 p s k ( t ) = s ( t ) c o s ( ω c t + θ ) , s ( t ) = 1 或 − 1 v_a=e_{2psk}(t)=s(t)cos(\omega_ct+\theta),s(t)=1或-1 va=e2psk(t)=s(t)cos(ωct+θ),s(t)=1或−1 | 输入的已调制信号 |
b | v b = c o s ( ω c t + φ ) v_b=cos(\omega_ct+\varphi) vb=cos(ωct+φ) | 压控振荡电压 |
c | v c = c o s ( ω c t + φ + π 2 ) = s i n ( ω c t + φ ) v_c=cos(\omega_ct+\varphi+{\pi \over 2} )=sin(\omega_ct+\varphi) vc=cos(ωct+φ+2π)=sin(ωct+φ) | b b b点电压移相 π 2 \pi \over 2 2π |
d | v d = v a v b = s ( t ) c o s ( ω c t + θ ) c o s ( ω c t + φ ) = 1 2 s ( t ) [ c o s ( φ − θ ) + c o s ( 2 ω c t + φ + θ ) ] v_d=v_av_b=s(t)cos(\omega_ct+\theta)cos(\omega_ct+\varphi)={1\over2}s(t)[cos(\varphi-\theta)+cos(2\omega_ct+\varphi+\theta)] vd=vavb=s(t)cos(ωct+θ)cos(ωct+φ)=21s(t)[cos(φ−θ)+cos(2ωct+φ+θ)] | a a a点电压乘以 b b b点电压 |
e | v e = v a v c = s ( t ) c o s ( ω c t + θ ) s i n ( ω c t + φ ) = 1 2 s ( t ) [ s i n ( φ − θ ) + s i n ( 2 ω c t + φ + θ ) ] v_e=v_av_c=s(t)cos(\omega_ct+\theta)sin(\omega_ct+\varphi)={1\over2}s(t)[sin(\varphi-\theta)+sin(2\omega_ct+\varphi+\theta)] ve=vavc=s(t)cos(ωct+θ)sin(ωct+φ)=21s(t)[sin(φ−θ)+sin(2ωct+φ+θ)] | a a a点电压乘以 c c c点电压 |
f | v f = 1 2 s ( t ) c o s ( φ − θ ) v_f={1\over2}s(t)cos(\varphi-\theta) vf=21s(t)cos(φ−θ) | 经过低通滤波器的 d d d点电压 |
g | v g = 1 2 s ( t ) s i n ( φ − θ ) v_g={1\over2}s(t)sin(\varphi-\theta) vg=21s(t)sin(φ−θ) | e e e点电压经过低通滤波器 |
h | v h = v f v g = 1 4 s 2 ( t ) s i n ( φ − θ ) c o s ( φ − θ ) = 1 8 s 2 ( t ) s i n ( 2 ( φ − θ ) ) = 1 8 s i n ( 2 ( φ − θ ) ) v_h=v_fv_g={1\over 4}s^2(t)sin(\varphi-\theta)cos(\varphi-\theta)={1\over 8}s^2(t)sin(2(\varphi-\theta))={1\over 8}sin(2(\varphi-\theta)) vh=vfvg=41s2(t)sin(φ−θ)cos(φ−θ)=81s2(t)sin(2(φ−θ))=81sin(2(φ−θ)) | f f f点电压乘以 g g g点电压 |
i | v i = 1 8 s i n ( 2 ( φ − θ ) ) ≈ 1 4 ( φ − θ ) ,当 φ 和 θ 很接近时 v_i={1\over 8}sin(2(\varphi-\theta)) \approx {1\over 4}(\varphi-\theta),当\varphi和\theta很接近时 vi=81sin(2(φ−θ))≈41(φ−θ),当φ和θ很接近时 | h h h点电压过滤噪声和干扰 |
其中 θ \theta θ为已调制信号载波相位, φ \varphi φ为压控振荡电压相位;
当 ( φ − θ ) ≠ 0 (\varphi-\theta) \neq 0 (φ−θ)=0时,即 v i ≠ 0 v_i \neq 0 vi=0说明输入信号和输出信号的相位差不为0,此时压控振荡器会接受到一个控制电压,对振荡器输出信号的频率实施控制,控制的结果是将压控振荡器的输出振荡频率 f o f_o fo拉向环路输入信号频率 f i f_i fi ,通过反馈机制最终二者相等时,环路被锁定 ,称为入锁。
当 φ − θ = 0 \varphi-\theta=0 φ−θ=0时, v i = 0 v_i=0 vi=0,压控振荡器的控制电压为0;此时压控振荡器的输出 v b = c o s ( ω c t + φ ) = c o s ( ω c t + θ ) v_b=cos(\omega_ct+\varphi)=cos(\omega_ct+\theta) vb=cos(ωct+φ)=cos(ωct+θ),和输入的已调制信号同频同相,即提取的相干载波;同时, v f = 1 2 s ( t ) c o s ( φ − θ ) = 1 2 s ( t ) v_f={1\over2}s(t)cos(\varphi-\theta)={1\over 2}s(t) vf=21s(t)cos(φ−θ)=21s(t),即可以作为相干解调信号的输出;
从相干解调的框图中也可以看出科斯塔斯环法中的 v f v_f vf对应于如下解调原理框图中的 v d v_d vd:
所以科斯塔斯环法同时具备提取相干载波和相干解调的功能;
科斯塔斯环的潜在问题:
锁相环(科斯塔斯环)在 ( φ − θ ) (\varphi - \theta) (φ−θ)稳定的点有两个,在 ( φ − θ ) = 0 和 π (\varphi - \theta)=0和\pi (φ−θ)=0和π处,所以科斯塔斯环提取的载频存在 18 0 o 180^o 180o的相位模糊性,即恢复的本地载波与所需的相干载波可能同相,也可能反相,这种相位关系的不确定性将会造成解调出的数字基带信号与发送的数字基带信号正好相反,即“1”变为“0”,“0”变为“1”,判决器输出数字信号全部出错。这种现象称为2PSK方式的“倒 π \pi π”现象。;
五、仿真系统组成和参数设置
整体电路图:
主要电路及参数:
模拟键控法调制:
模拟键控调制电路:
模拟键控调制电路参数设置:
元器件名称 | 参数 | 序号 | 功能 |
---|---|---|---|
PN序列产生器 | ![]() | 0 | 产生为0,1的随机序列码作为基带信号 |
余弦信号产生器 | ![]() | 2 | 产生余弦信号作为载波信号 |
乘法器 | ![]() | 22 | 将基带信号和载波信号相乘产生调制信号 |
噪声模拟:
噪声模拟电路:
噪声模拟电路参数设置:
科斯塔斯环实现相干载波提取和相干解调:
科斯塔斯环相干载波提取和相干解调电路:
科斯塔斯环参数设置:
元器件名称 | 参数 | 序号 | 功能 |
---|---|---|---|
乘法器 | ![]() | 12,13,14 | 信号相乘 |
低通滤波器 | ![]() | 9 | 去除高频分量 |
低通滤波器 | ![]() | 11 | 去除高频分量 |
低通滤波器 | 参数同序号11的低通滤波器 | 27 | 去除高频分量 |
抽样器 | ![]() | 16 | 去除毛刺信号 |
缓冲器 | ![]() | 20 | 变为幅度只取-1V,1V的信号(通过设置True Output和Flase Output实现) |
保持器 | ![]() | 17 | 实现相干解调信号放大 |
VCO压控振荡器 | ![]() | 8 | 生成相干载波 |
低通滤波器 | ![]() | 10 | 参与实现环路滤波器,滤去高频分量和噪声 |
放大器 | ![]() | 15 | 参与实现环路滤波器 |
六、实验结果分析:
时钟频率设置
根据奈奎斯特采样定理,当采样频率大于信号中最高频率的2倍时,采样之后的数字信号才能完整地保留了原始信号中的信息;
由于原信号频率为(PN 序列 20kHz, 载波信号100kHz)同时为了尽可能才更多的点以保证图像分辨率,选取 1000kHz 的采样频率。采样时长选取 5ms,得到总采样点数为 5001 个。
系统各点时域波形分析
总瀑布图分析
从上到下分别为:基带信号、载波信号、调制信号、相干载波、相干解调、解调信号、带高频分量的相干解调信号;
可以看出:
基带信号波形是矩形脉冲序列,载波信号波形为余弦信号;
调制信号在基带信号为-1时,相位取 π \pi π;基带信号为1时,相位取 0 0 0,符合2PSK调制原理;
载波信号和提取出的相干载波的波形完全相同;
相干解调信号中电平变化不是理想的阶跃型,但是仍然可以粗略看出与基带信号波形变化趋势相似;
最终解调信号波形和基带信号波形只有轻微的相位不同,可能是电路延时造成的,说明解调成功;
带高频分量的相干解调信号也可以看出变化趋势和基带信号基本一致;
基带信号波形
载波信号波形
调制信号波形
相干载波波形
相干解调信号波形
最终解调信号波形
滤波器的单位冲击响应和幅频特性曲线(20KHz的低通滤波器)
单位冲击响应
频幅特性曲线
可以看出在截止频率20KHz之后,输入信号开始失真;
主要信号频谱分析
输入基带信号(20KHz的PN序列码频谱):
载波信号频谱(20KHz的余弦信号):
载波信号的频谱中心频率在 100 e + 3 100e+3 100e+3处,和我们开始时载波设置的信号频率 100 K H z 100KHz 100KHz相符;
已调制信号的频谱:
上下边带为80KHz,120KHz,中心频率为100KHz,中心频率符合我们设置的载波频率100KHz;
解调信号频谱:
解调信号频率分布和基带信号大致相同;
眼图分析
眼图原理:
眼图分析法是用示波器观察接收到的信号波形的方法来分析码间串扰和噪声对系统性能的影响。如果将输入波形输入示波器的Y轴,并且当示波器的水平扫描周期和码元定时同步时,在示波器上显示的图形很像人的眼睛,故被称为眼图。
在无码间串扰和噪声的理想情况下,波形无失真,“眼”开启得最大;当有码间串扰时,波形失真,引起“眼”的部分闭合;若加上噪声的干扰,眼图的线条变得模糊,“眼”开启得小了。所以眼图的“眼”张开大小表示了失真程度;
眼图实验结果:
不同情况 | 眼图 | 变化描述 |
---|---|---|
理想情况(无噪声) | ![]() | 眼图线条清晰 |
增加噪声 | ![]() | 眼图线条模糊 |
眼图时间切片设置(5个眼): 时间切片长度 = 采用周期 T ∗ 眼个数 n = n f = 1 20000 ∗ 5 = 2.5 ∗ 1 0 − 4 时间切片长度=采用周期T*眼个数n={n\over f}={1\over 20000}*5=2.5*10^{-4} 时间切片长度=采用周期T∗眼个数n=fn=200001∗5=2.5∗10−4
实验结果
科斯塔斯环法提取的相干载波和原始载波相同,使用科斯塔斯环法并未出现倒 π \pi π现象,最终解调信号和原始基带信号基本相同;