Datarate = 20e3;在赫兹％比特率m = 2;％调制阶数（BPSK）ChipValues = [1; 1; 1; 1; 0; 1; 0; 1; 1; 0; 0; 1; 0; 0; 0];％为位0芯片的值。％芯片值1是相反的。

numSymbols = 1000;需要一个EVM值的符号的数量％numFrames = 100;％帧数nSamps = 4;样品的％数表示一个符号filtSpan = 8;％过滤器跨度在符号中增益=长度（chipValues）;％展开增益（每个符号的芯片数量）芯片速率=增益* DATARATE;％片速率采样率= nSamps *芯片速率;％最终采样率加numBits =小区（（numSymbols）/增益）+1;需要一个EVM值的位的数量％SNR = 60;DB中的模拟信噪比％

％映像币值chipValues = 1  -  2个* chipValues;％设计与滚降系数1升余弦滤波器rctfilt = comm.raisedcosinetransmitfilter（'RolloffFactor'，1，......'OutputSamplesPerSymbol'，nSamps，......'filterspaninsymbols'，filtSpan）;rcrFilt = comm.RaisedCosineReceiveFilter（'RolloffFactor'，1，......'InputSamplesPerSymbol'，nSamps，......'filterspaninsymbols'，filtSpan，......'DecimationFactor'，nsamps）;

EVM = comm.EVM（'正常化'那“峰值功率星座”）

％Tx和Rx滤波器的延迟是相同的并且等于一半的过滤器跨度。％总延迟等于两个过滤器延迟的总和，这是过滤器一个过滤器的％跨度。refSigDelay = rctFilt.FilterSpanInSymbols;％模拟帧中的符号数simNumSymbols =加numBits *增益;％初始化峰值RMS EVMpeakRMSEVM = -INF;％创建comm.DifferentialEncoder对象进行差分编码数据diffenc = comm.DifferentialEncoder;％创建comm.awgnchannel系统对象并设置其Noisemethod属性％至“信噪比（SNR）”陈= comm.AWGNChannel（'NoiseMethod'那'信噪比（SNR）'那......'SNR'，SNR）;在阵阵％环为了p = 1：numframes％生成随机数据B =兰迪（[0 M-1]，加numBits，1）;％差分编码d = diffenc（B）;％调制X = 1-2 * d;％将符号转换为芯片（传播）C =重塑（chipValues * X”，simNumSymbols，1）;％脉冲形状杯= rctfilt（c）;％计算并设置信道对象的“SignalPower”属性chan.SignalPower =总和（杯^ 2）/长度（CUP）;％添加噪声R =瓒（CUP）;％缩减像素采样接收到的信号。占滤波器延迟。RD = Rcrfilt（R）;使用EVM系统对象的％测量rmsEVM = EVM（络合物（RD（refSigDelay +（1：numSymbols））），......络合物（C（1：numSymbols）））;％更新峰值RMS EVM计算如果（peakrmsevm 结尾结尾％显示结果fprintf（'最坏情况RMS EVM（%%）：％1.2F \ N'，peakRMSEVM）