2022年通信系统课群综合训练与设计 .pdf

武汉理工大学通信系统课群综合训练与设计报告1 目录目录 .1 摘要 .2 Abstract.3 1.设计任务.4 2.实验电路原理分析.5 2.1 增量调制原理.5 2.2 Miller码原理.6 2.3 循环码原理.7 2.4 PSK 原理.9 2.5 衰落信道原理.10 3.实验过程及结果分析.13 5.参考文献.20 6.附录 .21 名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告2 摘要从人类信息交流和通信的演化进程可以清楚地体会信息技术的不断发展性.现代信息技术具有强大的社会功能,已经成为21 世纪推动社会生产力发展和经济增长的重要因素.信息技术在改变社会的产业结构和生产的同时,也对人类的思想观念、思维方式和生活方式产生着重大而深远的影响.MATLAB 是美国 MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB 和 Simulink两大部分。MATLAB 和 Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB 可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。能够在连续时间域、离散时间域或者两者的混合时间域里进行建模，其中通信系统工具箱中包含了对通信系统进行分析和仿真所 需的信源编码、纠错编码、信道、调制解调以及其他所用的库函数和模块。关键词：增量调制，Miller码，循环码，PSK，衰落信道名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告3 Abstract In todays information society,you can clearly feel that the continuous development of information technology,from evolutionary process of human information exchange and communication of modern information technology with a strong social function,the 21st century has become an important factor in promoting the development of social productivity and economic growth.IT in a changing society industrial structure and production,but also on the concept of human thinking,ways of thinking and way of life are significant and far-reaching impact.MATLAB is a U.S.commercial mathematics software for algorithm development,data visualization,data analysis and numerical calculation senior technical computing language and interactive environment,including MATLAB and Simulink two parts.MATLAB and Mathematical,Maple and called three mathematical software.It is second to none in terms of numerical mathematics class technology application software.MATLAB matrix operations can be carried out,drawing functions and data,algorithm,creating the user interface,the connection procedure other programming languages,mainly used in engineering calculations,control design,signal processing and communications,image processing,signal detection,financial modeling design and analysis and other fields.Modeling can be performed in a continuous time domain,the discrete time domain,or a mixture of both in the time domain,wherein the communication system includes a toolbox for the analysis and simulation of the communication system required for the source coding,error correction coding,channel modulation demodulation and other library functions and modules used.Keywords:delta modulation,Miller codes,cyclic codes,PSK fading channel 名师资料总结-精品资料欢迎下载-名师精心整理-第 3 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告4 1.设计任务利用 Matlab 设计完成一个典型的通信系统（如图1 所示）。图 1 典型的通信系统框图其中，首先输入模拟信号，给出此模拟信号的时域波形。数字化：将模拟信号进行数字化，得到数字信号，选择增量调制。信道编码：实现简单的信道编译码Miller码。信道码为循环码。调制方式为 PSK。信道类型为衰落信道。名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告5 2.实验电路原理分析2.1 增量调制原理增量调制简称 M或增量脉码调制方式（DM），它是继 PCM 后出现的又一种模拟信号数字化的方法。1946年由法国工程师 De Loraine 提出，目的在于简化模拟信号的数字化方法。主要在军事通信和卫星通信中广泛使用，有时也作为高速大规模集成电路中的A/D 转换器使用。增量调制最主要的特点就是它所产生的二进制代码表示模拟信号前后两个抽样值的差别(增加、还是减少)而不是代表抽样值本身的大小，因此把它称为增量调制。在增量调制系统的发端调制后的二进制代码1 和 0 只表示信号这一个抽样时刻相对于前一个抽样时刻是增加(用 1 码)还是减少(用 0 码)。收端译码器每收到一个 1 码，译码器的输出相对于前一个时刻的值上升一个量化阶，而收到一个 0 码，译码器的输出相对于前一个时刻的值下降一个量化阶。增量调制(DM)是 DPCM 的一种简化形式。在增量调制方式下,采用 1 比特量化器,即用 1 位二进制码传输样值的增量信息,预测器是一个单位延迟器,延迟一个采样时间间隔。预测滤波器的分子系数向量是0,1,分母系数为 1。当前样值与预测器输出的前一样值相比较,如果其差值大于零,则发 1 码,如果小于零则发0码。增量调制系统框图如图2 所示,其中量化器是一个零值比较器,根据输入的电平极性,输出为,预测器是一个单位延迟器,其输出为前一个采样时刻的解码样值,编码器也是一个零值比较器,若其输入为负值,则编码输出为0,否则输出为 1。解码器将输入 1,0 符号转换为,然后与预测值相加后得出解码样值输出,同时也作为预测器的输入。输入样值enen=sgn(en)传输sgn(en)xn=xn-1+sgn(en)xn+-+预测输出+n-1+预测输出解码样值输出xn-1 预测输入xn=xn-1+sgn(en)图 2 增量调制原理框图解码编码二电平量化单 位延迟单 位延迟名师资料总结-精品资料欢迎下载-名师精心整理-第 5 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告6 2.2 Miller 码原理Miller 码又称延迟调制码，它是双相码的一种变形。其编码规则如下：“1”码用码元中心点出现跃变来表示，即用“10”或“01”表示。“0”码有两种情况：单个零时，在码元持续时间内不出现电平跃变，且与相邻码元的边界处也不跃变，连“0”时，在两个“0”码的边界处出现电平跃变，即“00”与“11”交替出现。密勒码可由双相码的下降沿去触发双稳电路产生。密勒码最初用于气象卫星和磁记录，现在也用于低速基带数传机。密勒码的电平跃变如图所示：图 3 密勒码的电平跃变名师资料总结-精品资料欢迎下载-名师精心整理-第 6 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告7 2.3 循环码原理图 4 总原理方框图循环码的编、译码系统由下列五部分组成：定时单元、信码发生器及显示部分、编码器和译码器。（1）定时单元本单元提供编码器及译码器所需的时序信号。其时钟重复频率（CP）为 2MHZ。（2）信码发生器本单元给编码器提供一个信号源，手控开关（板上CDIN）置于+5V时，发光二极管亮，代表输出“1”码元。若开关置于“0”，代表输出“0”码元。（3）循环码编码器编 码 器 是本 实 验 的 主要 部 分。根 据生 成 多 项 式1)(4569xxxxxxg，采用 5 个异或门和 D触发器组成编码器。在K1信号的控制下，输入6 位信息码元 CDIN，一方面串行输入信道，另一方面通过与门送入除法电路进行计算。第6 位输入码元结束时，K1信号也为零，在 CP脉冲作用下，移位寄存器将计算的结果（CDOUT）送往信道，即在6 位信息码元后附加了 9 位监督码，使码长（n=K+r）为 15。（4）译码方框图及原理介绍图 5 译码器方框图经过信道加错后的信码，在K1信号的作用下，进入6 位移存器，同时另一路进入除法电路进行伴随式计算，当6 个信码全存人移存器时，电子开关置于名师资料总结-精品资料欢迎下载-名师精心整理-第 7 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告8“0”，此时信码保存在移存器中，同时另一路已进入除法器的信码，在CP脉冲的作用下，进入除法电路及正交方程形成网络、大数逻辑判决电路。由于本实验最小码间距离d0=6，故最多能纠正两位错码，若错码个数在2 个以内，该系统能自行纠正，纠正后的信码通过电子开关进入移存器，并在显示信号K3的作用下，若发光二极管亮表示“1”码，不亮表示“0”码。此时译码信号是并行输出至显示部分的，它显示的信号应与信源显示的一一对应（注意此时信道干扰产生的错码只能是 1 个或 2 个）。假如信道中错误个数已超过该码纠错能力（即超过2），那么译码显示与编码显示不能对应。正交方程的定义是：假定最高错误码元为e14，其次 e13，此类推至 e0，即e14 在每个方程中均出现一次，而其它错误码元在4 个方程组中出现一次，正交方程组如下：。，7810140026141151314231112143eeeeSeeeeSeeeeSeeeeS（1）由正交形成网络输出至大数逻辑判决电路，由四个三与非门及四输入与非门组成。该电路输出信号通过与门在K2信号作用下，一方面进入除法电路进行伴随式复位，另一方面进入异或门，与 6 位移位寄存器中相应的信号相加（最高位）。已纠正的码经过与门又存在相应的移位寄存器中。6 位移位寄存器将纠正后的正确码字在显示脉冲K3作用下并行输出。在 2 个错误码元以内本电路能自行校正，译码器显示的码与编码器显示的码应一致。如果错误码元超过2个，译码电路会产生错纠现象，译码显示电路显示的码与编码器不一样，即使一样也是巧合。名师资料总结-精品资料欢迎下载-名师精心整理-第 8 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告9 2.4 PSK 原理相位偏移调制，又称移相键控（PSK，Phase Shift Keying）是一种利用相位差异的信号来传送资料的调制方式。该传送信号必须为正交信号，其基底更须为单位化信号。一个信号所代表的数学公式（2）一般调制信号的改变部份可分为幅度A（ASK用）、相位（PSK用）及频率（FSK用）三种。其中PSK即利用相位差异来产生的调制方式。MPSK 通用的传输符号之公式。（3）PSK又可称 M-PSK或 MPSK，目前有 BPSK、QPSK、16PSK、64PSK等等，常用的只有 QPSK。而 M是代表传送信号的符号（symbol）种类。符号越多，传送的比特数越多，自然在固定时间可传送越多的资料量（bps）。传输量公式。（4）图 6 BPSK、QPSK、8PSK及 16PSK的 BER对 SNR图假设各 MPSK 皆在同一能量下传送，PSK会因为符号种类（M）的提升使比特错误率（Bits Error Rate，BER）快速上升。所以在符号数M大于 16后都由 QAM名师资料总结-精品资料欢迎下载-名师精心整理-第 9 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告10 来执行调制工作。QPSK 如果用格雷码对映的方式，其BER 会和 BPSK 一样。所以目前常用的只有 QPSK。2.5 衰落信道原理瑞利衰落信道（Rayleigh fading channel）是一种无线电信号传播环境的“统计模型(statistical model)”。这种模型假设信号通过无线信道之后，其信号幅度(amplitude)是随机的，即“衰落(fading)”，并且其包络(envelope)服从瑞利分布（Rayleigh distribution）。这一信道模型能够描述由电离层和对流层反射的短波信道，以及建筑物密集的城市环境。瑞利衰落只适用于从发射机到接收机不存在直射信号（LoS，Line of Sight）的情况，否则应使用莱斯衰落信道（Ricean fading channel）作为信道模型。瑞利衰落能有效描述存在能够大量散射无线电信号的障碍物的无线传播环境。若传播环境中存在足够多的散射，则冲激信号到达接收机后表现为大量统计独立的随机变量的叠加，根据中心极限定理，则这一无线信道的冲激响应将是一个高斯过程。如果这一散射信道中不存在主要的信号分量，通常这一条件是指不存在直射信号（LoS），则这一过程的均值为0，且相位服从 0 到 2 的均匀分布。即，信道响应的包络服从瑞利分布。设随机变量R，于是其概率密度函数为：（5）其中。（6）若信道中存在一主要分量，例如直射信号（LoS），则信道响应的包络服从莱斯分布，对应的信道模型为莱斯衰落信道。图 7 最大多普勒频移为10Hz的瑞利衰落信道。名师资料总结-精品资料欢迎下载-名师精心整理-第 10 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告11 图 8 最大多普勒频移为100Hz的瑞利衰落信道。图 9 瑞利衰落信道的自相关函数，其多普勒频移为10Hz。瑞利衰落模型适用于描述建筑物密集的城镇中心地带的无线信道。密集的建筑和其他物体使得无线设备的发射机和接收机之间没有直射路径，而且使得无线信号被衰减、反射、折射、衍射。在曼哈顿的实验证明，当地的无线信道环境确实接近于瑞利衰落。通过电离层和对流层反射的无线电信道也可以用瑞利衰落来描述，因为大气中存在的各种粒子能够将无线信号大量散射。瑞利衰落属于小尺度的衰落效应，它总是叠加于如阴影、衰减等大尺度衰落效应上。信道衰落的快慢与发射端和接收端的相对运动速度的大小有关。相对运动导致接收信号的多普勒频移。图中所示即为一固定信号通过单径的瑞利衰落信道后，在 1 秒内的能量波动，这一瑞利衰落信道的多普勒频移最大分别为10Hz和100Hz，在 GSM1800MHz的载波频率上，其相应的移动速度分别为约6 千米每小时和 60 千米每小时。特别需要注意的是信号的“深衰落”现象，此时信号能量的衰减达到数千倍，即3040 分贝。名师资料总结-精品资料欢迎下载-名师精心整理-第 11 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告12 名师资料总结-精品资料欢迎下载-名师精心整理-第 12 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告13 3.实验过程及结果分析图 3.1 原信号图3.2 增量调制信号名师资料总结-精品资料欢迎下载-名师精心整理-第 13 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告14 图 3.3 miller 编码图 3.4 循环码编码名师资料总结-精品资料欢迎下载-名师精心整理-第 14 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告15 图 3.5 PSK 调制信号图 3.6 PSK 解调信号名师资料总结-精品资料欢迎下载-名师精心整理-第 15 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告16 图 3.7 循环码解码信号图 3.8 miller 码解码信号名师资料总结-精品资料欢迎下载-名师精心整理-第 16 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告17 图 3.9 输出信号与原信号图 3.10 改变参数后的输出信号由图 3.9 和 3.10 的对比可发现，后者还原的更理想，其中3.9 的参数设置为采样间隔 Ts=10-3，量化阶距=0.4，而 3.10 的参数设置为采样间隔Ts=10-4,量化阶距=0.06。结果失真度分别为0.0508 和 0.0024（采用均方误差表示失名师资料总结-精品资料欢迎下载-名师精心整理-第 17 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告18 真度），这也验证了直观上的感觉。由参数可知采样间隔越小，失真度越小，可是量化阶距并非越小越好。由增量调制原理可知，译码器恢复的信号与编码器输入模拟信号的波形近似，但是存在失真。这种失真可称为量化噪声。由于编码、译码时用阶梯波形去近似表示模拟信号波形，阶梯本身的电压突跳会产生失真，这是增量调制的一般量化噪声，不可避免的；若信号上升的斜率超过阶梯波的最大可能斜率，则阶梯波的上升速度赶不上信号的上升速度，就发生了过载量化噪声。图 3.11 出现过载量化噪声如图 3.11 所示，模拟波形的上升速度明显高于阶梯波的上升，产生过载量化噪声，此时参数仅将图3.10 的改为 0,04，失真度就增至为0.0246，这也验证了量化阶距并非越小越好，一个阶梯台阶的斜率k=/Ts，它也就是阶梯波的最大可能斜率，或称为译码器的最大跟踪斜率。所以，为了避免发生过载量化噪声，必须使/Ts 够大，满足信号的斜率不会超过这个值。另一方面，值直接和基本量化噪声的大小有关，若取值太大，势必增大基本量化噪声。所以，用增大 fs，即减小 Ts 的办法，才能保证基本量化噪声和过载量化噪声两者都不超过要求。对于衰落信道而言，多普勒频移为0.09 时，传输的数字信号第17 位出错。多普勒效应就是当发射源与接收体之间存在相对运动时，接收体接收的发射源发射信息的频率与发射源发射信息频率不相同，这种现象称为多普勒效应，接收频率与发射频率之差称为多普勒频移。因而每次出错的地方不固定，符合实际情况。名师资料总结-精品资料欢迎下载-名师精心整理-第 18 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告19 4.总结通过这次课程设计，使我对通信原理和数字信号处理的知识在原有的的基础上得到了很大的巩固和提高，同时开阔了我的知识范围，增强了自己动手的能力以及独立思考问题的本能，为以后适应社会发展做了充分的铺垫。在软件仿真调试阶段，出现很多概念上的错误，通过查阅资料和同学的帮助最终解决，在发现错误、改正错误的过程中，暴露出理论方面的欠缺，还有从书本转到实践上的不足，但是，在实验完成后，这些方面的能力都有了很大的提高。经过这次课程设计的学习，让我收获多多。但在这中间，我也发现了我存在的很多不足。我的动手能力还不够强，当有些实验需要很强的动手能力时我还不能从容应对；我的探索方式还有待改善，当面对一些复杂的实验时我还不能很快很好的完成；我的数据处理能力还得提高，当眼前摆着一大堆复杂数据时我处理的方式及能力还不足，不能用最佳的处理手段使实验误差减小到最小程度总之，课程设计让我收获颇丰，同时也让我发现了自身的不足。在实验课上学得的，我将发挥到其它中去，也将在今后的学习和工作中不断提高、完善；在此间发现的不足，我将努力改善，通过学习、实践等方式不断提高，克服那些不应成为学习、获得知识的障碍。在今后的学习、工作中有更大的收获，在不断地探索中、在无私的学习、奉献中实现自己的人身价值！名师资料总结-精品资料欢迎下载-名师精心整理-第 19 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告20 5.参考文献1 樊昌信.通信原理.国防工业出版社,2006.2 管致中.信号与线性系统.高等教育出版社,2004.3 刘泉.通信电子线路.武汉理工大学出版社,2004.4 陈怀琛.MATLAB 及在电子信息课程中的应用.电子工业出版社,2006.5 郑君里.信号与系统.高等教育出版社,2000.6 薛定宇.基于 MATLAB/Simulink 的系统仿真技术与应用(第 2 版).清华大学出版社,2011.7 沈连丰.通信系统(第 4 版).电子工业出版社,2012 名师资料总结-精品资料欢迎下载-名师精心整理-第 20 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告21 6.附录function y=DMcoding()Ts=1e-3;t=0:Ts:20*Ts;x=sin(2*pi*50*t);delta=0.4;%量化阶距D(1+length(t)=0;%预测器初始状态for k=1:length(t)e(k)=x(k)-D(k);%误差信号e_q(k)=delta*(2*(e(k)=0)-1);%量化器输出D(k+1)=e_q(k)+D(k);%延迟器状态更新y(k)=(e_q(k)0);%编码输出endsubplot(2,1,1);plot(t,x,-o);axis(0 20*Ts-2 2);subplot(2,1,2);stairs(t,y);axis(0 20*Ts-2 2);function Miller=Millercoding(x)m=length(x);Miller=zeros(1,2*m);f=0;d=1;for i=1:m if x(i)=1 Miller(2*i-1)=f;Miller(2*i)=not(f);f=not(f);d=1;elseif d=1 Miller(2*i-1)=f;Miller(2*i)=f;d=0;elsef=not(f);Miller(2*i-1)=f;Miller(2*i)=f;d=0;endendendTs=1e-3;t=0:Ts:41*Ts;stairs(t,Miller);axis(0 41*Ts-2 2);名师资料总结-精品资料欢迎下载-名师精心整理-第 21 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告22 function T=Cycliccoding(Input)syms x%求出所有的生成多项式g=1 0 1 1 1;%选择任意一个作为（7,4）循环码的生成多项式r=4;%监督为数m1=poly2sym(Input);if m1=0 T=0 0 0 0 0 0 0;elsem11=expand(xr*m1);%用xr乘以 m1，相当于对 m1进行左移 r位的操作m2=sym2poly(m11);%将多项式转化为矩阵表示形式len=length(m2);if len7 for i=1:len m2(8-i)=m2(len-i+1);endfor i=1:7-len m2(i)=0;endendQ,R=deconv(m2,g);%求m2除以 g所得的余数R=abs(R);for i=1:length(R)R(i)=mod(R(i),2);end%转化为模 2运算T=R+m2;%T为生成的一个循环码字endTs=1e-3;t=0:Ts:44*Ts;stairs(t,T);axis(0 44*Ts-2 2);function re=PSKcoding(Input)cp=;mod1=;f=4*pi;t=0:2*pi/199:2*pi;for n=1:length(Input)if Input(n)=0 A=zeros(1,200);%每个值 200个点elseif Input(n)=1 A=ones(1,200);endcp=cp A;%s(t),码元宽度 200c=cos(f*t);%载波信号mod1=mod1 c;%与s(t)等长的载波信号变为矩阵形式名师资料总结-精品资料欢迎下载-名师精心整理-第 22 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告23 endfigure(1);subplot(2,1,1);plot(cp);grid on;axis(0 4000-2 2);title(二进制信号序列);cm=;mod=;for n=1:length(Input)if Input(n)=0 B=ones(1,200);c=cos(f*t);%载波信号elseif Input(n)=1 B=ones(1,200);c=cos(f*t+pi);%载波信号endcm=cm B;%码元宽度 200mod=mod c;%与 s(t)等长的载波信号endre=mod;%调制subplot(2,1,2);plot(re);grid on;axis(0 4000-2 2);title(2PSK调制信号);figure(2);subplot(2,1,1);plot(abs(fft(cp);axis(0 4000 0 800);title(原始信号频谱);subplot(2,1,2);plot(abs(fft(re);axis(0 4000 0 800);title(2PSK信号频谱);function y_out=Rayleighchannel(SignalInput)LengthOfSignal=length(SignalInput);%信号长度(最好大于两倍 fc)fm=51;%最大多普勒频移fc=510;%载波频率delay=0 31 71 109 173 251;power=0-1-9-10-15-20;%dBy_in=zeros(1,delay(6)SignalInput;%为时移补零y_out=zeros(1,LengthOfSignal);%用于信号输出for i=1:6 f=1:2*fm-1;%通频带长度y=0.5./(1-(f-fm)/fm).2).(1/2)/pi;%多普勒功率谱（基带）Sf=zeros(1,LengthOfSignal);Sf1=y;%多普勒滤波器的频响Sf(fc-fm+1:fc+fm-1)=y;%基带映射到载波频率上名师资料总结-精品资料欢迎下载-名师精心整理-第 23 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告24 x1=randn(1,LengthOfSignal);x2=randn(1,LengthOfSignal);nc=ifft(fft(x1+i*x2).*sqrt(Sf);%同相分量x3=randn(1,LengthOfSignal);x4=randn(1,LengthOfSignal);ns=ifft(fft(x3+i*x4).*sqrt(Sf);%正交分量r0=(real(nc)+j*real(ns);%瑞利信号r=abs(r0);%瑞利信号幅值y_out=y_out+r.*y_in(delay(6)+1-delay(i):delay(6)+LengthOfSignal-delay(i)*10(power(i)/20);end;figure(3);subplot(2,1,1);plot(SignalInput(delay(6)+1:LengthOfSignal);%去除时延造成的空白信号title(Signal Input);subplot(2,1,2);plot(y_out(delay(6)+1:LengthOfSignal);title(Signal Output);%figure(4);%subplot(2,1,1);%hist(r,256);%title(Amplitude Distribution Of Rayleigh Signal)；瑞利信号的振幅%subplot(2,1,2);%hist(angle(r0);%title(Angle Distribution Of Rayleigh Signal);瑞利信号的相位%figure(5);%plot(Sf1);%title(The Frequency Response of Doppler Filter);多普勒滤波器的频率响应function output=PSKdecoding(Input)t=0:2*pi/199:2*pi;f=4*pi;mod1=;for n=1:length(Input)/200 c=cos(f*t);%载波信号mod1=mod1 c;%与s(t）等长的载波信号变为矩阵形式endre=mod1.*Input;%同步解调figure(4);名师资料总结-精品资料欢迎下载-名师精心整理-第 24 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告25 subplot(2,1,1);plot(re);grid on;axis(0 4000-2 2);title(相乘后信号波形)subplot(2,1,2);plot(abs(fft(re);axis(0 4000 0 400);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(5);freqz(b,a,1000,11025);subplot(2,1,1);axis(0 4000-100 3);title(LPF幅频相频图);jt=filter(b,a,re);figure(6);subplot(2,1,1);plot(jt);grid onaxis(0 4000-2 2);title(经低通滤波器后的信号波形)subplot(2,1,2);plot(abs(fft(jt);axis(0 4000 0 400);title(经低通滤波器后的信号频谱);%抽样判决for m=1:length(Input)if jt(m)=0 jt(m)=0;if mod(m,200)=0 output(fix(m/200)+1)=jt(m-100);endendendfigure(7);subplot(2,1,1);plot(jt);grid onaxis(0 4000-2 2);title(PSK解调后波形)subplot(2,1,2);plot(abs(fft(jt);axis(0 4000 0 400);title(PSk解调后频谱);function re=Cyclicdecoding(Input)er=0 0 0 0 0 0 0;g=1 0 1 1 1;q,r=deconv(Input(1:7),g);名师资料总结-精品资料欢迎下载-名师精心整理-第 25 页，共 26 页 -武汉理工大学通信系统课群综合训练与设计报告26 s(1:7)=abs(r);s(1:7)=mod(s(1:7),2);if s(1:7)=0 0 0 1 0 1 1 er(1:7)=1 0 0 0 0 0 0;elseif s(1:7)=0 0 0 1 1 1 0 er(1:7)=0 1 0 0 0 0 0;else if s(1:7)=0 0 0 0 1 1 1 er(1:7)=0 0 1 0 0 0 0;elseif s(1:7)=0 0 0 1 0 0 0 er(1:7)=0 0 0 1 0 0 0;else if s(1:7)=0 0 0 0 1 0 0 er(1:7)=0 0 0 0 1 0 0;elseif s(1:7)=0 0 0 0 0 1 0 er(1:7)=0 0 0 0 0 1 0;else if s(1:7)=0 0 0 0 0 0 1 er(1:7)=0 0 0 0 0 0 1;endendendendendendend名师资料总结-精品资料欢迎下载-名师精心整理-第 26 页，共 26 页 -