CN114244664A - 一种航空数字通信仿真***及其仿真方法 - Google Patents

一种航空数字通信仿真***及其仿真方法 Download PDF

Info

Publication number
CN114244664A
CN114244664A CN202111543861.0A CN202111543861A CN114244664A CN 114244664 A CN114244664 A CN 114244664A CN 202111543861 A CN202111543861 A CN 202111543861A CN 114244664 A CN114244664 A CN 114244664A
Authority
CN
China
Prior art keywords
signal
code
function
code element
module
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111543861.0A
Other languages
English (en)
Inventor
张伟业
汪坤
叶晨
侯静
王晨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Comac Software Co ltd
Original Assignee
Comac Software Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Comac Software Co ltd filed Critical Comac Software Co ltd
Priority to CN202111543861.0A priority Critical patent/CN114244664A/zh
Publication of CN114244664A publication Critical patent/CN114244664A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0008—Modulated-carrier systems arrangements for allowing a transmitter or receiver to use more than one type of modulation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/02—Amplitude-modulated carrier systems, e.g. using on-off keying; Single sideband or vestigial sideband modulation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/10—Frequency-modulated carrier systems, i.e. using frequency-shift keying
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/18—Phase-modulated carrier systems, i.e. using phase-shift keying
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/18—Phase-modulated carrier systems, i.e. using phase-shift keying
    • H04L27/22—Demodulator circuits; Receiver circuits
    • H04L27/227—Demodulator circuits; Receiver circuits using coherent demodulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明提供一种航空数字通信仿真***及其仿真方法，该航空数字通信仿真***，包括数字基带传输***和数字带通***；所述数字基带传输***包括基本编码模块、***间串扰模块、眼图模块；所述数字带通***包括2ASK调制模块、2FSK调制模块、2PSK调制模块。本发明航空数字通信仿真***的知识针对性强、仿真效果好、知识面覆盖全，人机交互界面友好，且运算快。本发明的航空数字通信仿真方法包括基本编码设计，***间串扰设计，眼图设计，2ASK设计与解调，2FSK设计与解调，2PSK设计与解调。本发明的航空数字通信仿真方法能针对航空电子企业的设计培训特点，适应于对企业新员工进行基础的数字通信子***培训。

Description

