通信系统课群综合训练与设计.doc

武汉理工大学学科基础课群综合训练报告通信系统课群综合应用课程设计任务书学生姓名： 专业班级： 指导教师： 工作单位： 信息学院 题 目: 初始条件：理论方面电路分析基础、模拟电子、高频电子线路、通信原理、数字信号处理等专业知识设备方面微型计算机软件方面matlab仿真平台场地方面鉴主15楼通信实验室1要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1. 理论设计与分析，包括理论的推导和计算过程。2. 制作相应的matlab编制仿真程序，并与理论结果对比，最后得出结论。3. 完成课程设计报告，包括必要的公式、文字/图表说明和仿真结果。时间安排：本课程设计自任务书发布之日开始，分散进行，答辩时间统一。指导教师签名： 2016年12月4日系主任（或责任教师）签名： 2016年12月4日目录摘要I1 原理分析与方案论证11.1 通信系统架构11.2 信源调制DM调制21.3 基带码AMI码21.3.1 AMI码编码规则21.3.2 AMI解调原理31.4 信道码汉明码31.5 噪声信道AWGN52 各模块的MATLAB实现52.1 信号源的实现52.2 信源编码DM增量调制52.2.1 DM编码简介52.2.2 DM编码的实现62.3 基带编码AMI编码72.4 信道编码汉明码编码82.4.1 汉明码编码分析82.4.2 汉明码的MATLAB实现92.5 PSK调制102.5.1 PSK调制原理102.5.2 PSK调制的MATLAB实现122.6 信道噪声AWGN132.7 PSK解调142.8 汉明码解码172.8.1 汉明码校验172.8.2 汉明码译码程序182.9 AMI译码192.10 .DM译码193 仿真结果分析203.1 源信号与接收信号波形对比203.2 误码率统计214 心得体会225 参考文献236 附录24本科生课程设计成绩评定表30摘要 通信是通过某种媒体进行的信息传递，目的是传输信息，通信系统是用以完成信息传输过程的技术系统的总称，作用是将信息从信源发送到一个或多个目的地。调制与解调在信息的传输过程中占据着重要的地位，是不可或缺的，因此研究系统的调制和解调过程就极为重要。MATLAB是集数值计算、图形绘制、图像处理及系统仿真等强大功能于一体的科学计算语言，它强大的矩阵运算和图形可视化的功能以及丰富的工具箱，为通信系统的调制和解调过程的分析提供了极大的方便。 调制在通信系统中有十分重要的作用。通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于传播的已调信号，而且它对系统的传输有效性和传输的可靠性有着很大的影响，调制方式往往决定了一个通信系统的性能。MATLAB软件广泛用于数字信号分析，系统识别，时序分析与建模，神经网络、动态仿真等方面有着广泛的应用。本课题利用MATLAB软件完成模拟信号进行DM编码、AMI编码、汉明码、PSK调制、经过AWGN信道，再解调、译码的完整通信系统仿真，并通过统计误码率和对比前后波形，对这个通信系统进行评估。关键字：MATLAB 编码 调制 解调 译码 通信系统Abstract Communication is through some of the information transmission media, the purpose is to transmit information, communication system is used to complete the process of information transmission technology system body, the function is the information from the source sent to one or more of the destination. Modulation and demodulation in information transmission process plays an important role, is indispensable, so the system of modulation and demodulation process is extremely important. MATLAB is set numerical calculation, graphics, image processing and system simulation and powerful function in the integration of scientific computing language, its powerful matrix operation and graphics visual function and rich toolbox, as the communication system of the modulation and demodulation process analysis provides a great convenience. Modulation in a communication system has a very important role. Through the modulation, can not only move the modulation signal spectrum, and the spectrum of the move to want the position, and the modulation signal into suitable for transmission of the modulated signal, and it to the system transmission efficiency and the reliability of the transmission often determines a communication system performance. MATLAB software is widely used in digital signal analysis, system identification, time series analysis and modeling, neural networks, dynamic simulation for such applications. This topic using MATLAB software simulation signal PCM coding, Manchester coding, hamming code and ASK modulation, through AWGN channel, and demodulation, decoding the integrity of the communication system simulation,.Key words：MATLAB，Coding, modulation and demodulation, decoding, communication systemI1 原理分析与方案论证1.1 通信系统架构通信的目的是传输信息。通信系统的作用就是将信息从信息源发送到一个或多个目的地。对于任何一个通信系统，均可视为由发送端、信道和接收端三大部分组成（如图1-1所示）。接收设备发送设备信息源编码设备信 道信息源噪声源发送端接收端信道调制设备译码码设备解调设备图1-1 通信系统一般模型1、信息源:把原始信息变换成原始电信号。2、信源编码：实现模拟信号的数字化传输即完成A/D变化。提高信号传输的有效性。即在保证一定传输质量的情况下，用尽可能少的数字脉冲来表示信源产生的信息。信源编码也称作频带压缩编码或数据压缩编码。3、信道编码：信道编码的目的：信道编码主要解决数字通信的可靠性问题。信道编码的原理：对传输的信息码元按一定的规则加入一些冗余码（监督码），形成新的码字，接收端按照约定好的规律进行检错甚至纠错。信道编码又称为差错控制编码、抗干扰编码、纠错编码 。4、数字调制数字调制技术的概念：把数字基带信号的频谱搬移到高频处，形成适合在信道中传输的频带信号。数字调制的主要作用：提高信号在信道上传输的效率，达到信号远距离传输的目的。基本的数字调制方式：振幅键控ASK、频移键控FSK、相移键控PSK。5、信道：信道是信号传输媒介的总称，传输信道的类型有有线信道（如电缆、光纤）和无线信道（如自由空间）两种。6、噪声源：通信系统中各种设备以及信道中所固有的，为了分析方便，把噪声源视为各处噪声的集中表现而抽象加入到信道。1.2 信源调制-DM调制对模拟信号采样，并用每个样值与它的预测值的差值对周期脉冲序列进行调制，简称M或DM。已调脉冲序列以脉冲的有、无来表征差值的正负号，也就是差值只编成一位二进制码。增量调制的基本原理是于1946年提出的，它是一种最简单的差值脉冲编码。早期的语言增量调制编码器是由分立元件组成的。随着模拟集成电路技术的发展，70年代末出现了音节压扩增量调制集成单片，80年代出现了瞬时压扩集成单片，单片内包括了开关电容滤波器与开关电容积分器，集成度不断提高，使增量调制的编码器的体积减小，功耗降低。1.3 基带码AMI码AMI(Alternative Mark Inversion)码的全称是信号交替反转码，是通信编码中的一种，为极性交替翻转码，1码通常称为传号，0码则叫空号，这是沿用了早期电报通信中的叫法。从形态上看，它已是三状态信号，所以AMI码是伪三进制码。1.3.1 AMI码编码规则 消息代码中的0 传输码中的0消息代码中的1 传输码中的+1、-1交替例如:消息代码:1 0 1 0 1 0 0 0 1 0 1 1 1AMI码: +1 0 -1 0 +1 0 0 0 -1 0 +1 -1 +11.3.2 AMI解调原理 1 由AMI码确定的基带信号中正负脉冲交替，而0电位保持不变;所以由AMI码确定的基带信号无直流分量，且只有很小的低频分量;2 不易提取定时信号，由于它可能出现长的连0串。3 解码规则从收到的符号序列中将所有的-1变换成+1后，就可以得到原消息代码1.4 信道码汉明码汉明码Hamming code：用于数据传送，能检测所有一位和双位差错并纠正所有一位差错的二进制代码。与其他的错误校验码类似，汉明码也利用了奇偶校验位的概念，通过在数据位后面增加一些比特，可以验证数据的有效性。利用一个以上的校验位，汉明码不仅可以验证数据是否有效，还能在数据出错的情况下指明错误位置。在接受端通过纠错译码自动纠正传输中的差错来实现码纠错功能，称为前向纠错FEC。在数据链路中存在大量噪音时，FEC可以增加数据吞吐量。通过在传输码列中加入冗余位(也称纠错位)可以实现前向纠错。但这种方法比简单重传协议的成本要高。汉明码利用奇偶块机制降低了前向纠错的成本。现以数据码1101为例讲讲汉明码的编码原理，此时D8=1、D4=1、D2=0、D1=1，在P1编码时，先将D8、D4、D1的二进制码相加，结果为奇数3，汉明码对奇数结果编码为1，偶数结果为0（奇数位。若奇数结果编码为0.偶数结果为1，则叫偶数位），因此P1值为1，D8+D2+D1=2，为偶数，那么P2值为0，D4+D2+D1=2，为偶数，P3值为0。这样，参照上文的位置表，汉明码处理的结果就是1010101。在这个4位数据码的例子中，我们可以发现每个汉明码都是以三个数据码为基准进行编码的。从编码形式上，我们可以发现汉明码是一个校验很严谨的编码方式。在这个例子中，通过对4个数据位的3个位的3次组合检测来达到具体码位的校验与修正目的（不过只允许一个位出错，两个出错就无法检查出来了，这从下面的纠错例子中就能体现出来）。在校验时则把每个汉明码与各自对应的数据位值相加，如果结果为偶数（纠错代码为0）就是正确，如果为奇数（纠错代码为1）则说明当前汉明码所对应的三个数据位中有错误，此时再通过其他两个汉明码各自的运算来确定具体是哪个位出了问题。还是刚才的1101的例子，正确的编码应该是1010101，如果第三个数据位在传输途中因干扰而变成了1，就成了1010111。检测时，P1+D8+D4+D1的结果是偶数4，第一位纠错代码为0，正确。P2+D8+D2+D1的结果是奇数3，第二位纠错代码为1，有错误。P3+D4+D2+D1的结果是奇数3，第三位纠错代码为1，有错误。那么具体是哪个位有错误呢？三个纠错代码从高到低排列为二进制编码110，换算成十进制就是6，也就是说第6位数据错了，而数据第三位在汉明码编码后的位置正好是第6位。针对4位数据的汉明码编码示意图汉明码是一个在原有数据中插入若干校验码来进行错误检查和纠正的编码技术。以典型的4位数据编码为例，汉明码将加入3个校验码，从而使实际传输的数据位达到7个（位），它们的位置如果把上图中的位置横过来就是：数据位1234567代码P1 P2 D8 P3 D4 D2 D1说明第1个汉明码第2个汉明码第1个数据码第3个汉明码第2个数据码第3个数据码第4个数据码1.5 噪声信道AWGN加性高斯白噪声AWGN(Additive White Gaussian Noise) 是最基本的噪声与干扰模型。加性噪声：叠加在信号上的一种噪声，通常记为n(t)，而且无论有无信号，噪声n(t)都是始终存在的。因此通常称它为加性噪声或者加性干扰。白噪声：噪声的功率谱密度在所有的频率上均为一常数，则称这样的噪声为白噪声。如果白噪声取值的概率分布服从高斯分布，则称这样的噪声为高斯白噪声。2 各模块的MATLAB实现2.1 信号源的实现产生随机信源的方式很多，这里选择在周期远远大于选择域的正余弦组合信源。%信源的产生ts = 1e-3;t = 0:ts:20*ts ;x = sin( 2*pi*50*t )+0.5*sin( 2*pi*150*t );figure(1)subplot(211);plot( t, x, -o );axis( 0 20*ts, -2 2 );title(随机信源 ); 图2-1 信号源时间函数波形2.2 信源编码DM增量调制2.2.1 DM编码简介增量调制DM是继PCM之后推出的一种简单的模拟信号数字化方法，1946年由法国工厂De Loraine提出，又称M。DM可以看成是特殊的DPCM，在DPCM中当对预测误差进行移位二进制码时，就是DM系统。1bit编码输出对应的量化电平是+s或-s，分别表示误差的极性，是正还是负，当抽样频率足够高时，相邻的样值电平差很小，用1位码有可能表示相邻样值的变化规律。一般分为：一般量化和过载量化。过载量化时，信号变化速率大于阶梯波的最大斜率，产生较大的量化失真，所以应该避免。2.2.2 DM编码的实现设计思路是：设定一个差值e（k），表示两个样值之间的差，大于零则返1，反之返0，因为是正余弦组合函数，不存在单调为0的情况。程序如下：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); codeout(k) = ( e_q(k)0 ); 图2-3 DM编码结果2.3 基带编码AMI编码编写实现AMI编码子程序，即对DM编码后的信号进行处理，0保持不变，1用+-1表示交替出现 %AMI编码xn = codeout*1 ; yn = xn;num = 0;for k=1:length(yn) if xn(k) = 1 num = num + 1; if num/2 = fix(num/2) yn(k) = 1; else yn(k) = -1; end endendfigure(2)stairs( 0:length(yn)-1 , yn );axis( 0 length(yn),-2 2 );title(AMI编码); 图2-4 AMI编码2.4 信道编码汉明码编码2.4.1 汉明码编码分析现在以(7,4)分组码为例来说明汉明码的特点。其主要参数如下：码长：信息位：校验位：，且最小距离：其生成矩阵G（前四位为信息位，后三位为冗余位）如下： （式2-1）系统码可分为消息部分和冗余部分两部分，根据生成矩阵，输出码字可按下 式计算： （式2-2）所以有信息位冗余位由以上关系可以得到(7,4)汉明码的全部码字如表2-2所示。表2-2 (7,4)汉明码的全部码字序号信息码元冗余元序号信息码元冗余元000000008100011110001011910011002001010110101001030011110111011001401001101211000015010110113110101060110011141110100701110001511111112.4.2 汉明码的MATLAB实现%汉明码编码HanMing = encode( xn, 7, 4, hamming );figure(3)stairs( 0 : 41 , HanMing );axis(0 41, 0 1.5);title(汉明码编码);汉明码的编码和译码都能由matlab自带的函数库提供，以（7,4）汉明码为例，encode( xn, 7, 4, hamming );图2-5 汉明码编码2.5 PSK调制2.5.1 PSK调制原理在PSK调制时，载波的相位随调制信号状态不同而改变。如果两个频率相同的载波同时开始振荡，这两个频率同时达到正最大值，同时达到零值，同时达到负最大值，此时它们就处于“同相”状态；如果一个达到正最大值时，另一个达到负最大值，则称为“反相”。把信号振荡一次（一周）作为360度。如果一个波比另一个波相差半个周期，两个波的相位差180度，也就是反相。当传输数字信号时，“1”码控制发0度相位，“0”码控制发180度相位。PSK相移键控调制技术在数据传输中，尤其是在中速和中高速的数传机（2400bit/s4800bit/s）中得到了广泛的应用。相移键控有很好的抗干扰性,在有衰落的信道中也能获得很好的效果。主要讨论二相和四相调相，在实际应用中还有八相及十六相调相。PSK也可分为二进制PSK（2PSK或BIT/SK）和多进制PSK（MPSK）。在这种调制技术中，载波相位只有0和两种取值，分别对应于调制信号的“0”和“1”。传“1“信号时，发起始相位为的载波；当传“0”信号时，发起始相位为0的载波。由“0”和“1”表示的二进制调制信号通过电平转换后，变成由“1”和“1”表示的双极性NRZ（不归零）信号，然后与载波相乘，即可形成2PSK信号，在MPSK中，最常用的是四相相移键控，即QPSK（Quadrature Phase Shift Keying），在卫星信道中传送数字电视信号时采用的就是QPSK调制方式。可以看成是由两个2PSK调制器构成的。输入的串行二进制信息序列经串并变换后分成两路速率减半的序列，由电平转换器分别产生双极性二电平信号I（t）和Q（t），然后对载波Acos2fct和Asin2fct进行调制，相加后即可得到QPSK信号。PSK信号也可以用矢量图表示，矢量图中通常以零度载波相位作为参考相位。四相相移调制是利用载波的四种不同相位差来表征输入的数字信息，是四进制移相键控。QPSK是在M=4时的调相技术，它规定了四种载波相位，分别为45，135，225，315。调制器输入的数据是二进制数字序列，为了能和四进制的载波相位配合起来，则需要把二进制数据变换为四进制数据，这就是说需要把二进制数字序列中每两比特分成一组，共有四种组合，即00，01，10，11，其中每一组称为双比特码元。每一个双比特码元是由两位二进制信息比特组成的，它们分别代表四进制四个符号中的一个符号。QPSK中每次调制可传输2个信息比特，这些信息比特是通过载波的四种相位来传递的。解调器根据星座图及接收到的载波信号的相位来判断发送端发送的信息比特。2.5.2 PSK调制的MATLAB实现程序选用2PSK，如下：g=(HanMing); f=2;t=0:0.01:0.99;cp=;sp=;mod=;bit=; for n=1:length(g);if g(n)=0;c1=-ones(1,100);c2=zeros(1,100);else g(n)=1;c1=ones(1,100);c2=ones(1,100); endc=sin(2*pi*f*t);cp=cp c1;mod=mod c; bit=bit c2;end psk=cp.*mod;figure(4)subplot(2,1,1);stairs(bit,LineWidth,1.5);grid on;title(Binary Signal);axis(0 500 -2.5 2.5); subplot(2,1,2);plot(psk,LineWidth,1.5);grid on;title(PSK modulation);axis(0 500 -2.5 2.5); 图2-7 2PSK调制信号2.6 信道噪声AWGNMatlab本身自带了于某信号中加入高斯白噪声的函数：AWGNy = awgn(x,SNR) 在信号x中加入高斯白噪声。SNR为信噪比，以dB为单位。图2-8 加入高斯白噪声的调制波2.7 PSK解调本设计采用相干解调，即采用与载波信号同步的信号与已调信号相乘，再通过FIR低通滤波器，即可解调出调制信号。解调框图如图2-9所示：FIR滤波器乘法器rpskz频率为fc的相干波car 图2-9 PSK解调框图先和载波相乘：repsk=pskz.*mod;%pskz是加噪声后的psk输出，repsk为称值figure(6)subplot(2,2,1);plot(repsk);title(repsk)subplot(2,2,2);plot(pskz);title(pskz) %对比pskz与repsk设计滤波器：wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;N,wc=buttord(wp,ws,Rp,As,s) B,A=butter(N,wc,s); fk=0:800/512:8000;wk=2*pi*fk;Hk=freqs(B,A,wk);figureplot(fk/1000,20*log10(pbs(Hk);axis(0,4,-35,5)滤波器图像： 图2-10 滤波器频谱由于载波频率设置为fc=3khz，所以解调时保留与m（t）相关的低频波，滤除2fc大于4的其他频波。再通过滤波器：fft1=fft(repsk);subplot(2,2,3);plot(abs(fft1)fft1=fft1 zeros(1,921); Ak=fft1.*Hk;at=ifft(Ak);求得的at即是未经判决的粗汉明码，再经过判决程序：hn=;ren=;for k=1:100:length(at) hn(k)=at(k)-1.5ren(k)=(hn(k)0); endfigure(7)plot(ren); 图2-10 PSK解调后波形2.8 汉明码解码2.8.1 汉明码校验 在发送端信息码元M利用式（8-16），实现信道编码，产生线性分组码A；在传输过程中有可能出现误码，设接收到的码组为B。则收发码组之差为：（式2-3） 这里，表示i位有错，表示i位无错。基于这样的原则接收端利用接收到的码组B计算校正子：（式2-4） 因此，校正子仅与E有关，即错误图样与校正子之间有确定的关系。对于上述（7，4）码，校正子S与错误图样的对应关系可由式（2-4）求得，其计算结果见表2-3所示。在接收端的译码器中有专门的校正子计算电路，从而实现检错和纠错。表2-3 （7，4）码校正子与错误图样的对应关系序号错误码位ESe6e5e4e3e2e1e0S3S2S101234567/b0b1b2b3b4b5b6000000000000010000010000010000010000010000010000010000000000010101000111011101112.8.2 汉明码译码程序 汉明码译码可由matlab自带函数库函数实现：%汉明码译码HanMing_decode = decode( HanMing,7 ,4, hamming );figure(8)stairs( 0 : 23 , HanMing_decode );axis(0 23, 0 1.5);title(汉明译码); 图2-11 汉明码解码2.9 AMI译码AMI译码即为其编码的逆过程，取AMI码的绝对值，即可把-1翻1,1和0保留 实现译码。%AMI译码decodeAMI = abs(yn);figure(9)stairs( 0 : 20 , decodeAMI );axis(0 20, 0 1.5);title(AMI);解码后的结果如图2-12所示: 图2-12 AMI解码2.10 .DM译码%DM译码解调ts = 1e-3;t = 0:ts:20*ts ;Dr( 1+length(t) ) = 0;for k=1:length(t) eq(k) = delta * ( 2*decodeAMI(k) - 1 ); xr(k) = eq(k) + Dr(k);endfigure(10)subplot(2,1,1);stairs(t, xr);hold on;subplot(2,1,1);plot(t, x);title(译码源码对比 ); 图2-13 DM译码后的波形3 仿真结果分析 3.1 源信号与接收信号波形对比对比输入的源信号与经过噪声信道后接收到的信号波形，分别如图3-1a和图3-1b所示： 图3-1a 信源信号 图3-1b 接收信号如图可见，复原后的信号大致与源信号相同，但还是出现了变形，这是因为PSK调制方式的抗噪声性能并不是很好，当高斯白噪声同时干扰了一组码中的二个以上的码元时，便会出现误码。3.2 误码率的计算 通过比较DM编码后与AMI译码后的两组数据，统计出不一致的数据的个数，除以总的DM编码个数，即可得到误码率。通过分析误码率的大小可以判断此通信系统的好坏。len=length(codeout);errornum=0;for i=1:len if(codeout(i)=decodeAMI(i) errornum=errornum+1; endend disp(误码率为 );errornum/len4 心得体会 课程设计开始的时候是在准备考研那段期之后，没仔细认真的看设计课题，当时在做MATLAB综合设计性实验的时候比较吃力，觉得使用MATLAB编写程序好难。总是要开好多Mfile文件进行编写函数，最麻烦的是写完一个子函数在单独对其进行验证就存在一些函数输入的问题。当时间过了2天之后，才拿起设计要求认真看起来，想起老师在课堂上进行MATLAB讲解时还是很感兴趣的，于是按照老师的要求开始完成本次课程设计。对通信系统的理解：因为之前有过课程设计的经验，于是刚开始的时候没有急于编写代码，而是仔细认真研究书本上的知识，进行消化和理解。定性了分析了整个通信系统；了解了通信系统研究的意义；通信系统课程设计是培养我们综合运用所学知识,发现、提出、分析和解决实际问题,锻炼实践能力的重要环节,是对学生实际工作能力的具体训练和考察过程。课程设计的收获：本次课程设计结束了，受益良多。它带来不仅是知识上的收获，而且精神上的收获也是无尽的。本次课程设计是我在大学寥寥可数的几次设计性实验之一，是和同学们一起奋斗的一场经历。看到同学们的努力，看到老师的辛勤，自己也都不敢懈怠了。实验开始的那一刻，告诫自己没有失败，老师既然出了设计项目，既然已经选择的这个设计项目，就不许中途退出，这不是游戏，我想那一刻要的是勇气，是决心吧！实验成功的那一刻，收获的是信心，是对自己本学期以来学习通信专业课的一个肯定，是对以后的学习的一种督促，也是对自己的未来大学校门外的生活的一个鼓励。还有就是在设计肯定会有不甚理解的地方，这就需要我们自己动手查找资料，增加的不仅是查找的资料上面的知识，更多的搜索的能力。编程体会：在编写代码的时候，不急于刚开始就对着电脑编写程序，一定要理解和整理出来你所想要做什么；从前都没有画流程图的习惯，不过此次试验之前都有仔细分析，画流程图，发现在代码编写的过程中这样不仅效率高，而且成功率还快；还有就是因为在程序中过多的嵌套过多for循环语句和if选择语句，所以很容易不小心就会导致代码错误，所以程序中应尽量避免。还有就是一定要仔细，仔细，再仔细。5 参考文献1樊昌信，曹丽娜. 通信原理（第六版）M,北京：国防工业出版社，20092程相君，陈生潭 信号与系统 M,西安： 西安电子科技大学出版社,19903曹志刚，钱压生 现代通信原理M,北京： 清华大学出版社，19924刘敏，魏玲 MATLAB通信仿真与应用M，北京：国防工业出版社，20015唐向.岳恒立.学峰 ATLAB及在电子信息类课程中的应用M，北京： 电子工业出版社，2006 6陈杰 MATLAB宝典M，北京：电子工业出版社7张德丰 MATLAB/Simulink建模与仿真M，北京 电子工业出版社 6 附录ts = 1e-3;t = 0:ts:20*ts ;x = sin( 2*pi*50*t )+0.5*sin( 2*pi*150*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);codeout(k) = ( e_q(k)0 );end figure(1)subplot(211);plot( t, x, -o );axis( 0 20*ts, -2 2 );title(随机信源);subplot(212);stairs( t, codeout );axis( 0 20*ts, -2 2 );title(增量编码); %AMI调制xn = codeout*1 ;yn = xn;num = 0;for k=1:length(yn)if xn(k) = 1num = num + 1;if num/2 = fix(num/2)yn(k) = 1;elseyn(k) = -1;endendendfigure(2)stairs( ( 0:length(yn)-1 ), yn );axis( 0 length(yn),-2 2 );title(AMI); %汉明码编码HanMing = encode( xn, 7, 4, hamming );figure(3)stairs( ( 0 : 41 ), HanMing );axis(0 41, 0 1.5);title(汉明编码); %pskg=(HanMing); f=2;t=0:0.01:0.99;cp=;sp=;mod=;bit=; for n=1:length(g);if g(n)=0;c1=-ones(1,100);c2=zeros(1,100);else g(n)=1;c1=ones(1,100);c2=ones(1,100); endc=sin(2*pi*f*t);cp=cp c1;mod=mod c; bit=bit c2;end psk=cp.*mod;figure(4)subplot(2,1,1);stairs(bit,LineWidth,1.5);grid on;title(Binary Signal);axis(0 500 -2.5 2.5); subplot(2,1,2);p