正弦合成法实现的瑞利多径信道MATLAB源码

用MATLAB实现的瑞利多径信道，基于正弦合成法

function [r]=ray_doppler（fm M dt N）

% Define variables:
%   fm       -- the maximum Doppler frequency in Hz
%   M        -- the number of sinusoids to generate the Rayleigh process M > 16
%   dt       -- one symbol duration in second
%   N        -- the length of the fading sequence in symbol
%
%   Reference： N. Kostov ‘‘Mobile Radio Channels Modeling in MATLAB‘‘
%   Radioengineering Vol 12 Issue 4 Dec 2003
%
%   Last updated and verified by hangasyougo 2009Dec01

T=N*dt-dt;
t=0:dt:T;
c=sqrt（2/M）;     % scaling factor of power
w=2*pi*fm;     % maximum Doppler frequency in rad
x=0;
y=0;

for n=1:M

    alpha=（2*pi*n-pi+（2*pi*rand-pi））/（4*M）;
    ph1=2*pi*rand-pi;
    ph2=2*pi*rand-pi;
    x=x+c*cos（w*t*cos（alpha）+ph1）; % x-axis Gaussian process power is 1
    y=y+c*cos（w*t*sin（alpha）+ph2）; % y-axis Gaussian process power is 1

end

% generate a complex-valued sequence
% its amplitude is Rayleigh distributed
% its angle is uniformly distributed

r=（x+sqrt（-1）*y）/sqrt（2）;  % normalized to the unit power power is -2

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

     文件       1609  2010-04-04 09:51  test_ray_doppler.m

     文件       1075  2010-03-16 16:29  ray_doppler.m

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

                 2684                    2