一种航空数字通信仿真***及其仿真方法
技术领域
本发明涉及民用航空技术领域，尤其是涉及一种航空数字通信仿真***及其仿真方法。
背景技术
通信***是民用飞机航空电子***的重要组成之一，也是航空电子***设计公司重点培训基础课程之一，通信***原理培训有着理论抽象、数学运算推导复杂、培训课件呆板、波形图手动绘制、实践操作性强、新理论新概念频出等特点，导致培训枯燥并难以理解，新员工入职培训效果大打折扣，使后续通信***设计类项目开展困难。
当前航空领域通信***原理培训主要以硬件设备及专业示波器为主，需要特定的培训室和通信***专业仪器仪表设备，存在投入成本高、仪表设备损坏率高、教员专业性强等问题，而纯软件人机界面交互式培训正好能避免上述问题。
市场急需一种知识针对性强、仿真效果好、知识面覆盖全，人机交互界面友好，且运算量不太大，能针对航空电子企业的设计培训特点，适应于对企业新员工进行基础的数字通信子***培训仿真***。
发明内容
为解决现有技术中的问题，本发明提供一种知识针对性强、仿真效果好、知识面覆盖全，人机交互界面友好，且运算快的航空数字通信仿真***。该航空数字通信仿真***能针对航空电子企业的设计培训特点，适应于对企业新员工进行基础的数字通信子***培训。
为实现上述目的，本发明提供如下技术方案：一种航空数字通信仿真***，包括数字基带传输***和数字带通***；所述数字基带传输***包括基本编码模块、***间串扰模块、眼图模块；所述数字带通***包括2ASK调制模块、2FSK调制模块、2PSK调制模块；
所述基本编码模块用于对单极性归零码、单极性不归零码、双极性归零码、双极性不归零码、AMI码、HDB3码的编码设计，并仿真出对应编码的波形；所述***间串扰模块用于对***间串扰设计，仿真出对应的时域和频率波形；所述眼图模块用于眼图设计，仿真出不同传输速率信号的时域波形和对应的眼图波形；
所述2ASK调制模块包括2ASK设计模块和2ASK解调模块，所述2ASK设计模块用于仿真生成已调2ASK信号，所述2ASK解调模块用于根据已调2ASK信号，恢复输出对应的原始波形；所述2FSK调制模块包括2FSK设计模块和2FSK解调模块，所述2FSK设计模块用于仿真生成已调2FSK信号，所述2FSk解调模块用于根据已调2FSK信号，恢复输出对应的原始波形；所述2PSK调制模块包括2PSK设计模块和2PSK解调模块，所述2PSK设计模块用于仿真生成已调2PSK信号，所述2PSK解调模块用于根据已调2PSK信号恢复输出对应的原始波形。
优选的技术方案，所述基本编码模块采用rand函数随机生成随机序列，并用round函数四舍五入取值，再用for循环对随机序列进行读取，编写出单极性或双极性归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则前t0/2的时间间隔里设置为1，后t0/2的时间间隔里设置为0；如果码元为0，则前t0/2的时间间隔里设置为0(-1)，后t0/2的时间间隔里设置为0；再编写单极性或双极性非归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则整个t0时间间隔里设置为1；如果码元为0，则整个t0时间间隔里设置为0(-1)；编写AMI码，读取随机序列，用for循环分别记录原码中1的数量和0的数量，并对1码进行交替赋值为+1和-1,0则保持不变；在AMI码的基础上，继续采用for循环0的数目判断，连续的0码数目不能超过3个，采用plot函数仿真出单极性归零码、单极性不归零码、双极性归零码、双极性不归零码、AMI码、HDB3码对应的波形。
优选的技术方案，所述***间串扰模块采用for循环，生成升余弦滚降函数波形，利用fftshift函数将升余弦滚降函数进行FFT变换，再采用plot函数将升余弦滚降函数仿真出对应的时域和频域波形。
优选的技术方案，所述眼图模块采用randn函数产生随机序列，再用sign函数生成符号函数，然后利用conv函数求出符号函数与两个不同传输速率信号的向量卷积，进而分别仿真出不同传输速率信号的时域波形，以及对应的眼图波形。
优选的技术方案，所述2ASK设计模块采用randn函数随机产生一组二进制序列st1，并四舍五入取值，依旧采用for循环读取二进制序列，当码元小于1时，读取码元为0，否则读取码元为1；载波信号的表达式为S1，f0为S1的频率，其中f0为可变参数，根据S1和f0生成载波，并将载波与序列相乘得到仿真出的已调2ASK信号。
优选的技术方案，所述2ASK解调模块将已调的2ASK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过设置门限高度为0.25的判决器，得到判决之后的波形，最终输出恢复之后的原始波形。
优选的技术方案，所述2FSK设计模块采用randn函数随机产生一组二进制序列，采用for循环读取二进制序列st1，即当码元小于1时，读取码元为0，否则读取码元为1；再用另一个for循环读取二进制序列st1，当码元小于1时，读取码元为1，否则读取码元为0，得到st1序列的反码st2；用st1与其对应的载波信号1相乘得到信号F1，信号F1＝st1*cosω1t，cosω1t为与st1序列对应的载波信号1；用st2与其对应的载波信号2相乘得到信号F2，信号F2＝st2*cosω2t，cosω2t为与st2序列对应的载波信号2；最后将信号F1和信号F2相加，仿真出已调2FSK信号。
优选的技术方案，所述2FSK解调模块将已调的2FSK信号分别送入带通滤波器H1和带通滤波器H2，得到信号h1和h2，将信号h1和h2分别经过相乘器SW1和相乘器SW2，得到信号sw1和sw2，将信号sw1和sw2分别经过低通滤波器ST1和ST2，得到信号stt1和stt2，将信号stt1和stt2相加，得到判决之后的波形，即可恢复输出对应的原始波形。
优选的技术方案，所述2PSK设计模块采用randn函数随机产生一组二进制序列st1，并用round函数四舍五入取值，采用for循环读取二进制序列，当码元小于1时，读取码元为-1，否则读取码元为1；输入频率f0，仿真生成载波，将载波与序列相乘st1，即得到已调2PSK信号；所述2PSK解调模块将已调的2PSK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过门限高度为0.25的判决器，得到判决之后的波形，即可恢复输出对应的原始波形。
本发明的再一目的是提供一种航空数字通信仿真方法，其包括以下步骤：
基本编码设计，***间串扰设计，眼图设计，2ASK设计与解调，2FSK设计与解调，2PSK设计与解调；
所述基本编码设计的过程如下：
采用rand函数随机生成随机序列，并用round函数四舍五入取值，再用for循环对随机序列进行读取，编写出单极性或双极性归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则前t0/2的时间间隔里设置为1，后t0/2的时间间隔里设置为0；如果码元为0，则前t0/2的时间间隔里设置为0(-1)，后t0/2的时间间隔里设置为0；再编写单极性或双极性非归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则整个t0时间间隔里设置为1；如果码元为0，则整个t0时间间隔里设置为0(-1)；编写AMI码，读取随机序列，用for循环分别记录原码中1的数量和0的数量，并对1码进行交替赋值为+1和-1,0则保持不变；在AMI码的基础上，继续采用for循环0的数目判断，连续的0码数目不能超过3个，采用plot函数仿真出单极性归零码、单极性不归零码、双极性归零码、双极性不归零码、AMI码、HDB3码对应的波形；具体公式如下：
基带信号是由rand函数随机生成随机序列，其表达式为：
st1＝10111010100101001011
单极性归零码的编码规则：其中单极性是指用0电平表示“0”码，用正电平表示“1”码，归零码则是指表示每个码元的电波形持续时间小于一个码元周期。其表达式为：st2＝1000101010001000100000100010000010001010
单极性不归零码的编码规则：其中单极性是指用0电平表示“0”码，用正电平表示“1”码，不归零码则是指表示每个码元的电波形持续时间占满一个码元周期。其表达式为：st3＝10111010100101001011
双极性归零码的编码规则：其中双极性是指用负电平表示“0”码，用正电平表示“1”码，归零码则是指表示每个码元的电波形持续时间小于一个码元周期。其表达式为：
st4＝+10-10+10+10+10-10+10-10+10-10-10+10-10+10-10-10+10-10+10+10
双极性不归零码的编码规则：其中双极性是指用负电平表示“0”码，用正电平表示“1”码，不归零码则是指表示每个码元的电波形持续时间占满一个码元周期。其表达式为：
st5＝+1-1+1+1+1-1+1-1+1-1-1+1-1+1-1-1+1-1+1+1
AMI码的编码规则，其中AMI码对应的波形是具有正、负、零三种电平的脉冲序列，消息码的“1”码(传号)交替地变换为“+1”和“-1”，消息码的“0”码(空号)保持不变。其表达式为：
st6＝-10+1-1+10-10+100-10+100-10+1-1
HDB3码为AMI码的改进型，使连“0”个数不超过3个。“0”数目小于等于3时，HDB3码与AMI码一样，+1与-1交替；连“0”数目超过3时，将每4个连“0”化为一节，定义为B00V，称为破坏节；V称为破坏脉冲，而B称为调节脉冲，V与前一个相邻的非“0”脉冲的极性相同(这破坏了极***替的规则，所以V称为破坏脉冲)；相邻的V码之间极性必须时交替，V的取值为+1或-1，V码后面的传号码极性也要交替；B的取值可选0、+1或-1，以使V同时满足要求。其表达式为：
st6＝-10+1-1+10-10+100-10+100-10+1-1
所述***间串扰设计的过程如下：
采用for循环，生成升余弦滚降函数波形，利用fftshift函数将升余弦滚降函数进行FFT变换，再采用plot函数将升余弦滚降函数仿真出对应的时域和频域波形；具体公式如下：
Figure BDA0003415159580000061
Figure BDA0003415159580000062
式中，Ts为发送序列的时间间隔，α为滚降系数，用于描述滚降程度。
其中α＝0时，为理想低通***，α＝1，即为升余弦滚降特性。
所述眼图设计的过程如下：
采用randn函数产生随机序列，再用sign函数生成符号函数，然后利用conv函数求出符号函数与两个不同传输速率信号的向量卷积，进而分别仿真出不同传输速率信号的时域波形，以及对应的眼图波形；其中眼图波形采用eyediagram(x,n)函数，其中，x代表被测信号，每个轨迹包括n个采样点，即n为x信号的周期。
所述2ASK设计与解调的过程如下：
采用randn函数随机产生一组二进制序列st1，并四舍五入取值，依旧采用for循环读取二进制序列，当码元小于1时，读取码元为0，否则读取码元为1；载波信号的表达式为S1，f0为S1的频率，其中f0为可变参数，根据S1和f0生成载波，并将载波与序列相乘得到仿真出的已调2ASK信号；将已调的2ASK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过设置门限高度为0.25的判决器，得到判决之后的波形，最终输出恢复之后的原始波形。具体公式如下：
2ASK的产生：
基带信号形成器
Figure BDA0003415159580000071
式中，g(t)为高度为1，宽度等于Ts的基带矩形脉冲，an为第N个符号的电平，
Figure BDA0003415159580000072
再与载波cosωct相乘，得到2ASK的表达式为e2ASK(t)＝s(t)cosωct
2ASK的解调：
2ASK的解调有两种方法，即相干解调与非相干解调，此处用到的是相干解调：首先经过带通滤波器，得到我们所需要的频段信号，然后与载波cosωct相乘，即
Figure BDA0003415159580000073
由于式中的第二项的频率远大于基带信号的频率，可以用低通滤波器滤掉，即得到
Figure BDA0003415159580000081
最后进行抽样判决器后，即通过设置门限高度为0.25的判决器，得到判决之后的波形s(t)。
所述2FSK设计与解调的过程如下：
采用randn函数随机产生一组二进制序列，采用for循环读取二进制序列st1，即当码元小于1时，读取码元为0，否则读取码元为1；再用另一个for循环读取二进制序列st1，当码元小于1时，读取码元为1，否则读取码元为0，得到st1序列的反码st2；用st1与其对应的载波信号1即cosω1t相乘得到F1，即F1＝st1*cosω1t；用st2与其对应的载波信号2即cosω2t相乘得到F2，即F2＝st2*cosω2t，最后将F1和F2相加，仿真出已调2FSK信号；将已调的2FSK信号分别送入带通滤波器H1和带通滤波器H2，得到信号h1和h2，将信号h1和h2分别经过相乘器SW1和相乘器SW2，得到信号sw1和sw2，将信号sw1和sw2分别经过低通滤波器ST1和ST2，得到信号stt1和stt2，将信号stt1和stt2相加，得到判决之后的波形，即可恢复输出对应的原始波形；具体公式如下：
2FSK的产生：
e2FSK(t)＝s1(t)cosω1t+s2(t)cosω2t
其中
Figure BDA0003415159580000082
g(t)为单个矩形脉冲，Ts为脉冲持续时间，an为第N个符号的电平，
Figure BDA0003415159580000083
为电平取反，即
Figure BDA0003415159580000084
2FSK的解调：
2FSK的解调有两种方法，即相干解调与非相干解调，此处用到的是相干解调：首先经过带通滤波器，得到我们所需要的频段信号，然后分别与载波cosω1t和cosω2t相乘，即
Figure BDA0003415159580000091
和
Figure BDA0003415159580000092
由于式中的第二项的频率远大于基带信号的频率，可以用低通滤波器滤掉，即得到
Figure BDA0003415159580000093
和
Figure BDA0003415159580000094
最后将其相加进行抽样判决器后，得到判决之后的波形s(t)。
所述2PSK设计与解调的过程如下：
采用randn函数随机产生一组二进制序列st1，并用round函数四舍五入取值，采用for循环读取二进制序列，当码元小于1时，读取码元为-1，否则读取码元为1；输入频率f0，仿真生成载波，将载波与序列相乘st1，即得到已调2PSK信号；将已调的2PSK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过门限高度为0.25的判决器，得到判决之后的波形，即可恢复输出对应的原始波形；具体公式如下：
2PSK信号的产生：
基带信号形成器
Figure BDA0003415159580000095
式中，g(t)为高度为1，宽度等于Ts的基带矩形脉冲，an为第N个符号的电平，
Figure BDA0003415159580000096
再与载波cosωct相乘，得到2PSK的表达式为e2PSK(t)＝s(t)cosωct
其中，发送符号“0”时，e2PSK(t)取0相位(an取+1)；
发送符号“1”时，e2PSK(t)取π相位(an取－1)。
2PSK的解调：
2PSK的解调此处用到的是相干解调：首先经过带通滤波器，得到我们所需要的频段信号，然后与载波cosωct相乘，即
Figure BDA0003415159580000101
由于式中的第二项的频率远大于基带信号的频率，可以用低通滤波器滤掉，即得到
Figure BDA0003415159580000102
最后进行抽样判决器后，即通过设置门限高度为0.25的判决器，得到判决之后的波形s(t)。
与现有技术相比，本发明的有益效果是：
本发明知识针对性强、仿真效果好、知识面覆盖全，人机交互界面友好，且运算快的航空数字通信仿真***。该航空数字通信仿真***能针对航空电子企业的设计培训特点，适应于对企业新员工进行基础的数字通信子***培训。本发明替代了实验室的硬件和仪器仪表设备，并完整实现了对航空数字通信***进行了仿真，仿真***完整，具有演示性强，操作简单实用，适用于航空企业新员工培训，并能够保证企业培训效果，为后续项目设计奠定基础。同时，相对于现有技术中的硬件演示，大大节约了成本。
数字通信子***里从数字基带传输***和数字带通传输***两个方面出发，可以实现NRZ/RZ/AMI/HDB3的四种基本的编码方式，其中每一种编码方式的规则都大不相同，在仿真模型的实际传输中，可以根据各种编码的不同优缺点来选择合适的编码方式。本仿真***还可以实现2ASK/2FSK/2PSK的三种方法的数字调制与解调的过程的软件仿真。其中，2ASK信号是利用幅度传达信息，2FSK信号是通过不同的频率传达信息，2PSK信号则是通过相位的不同来传达信息。
附图说明
图1为基本编码设计设计过程示意图；
图2为***间串扰设计过程示意图；
图3为眼图设计过程示意图；
图4为2ASK设计与解调过程示意图；
图5为2FSK设计过程示意图；
图6为2FSK解调过程示意图；
图7为2PSK设计与解调过程示意图。
具体实施方式
本发明一种航空数字通信仿真***，包括数字基带传输***和数字带通***；所述数字基带传输***包括基本编码模块、***间串扰模块、眼图模块；所述数字带通***包括2ASK调制模块、2FSK调制模块、2PSK调制模块；
所述基本编码模块用于对单极性归零码、单极性不归零码、双极性归零码、双极性不归零码、AMI码、HDB3码的编码设计，并仿真出对应编码的波形。基本编码模块采用rand函数随机生成随机序列，并用round函数四舍五入取值，再用for循环对随机序列进行读取，编写出单极性或双极性归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则前t0/2的时间间隔里设置为1，后t0/2的时间间隔里设置为0；如果码元为0，则前t0/2的时间间隔里设置为0(-1)，后t0/2的时间间隔里设置为0；再编写单极性或双极性非归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则整个t0时间间隔里设置为1；如果码元为0，则整个t0时间间隔里设置为0(-1)；编写AMI码，读取随机序列，用for循环分别记录原码中1的数量和0的数量，并对1码进行交替赋值为+1和-1,0则保持不变；在AMI码的基础上，继续采用for循环0的数目判断，连续的0码数目不能超过3个，采用plot函数仿真出单极性归零码、单极性不归零码、双极性归零码、双极性不归零码、AMI码、HDB3码对应的波形。
所述***间串扰模块用于对***间串扰设计，仿真出对应的时域和频率波形。所述***间串扰模块采用for循环，生成升余弦滚降函数波形，利用fftshift函数将升余弦滚降函数进行FFT变换，再采用plot函数将升余弦滚降函数仿真出对应的时域和频域波形。
所述眼图模块用于眼图设计，仿真出不同传输速率信号的时域波形和对应的眼图波形。所述眼图模块采用randn函数产生随机序列，再用sign函数生成符号函数，然后利用conv函数求出符号函数与两个不同传输速率信号的向量卷积，进而分别仿真出不同传输速率信号的时域波形，以及对应的眼图波形。
所述2ASK调制模块包括2ASK设计模块和2ASK解调模块，所述2ASK设计模块用于仿真生成已调2ASK信号，所述2ASK解调模块用于根据已调2ASK信号，恢复输出对应的原始波形。所述2ASK设计模块采用randn函数随机产生一组二进制序列st1，并四舍五入取值，依旧采用for循环读取二进制序列，当码元小于1时，读取码元为0，否则读取码元为1；载波信号的表达式为S1，f0为S1的频率，其中f0为可变参数，根据S1和f0生成载波，并将载波与序列相乘得到仿真出的已调2ASK信号。所述2ASK解调模块将已调的2ASK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过设置门限高度为0.25的判决器，得到判决之后的波形，最终输出恢复之后的原始波形。
所述2FSK调制模块包括2FSK设计模块和2FSK解调模块，所述2FSK设计模块用于仿真生成已调2FSK信号，所述2FSk解调模块用于根据已调2FSK信号，恢复输出对应的原始波形。所述2FSK设计模块采用randn函数随机产生一组二进制序列，采用for循环读取二进制序列st1，即当码元小于1时，读取码元为0，否则读取码元为1；再用另一个for循环读取二进制序列st1，当码元小于1时，读取码元为1，否则读取码元为0，得到st1序列的反码st2；用st1与其对应的载波信号1相乘得到信号F1，信号F1＝st1*cosω1t，cosω1t为与st1序列对应的载波信号1；用st2与其对应的载波信号2相乘得到信号F2，信号F2＝st2*cosω2t，cosω2t为与st2序列对应的载波信号2；最后将信号F1和信号F2相加，仿真出已调2FSK信号。所述2FSK解调模块将已调的2FSK信号分别送入带通滤波器H1和带通滤波器H2，得到信号h1和h2，将信号h1和h2分别经过相乘器SW1和相乘器SW2，得到信号sw1和sw2，将信号sw1和sw2分别经过低通滤波器ST1和ST2，得到信号stt1和stt2，将信号stt1和stt2相加，得到判决之后的波形，即可恢复输出对应的原始波形。
所述2PSK调制模块包括2PSK设计模块和2PSK解调模块，所述2PSK设计模块用于仿真生成已调2PSK信号，所述2PSK解调模块用于根据已调2PSK信号恢复输出对应的原始波形。所述2PSK设计模块采用randn函数随机产生一组二进制序列st1，并用round函数四舍五入取值，采用for循环读取二进制序列，当码元小于1时，读取码元为-1，否则读取码元为1；输入频率f0，仿真生成载波，将载波与序列相乘st1，即得到已调2PSK信号；所述2PSK解调模块将已调的2PSK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过门限高度为0.25的判决器，得到判决之后的波形，即可恢复输出对应的原始波形。
上述航空数字通信仿真***采用以下仿真方法得以实现，如图1至图7所示，其包括以下仿真步骤：
基本编码设计，***间串扰设计，眼图设计，2ASK设计与解调，2FSK设计与解调，2PSK设计与解调；
所述基本编码设计的过程如下：
采用rand函数随机生成随机序列，并用round函数四舍五入取值，再用for循环对随机序列进行读取，编写出单极性或双极性归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则前t0/2的时间间隔里设置为1，后t0/2的时间间隔里设置为0；如果码元为0，则前t0/2的时间间隔里设置为0(-1)，后t0/2的时间间隔里设置为0；再编写单极性或双极性非归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则整个t0时间间隔里设置为1；如果码元为0，则整个t0时间间隔里设置为0(-1)；编写AMI码，读取随机序列，用for循环分别记录原码中1的数量和0的数量，并对1码进行交替赋值为+1和-1,0则保持不变；在AMI码的基础上，继续采用for循环0的数目判断，连续的0码数目不能超过3个，采用plot函数仿真出单极性归零码、单极性不归零码、双极性归零码、双极性不归零码、AMI码、HDB3码对应的波形；具体公式如下：
基带信号是由rand函数随机生成随机序列，其表达式为：
st1＝10111010100101001011
单极性归零码的编码规则：其中单极性是指用0电平表示“0”码，用正电平表示“1”码，归零码则是指表示每个码元的电波形持续时间小于一个码元周期。其表达式为：st2＝1000101010001000100000100010000010001010
单极性不归零码的编码规则：其中单极性是指用0电平表示“0”码，用正电平表示“1”码，不归零码则是指表示每个码元的电波形持续时间占满一个码元周期。其表达式为：st3＝10111010100101001011
双极性归零码的编码规则：其中双极性是指用负电平表示“0”码，用正电平表示“1”码，归零码则是指表示每个码元的电波形持续时间小于一个码元周期。其表达式为：
st4＝+10-10+10+10+10-10+10-10+10-10-10+10-10+10-10-10+10-10+10+10
双极性不归零码的编码规则：其中双极性是指用负电平表示“0”码，用正电平表示“1”码，不归零码则是指表示每个码元的电波形持续时间占满一个码元周期。其表达式为：
st5＝+1-1+1+1+1-1+1-1+1-1-1+1-1+1-1-1+1-1+1+1
AMI码的编码规则，其中AMI码对应的波形是具有正、负、零三种电平的脉冲序列，消息码的“1”码(传号)交替地变换为“+1”和“-1”，消息码的“0”码(空号)保持不变。其表达式为：
st6＝-10+1-1+10-10+100-10+100-10+1-1
HDB3码为AMI码的改进型，使连“0”个数不超过3个。“0”数目小于等于3时，HDB3码与AMI码一样，+1与-1交替；连“0”数目超过3时，将每4个连“0”化为一节，定义为B00V，称为破坏节；V称为破坏脉冲，而B称为调节脉冲，V与前一个相邻的非“0”脉冲的极性相同(这破坏了极***替的规则，所以V称为破坏脉冲)；相邻的V码之间极性必须时交替，V的取值为+1或-1，V码后面的传号码极性也要交替；B的取值可选0、+1或-1，以使V同时满足要求。其表达式为：
st6＝-10+1-1+10-10+100-10+100-10+1-1
所述***间串扰设计的过程如下：
采用for循环，生成升余弦滚降函数波形，利用fftshift函数将升余弦滚降函数进行FFT变换，再采用plot函数将升余弦滚降函数仿真出对应的时域和频域波形；具体公式如下：
Figure BDA0003415159580000161
Figure BDA0003415159580000162
式中，Ts为发送序列的时间间隔，α为滚降系数，用于描述滚降程度。
其中α＝0时，为理想低通***，α＝1，即为升余弦滚降特性。
所述眼图设计的过程如下：
采用randn函数产生随机序列，再用sign函数生成符号函数，然后利用conv函数求出符号函数与两个不同传输速率信号的向量卷积，进而分别仿真出不同传输速率信号的时域波形，以及对应的眼图波形；其中眼图波形采用eyediagram(x,n)函数，其中，x代表被测信号，每个轨迹包括n个采样点，即n为x信号的周期。
所述2ASK设计与解调的过程如下：
采用randn函数随机产生一组二进制序列st1，并四舍五入取值，依旧采用for循环读取二进制序列，当码元小于1时，读取码元为0，否则读取码元为1；载波信号的表达式为S1，f0为S1的频率，其中f0为可变参数，根据S1和f0生成载波，并将载波与序列相乘得到仿真出的已调2ASK信号；将已调的2ASK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过设置门限高度为0.25的判决器，得到判决之后的波形，最终输出恢复之后的原始波形。具体公式如下：
2ASK的产生：
基带信号形成器
Figure BDA0003415159580000171
式中，g(t)为高度为1，宽度等于Ts的基带矩形脉冲，an为第N个符号的电平，
Figure BDA0003415159580000172
再与载波cosωct相乘，得到2ASK的表达式为e2ASK(t)＝s(t)cosωct
2ASK的解调：
2ASK的解调有两种方法，即相干解调与非相干解调，此处用到的是相干解调：首先经过带通滤波器，得到我们所需要的频段信号，然后与载波cosωct相乘，即
Figure BDA0003415159580000173
由于式中的第二项的频率远大于基带信号的频率，可以用低通滤波器滤掉，即得到
Figure BDA0003415159580000174
最后进行抽样判决器后，即通过设置门限高度为0.25的判决器，得到判决之后的波形s(t)。
所述2FSK设计与解调的过程如下：
采用randn函数随机产生一组二进制序列，采用for循环读取二进制序列st1，即当码元小于1时，读取码元为0，否则读取码元为1；再用另一个for循环读取二进制序列st1，当码元小于1时，读取码元为1，否则读取码元为0，得到st1序列的反码st2；用st1与其对应的载波信号1即cosω1t相乘得到F1，即F1＝st1*cosω1t；用st2与其对应的载波信号2即cosω2t相乘得到F2，即F2＝st2*cosω2t，最后将F1和F2相加，仿真出已调2FSK信号；将已调的2FSK信号分别送入带通滤波器H1和带通滤波器H2，得到信号h1和h2，将信号h1和h2分别经过相乘器SW1和相乘器SW2，得到信号sw1和sw2，将信号sw1和sw2分别经过低通滤波器ST1和ST2，得到信号stt1和stt2，将信号stt1和stt2相加，得到判决之后的波形，即可恢复输出对应的原始波形；具体公式如下：
2FSK的产生：
e2FSK(t)＝s1(t)cosω1t+s2(t)cosω2t
其中
Figure BDA0003415159580000181
g(t)为单个矩形脉冲，Ts为脉冲持续时间，an为第N个符号的电平，
Figure BDA0003415159580000182
为电平取反，即
Figure BDA0003415159580000183
2FSK的解调：
2FSK的解调有两种方法，即相干解调与非相干解调，此处用到的是相干解调：首先经过带通滤波器，得到我们所需要的频段信号，然后分别与载波cosω1t和cosω2t相乘，即
Figure BDA0003415159580000184
和
Figure BDA0003415159580000185
由于式中的第二项的频率远大于基带信号的频率，可以用低通滤波器滤掉，即得到
Figure BDA0003415159580000186
和
Figure BDA0003415159580000187
最后将其相加进行抽样判决器后，得到判决之后的波形s(t)。
所述2PSK设计与解调的过程如下：
采用randn函数随机产生一组二进制序列st1，并用round函数四舍五入取值，采用for循环读取二进制序列，当码元小于1时，读取码元为-1，否则读取码元为1；输入频率f0，仿真生成载波，将载波与序列相乘st1，即得到已调2PSK信号；将已调的2PSK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过门限高度为0.25的判决器，得到判决之后的波形，即可恢复输出对应的原始波形；具体公式如下：
2PSK信号的产生：
基带信号形成器
Figure BDA0003415159580000191
式中，g(t)为高度为1，宽度等于Ts的基带矩形脉冲，an为第N个符号的电平，
Figure BDA0003415159580000192
再与载波cosωct相乘，得到2PSK的表达式为e2PSK(t)＝s(t)cosωct
其中，发送符号“0”时，e2PSK(t)取0相位(an取+1)；
发送符号“1”时，e2PSK(t)取π相位(an取－1)。
2PSK的解调：
2PSK的解调此处用到的是相干解调：首先经过带通滤波器，得到我们所需要的频段信号，然后与载波cosωct相乘，即
Figure BDA0003415159580000193
由于式中的第二项的频率远大于基带信号的频率，可以用低通滤波器滤掉，即得到
Figure BDA0003415159580000194
最后进行抽样判决器后，即通过设置门限高度为0.25的判决器，得到判决之后的波形s(t)。
本发明的航空数字通信仿真***的构成和工作原理如下：该航空数字通信仿真***包括数字基带传输***和数字带通传输***。数字基带传输***包括几种基本编码、***间串扰、眼图等。如图1所示，几种基本编码的设计过程：首先设置基带信号码元数为20，采用rand函数随机生成1-20个随机序列，并用round函数四舍五入的取值，再用for循环对随机序列进行读取。先编写单极性(双极性)归零码，设置每一个间隔周期为t0，用for循环，对其从第一个码元开始判断。如果码元为1，则前t0/2的时间间隔里设置为1，后t0/2的时间间隔里设置为0；如果码元为0，则前t0/2的时间间隔里设置为0(-1)，后t0/2的时间间隔里设置为0。再编写单极性(双极性)非归零码，设置每一个间隔周期为t0，用for循环，对其从第一个码元开始判断。如果码元为1，则整个t0时间间隔里设置为1；如果码元为0，则整个t0时间间隔里设置为0(-1)。再编写AMI码，读取随机序列，用for循环分别记录原码中1的数量和0的数量。并对1码进行交替赋值为+1和-1,0则保持不变。在AMI码的基础上，继续采用for循环0的数目判断，其连续的0码数目不能超过3个。最后采用plot函数将以上所有波形画出，并显示在相应坐标轴上。如图2所示，***间串扰设计过程：采用for循环，将三种情况分别设定符合其规定的判断条件，最终生成升余弦滚降函数波形。利用fftshift函数将升余弦滚降函数进行FFT变换。最后，采用plot函数将升余弦滚降函数的时域和频域波形画出。如图3所示，眼图设计过程：首先用randn函数随机产生序列，在用sign函数生成符号函数，再设置两个不同传输速率的信号。然后利用conv函数实现符号函数分别与两个不同传输速率信号的向量卷积。最后分别画出不同传输速率的信号的时域波形和所对应的眼图的波形。如图4所示，2ASK调制设计过程：首先设置基带信号码元数，采用randn函数随机产生一组二进制序列st1，并四舍五入取值，依旧采用for循环读取二进制序列，当码元小于1时，读取码元为0，否则码元为1。设载波信号的表达式为S1，f0则为S1的频率，其中f0设置为可变参数，输入参数f0，***读取参数之后生成载波，并将载波与序列相乘，即得到已调信号。随后将已调的2ASK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过设置门限高度为0.25的判决器，得到判决之后的波形，最终输出恢复之后的原始波形。最后分别用plot函数画出各个阶段的波形图。如图5所示，2FSK调制设计过程：首先设置基带信号码元数，采用randn函数随机产生一组二进制序列，依旧采用for循环读取二进制序列st1，即当码元小于1时，读取码元为0，否则码元为1。紧接着再用一个for循环，采用和上一个for循环相反的判断条件读取二进制序列st1，即当码元小于1时，读取码元为1，否则码元为0，此时得到的便是得到st1序列的反码st2。其中f1和f2设置为可变参数，输入参数f1和f2，***读取参数之后分别生成载波信号1和载波信号2，此时再用st1与载波信号1相乘得到F1，st2与载波信号2相乘得到F2，最后将F1和F2相加，即得到最终的已调信号。如图6所示，2FSK解调设计过程：随后将已调的2FSK信号分别送入带通滤波器H1和带通滤波器H2，得到信号h1和h2，再分别经过相乘器SW1和相乘器SW2，得到信号sw1和sw2，然后分别经过低通滤波器ST1和ST2，得到信号stt1和stt2，最后将stt1和stt2相加，得到判决之后的波形，此时不需要设置判决门限，最终恢复输出原始波形。最后分别用plot函数画出各个阶段的波形图。如图7所示，2PSK调制设计过程：首先设置基带信号码元数，采用randn函数随机产生一组二进制序列st1，并用round函数四舍五入取值，依旧采用for循环读取二进制序列，此时for循环读取的方式和上文2ASK和2FSK的有所不同，当码元小于1时，读取码元为-1，否则码元为1。输入参数f0，***读取参数之后生成载波，并将载波与序列相乘，即得到已调信号。随后将已调的2PSK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过设置门限高度为0.25的判决器，得到判决之后的波形，最终输出恢复之后的原始波形。最后分别用plot函数画出各个阶段的波形图。
以上仅为本发明的较佳实施例而已，并不用以限制本发明，应当指出的是，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种航空数字通信仿真***，其特征在于：其包括数字基带传输***和数字带通***；所述数字基带传输***包括基本编码模块、***间串扰模块、眼图模块；所述数字带通***包括2ASK调制模块、2FSK调制模块、2PSK调制模块；
所述基本编码模块用于对单极性归零码、单极性不归零码、双极性归零码、双极性不归零码、AMI码、HDB3码的编码设计，并仿真出对应编码的波形；所述***间串扰模块用于对***间串扰设计，仿真出对应的时域和频率波形；所述眼图模块用于眼图设计，仿真出不同传输速率信号的时域波形和对应的眼图波形；
所述2ASK调制模块包括2ASK设计模块和2ASK解调模块，所述2ASK设计模块用于仿真生成已调2ASK信号，所述2ASK解调模块用于根据已调2ASK信号，恢复输出对应的原始波形；所述2FSK调制模块包括2FSK设计模块和2FSK解调模块，所述2FSK设计模块用于仿真生成已调2FSK信号，所述2FSk解调模块用于根据已调2FSK信号，恢复输出对应的原始波形；所述2PSK调制模块包括2PSK设计模块和2PSK解调模块，所述2PSK设计模块用于仿真生成已调2PSK信号，所述2PSK解调模块用于根据已调2PSK信号恢复输出对应的原始波形。
2.根据权利要求1所述的一种航空数字通信仿真***，其特征在于：所述基本编码模块采用rand函数随机生成随机序列，并用round函数四舍五入取值，再用for循环对随机序列进行读取，编写出单极性或双极性归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则前t0/2的时间间隔里设置为1，后t0/2的时间间隔里设置为0；如果码元为0，则前t0/2的时间间隔里设置为0(-1)，后t0/2的时间间隔里设置为0；再编写单极性或双极性非归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则整个t0时间间隔里设置为1；如果码元为0，则整个t0时间间隔里设置为0(-1)；编写AMI码，读取随机序列，用for循环分别记录原码中1的数量和0的数量，并对1码进行交替赋值为+1和-1,0则保持不变；在AMI码的基础上，继续采用for循环0的数目判断，连续的0码数目不能超过3个，采用plot函数仿真出单极性归零码、单极性不归零码、双极性归零码、双极性不归零码、AMI码、HDB3码对应的波形。
3.根据权利要求1所述的一种航空数字通信仿真***，其特征在于：所述***间串扰模块采用for循环，生成升余弦滚降函数波形，利用fftshift函数将升余弦滚降函数进行FFT变换，再采用plot函数将升余弦滚降函数仿真出对应的时域和频域波形。
4.根据权利要求1所述的一种航空数字通信仿真***，其特征在于：所述眼图模块采用randn函数产生随机序列，再用sign函数生成符号函数，然后利用conv函数求出符号函数与两个不同传输速率信号的向量卷积，进而分别仿真出不同传输速率信号的时域波形，以及对应的眼图波形。
5.根据权利要求1所述的一种航空数字通信仿真***，其特征在于：所述2ASK设计模块采用randn函数随机产生一组二进制序列st1，并四舍五入取值，依旧采用for循环读取二进制序列，当码元小于1时，读取码元为0，否则读取码元为1；载波信号的表达式为S1，f0为S1的频率，其中f0为可变参数，根据S1和f0生成载波，并将载波与序列相乘得到仿真出的已调2ASK信号。
6.根据权利要求5所述的一种航空数字通信仿真***，其特征在于：所述2ASK解调模块将已调的2ASK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过设置门限高度为0.25的判决器，得到判决之后的波形，最终输出恢复之后的原始波形。
7.根据权利要求1所述的一种航空数字通信仿真***，其特征在于：所述2FSK设计模块采用randn函数随机产生一组二进制序列，采用for循环读取二进制序列st1，即当码元小于1时，读取码元为0，否则读取码元为1；再用另一个for循环读取二进制序列st1，当码元小于1时，读取码元为1，否则读取码元为0，得到st1序列的反码st2；用st1与其对应的载波信号1相乘得到信号F1，即信号F1＝st1*cosω1t，cosω1t为与st1序列对应的载波信号1；用st2与其对应的载波信号2相乘得到信号F2，即信号F2＝st2*cosω2t，cosω2t为与st2序列与对应的载波信号2；最后将信号F1和信号F2相加，仿真出已调2FSK信号。
8.根据权利要求7所述的一种航空数字通信仿真***，其特征在于：所述2FSK解调模块将已调的2FSK信号分别送入带通滤波器H1和带通滤波器H2，得到信号h1和h2，将信号h1和h2分别经过相乘器SW1和相乘器SW2，得到信号sw1和sw2，将信号sw1和sw2分别经过低通滤波器ST1和ST2，得到信号stt1和stt2，将信号stt1和stt2相加，得到判决之后的波形，即可恢复输出对应的原始波形。
9.根据权利要求1所述的一种航空数字通信仿真***，其特征在于：所述2PSK设计模块采用randn函数随机产生一组二进制序列st1，并用round函数四舍五入取值，采用for循环读取二进制序列，当码元小于1时，读取码元为-1，否则读取码元为1；输入频率f0，仿真生成载波，将载波与序列相乘st1，即得到已调2PSK信号；所述2PSK解调模块将已调的2PSK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过门限高度为0.25的判决器，得到判决之后的波形，即可恢复输出对应的原始波形。
10.一种航空数字通信仿真方法，其特征在于：其包括以下步骤：
基本编码设计，***间串扰设计，眼图设计，2ASK设计与解调，2FSK设计与解调，2PSK设计与解调；
所述基本编码设计的过程如下：
采用rand函数随机生成随机序列，并用round函数四舍五入取值，再用for循环对随机序列进行读取，编写出单极性或双极性归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则前t0/2的时间间隔里设置为1，后t0/2的时间间隔里设置为0；如果码元为0，则前t0/2的时间间隔里设置为0(-1)，后t0/2的时间间隔里设置为0；再编写单极性或双极性非归零码，设置每一个间隔周期为t0，用for循环，对随机序列从第一个码元开始判断，如果码元为1，则整个t0时间间隔里设置为1；如果码元为0，则整个t0时间间隔里设置为0(-1)；编写AMI码，读取随机序列，用for循环分别记录原码中1的数量和0的数量，并对1码进行交替赋值为+1和-1,0则保持不变；在AMI码的基础上，继续采用for循环0的数目判断，连续的0码数目不能超过3个，采用plot函数仿真出单极性归零码、单极性不归零码、双极性归零码、双极性不归零码、AMI码、HDB3码对应的波形；具体公式如下：
基带信号是由rand函数随机生成随机序列，其表达式为：
st1＝10111010100101001011
单极性归零码的编码规则：用0电平表示“0”码，用正电平表示“1”码，每个码元的电波形持续时间小于一个码元周期，单极性归零码表达式为：st2＝1000101010001000100000100010000010001010
单极性不归零码的编码规则：用0电平表示“0”码，用正电平表示“1”码，每个码元的电波形持续时间占满一个码元周期，其表达式为：st3＝10111010100101001011
双极性归零码的编码规则：用负电平表示“0”码，用正电平表示“1”码，每个码元的电波形持续时间小于一个码元周期，其表达式为：
st4＝+10-10+10+10+10-10+10-10+10-10-10+10-10+10-10-10+10-10+10+10
双极性不归零码的编码规则：用负电平表示“0”码，用正电平表示“1”码，每个码元的电波形持续时间占满一个码元周期，其表达式为：
st5＝+1-1+1+1+1-1+1-1+1-1-1+1-1+1-1-1+1-1+1+1
AMI码的编码规则：AMI码对应的波形具有正、负、零三种电平的脉冲序列，消息码的“1”码交替地变换为“+1”和“-1”，消息码的“0”码保持不变，其表达式为：
st6＝-10+1-1+10-10+100-10+100-10+1-1
HDB3码的编码规则：HDB3码对应的波形具有正、负、零三种电平的脉冲序列，消息码的“1”码交替地变换为“+1”和“-1”，消息码的“0”码保持不变，且HDB3码连“0”个数不超过3个；“0”数目小于等于3时，+1与-1交替；连“0”数目超过3时，将每4个连“0”化为一节，定义为B00V，称为破坏节；V称为破坏脉冲，B称为调节脉冲，V与前一个相邻的非“0”脉冲的极性相同；相邻的V码之间极性必须交替，V的取值为+1或-1，V码后面的传号码极性也交替；B的取值为0、+1或-1，其表达式为：
st6＝-10+1-1+10-10+100-10+100-10+1-1
所述***间串扰设计的过程如下：
采用for循环，生成升余弦滚降函数波形，利用fftshift函数将升余弦滚降函数进行FFT变换，再采用plot函数将升余弦滚降函数仿真出对应的时域和频域波形；具体公式如下：
Figure FDA0003415159570000061
Figure FDA0003415159570000062
式中，Ts为发送序列的时间间隔，α为滚降系数，用于描述滚降程度；
其中α＝0时，为理想低通***，α＝1，为升余弦滚降特性；
所述眼图设计的过程如下：
采用randn函数产生随机序列，再用sign函数生成符号函数，然后利用conv函数求出符号函数与两个不同传输速率信号的向量卷积，进而分别仿真出不同传输速率信号的时域波形，以及对应的眼图波形；眼图波形采用eyediagram(x,n)函数，其中，x代表被测信号，每个轨迹包括n个采样点，n为被测信号的周期；
所述2ASK设计与解调的过程如下：
采用randn函数随机产生一组二进制序列st1，并四舍五入取值，依旧采用for循环读取二进制序列，当码元小于1时，读取码元为0，否则读取码元为1；载波信号的表达式为S1，f0为S1的频率，其中f0为可变参数，根据S1和f0生成载波，并将载波与序列相乘得到仿真出的已调2ASK信号；将已调的2ASK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过设置门限高度为0.25的判决器，得到判决之后的波形，最终输出恢复之后的原始波形。具体公式如下：
2ASK信号的产生：
基带信号形成器
Figure FDA0003415159570000071
式中，g(t)为高度为1，宽度等于Ts的基带矩形脉冲，an为第N个符号的电平，
Figure FDA0003415159570000072
再与载波cosωct相乘，得到2ASK信号的表达式为e2ASK(t)＝s(t)cosωct
2ASK信号的解调：
首先经过带通滤波器，得到所需要的频段信号，然后与载波cosωct相乘，即
Figure FDA0003415159570000073
通过低通滤波器过滤，即得到
Figure FDA0003415159570000074
最后通过设置门限高度为0.25的判决器进行抽样判决器后，得到判决之后的波形s(t)；
所述2FSK设计与解调的过程如下：
采用randn函数随机产生一组二进制序列，采用for循环读取二进制序列st1，即当码元小于1时，读取码元为0，否则读取码元为1；再用另一个for循环读取二进制序列st1，当码元小于1时，读取码元为1，否则读取码元为0，得到st1序列的反码st2；用st1与其对应的载波信号1相乘得到信号F1，F1＝st1*cosω1t，cosω1t为与st1序列对应的载波信号1；用st2与其对应的载波信号2相乘得到信号F2，F2＝st2*cosω2t，cosω2t为与st2序列对应的载波信号2；最后将信号F1和信号F2相加，仿真出已调2FSK信号；将已调的2FSK信号分别送入带通滤波器H1和带通滤波器H2，得到信号h1和h2，将信号h1和h2分别经过相乘器SW1和相乘器SW2，得到信号sw1和sw2，将信号sw1和sw2分别经过低通滤波器ST1和ST2，得到信号stt1和stt2，将信号stt1和stt2相加，得到判决之后的波形，即可恢复输出对应的原始波形；具体公式如下：
2FSK信号的产生：
e2FSK(t)＝s1(t)cosω1t+s2(t)cosω2t
其中
Figure FDA0003415159570000081
g(t)为单个矩形脉冲，Ts为脉冲持续时间，an为第N个符号的电平，
Figure FDA0003415159570000082
为电平取反，即
Figure FDA0003415159570000083
2FSK信号的解调：
首先经过带通滤波器，得到所需要的频段信号，然后分别与载波cosω1t和cosω2t相乘，即
Figure FDA0003415159570000084
和
Figure FDA0003415159570000085
通过低通滤波器过滤，得到
Figure FDA0003415159570000086
和
Figure FDA0003415159570000087
最后通过判决器进行抽样判决器后，得到判决之后的波形s(t)；所述2PSK设计与解调的过程如下：
采用randn函数随机产生一组二进制序列st1，并用round函数四舍五入取值，采用for循环读取二进制序列，当码元小于1时，读取码元为-1，否则读取码元为1；输入频率f0，仿真生成载波，将载波与序列相乘st1，得到已调2PSK信号；将已调的2PSK信号送入带通滤波器H1得到信号h1，再经过相乘器SW1得到信号sw1，然后经过低通滤波器ST1，得到信号stt1，最后通过门限高度为0.25的判决器，得到判决之后的波形，即可恢复输出对应的原始波形；具体公式如下：
2PSK信号的产生：
基带信号形成器
Figure FDA0003415159570000088
式中，g(t)为高度为1，宽度等于Ts的基带矩形脉冲，an为第N个符号的电平，
Figure FDA0003415159570000089
再与载波cosωct相乘，得到2PSK的表达式为e2PSK(t)＝s(t)cosωct
其中，发送符号“0”时，e2PSK(t)取0相位，an取+1；
发送符号“1”时，e2PSK(t)取π相位，an取－1；
2PSK信号的解调：
首先经过带通滤波器，得到所需要的频段信号，然后与载波cosωct相乘，
Figure FDA0003415159570000091
通过低通滤波器过滤，得到
Figure FDA0003415159570000092
通过设置门限高度为0.25的判决器进行抽样判决后，得到判决之后的波形s(t)。
CN202111543861.0A 2021-12-16 2021-12-16 一种航空数字通信仿真***及其仿真方法 Pending CN114244664A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111543861.0A CN114244664A (zh) 2021-12-16 2021-12-16 一种航空数字通信仿真***及其仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111543861.0A CN114244664A (zh) 2021-12-16 2021-12-16 一种航空数字通信仿真***及其仿真方法

