LTE链路性能仿真

clear all;clc;close all;
N = 2048;
% 各种参数设定

B = 2e7;%信号带宽Hz
N = 2048;      %子载波数目
delta_fp = B/N;%子载波间隔
T_ofdm = 1/delta_fp;%不考虑循环前缀的OFDM符号周期
Ts = T_ofdm/N;%采样时间

V = 150/9;%移动速度    m/s
c = 3e8;%光速c   m/s
fc = 2e9;%Hz
fd = V*fc/c;%最大多普勒频移Hz
fp = 800;%频偏800Hz

DelayVec=[0 3e-8 1.5e-7 3.1e-7 3.7e-7 7.1e-7 1.09e-6 1.73e-6 2.51e-6];%相对时延s
GainVec=[0 -1.5 -1.4 -3.6 -0.6 -9.1 -7.0 -12.0 -16.9];%相对功率dB


chan = rayleighchan（TsfdDelayVecGainVec）;%瑞利衰落信道    
legacychannelsim（true）;
channel_state = [82140 44470];

snr_dB = 17;

for I = 1:length（snr_dB）
    m1 = 0;
    D = 100000;
    for r = 1:D
    reset（chanchannel_state（1:2）.‘）;%瑞利信道的初始化
    
 	SNR = 10.^（snr_dB（I）./10）;%将信噪比转化为线性值
    
    %生成m序列
    h1 = commsrc.pn（‘GenPoly‘[1 1 0 0 1]‘InitialStates‘[1 1 1 1]‘NumBitsOut‘256）;
    mseq = generate（h1）;
    K = length（mseq）;%PN序列的长度
    G = zeros（120）;%长度为20的全零保护间隔
    mseq_t = [G mseq‘ G];%时间同步序列
    mseq_f = [G mseq‘ mseq‘ mseq‘ mseq‘ mseq‘ mseq‘ G];%频率同步序列

    t_f_seq = [mseq_tmseq_f];
    t_f_bpsk = 1-2*t_f_seq;%BPSK调制
    
    n = [1:length（t_f_bpsk）];
    fp_out = t_f_bpsk.*exp（2*j*pi*fp*n*Ts）;%加800Hz频偏
    
    ray_out = filter（chanfp_out）;%过瑞利信道
    awgn_out = awgn（ray_outSNR‘measured‘‘linear‘）;%过高斯信道
    
%     时间同步
    for i = 1:40
        gama（i） = 0;
        belta（i） = 0;