clear?all;clc;close?all;
N?=?2048;
%?各種參數(shù)設(shè)定

B?=?2e7;%信號帶寬Hz
N?=?2048;??????%子載波數(shù)目
delta_fp?=?B/N;%子載波間隔
T_ofdm?=?1/delta_fp;%不考慮循環(huán)前綴的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）;%將信噪比轉(zhuǎn)化為線性值
????
????%生成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的全零保護(hù)間隔
????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調(diào)制
????
????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;