流(．．.“DVB-S.2数字视频广播标准误码率仿真\n\n'）;

流(．．.CPU加速解码器和GPU加速解码器的性能比较）;

试一试%查询图形处理器dev = parallel.gpu.GPUDevice.current;打印找到的图形处理器的信息流(．．.'检测到GPU (%s， %d多处理器，计算能力%s)\n',．．.dev.Name、dev.MultiprocessorCount dev.ComputeCapability);%包含一个基于gpu的模拟。doGPU = true;抓% #好< CTCH >% GPU不支持或不存在，或并行计金宝app算%Toolbox不存在且未获得许可。考虑一个只使用cpu的模拟。输入=输入([“***注意:未检测到GPU。”,．．.“继续只用cpu进行模拟?”[Y] / N ']，“年代”）;如果strcmpi (inp,“y”) || is空(inp) doGPU = false;其他的返回；结束结束

% DVB-S.2系统参数subsystemType =“QPSK 1/2”；星座和LDPC码率EsNodB = 0.75;每符号的能量与噪声PSD比，单位为dBnumFrames = 10;%要模拟的帧数maxNumLDPCIterations = 50;% LDPC解码器迭代dvb = getParamsDVBS2Demo (subsystemType EsNodB maxNumLDPCIterations);%创建和配置BCH编码器和解码器，调制器，解调器，% AWGN信道。createSimObjDVBS2Demo;%构造LDPC编码器encldpc = comm.LDPCEncoder (dvb.LDPCParityCheckMatrix);%构造LDPC解码器系统对象LDPC解码器配置ldpcPropertyValuePairs = {．．.“MaximumIterationCount”, dvb。LDPCNumIterations,．．.“ParityCheckMatrix”, dvb。LDPCParityCheckMatrix,．．.“DecisionMethod”,“艰难的决定”,．．.“IterationTerminationCondition”,“最大迭代数”,．．.“OutputValue”,“信息部分”}；构建LDPC解码器系统对象(CPU和GPU版本)decldpc = comm.LDPCDecoder (ldpcPropertyValuePairs {:});如果doGPU gpuLDPCDecoder = comm.gpu.LDPCDecoder(ldpcPropertyValuePairs{:});结束%创建ErrorRate对象，分析误码率的差异CPU与GPU之间的百分比。数量= comm.ErrorRate;

如果doGPU架构= 2;其他的架构= 1;结束初始化运行时结果向量运行时= 0(架构、numFrames);avgtime = 0(1、架构);% Seed用于通道和消息的随机数生成器%的创造。这将允许CPU和GPU之间公平的误码率比较。%缓存原始随机流以便稍后恢复。original_rs = RandStream.getGlobalStream;rs = RandStream.create (“mrg32k3a”,“种子”25);RandStream.setGlobalStream (rs);%循环每个处理单元- CPU和GPU为2 = 1:架构为执行循环做一些初始设置如果(ii == 1“CPU”；解码器= decldpc;%使用CPU LDPC解码器其他的拱=“图形”；解码器= gpuLDPCDecoder;%使用GPU LDPC解码器结束%重置错误率对象重置(BER);重置随机流重置(rs);%通知用户DVB-S。模拟开始。流([“\新加坡国立大学”拱“基于LDPC译码器:\ n”]);德尔= repmat (“\ b”流(“初始化……”))；%主模拟循环。运行numFrames+1次并忽略第一个%帧(具有初始化开销)%计算。使用第一次运行进行误码率计算。为rr = 1:(numFrames + 1)%开始计时ts =抽搐;% ***创建输入消息*** %味精= 0 (encbch。MessageLength, 1);味精(1:dvb.NumInfoBitsPerCodeword) =．．.逻辑(randi ([0 1] dvb.NumInfoBitsPerCodeword 1));% * * *发送* * * %bchencOut = encbch(味精);ldpcencOut = encldpc (bchencOut);xlvrOut = intrlv (ldpcencOut dvb.InterleaveOrder);modOut = pskModulator (xlvrOut);*** *噪音*** %陈chanOut = (modOut);% * * *获得* * * %demodOut = pskDemodulator (chanOut);dexlvrOut = deintrlv (demodOut dvb.InterleaveOrder);%使用上面设置的LDPC解码器。ldpcdecOut =解码器(dexlvrOut);bchdecOut = decbch (ldpcdecOut);% ***计算误码率*** %计算LDPC输出的误码率，而不是BCH。数量=误码率(逻辑(bchencOut), ldpcdecOut);%停止计时器Runtime (ii, rr) = toc(ts);%不报告带有初始化开销的第一帧。如果(rr > 1) fprintf(dels);newCharsToDelete =流('帧%d解码:%。2 f交会”,．．.rr-1,运行时(ii, rr));德尔= repmat (“\ b”1、newCharsToDelete);结束结束%结束通过DVB-S运行一帧。2系统。将运行时结果报告到命令窗口。流(del);删除打印出来的最后一行。%计算平均运行时间。不包括坐标系1，因为它%包含一些系统对象初始化时间。avgtime (ii) =意味着(运行时(ii, 2:结束));流(' %d帧解码，%。2 f秒/框架\ n '、numFrames avgtime (ii));流('误码率:%g \n'ber (1));结束%建筑环

将随机流重置为缓存对象RandStream.setGlobalStream (original_rs);

如果~ doGPU流('\n*** GPU不存在***\n '）;其他的%计算系统加速流(['\nFull系统模拟运行%。使用'快2f倍．．.基于gpu的LDPC解码器。\n\n'), avgtime (1) / avgtime (2));结束