1、1无线衰落信道及仿真无线衰落信道在无线移动信道中，信号从发射天线经过复杂的传播环境到达接收天线，接收信号为各反射、衍射和散射分量以及信道噪声的复合，因而会产生严重的失真。另外，移动信道中由于移动台运动或信道环境的改变会使信道特性随时间随机变化，接收到的信号由于 Doppler 效应会产生更为严重的失真。信号在无线移动信道中传播，除了自由空间固有的传输损耗之外，还会由于受到建筑物、地形等的阻挡而引起信号功率的衰减。除了这些衰减作用之外，影响接收信号的主要因素包括：1 多径传播无线移动信道中，由于反射、散射等的影响，实际到达接收机的信号为发射信号经过多个传播路径之后各分量的叠加。不同路径分量的幅度
2、、相位、到达时间和入射角各不相同，使接收到的复合信号在幅度和相位上都产生了严重的失真。多径传播会引起信号在时间上的展宽，从而带来符号间的干扰（ISI ） 。2 移动台的运动速度在无线移动系统中，需要使用很高的载波频率进行信号传送。如果移动台相对于基站运动，由于各入射信号的入射角不相同，各路径分量受到不同的Doppler 频率调制，使接收到的复合信号产生非线性失真。若所使用的载波频率一定，移动台的移动速度越高，Doppler 频移对接收信号的影响就越严重。3 传播环境中物体的运动如果无线信道环境中存在运动的物体，会使到达接收天线的某些多径分量随时间变化。如果移动物体处于发射或接收天线附近且具有较
3、高的速度，这时，移动环境中运动物体引起的 Doppler 频移对信号的影响就必须加以考虑。24 信号的物理带宽宽带信号和窄带信号在多径信道中的表现出不同的衰落特性。如果传送信号的物理带宽比“信道带宽” （相干带宽）更宽，接收信号将产生失真。但如果信号带宽比 Doppler 带宽大很多，信号对 Doppler 频移引起的失真将不敏感。如果传送信号的物理带宽比信道带宽窄，则接收信号波形在时间上不会引起明显的失真。但如果信号带宽窄到可以与 Doppler 带宽相比拟时，信号对 Doppler 频移引起的失真将较为敏感。平衰落信道（多径传播，无附加时间扩展）传播特性与仿真模型(1) 连续导频正弦波通过
4、存在 Doppler 频移的快衰落信道假设基站使用一固定的垂直极化天线，移动台垂直于地面的接收天线接收了来自 个入射方位角的电磁波，每一入射电磁波的方向、相位是随机的，并N假设每一入射波能量相同。图 1 给出了移动台入射角为 的一束平面波示意图。实际上，如果移动台与基站之间不存在可视路径（LOS） ，则接收天线的每一入射分量所经历的散射衰耗基本相同，等能量假设是合理的。设发射机发送一角频率为 的连续单频正弦信号0（1）tAts0cos)(设移动台相对于基站的运动速度为 ，第 条入射电磁波与移动台运动方vn向夹角为 ，则其 Doppler 频移为n3xyzvx-y平 面 入 射 波图 1 移动台
5、随机入射角平面波示意图（2） ndndnnffcvcoscs2os20式中 为真空中光速， 为最大 Doppler 频率偏移。该单频信号经过多径信道到达移动台的接收收信号为（3）ttTttCAtr scNnn 0010 in)(o)(cos()( 式中 和 分别为tTc)ts（4）)o()(1nNnctt（5）)si()(1nns tCAtT和 为二不相关的零均值高斯过程，并具有相同的功率谱密度。tc)ts现考察接收信号 的功率谱密度。考虑 的相关函数(r)(tr)(sin)(cossin)(cos),( 0000 ttTtTttTEtRr i)( 0Esc（6）0s)cT式中假定 和 都为平
6、稳随机过程并具有相同的自相关函数。设随机信号(t)ts和 的功率谱密度分别为 和 。对式（6）作傅立叶变换，)(tcs )(cTP)(sT可得到信号 的功率谱密度)(tr（7）2/)()00ccTTrP4设接收天线的微小入射角 内入射功率为 ，天线方向增益为ddp)(,当入射波数 时，接收功率可以表示为)(GN（8）20)(pPr而入射角频率（9）cos)(0d从上式可见，从 和 方向的入射波将产生相同的频移。另外，考虑信号为实信号，其双边功率谱的正负边带关于纵轴对称，即 ，于)()rrP是有（10）|)(|)(|)(|)(210 dpGdpGdPdrr 由式（9）得到 ，采用 天线（ ）并假
7、|sin| 4/5.1定各方向入射功率相等（ ） 。令 ，并注意到功率谱正负频域2/1)(p10的对称性，可得（11） 其 他， ，0|)(1/5. |.)( 020 dddrP比较式（7）和（11）可以得到 和 的功率谱)(tTcts（12）其 他， ，0|1/.3)()( 2dddTTscP(2)复模型（基带等效模型）的构造I(t)和 Q(t)分别为发送信号的同行和正交分量。对于 twtQtIccsin)(os)(iRe tjj c 根据： )(2Re)(tjwcftf则 为包络5)(21*RetjwcejQtI其中包络为 ，t对于 )(tI )(cos(1 tItwCAtrNn ncI)
8、(sin)si(s)( 11 tIwtACtt cNnncnn )(icos)( tIwtTttTcs对于 )(tQ)(sin(1 tQtwCAtrNn nc)(cos)si(si)o( 11 tQwtACttNnncnn i(c)( twtTttTQccs)(sin)(cos)()( tIwtTttTtrt cQI Qcs )(sin)(os)()( 222 tIwtTttwtTttTt cccscccsc 其中： )()(1Nnnc tACt)sin()(1Nnns twACt其中： )()(1Nnnc tACt)sin()(1Nnns tACt6)(sin)(cos)()( 222 tQ
9、wtTttwtTttTt ccccssc )(os)()(22 tItttcc )(inQwTs)(cos)()(22 twtITtsc inQcs)(sin)(s)()(Re2 twjttjtIt ccsc 1)(22 tjtjsc cejtTt 包络为 )(22)()( tjsc ejQtItT)()(sinco12 tjQtIjtttsc )()(Ijs可用下图描述： )()(21tjQtI)()(21tjTtTsc单延迟、多径衰落和 均为实信号，二者功率谱相同， 为多普勒滤波器)(tTcts )(thD其中： )()(cos2tttcc )()(i2tTtccs7SUI 信道Stand
10、frod UniversityInterim (SUI) channel modelsA 类地形：中高密林山地（高损耗）B 类地形：低密林山地（中损耗）C 类地形：低密林（近乎）平地 （低损耗）地形类型 SUI 模型C SUI-1，SUI-2B SUI-3，SUI-4A SUI-5，SUI-6SUI-1,2 多普勒频移低 时间扩展小 C 地形（低密林近平地）SUI-3 多普勒频移低 时间扩展小 B 地形（低密林山地）SUI-4 多普勒频移高 时间扩展中等 B 地形（低密林山地）SUI-5 多普勒频移低 时间扩展大 A 地形（中高密林山地）SUI-6 多普勒频移高 时间扩展大 A 地形（中高密林
11、山地）8信道仿真独 立 复 高 斯白 噪 声 源独 立 复 高 斯白 噪 声 源独 立 复 高 斯白 噪 声 源 DoplerfitDoplerfitDoplerfit 120CN（ ， ）21CN（ 0， ）2CN（ 0， ） ()S 20m1m2m201normFFFF经大尺度衰落 )(r)(1nh)(2nh)(3nhSUI-3 信道9附：SUI-3 信道仿真的 MATLAB 程序close all;clear all;N=10000; %number of independent realizationOR=4; %observation rate in HzM=256; %number
12、 of taps of the Doppler filterDop_res=0.1; %Doppler resolution of SUI parameter in Hz(used in resampling-process)res_accu=20; %accuracy of resamplingP=0 -5 -10; %power in each tap in dBK=1 0 0; %Reciean K-factor in linear scaletau=0.0 0.5 1.0; %tap delay in usDop=0.4 0.4 0.4; %Doppler maximal freque
13、ncy parameter in Hzant_corr=0.4; %antenna correlation (envelope correlation coefficient)Fnorm=-1.5113; %gain normalizationP=10.(P/10); %calculate linear powers2=P./(K+1); %calculate variancem2=P.*(K./(K+1); %calculate constant powerm=sqrt(m2);%Addition info:RMS delay spreadrmsdel= sqrt( sum(P.*(tau.
14、2)/sum(P) -( sum(P.*tau)/sum(P)2 );fprintf(rms delay spread %6.3f usn,rmsdel);10%Now we create the Ricean channel coefficients with the specified powersL=length(P);paths_r=sqrt(1/2) *(randn(L,N)+j*randn(L,N) .*(sqrt(s2)*ones(1,N);paths_c=m*ones(1,N);for p=1:LD=Dop(p)/max(Dop)/2; %normalize to highes
15、t Dopplerf0=0:M*D/(M*D); %frequency vectorPSD=0.785*f0.4-1.72*f0.2+1.0;filt= PSD(1:end-1) zeros(1,M-2*M*D) PSD(end:-1:2); %S(f)filt=sqrt(filt); %from S(f) to |H(f)|filt=ifftshift(ifft(filt); %get impulse responsefilt=real(filt); %want a real-valued filterfilt=filt/sqrt(sum(filt.2); %normalize filter
16、path=fftfilt(filt, paths_r(p,:) zeros(1,M) );paths_r(p,:)=path(1+M/2:end-M/2);end;paths=paths_r+paths_c;paths=paths*10(Fnorm/20); %Multiply all coefficients with F%Additional Info:average total tap powerPest=mean(abs(paths).2,2);fprintf(tap mean power level:%0.2f dBn,10*log10(Pest);%Additional Info:
17、spectral power distribution%figure,psd(paths(1,:),512,max(Dop) );figure,psd(paths(1,:),512,2*max(Dop) ); %should be 2 times max(Dop)SR=max(Dop)*2; %implicit sample ratem=lcm(SR/Dop_res,OR/Dop_res);11P=m/SR*Dop_res; %find nominatorQ=m/OR*Dop_res; %find denominatorpath_OR=zeros(L,ceil(N*P/Q); %create new arrayfor p=1:Lpaths_OR(p,:)=resample(paths(p,:),P,Q,res_accu);end;
