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

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

试一试%查询GPUdev = parallel.gpu.GPUDevice.current;打印找到的GPU的信息。%流(.．.'已检测到GPU (%s， %d多处理器，计算能力%s)\n'，.．.dev.Name、dev.MultiprocessorCount dev.ComputeCapability);包含一个基于gpu的模拟。doGPU = true;抓% #好< CTCH >% GPU不支持或不存在，或并行计金宝app算%工具箱不存在并已获得许可。考虑一个仅cpu的模拟。Inp =输入(['***注意:未检测到GPU。'，.．.“继续只使用cpu的模拟?”[Y] / N ']，“年代”）;如果strcmpi (inp,“y”) ||为空(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编码器配置对象encoderCfg = ldpcEncoderConfig(dvb.LDPCParityCheckMatrix);% LDPC解码器配置ldpcPropertyValuePairs = {.．.“MaximumIterationCount”, dvb。LDPCNumIterations,.．.“ParityCheckMatrix”, dvb。LDPCParityCheckMatrix,.．.“DecisionMethod”，“艰难的决定”，.．.“IterationTerminationCondition”，“最大迭代次数”，.．.“OutputValue”，“信息部分”};构造LDPC解码器配置对象decoderCfg = ldpcDecoderConfig(dvb.LDPCParityCheckMatrix);如果doGPU构造一个基于gpu的LDPC解码器系统对象gpuLDPCDecoder = com .gpu. ldpcdecoder (ldpcPropertyValuePairs{:});结束创建一个ErrorRate对象，用于分析误码率的差异CPU和GPU之间的%。BER = com . errorrate;

如果doGPU架构= 2;其他的Architectures = 1;结束初始化运行时结果向量运行时= 0(架构，numFrames);Avgtime = 0(1，架构);为通道和消息使用的随机数生成器提供种子%的创造。这将允许CPU和GPU之间进行公平的误码率比较。缓存原始的随机流，以便稍后恢复。original_rs = RandStream.getGlobalStream;rs = RandStream.create(“mrg32k3a”，“种子”25);RandStream.setGlobalStream (rs);%循环每个处理单元- CPU和GPU为Ii = 1:架构为执行循环做一些初始设置如果(ii == 1) arch =“CPU”；%使用CPU LDPC解码器其他的拱=“图形”；解码器= gpuLDPCDecoder;%使用GPU LDPC解码器结束重置“错误率”对象重置(BER);重置随机流重置(rs);%通知用户DVB-S。2模拟开始。流([“\新加坡国立大学”拱基于LDPC的解码器:\n']);Dels = repmat(“\ b”流(“初始化…”）);主要模拟循环。运行numFrames+1次，忽略第一次% frame(有初始化开销)%计算。使用第一次运行进行BER计算。为rr = 1:(numFrames+1)%开始计时器Ts = tic;% ***创建输入消息*** %MSG = 0 (encbch。MessageLength, 1);味精(1:dvb.NumInfoBitsPerCodeword) =.．.逻辑(randi ([0 1] dvb.NumInfoBitsPerCodeword 1));% ***发送*** %bchencOut = encbch(msg);ldpcencOut = ldpcEncode(bchencOut,encoderCfg);xlvrOut = intrlv(ldpcencOut,dvb.InterleaveOrder);modOut = pskModulator(xlvrOut);% ***腐败噪音*** %chanOut = chan(modOut);% ***接收*** %ydemodOut = pskDemodulator(chanOut);dexlvrOut = deintrlv(demodOut,dvb.InterleaveOrder);使用合适的LDPC解码器。如果比较字符串(拱,“CPU”) ldpcdecOut = logical(ldpcDecode(dexlvrOut,decoderCfg,dvb. cfg)LDPCNumIterations,“DecisionType”，“硬”，“终止”，“马克斯”，“OutputFormat”，“信息”）);其他的ldpcdecOut =解码器(dexlvrOut);结束bchdecOut = decbch(ldpcdecOut);% ***计算误码率*** %计算LDPC输出的误码率，而不是BCH。ber = ber (logical(bchencOut)，ldpcdecOut);%停止计时器Runtime (ii, rr) = toc(ts);不要报告带有初始化开销的第一帧。如果(rr > 1) fprintf(dels);newCharsToDelete = fprintf('帧%d解码:%。2 f交会”，.．.rr-1,运行时(ii, rr));Dels = repmat(“\ b”1、newCharsToDelete);结束结束通过DVB-S运行帧的%结束。2系统。%向命令窗口报告运行时结果。流(del);删除打印出的最后一行。计算平均运行时间。不包括坐标系1，因为它%包括一些系统对象初始化时间。Avgtime (ii) = mean(runtime(ii,2:end));流(' %d帧解码，%。2 f秒/框架\ n '、numFrames avgtime (ii));流(误码率:%g \nber (1));结束%体系结构循环

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

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