Publications (1)

Publication Number Publication Date
CN114244664A true CN114244664A (zh) 2022-03-25

Family

ID=80757445

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111543861.0A Pending CN114244664A (zh) 2021-12-16 2021-12-16 一种航空数字通信仿真***及其仿真方法

Country Status (1)

Country Link
CN (1) CN114244664A (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN203659297U (zh) * 2013-11-20 2014-06-18 武汉大学 一种教学用通信原理实验箱装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN203659297U (zh) * 2013-11-20 2014-06-18 武汉大学 一种教学用通信原理实验箱装置

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
付国兰;刘晓山;刘正奇;: "基于虚拟仪器的2ASK数字通信***仿真设计", 实验科学与技术, no. 06, pages 1 - 2 *
俞兴明;张建华;: "数字基带脉冲发送滤波器的DSP实现", 苏州市职业大学学报, no. 02, pages 1 - 3 *
周灼荣;聂涛;靳小超;代忠;: "基于FPGA的基带码发生器设计", 现代电子技术, no. 11, pages 1 *
王战强;王大庆;: "高速通信QPSK数字调制器设计与实现", 空间电子技术, no. 06 *
王晓玲;: "曼彻斯特码传输***仿真设计与性能分析", 信息化研究, no. 05 *
王晓玲;胡剑凌;: "数字基带传输***可靠性研究及相关因素分析", 现代电子技术, no. 20, pages 1 - 2 *
赵琳;邵敏敏;: "基于Simulink的数字通信***的建模与仿真", 电脑知识与技术, no. 27 *
邵玉斌: "Matlab／Simulink通信***建模与仿真实例分析", 30 June 2008, 清华大学出版社 *
陈冬英;: "基于综合能力培养的通信原理课程教学改革", 高师理科学刊, no. 08, pages 2 *

Similar Documents

Publication Publication Date Title
Pawar et al. Modulation recognition in continuous phase modulation using approximate entropy
CN107102256A (zh) 用于对pam 发射机中的抖动进行限制的失真测量
CN107231228A (zh) 混合混沌的低截获通信***设计方法
Uluışık et al. A LabVIEW-based analog modulation tool for virtual and real experimentation
CN114244664A (zh) 一种航空数字通信仿真***及其仿真方法
CN108123724A (zh) 一种基于短波窄带波形的通信***
Abdullah et al. FPAA implementation of chaotic modulation based on Nahrain map
Tanyel et al. Preparation of a virtual toolkit for Communication systems
CN108418672A (zh) 一种基于逆时混沌匹配滤波器的语音通信加密方法
Welch et al. commDSK: a tool for teaching modem design and analysis
CN108243134A (zh) 一种qpsk通信中基于方波节拍的超简spwm调制***
Kamil et al. Lorenz-Based chaotic secure communication schemes
CN111385234A (zh) Ofdm雷达通信一体化波形处理方法、装置、终端及存储介质
Gunther et al. Synchronization and demodulation programming projects to accompany a first course on digital communications
CN117728908A (zh) 一种针对于高速信号传输特性的数学建模分析方法
Ainsleigh et al. Modeling exponential signals in a dispersive multipath environment
CN112804179A (zh) 一种多载波调制装置及其调制方法
Gupta et al. Implementation and Analysis of Different Line Coding Schemes using Verilog
Laskov et al. Analysis of Amplitude Modulation and Demodulation in MATLAB Simulink Environment
Radjapov ORGANIZING AND TESTING QPSK (QUADRATURE PHASE SHIFT KEYING) SIGNALS IN TELECOMMUNICATION TECHNOLOGIES SPHERE IN MATLAB PROGRAMME
Wright et al. An inexpensive method to teach hands-on digital communications
Lazović et al. Performance Analysis of an Underwater Acoustic Communication System Based on DCSK Modulation
Khan et al. Enhancing Students’ Lab Experiences using Simulink-based Pre-Labs of Corresponding Hardware-based Labs
Davarian et al. A multiple digital modulator
Kolumban University Ivory Tower: To Stay in or Break Out SDE: A Direct Jump from Scientific Research to Industrial Applications

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination