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

%----------------三徑信道的參數----------------%
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）];????????%?第二徑延時4個點
h3=[zeros（18）?h2（1:end-8）];????????%?第三徑延時8個點


%-----------------產生訓練序列-----------------%
%-訓練符號頻域數據采用802.11a中的長訓練符號數據-%
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）;??????????????????????????%?訓練符號時域數據
preamble1=[preamble1（Nfft-Ncp+1:end）?preamble1];????%?加入16個點的循環前綴

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

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

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????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