瑞利信道仿真
瑞利衰落信道的matlab 仿真
一 瑞利衰落信道简介
瑞利衰落信道（Rayleigh fading channel）是一种无线电信号传播环境的统计模型。这种模型假设信号通过无线信道之后，其信号幅度是随机的，即“衰落”，并且其包络服从瑞利分布。
二 模型的适用
瑞利衰落模型适用于描述建筑物密集的城镇中心地带的无线信道。密集的建筑和其他物体使得无线设备的发射机和接收机之间没有直射路径，而且使得无线信号被衰减、反射、折射、衍射。在曼哈顿的实验证明，当地的无线信道环境确实接近于瑞利衰落。[3] 通过电离层和对流层反射的无线电信道也可以用瑞利衰落来描述，因为大气中存在的各种粒子能够将无线信号大量散射。瑞利衰落属于小尺度的衰落效应，它总是叠加于如阴影、衰减等大尺度衰落效应上。信道衰落的快慢与发射端和接收端的相对运动速度的大小有关。相对运对导致接收信号的多普勒频移。图中所示即为一固定信号通过单径的瑞利衰落信道后，在1秒内的能量波动，这一瑞利衰落信道的多普勒频移最大分别为10Hz 和100Hz ，在GSM1800MHz 的载波频率上，其相应的移动速度分别为约6千米每小时和60千米每小时。特别需要注意的是信号的“深衰落”现象，此时信号能量的衰减达到数千倍，即30～40分贝。
三 性质
多普勒功率普密度
,
四 瑞利衰落信道的仿真
根据上文所 述，瑞利衰落信道可以通过发生实部和虚部都服从独立的高斯分布变量来仿真生成。不过，在有些情况下，研究者只对幅度的波动感兴趣。针对这种情况，有两种方 法可以仿真产生瑞利衰落信道。这两种方法的目的是产生一个信号，有着上文所示的多普勒功率谱或者等效的自相关函数。这个信号就是瑞利衰落信道的冲激响应。
Jakes 模型和clark 模型
本次只以下图所示的模型来仿真单路信号的产生。课本上也有相关的分析。
仿真结果如下：当终端移动速度为30km/h时，瑞利分布的包络为：
当终端移动速度为120km/h时，瑞利分布的包络为：
五 源程序
function [h]=rayleigh(fd,t) %产生瑞利衰落信道
fc=900*10^6; %选取载波频率
v1=30*1000/3600; %移动速度v1=30km/h
c=3*10^8; %定义光速
fd=v1*fc/c; %多普勒频移
ts=1/10000; %信道抽样时间间隔
t=0:ts:1; %生成时间序列
h1=rayleigh(fd,t); %产生信道数据
v2=120*1000/3600; %移动速度v2=120km/h
fd=v2*fc/c; %多普勒频移
h2=rayleigh(fd,t); %产生信道数据
plot(20*log10(abs(h1(1:10000))))
title('v=30km/h时的信道曲线' )
xlabel(' 时间' );ylabel(' 功率' )
plot(20*log10(abs(h2(1:10000))))
title('v=120km/h时的信道曲线' )
xlabel(' 时间' );ylabel(' 功率' )
function [h]=rayleigh(fd,t)
%该程序利用改进的jakes 模型来产生单径的平坦型瑞利衰落信道
%输入变量说明：
% fd：信道的最大多普勒频移 单位Hz
% t :信号的抽样时间序列，抽样间隔单位s
% h为输出的瑞利信道函数，是一个时间函数复序列
N=40; %假设的入射波数目
wm=2*pi*fd;
M=N/4; %每象限的入射波数目即振荡器数目
Tc=zeros(1,length(t)); %信道函数的实部
Ts=zeros(1,length(t)); %信道函数的虚部
P_nor=sqrt(1/M); %归一化功率系
theta=2*pi*rand(1,1)-pi; %区别个条路径的均匀分布随机相位
for n=1:M
%第i 条入射波的入射角
alfa(n)=(2*pi*n-pi+theta)/N;
fi_tc=2*pi*rand(1,1)-pi; %对每个子载波而言在(-pi,pi)之间均匀分布的随机相位
fi_ts=2*pi*rand(1,1)-pi;
Tc=Tc+2*cos(wm*t*cos(alfa(n))+fi_tc);
Ts=Ts+2*cos(wm*t*sin(alfa(n))+fi_ts); %计算冲激响应函数
end ;
h= P_nor*(Tc+j*Ts); %乘归一化功率系数得到传输函数