[信息与通信]第7章Matlab在通信仿真中的应用.pptVIP

第七章MATLAB在通信仿真中的应用 7.1数字基带传输仿真模型 SNRdb=5; out_channel=awgn(out_info,SNRdb,measured); * * Data_in= randint(1,Len_B); %产生等概分布的0,1随机序列 alpha=randperm(size(Data_in,2)); InterCode =Data_in(alpha); %输出交织后的编码序列 发送端的交织 接收端的解交织 function Data_derandperm=derandperm(Data_in,alpha) %alpha为发送端产生的交织序列号 Data_derandperm(alpha)=Data_in; %输出交织后的编码序列 一.数字信源信号的产生 Len_frame=1024; BinSignal = randint(1,Len_frame); 二.差错控制编码 1.信道编码 2.交织 图7-2 BPSK、QPSK、16QAM比特编码映射图 三.调制映射 调制映射的基本思想就是将二进制序列，分成长度为n=1、2、4等?的组，然后将这些二进制?序列组分别映射为如图7-2所示BPSK、QPSK、16QAM星座映射图中对应的复数点。 其实就是一种查表的方法。 如图7-2所示， BPSK调制将0映射成1，1映射成-1+0j； QPSK将每2位比特流映射成星座坐标图中的4个点。其中QPSK有GRAY mapping 和 Natural mapping两种映射方式，Gray mpping表示按格雷码方式调制，Natural mapping 表示按自然码方式调制。 16QAM正交振幅调制将每4位比特流映射映射成星座坐标图中的16个点； 在星座图中，每个点的横坐标对应于该点的实部，纵坐标对应该点的虚部，这样通过星座映射将0，1比特流转换成了复数信号流。 如若特流进行Gray mapping 的QPSK调制，则变成1+j，-1-j，1-j，-1＋j?。 图7-2 MPSK星座图 由awgn函数实现加噪 一.信号加噪 7.2信道的仿真 二.信号加多径 Multipath=[1 0 0 0 0 0.707]; %有延迟5个码片,衰落系数为0.707的多径信号 1.莱斯多径的表示方法 2.信号加多径的方法:用滤波函数 BaseSignal_filter = filter(Multipath,1, out_awgninfo); %信号加多径，用滤波 基本思路： 对某一信噪比SNR，测试多帧数据，从而计算在该信噪比下系统总的误比特率（误码率） 给出从SNRdb=0:10,计算出11个BER值，画出SNR－BER关系图 7.3 系统SNR-BER性能仿真 Error = biterr( Data_in,Data_out); %Error为一帧信号中的错误码元个数 Totalerror = Totalerror+Error; %统计总错误码元个数 Totalber=Totalerror/(Nframe*Len_frame); %Totalber统计总误码率 ％Len_frame为一帧数据长度 ％Nframe帧数 ％Totalerror总误比特数 绘图语句（ BPSK抗噪与抗多径性能比较绘图程序） clc, clear all load dataAWGN load dataMultipath SNR0db=0:1:10; semilogy(SNR0db,Totalber_AWGN,r p -, … SNR0db,Totalber_Multipath,b s -); xlabel(SNR); ylabel(BER); title(‘随机交织BPSK-抗噪及抗多径性能) legend(抗噪性能,抗多径性能); grid on; function BPSKCode = mappingBPSK(Xframe) Xframe = Xframe*pi;%将数字转换成相位 [Xp,Yp] = pol2cart(Xframe,ones(1,length(Xframe))); %极坐标转换成直角坐标，极半径为1 BPSKCode = complex(Xp,Yp); %映射成复数 BPSK逆映射目标：相位相比，离0相位近判为0相位，离pi相位近判为pi相位 Conang =[0,pi]; %生成参考相位 Gray=[0 1];