clear?all;
close?all;
clc;
%----------------?參數(shù)設(shè)置部分----------------%
Nsp=52;?????????????%系統(tǒng)子載波數(shù)（不包括直流載波）
Nfft=64;????????????%?FFT?長(zhǎng)度
Ncp=16;?????????????%?循環(huán)前綴長(zhǎng)度
Ns=Nfft+Ncp;????????%?1個(gè)完整OFDM符號(hào)長(zhǎng)度
noc=53;?????????????%?包含直流載波的總的子載波數(shù)
Nd=6;???????????????%?每幀包含的OFDM符號(hào)數(shù)（不包括訓(xùn)練符號(hào)）
M1=4;???????????????%?QPSK調(diào)制
sr=250000;??????????%?OFDM符號(hào)速率
SNR=20;?????????	%?信噪比
ts=1/sr/Ns;????????????????????%?OFDM符號(hào)抽樣時(shí)間間隔
t=0:ts:（Ns*（Nd+1）-1）*ts;???????%?抽樣時(shí)刻
fd=100;????????????????????????%?最大多普勒頻移

%----------------三徑信道的參數(shù)----------------%
h=rayleigh（fdt）;???????????????????%?生成單徑Rayleigh衰落信道
h1=sqrt（1/2）*h;?????????????????????%?第一徑的功率是總功率的1/2
h2=sqrt（1/3）*h;?????????????????????%?第二徑的功率是總功率的1/3
h3=sqrt（1/6）*h;?????????????????????%?第三徑的功率是總功率的1/6
h2=[zeros（14）?h2（1:end-4）];????????%?第二徑延時(shí)4個(gè)點(diǎn)
h3=[zeros（18）?h2（1:end-8）];????????%?第三徑延時(shí)8個(gè)點(diǎn)


%-----------------產(chǎn)生訓(xùn)練序列-----------------%
%-訓(xùn)練符號(hào)頻域數(shù)據(jù)采用802.11a中的長(zhǎng)訓(xùn)練符號(hào)數(shù)據(jù)-%
Preamble=[1?1?-1?-1?1?1?-1?1?-1?1?1?1?1?1?1?-1?-1?1?1?-1?1?-1?1?1?1?1?...
????1?-1?-1?1?1?-1?1?-1?1?-1?-1?-1?-1?-1?1?1?-1?-1?1?-1?1?-1?1?1?1?1];
Preamble1=zeros（1Nfft）;????????????????????????????%?1X64全零矩陣
Preamble1（2:27）=Preamble（27:end）;???????????????????%?交織
Preamble1（39:end）=Preamble（1:26）;
preamble1=ifft（Preamble1）;??????????????????????????%?訓(xùn)練符號(hào)時(shí)域數(shù)據(jù)
preamble1=[preamble1（Nfft-Ncp+1:end）?preamble1];????%?加入16個(gè)點(diǎn)的循環(huán)前綴

?%-------------------發(fā)射機(jī)部分-------------------%
?msg1=randsrc（NspNd[0:M1-1]）;?????????%?QPSK信息數(shù)據(jù)（52X6的矩陣）
?
?data1=qammod（msg1M1）/sqrt（2）;?????????%?QPSK調(diào)制并歸一化
?
?data2=zeros（NfftNd）;??????????????????%?根據(jù)FFT要求，對(duì)數(shù)據(jù)重排（64X6的矩陣）
?
?data2（2:27:）=data1（27:end:）;
?data2（39:end:）=data1（1:26:）;
?
?data2=ifft（data2）;??????????????????????????%?IFFT變換
?
?data2=[data2（Nfft-Ncp+1:end:）;data2];??????%?加入循環(huán)前綴
?
?spow1=norm（data2‘fro‘）.^2/（Nsp*Nd）;????????%?計(jì)算符號(hào)能量
?
?%下面進(jìn)行的是加入導(dǎo)頻的工作，加入導(dǎo)頻后，每一幀含有7個(gè)符號(hào)，每個(gè)符號(hào)中由于
?%加入了循環(huán)前綴，因此含有80個(gè)點(diǎn)
?data3=zeros（Ns（Nd+1））;????????????????%?加入訓(xùn)練符號(hào)（80X7的矩陣）

data3（:1）=preamble1.‘;?%?在每一幀的開頭加入導(dǎo)頻序列
data3（:2:（Nd+1））=data2（:1:Nd）;
%?在導(dǎo)頻序列之后加入6個(gè)符號(hào)構(gòu)成一幀
?
?data3=reshape（data31Ns*（Nd+1））;??????%?并串變換（1X560的矩陣）
?
?data31=zeros（1length（data3））;
?data32=zeros（1length（data3））;
?data31（5:end）=data3（1:end-4）;???????????????%?第二徑接收到的符號(hào)
?data32（9:end）=data3（1:end-8）;???????????????%?第三徑接收到的符號(hào)
?
?sigma1=sqrt（1/2*spow1/log2（M1）*10.^（-SNR/10））;?%?根據(jù)SNR計(jì)算噪聲標(biāo)準(zhǔn)差
?
?dd1=data3（1:Ns*（Nd+1））;?%取出第一徑每一幀的數(shù)據(jù)
?dd2=data31（1:Ns*（Nd+1））;%取出第二徑每一幀的數(shù)據(jù)
?dd3=data32（1:Ns*（Nd+1））;%取出第三徑每一幀的數(shù)據(jù)
?
?hh1=h1（1:Ns*（Nd+1））;????%?當(dāng)前幀的3徑信道參數(shù)
?hh2=h2（1:Ns*（Nd+1））;
?hh3=h3（1:Ns*（Nd+1））;
?
?%?信號(hào)通過3徑衰落信道，并加入高斯白噪聲
?r1=hh1.*dd1+hh2.*dd2+hh3.*dd3+sigma1*（randn（1length（dd1））+j*randn（1length（dd1）））;
?
?%-------------------接收機(jī)部分-------------------%
?r1=reshape（r1NsNd+1）;?????????????%?串并變換（80X7的矩陣）
?
?r1=r1（Ncp+1:end:）;?????????????????%?移除循環(huán)前綴
?
?R1=fft（r1）;?????????????????????????%?fft運(yùn)算
?
?R1=[R1（39:end:）;R1（2:27:）];???????%?數(shù)據(jù)重排（解交織）
?
?HH1=（Preamble.‘）./R1（:1）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-06-15?15:34??代碼\
?????文件??????705654??2018-06-15?15:34??代碼\1.bmp
?????文件??????705654??2018-06-15?15:34??代碼\2.bmp
?????文件??????705654??2018-06-15?15:34??代碼\3.bmp
?????文件????????4461??2018-06-15?15:33??代碼\OFDM.m
?????文件????????1009??2017-05-11?10:07??代碼\rayleigh.m
?????文件??????????91??2018-06-15?15:36??代碼\說明.txt