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

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

尝试％查询GPUdev = parallel.gpu.gpudevice.current;％打印有关找到的GPU的信息流(．．．'检测到GPU (%s， %d多处理器，计算能力%s)\n'，．．．开发人员名称、开发人员多处理器计数、开发人员计算能力）；%包括基于GPU的模拟。dogpu = true;抓住%#好的% GPU不支持或不存在，或并行计万博1manbetx算%工具箱不存在且未获得许可。考虑只CPU的模拟。inp=输入([“***注意：未检测到GPU。”，．．．“继续只用cpu进行模拟?”[Y] / N ']，“年代”）;如果strcmpi (inp,'是') || is空(inp) doGPU = false;别的返回；结尾结尾

% DVB-S.2系统参数subsystemtype =“QPSK 1/2”；％星座和LDPC码率esnodb = 0.75;%每符号能量与噪声PSD之比（dB）numframes = 10;百分比要模拟的帧数maxnumldpciterations = 50;% LDPC解码器迭代DVB = GetParamsdvbs2demo（子系统rtype，esnodb，maxnumldpciterations）;%创建并配置BCH编码器和解码器、调制器、解调器、，％AWGN频道。createSimObjDVBS2Demo;％构造LDPC编码器配置对象encoderCfg = ldpcEncoderConfig (dvb.LDPCParityCheckMatrix);％LDPC解码器配置ldpcPropertyValuePairs = {．．．“MaximumIterationCount”，dvb.ldpcnumiterations，．．．“ParityCheckMatrix”，dvb.ldpcparityCheckmatrix，．．．“决策方法”，“艰难的决定”，．．．“IterationTerminationCondition”，'最大迭代计数'，．．．'opportValue'，“信息部分”};%构造LDPC解码器配置对象decodercfg = ldpcdecoderconfig（dvb.ldpcparitycheckmatrix）;如果DOGPU.％构造基于GPU的LDPC解码器系统对象gpuLDPCDecoder=comm.gpu.LDPCDecoder（ldpcPropertyValuePairs{：}）；结尾％创建错误对象以分析误码率的差异CPU与GPU之间的百分比。数量= comm.ErrorRate;

如果Dogpu架构= 2;别的体系结构= 1;结尾％初始化运行时间结果向量运行时=零（架构，Numframes）;AvgTime = Zeros（1，架构）;%种子用于通道和消息的随机数生成器%这将允许CPU和GPU之间进行公平的误码率比较。％缓存原始随机流以稍后还原。原始\u rs=RandStream.getGlobalStream；rs=RandStream.create('MRG32K3A'，'种子'，25）;randstream.setGlobalstream（RS）;%循环每个处理单元- CPU和GPU为了II = 1：架构%对执行循环进行一些初始设置如果（ii == 1）拱=“CPU”；％使用CPU LDPC解码器别的拱门='GPU'；解码器= gpuLDPCDecoder;%使用GPU LDPC解码器结尾％重置错误率对象重置（BER）;重置随机流重置(rs);%通知用户DVB-S。模拟开始。fprintf（[“\新加坡国立大学”拱'基于LDPC解码器：\ n']）;dels = repmat（“\b”，1，fprintf（'初始化......'））;%主模拟循环。运行numFrames+1次，忽略第一次运行时％帧（具有初始化开销）％ 计算。使用第一次运行进行BER计算。为了rr=1：（numFrames+1）%开始计时ts =抽搐;% ***创建输入消息*** %味精= 0 (encbch。MessageLength, 1);味精(1:dvb.NumInfoBitsPerCodeword) =．．．逻辑（randi（[0 1]，DVB.NuminfobitSperCodeword，1））;％***传输***％bchencOut=encbch（msg）；ldpcentcout=ldpcEncode（bchencOut，encoderCfg）；xlvrOut=intrlv（ldpcentcout，dvb.InterleaveOrder）；modOut=pskModulator（xlvrOut）；*** *噪音*** %chanout = chan（modout）;％***收到***％ydemodOut = pskDemodulator (chanOut);dexlvrOut = deintrlv (demodOut dvb.InterleaveOrder);%使用合适的LDPC解码器。如果比较字符串(拱,“CPU”)ldpcdecOut=逻辑（ldpcDecode）（dexlvrOut、decoderCfg、dvb.LDPCNUMIATIONS、，'决定'，“硬”，“终止”，“马克斯”，'输出格式'，“信息”））;别的ldpcdecOut=解码器（dexlvrOut）；结尾bchdecout = decbch（ldpcdecout）;％*** Compute Ber ***％计算LDPC输出的BER，而不是BCH。ber=ber（逻辑（bchencOut），ldpcdecOut）；%停止计时器Runtime (ii, rr) = toc(ts);%不报告带有初始化开销的第一帧。如果（rr> 1）fprintf（dels）;newcharstodelete = fprintf（'帧％d解码：％.2f sec'，．．．RR-1，运行时（II，RR））;dels = repmat（“\b”，1，newCharsToDelete）；结尾结尾%通过DVB-S.2系统运行帧的结束。将运行时结果报告到命令窗口。流(del);%删除打印出来的最后一行。%计算平均运行时间。不包括坐标系1，因为它％包括一些系统对象初始化时间。avgtime (ii) =意味着(运行时(ii, 2:结束));流('％d帧解码，％.2f sec /框架\ n'，numframes，avgtime（ii））;流('塔错误率：％g \ n'，BER（1））;结尾%体系结构循环

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

如果〜dogpu fprintf（'\ n *** gpu不呈现*** \ n \ n'）;别的%计算系统加速比fprintf（['\ nfull系统仿真使用的速度速度为0.2f次．．．基于gpu的LDPC解码器。\n\n'), avgtime (1) / avgtime (2));结尾