F-OFDM和OFDM调制

s = rng (211);%设置RNG状态可重复性

numFFT = 1024;% FFT点数numRBs = 50;%资源块的数量rbSize = 12;%每个资源块的副载波数量cpLen = 72;%循环前缀长度的样品bitsPerSubCarrier = 6;% 2:正交相移编码,4:16 qam, 6: 64 qam, 8: 256 qamsnrdB = 18;在dB %信噪比toneOffset = 2.5;%的语气抵消或多余的带宽(副载波)L = 513;%滤波器长度(= filterOrder + 1),奇怪

numDataCarriers = numRBs * rbSize;%副环带副载波的数据的数量halfFilt =地板(L / 2);n = -halfFilt: halfFilt;% Sinc函数原型滤波器pb = sinc ((numDataCarriers + 2 * toneOffset)。* n / numFFT);% Sinc截断窗口w = (0.5 * (1 + cos(2 *π。* n / (l - 1))))。^ 0.6;%归一化低通滤波器系数fnum = (pb。* w) / (pb。* w)总和;%滤波器脉冲响应h = fvtool (fnum,“分析”,“冲动”,“Fs”15.36 e6);

%使用dsp过滤器过滤的对象filtTx = dsp.FIRFilter (“结构”,直接形成对称的,…“分子”,fnum);filtRx =克隆(filtTx);Rx %匹配滤波器

%建立图的谱图hFig =图(“位置”figposition ([46 50 30 30]),菜单条的,“没有”);轴(0.5 [-0.5 -200 -20]);持有在;网格在包含(归一化频率的);ylabel (“PSD(瓦分贝/ Hz)”)标题(“F-OFDM”,num2str (numRBs)“资源块”,…num2str (rbSize)每副载波的])%生成数据符号bitsIn =兰迪([0,1],bitsPerSubCarrier * numDataCarriers, 1);% QAM符号映射器symbolsIn = qammod (bitsIn 2 ^ bitsPerSubCarrier,“InputType”,“一点”,…“UnitAveragePower”,真正的);%包数据到一个OFDM符号偏移量= (numFFT-numDataCarriers) / 2;%带中心symbolsInOFDM =[0(抵消,1);symbolsIn;…0 (numFFT-offset-numDataCarriers 1)];ifftOut =传输线(ifftshift (symbolsInOFDM));%预谋循环前缀txSigOFDM = [ifftOut (end-cpLen + 1:结束);ifftOut];%过滤器,平尾补零。得到传输信号txSigFOFDM = filtTx ([txSigOFDM;0 (l - 1, 1)]);%绘制功率谱密度(PSD)[psd f] =周期图(txSigFOFDM rectwin(长度(txSigFOFDM)),…numFFT * 2, 1“中心”);情节(f, 10 * log10 (psd));

%计算peak-to-average-power比(地表铺面)点=功率计(测量=“Peak-to-average功率比”,ComputeCCDF = true);下午paprFOFDM = (txSigFOFDM);disp ([“Peak-to-Average-Power-Ratio F-OFDM = 'num2str (paprFOFDM)“数据库”]);

%绘制功率谱密度(PSD)的OFDM信号[psd f] =周期图(txSigOFDM rectwin(长度(txSigOFDM)), numFFT * 2,…1,“中心”);hFig1 =图(“位置”figposition ([46 15 30 30]));情节(f, 10 * log10 (psd));网格在轴(0.5 [-0.5 -100 -20]);包含(归一化频率的);ylabel (“PSD(瓦分贝/ Hz)”)标题(OFDM,num2str (numRBs * rbSize)副载波的])

%计算peak-to-average-power比(地表铺面)点=功率计(测量=“Peak-to-average功率比”,ComputeCCDF = true);下午paprOFDM = (txSigOFDM);disp ([对OFDM的Peak-to-Average-Power-Ratio = 'num2str (paprOFDM)“数据库”]);

%添加WGNrxSig = awgn (txSigFOFDM snrdB,“测量”);

%接收匹配滤波器rxSigFilt = filtRx (rxSig);%占过滤器延迟rxSigFiltSync = rxSigFilt (L:结束);%去除循环前缀rxSymbol = rxSigFiltSync (cpLen + 1:结束);%执行FFTRxSymbols = fftshift (fft (rxSymbol));%选择数据副载波dataRxSymbols = RxSymbols(抵消+ (1:numDataCarriers));%的阴谋收到星座符号开关bitsPerSubCarrier情况下2%正交相移编码refConst = qammod ((0:3)。“4“UnitAveragePower”,真正的);情况下4% 16 qamrefConst = qammod ((0:15)。”,16日“UnitAveragePower”,真正的);情况下6% 64 qamrefConst = qammod ((0:63)。”,64年,“UnitAveragePower”,真正的);情况下8% 256 qamrefConst = qammod ((0:255)。”,256年,“UnitAveragePower”,真正的);结束constDiagRx = comm.ConstellationDiagram (…“ShowReferenceConstellation”,真的,…“ReferenceConstellation”refConst,…“位置”figposition ([20 15 30 40]),…“EnableMeasurements”,真的,…“MeasurementInterval”、长度(dataRxSymbols),…“标题”,“F-OFDM解调的符号”,…“名字”,“F-OFDM接待”,…“XLimits”(-1.5 - 1.5),“YLimits”[-1.5 - 1.5]);constDiagRx (dataRxSymbols);

%通道均衡没有必要在这里没有信道建模%的误码率计算数量= comm.ErrorRate;%执行困难的决定和测量错误rxBits = qamdemod (dataRxSymbols 2 ^ bitsPerSubCarrier,“OutputType”,“一点”,…“UnitAveragePower”,真正的);数量=误码率(bitsIn rxBits);disp ([“F-OFDM接待,误码率= 'num2str (ber) (1)“在信噪比= '…num2str (snrdB)“数据库”]);

%恢复RNG状态rng(年代);