基于MATLAB的无线信道仿真及研究

本设计课题任务的内容为：对OFDM系统无线信道进行研究,利用仿真器进行仿真,研究分析电磁波在该无线信道中的传播和变化规律。 具体要求：(1)在研究无线信道传播理论基础上,分析无线信道传播特性,建立各种衰落信道的结构模型,设计无线信道抽头延迟线模型和Jakes仿真模型。 (2)对路径损耗信道模型进行分析,比较各模型的特点,仿真分析模型误差,提出各种模型的适用环境。 (3)利用Jakes仿真器,对小尺度衰落信道进行计算机仿真,验证平坦衰落和频率选择性衰落信道特性,分析小尺度衰落的各种性能参数。 (4)对OFDM系统进行仿真,通过比较加保护间隔和不加保护间隔系统的误码率,给出OFDM具有独特的抗多径衰落特性。 (5)通过分析移动台移动速度和周围环境对系统误码率、信号包络、多普勒功率谱和传递函数等系统参数的影响,给出小尺度衰落随移动台移动速度和周围环境的变化关系。

para=52;    
fftlen=64; 
noc=53;     
nd=6;       
knd=1;     
ml=2;       
sr=250000; 
br=sr.*ml;  
gilen=16;   


tstp=1/sr/（fftlen+gilen）; % Time resolution
itau=[0];       
dlvl1=[0];     
n0=[6];	       
th1=[0.0];     
itnd1=[1000];         	
now1=1;         
fd=150;         
flat=0;         
itnd0=nd*（fftlen+gilen）*20; 
nloop=1000;  
noe = 0;   
nod = 0;   
eop=0;      
nop=0;     
%画误码的初始、循环 
ebn0_min=1;
ebn0_max=10;
ebn0=ebn0_min:1:ebn0_max
len_ebn=length（ebn0）; 
for j=1:len_ebn;
    ebn0=j;

for iii=1:nloop
   
seridata=rand（1para*nd*ml）>0.5; 

paradata=reshape（seridataparand*ml）; 

[ichqch]=qpskmod（paradataparandml）;
kmod=1/sqrt（2）;
ich=ich.*kmod;
qch=qch.*kmod;

% CE data generation
kndata=zeros（1fftlen）;
kndata0=2.*（rand（152）>0.5）-1;
kndata（2:27）=kndata0（1:26）;
kndata（39:64）=kndata0（27:52）;
ceich=kndata; % CE:BPSK
ceqch=zeros（164）;

%------------- data mapping （DC=0） -----------

[ich1qch1]=crmapping（ichqchfftlennd）;

ich2=[ceich.‘ ich1]; % I-channel transmission data
qch2=[ceqch.‘ qch1]; % Q-channel transmission data

%------------------- IFFT  -------------------

x=ich2+qch2.*i;
y=ifft（x）;
ich3=real（y）;
qch3=imag（y）;

%---------- Gurad interval insertion ---------

fftlen2=fftlen+gilen;
[ich4qch4]= giins（ich3qch3fftlengilennd+1）;

%---------- Attenuation Calculation ----------

spow=sum（ich4.^2+qch4.^2）/nd./para;
attn=0.5*spow*sr/br*10.^（-ebn0/10）;
attn=sqrt（attn）;


%********************** fading channel ****************************** 
%If you would like to simulate performance under fading please remove “*“
%from the following four sentenses
[ifadeqfaderamprcosrsin]=sefade（ich4qch4itaudlvl1th1n0itnd1now1length（ich4）tstpfdflat）;
itnd1 = itnd1+itnd0;  % Updata fading counter
ich4=ifade;
qch4=qfade;

%***************************  Receiver  *****************************
%--------------- AWGN addition --------------- 
[ich5qch5]=comb（ich4qch4attn）;

%----Perfect fading compensation for one path fading ----
%If you would like to simulate performance under perfect compensation please remove “*“
%from the following four sentenses
ifade2=1./ramp.*（rcos（1:）.*ich5+rsin（1:）.*qch5）;
qfade2=1./ramp.*（-rsin（1:）.*ich5+rcos（1:）.*qch5）;
ich5=ifade2;
qch5=qfade2;

%----------- Guard interval removal ----------
[ich6qch6]= girem（ich5qch5fftlen2gilennd+1）;

%------------------  FFT  --------------------
rx=ich6+qch6.*i;
ry=fft（rx）;
ich7=real（ry）;
qch7=imag（ry）;

%-------------- Fading compensation by CE symbol --------------
%
%If you would like to simulate performance under CE-based compensation please remove “*“
%in this area
%

% preparation known CE data
ce=1;
ice0=ich2（:ce）;
qce0=qch2（:ce）;

% taking CE data out of received data
ice1=ich7（:ce）;
qce1=qch7（:ce）;

% calculating reverse rotation 
iv=real（（1./（ice1.^2+qce1.^2））.*（ice0+i.*qce0）.*（ice1-i.*qce1））;
qv=imag（（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        725  2013-03-21 15:19  OFDM程序\AC_PSD.m

     文件        333  2013-03-21 15:20  OFDM程序\delay.m

     文件       1764  2013-03-21 15:19  OFDM程序\fade.m

     文件        621  2013-03-21 15:20  OFDM程序\fade2.m

     文件        648  2013-03-21 15:20  OFDM程序\FadingEnv.m

     文件        515  2013-03-21 15:20  OFDM程序\IR_TF.m

     文件       1042  2013-03-21 15:20  OFDM程序\Jake.m

     文件       2321  2013-03-21 15:21  OFDM程序\main.m

     文件       1073  2013-03-21 15:21  OFDM程序\phase.m

     文件       2743  2013-03-21 15:21  OFDM程序\setPara.m

     目录          0  2008-03-29 14:17  OFDM程序

----------- ---------  ---------- -----  ----

                11785                    11