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