📄 rayleigh_fading.mht
字号:
From: <Saved by Windows Internet Explorer 7>
Subject:
Date: Sat, 21 Jun 2008 15:54:18 +0330
MIME-Version: 1.0
Content-Type: text/html;
charset="Windows-1252"
Content-Transfer-Encoding: quoted-printable
Content-Location: http://www.proxoholic.com/index.php?q=aHR0cDovL3d3dy5tYXRod29ya3MuY29tL21hdGxhYmNlbnRyYWwvZmlsZXMvOTIzOC9SYXlsZWlnaF9mYWRpbmcubQ%3D%3D
X-MimeOLE: Produced By Microsoft MimeOLE V6.00.2900.3138
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=3DContent-Type content=3D"text/html; =
charset=3Dwindows-1252">
<META content=3D"MSHTML 6.00.6000.20591" name=3DGENERATOR></HEAD>
<BODY><PRE>% Program to simulate Rayleigh fading using a p-th order =
autoregressive model AR(p) according to=20
% Baddour's work: "Autoregressive modeling for fading channel =
simulation", IEEE Transaction on Wireless Communications, July 2005.
function [chann]=3DRayleigh_fading(P,M,fm,fs,epselonn)
% P: AR model order
% M: number of samples
% fd: maximum doppler frequency in Hz
% fs: Symbol frequency in ksps
% epselonn: added bias, depends on the Doppler rate, see the paper =
"Autoregressive modeling for fading channel simulation".=20
=20
% usage: Rayleigh_fading(100,10000,150,3,0.00000001)
=20
=20
=
%------------------------------------------------------------------------=
------------------------------------------------- =20
=20
for p=3D1:P+1
vector_corr(p)=3Dbesselj(0,2*pi*fm*(p-1)/(fs*1000)); % Bessel =
autocorrelation function according to Jakes' model
end
=20
=
auto_correaltion_matrix=3Dtoeplitz(vector_corr(1:P))+eye(P)*epselonn; % =
adding a small bias, epselonn, to the autocorrelation matrix to overcome =
the ill conditioning of Yule-Walker equations
=
AR_parameters=3D-inv(auto_correaltion_matrix)*vector_corr(2:P+1)'; % =
Solving the Yule-Walker equations to obtain the model parameters
=
segma_u=3Dauto_correaltion_matrix(1,1)+vector_corr(2:P+1)*AR_parameters;
=20
KKK=3D2000;
=20
h=3Dfilter(1,[1 =
AR_parameters.'],wgn(M+KKK,1,10*log10(segma_u),'complex')); % Use the =
function Filter to generate the channel coefficients
chann=3Dh(KKK+1:end,:); % Ignore the first KKK samples
=20
=
%------------------------------------------------------------------------=
-------------------------------------------------=20
=20
</PRE></BODY></HTML>
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -