通信系统仿真课设——多径时变信道

1、*实践教学*兰州理工大学计算机与通信学院2013年春季学期通信系统仿真训练课程设计题 目：多径时变信道模型的仿真与性能分析专业班级： 通信工程（1）班 姓 名： 学 号： 指导教师： 成 绩： 摘 要本课程设计从信道的随机过程特征出发，分析了多径衰落信道的时间选择性和频率选择性，讨论信道的衰落特征，并重点针对服从瑞利分布的多径衰落信道模型在matlab环境下进行模拟仿真，然后观察单频信号、数字信号经过多径时变信道后的时域波形和频谱，通过仿真分析了多径效应和多普勒效应对信道以及信号的影响,有助于对信道的进一步理解和掌握以及工程实践提供一定的借鉴。关键词：多径效应；时变信道；瑞利衰落；matlab

2、a目 录前 言1一、多径时变信道概述21.1 多径衰落信道的分类21.1. 1 平坦衰落和频率选择性衰落21.1. 2 快衰落和慢衰落31.1.3 时不变多径衰落信道31.1.4时变多径衰落信道31.2多径衰落信道基本特性41.2.1 时间选择性与频率选择性41.2.2多径效应与频率选择性41.3瑞利衰落信道51.3.1瑞利衰落信道简介51.3.2瑞利衰落信道基本模型61.3.3产生服从瑞利分布的路径衰落r(t)61.3.4产生多径延时71.3.5仿真框架8二、时变多径信道仿真程序设计92.1时变多径信道仿真模块92.2瑞利信道仿真模块92.3单频信号经过多径信道仿真模块102.4数字信号经过

3、多径信道仿真模块10三、设计结果与测试113.1时变多径信道仿真113.1.1仿真移动台在不同位置的多径信号113.1.2仿真不同频率的信号的多径信号改变123.1.3仿真移动台不同速度的信号的多径信号改变133.2瑞利信道仿真133.3多径信道对单音频信号的影响仿真143.4 多径信道对数字信号影响15总结17参考文献18附录19致谢252前 言 由于信道中电磁波受到反射、绕射、散射、多径传播等因素的影响，接收端所接收到的信号是各个方向到达电磁波的叠加，使信号在小范围内引起剧烈的波动，称之为多径衰落，亦称为小尺度衰落。小尺度衰落直接体现了无线信道的复杂性和随机性，是决定无线通信系统性能的基本

4、问题。小尺度衰落信道按相干带宽分为平坦衰落信道和频率选择性衰落信道，按相干时间与信号周期的关系分为慢衰落信道和快衰落信道。由于频率选择性信道可以通过平坦衰落信道实现，因此，平坦衰落信道的计算机仿真是研究多径信道、MIMO信道的重要方法，更是进一步研究信道容量、信道编码及均衡器等模型的前提。一般而言，在市区内，平坦衰落信道响应的包络服从瑞利分布，而相位服从均匀分布。因此，无线信道建模的核心与关键即是如何简单有效地产生瑞利衰落包络，这也是研究与分析各种移动通信系统的首要任务。当信号的多径发生在发送信号经由传播路径以不同的延迟到达接收机的时候,一般会引起数字通讯系统中的符号间干扰。而且,由不同传播路

5、径到达的各信号分量会相互削弱,导致信号能量衰减,造成信噪比降低。多径传播导致信号在不同域中的扩展,包括角度(或空间)扩展、时延(或时间)扩展和多普勒(或频率)扩展,这些扩展对信号有显著影响。研究多径衰落信道的特性有利于找到适合不同通讯要求的数字通讯方式,特别是针对军事领域电子对抗的通讯设计和有限信道资源的合理使用。25一、多径时变信道概述多径时变信道是指信道参数随时间变化，它对信号传输的影响是输入信号的频率弥散。当输入信号为单频信号，经过时变信道后的输出不再是单频信号，而是一个窄带信号，带宽大小视时变因素的快慢而定。多径时变信道是信号传输的路径不止一条，接收端同时受到来自多条传输路径的信号，这

6、些信号可能同相相加或反向相消。信号经多径时变信道后，会产生码间干扰和衰落。其中衰落快慢取决于信道随时间变化的快慢，码间干扰的严重程度取决于码元间隔和多径间的时延差的相对关系。1.1 多径衰落信道的分类移动通信中的时间色散和频率色散可能产生4种衰落效应1,这是由信号、信道以及相对运动的特性引起的。根据信号带宽和信道带宽的比较,可将信道分为平坦衰落信道和频率选择性衰落信道;而根据发送信号与信道变化快慢程度的比较,可以将信道分为快衰落和慢衰落信道。1.1. 1 平坦衰落和频率选择性衰落如果信道带宽大于发送信号带宽,且在带宽范围内有恒定增益,且线性相关,则接收机信号就会经历平坦衰落过程。在平坦衰落情况

7、下,信道的多径结构使发送信号的频谱特性在接收机内仍能保持不变,所以平坦衰落也称为频率非选择性衰落2。平坦衰落信道的条件可以概括为:Bs << BC ,Ts >> 其中, Ts 为信号周期(信号带宽Bs 的倒数) ;t 为时延扩展; BC 为相关带宽。如果信道具有恒定增益和线性相位的带宽范围小于发送信号带宽,则该信道特性会导致接收信号产生频率选择性衰落。此时,信道冲激响应具有多径时延扩展,其值大于发送信号带宽的倒数。在这种情况下,接收信号中包含经历了衰减和时延的发送信号的多径波,因而产生接收信号失真。频率选择性衰落是由信道中发送信号的时间色散引起的,这种色散会引起符号间干

8、扰。对于频率选择性衰落而言,发送信号的带宽大于信道的相关带宽,有频域可以看出,不同频率信号获得不同增益,产生频率选择性衰落,其条件是:Bs > BC ,Ts < 1.1. 2 快衰落和慢衰落当信道的相关时间比发送信号周期短,且信号的带宽Bs 小于多普勒扩展BD 时,信道冲激响应在符号周期内变化很快,从而导致信号失真,产生衰落,此衰落为快衰落1。当信号的相关时间远大于发送信号的周期,且信号的带宽Bs 远远大于多普勒扩展BD 时,信道冲激响应变化比要传送的信号码元的周期低的多,可以认为该信道是慢衰落信道1。1.1.3 时不变多径衰落信道在没有多普勒效应的情况下(收发两端不存在相对运动)

9、 ,信道为时不变多径信道2,信道的转移函数只是频率的函数, 对于时不变多径信道,信道转移函数随着频率的变化而改变。因为多径波幅度和到达时间取决于收发两端的距离,所以信号接收强度同样取决于收发两端的距离。转移函数随着距离增加衰落加剧,它同样是距离的函数。可以看出对于时不变多径衰落信道而言,如果信道的时延扩展越大,则输入信号衰落越严重。1.1.4时变多径衰落信道在此种情况下,由于存在多条到达路径,所以每条路径的到达角度和到达时间都不相同,并且各自相互独立。当收发两端存在相对运动时,则信道变为了时变信道2。下面以正弦信号举例说明时不变信道与时变信道对信号的影响。信号经过时变信道后产生很大的衰落,其幅

10、度随着时间的变化而变化,而时不变信道对信号基本没有改变。时变信道对信号的衰落随着多普勒扩展的增大而加剧。1.2多径衰落信道基本特性1.2.1 时间选择性与频率选择性由于移动台随机运动的多普勒2 (Doppler) 效应引起的信道频率色散,造成了传播信道的时间选择性衰落;而微区环境使信号传播多径时延引起的信道时间色散,造成了传播信道的频率选择性衰落。多普勒效应2与时间选择性信道的时变特性是通过频率色散参数来描述的。这种时变特性或是由移动台和基站之间的相对运动引起的,或是由信道中物体运动引起的。当信道是时变的,则这种信道具有时间选择性衰落。时间选择性衰落会造成信号的失真,这是由于信号的接收端和发送

11、端的信道特性改变所致,信号在时刻改变的信道中传输自然会导致信号的失真。由于移动台的移动,产生了多普勒频移,也就是频率色散,使得信道是时变的。如果信号持续的时间比较短,以至于信道在信持续的时间里没有显著的变化,那么信号受信道时间选择性衰落较小;但随着信号持续时间的增加,信道的特性产生较为显著的变化,那么信号就会产生时间选择性衰落,产生失真。信号的失真会随着持续时间的增加而加剧。1.2.2多径效应与频率选择性时间色散和频率选择性是同时出现的,只是表现的形式不同,都是依赖与所处的环境中的因素变化,具体来说就是因为不同时延的多径波叠加的结果。从时域上说,由于经历了多径效应的影响,导致在接收端不同分量的

12、叠加,展宽了时域。从频域上说,由于信道也就是起到了滤波的作用,当发送信号的带宽极窄时,那么接收端的信号的所有频率成分都经历了相同的衰减,也就是频率非选择性衰落(平坦衰落) 。随着发送信号的带宽增加,那么接收信号在传输过程中,频率接近的分量经历了相似的衰减,而频率相差很远的分量就会经历不同的衰减。也就是频率选择性衰落。当发送信号的带宽非常大的时候,那么在接收端接收到的信号会与发送信号大相径庭,接收机会受到时间色散的影响,产生严重的衰落。在数字通信系统中,会导致严重的符号间干扰。时延扩展:用来描述时延扩展的参数有:平均附加时延, rms 时延扩展和附加时延扩展。相关带宽:相关带宽BC 表示信号包络

13、相关度为某一特定值时的信号带宽。也就是说,当频率间隔小于BC 时,信号包络(幅度) 具有很强的相关性,当频率间隔大于BC 时,信号包络(幅度) 相关性很小。1.3瑞利衰落信道1.3.1瑞利衰落信道简介瑞利衰落信道3（Rayleigh fading channel）是一种无线电信号传播环境的统计模型。这种模型假设信号通过无线信道之后，其信号幅度是随机的，表现为“衰落”特性，并且多径衰落的信号包络服从瑞利分布。由此，这种多径衰落也称为瑞利衰落。 （1）环境条件：通常在离基站较远、反射物较多的地区，发射机和接收机之间没有直射波路径，存在大量反射波；到达接收天线的方向角随机且在（02）均匀分布；各反射

14、波的幅度和相位都统计独立。（2）幅度、相位的分布特性：包络 r 服从瑞利分布，在02内服从均匀分布。瑞利分布的概率分布密度如图1-1所示：图1-1 瑞利分布的概率分布密度1.3.2瑞利衰落信道基本模型根据ITU-RM.1125标准，离散多径衰落信道模型3为 (1-1)其中复路径衰落，服从瑞利分布； 是多径时延。多径衰落信道模型框图如图1-2所示：图1-2 多径衰落信道模型框图假设经反射（或散射）到达接收天线的信号为N个幅值和相位均随机的且统计独立的信号之和。信号振幅为r,相位为,则信号经过瑞利衰落信道其包络概率密度函数为P(r)= (r0) （1-2）相位概率密度函数为：P()=1/2 （）

15、（1-3）1.3.3产生服从瑞利分布的路径衰落r(t)利用窄带高斯过程的特性，其振幅服从瑞利分布4，即 （1-4）上式中，、分别为窄带高斯过程的同相和正交支路的基带信号。首先产生独立的复高斯噪声的样本，并经过FFT后形成频域的样本，然后与S（f）开方后的值相乘，以获得满足多普勒频谱特性要求的信号，经IFFT后变换成时域波形，再经过平方，将两路的信号相加并进行开方运算后，形成瑞利衰落的信号r(t)。如下图1-3所示:图1-3 瑞利衰落的产生示意图其中， （1-5）1.3.4产生多径延时多径/延时参数5如表1-1所示：表1-1 多径延时参数TapRelative delay (ns)Average

16、 power (dB)1002310-1.03710-9.041 090-10.051 730-15.062 510-20.01.3.5仿真框架根据多径衰落信道模型5（见图1-2），利用瑞利分布的路径衰落（见图1-3）和多径延时参数（见表1-1），我们可以得到多径信道的仿真框图，如图1-4所示：图1-4 多径信道的仿真框图二、时变多径信道仿真程序设计本设计分为三大模块：时变多径信道仿真模块，瑞利信道仿真，单频信号经过多径信道模块，数字信号经过多径信道模块移动台与基站距离不同频率的信号不同速度绘图比较2.1时变多径信道仿真模块图2-1 时变多径信道仿真模块设计流程图2.2瑞利信道仿真模块设计瑞利

17、信道绘出瑞利信道的概率分布对比理论和实际的瑞利信道分布 图2-2 瑞利信道仿真模块设计流程图单频信号20径衰落信道对接收信号进行时域频域分析2.3单频信号经过多径信道仿真模块图2-3 单频信号经过多径信道仿真模块设计流程图数字信号3径衰落信道对接收信号进行时域频域分析对接收信号进行误码分析2.4数字信号经过多径信道仿真模块图2-4 数字信号经过多径信道仿真模块设计流程图三、设计结果与测试3.1时变多径信道仿真3.1.1仿真移动台在不同位置的多径信号r0=3000、9000时的直射信号、反射信号与合成信号如图3-1、3-2图3-1 r0=3000时的直射信号、反射信号与合成信号图3-2 r0=9

18、000时的直射信号、反射信号与合成信号从图3-1、3-2中可以看出，即使移动台是静止的，由于反射径的存在，使得接收到的合成信号最大值要小于直射径的信号。r0=9000由于靠近反射墙的位置，直射信号要比r0=3000处要弱一些，反射信号要比r0=3000位置处的信号要强一些，但移动台接收到的合成信号更弱了。不仅要小于直射径的信号更小于反射径的信号。结论，即使移动台静止，由于反射径的存在，使接收信号要比没有反射径时的信号弱，衰落由此产生7。3.1.2仿真不同频率的信号的多径信号改变分别让f=1和f=5，画出的图形分别如图3-3、图3-4所示：图3-3 f=1，c=10时的直射信号、反射信号与合成信

19、号图3-4 f=5、c=10时的直射信号、反射信号与合成信号从结果我们明显可以看出，这次是f=1的信号得到增强，而f=5的信号得到削弱。这就是频率选择性衰落。在同一位置，由于反射径信号的存在，发射不同频率的信号时，在接收机处接收到信号有的频率是被增强了，有的频率是被削弱了。频率选择性由此产生7。3.1.3仿真移动台不同速度的信号的多径信号改变移动台向反射墙运动，速度为v图3-5 移动台运动时的直射信号、反射信号与合成信号当移动台运动起来以后，我们发现即使同一频率，在不同的时间点，合成信号的强度也是不一样的。在图3-5中，我们可以看到在t=2，4.5，7，9.5s时，接收信号的强度相对处于波谷位

20、置，特别是在t=9.5s时，接收的合成信号几乎为0，而我们对照一下t=9.5s时的直射信号和反射信号，它们都比合成信号大很多。而在t=3，5.5，8，10.5位置，接收信号的强度相对处于波峰位置。这种由于移动台运动而导致的信号增强或削弱的情况就是时间选择性衰落8。3.2瑞利信道仿真瑞丽分布的概率密度函数如图所示： 图3-6 瑞利分布的概率密度函数图分析：包络 r 服从瑞利分布，在02内服从均匀分布。瑞利分布的概率分布密度8如图3-6，经过分析对照可以发现，输出的的统计特性完全符合要求。3.3多径信道对单音频信号的影响仿真多径效应总的来说有三点，即对单一正弦波产生频域弥散，对宽带信号频率选择性衰

21、落以及对数字信号产生时域弥散。程序见附录2（1）输入的单频正弦信号时域及频域图像如下图3-7图3-7 单频正弦信号时域及频域图像分析：由图可知，原信号为单一正弦信号，时域标准正弦，频域单一冲激。（2）信号经多径传播后在接收端所得信号时域及频域图像如下图3-8图3-8 单频正弦信号经多径传播后在接收端所得信号时域及频域图像分析：由图可知，过多径后，不同路径时延不同。不同时延的信号叠加，导致时域图形不再是单一正弦。（由于采样频率的关系，看起来还是光滑曲线）频域出现了毛刺，即频域扩散9。3.4 多径信道对数字信号影响 程序见附录3 信道幅频和相频响应如图3-9所示： 图3-9信道幅频和相频响应 分析

22、：通过该信道的幅频响应和相频响应可以看出，信道的幅度和频率响应都不理想，输入信号过信道会造成幅度频率失真。 信号进入多径信道前后的时域频域图如图3-10所示： 图3-10 信号进入多径信道前后的时域频域图分析：如图3-10所示，因为信道幅频特性不理想，会造成输入信号失真，针对模拟信号体现在波形失真，针对数字信号体现为码间干扰。由于各径时延差不同，通过各路径的信号衰减不同，信号经过多径后形成码间干扰90。总结本课程设计从信道的随机过程特征出发，研究信道的平稳过程,讨论信道衰落特征，并重点针对服从瑞利分布的多径衰落信道模型在matlab环境下进行模拟仿真，观察了单频信号、数字信号经过多径时变信道后

23、的时域波形和频谱。首先我们查阅了许多有关多径时变信道、瑞利信道的资料，了解了多径时变信道的特点，掌握独立瑞利衰落的多径信道的模型，然后我们用matlab软件仿真出瑞利衰落衰落的多径信道，看到了这种信道的瑞利分布图，接着我们在matlab环境下，将单频信号通过我们所设计的多径时变信道，然后在接收端观察接收信号的时域及频域波形，看到了时域波形出现弥散的结果，最后我们将数字信号通过多径时变信道，并观察起输出信号的波形和频谱，发现它产生了码间干扰并造成误码，这些实验结果符合我们在理论上分析的结果。在这次课程设计过程中我锻炼了自己的思考能力和动手能力。通过对课题的分析和仿真设计，加强了我思考问题的完整性

24、和实际生活联系的可行性。在方案设计选择和衰落的选择上，培养了我们对各种衰落的认识。也进一步锻炼了我们个人查阅技术资料的能力，动手能力，发现问题，解决问题的能力。总之，这次课程设计中我学到了很多知识，拓展了我的思维，而且让我明白了怎样理论联系实际，在实践过程中遇到的问题怎样用理论来解决，为以后的工作积累了经验，增强了信心。参考文献 1 BELLO. P. A. Characterization of randomly time variant lin2ear channelsJ . IEEE Trans. Communications Systems , 1963 ,11 (4) :360 -

25、393.2 CLARKE ,R. H.A statistical theory of mobile radio receptionJ . Bell Systems Technical Joumal ,1968 ,47 (6) : 957 - 1000.3樊昌信，曹丽娜编著，通信原理，北京：国防工业出版社，2010：73-844樊昌信，通信原理，北京：国防工业出版社，2002:70-80:67-905曹志刚等著，现代通信原理，北京：清华大学出版社，2001.5。：65-786吴伟陵等著，移动通信原理，北京：电子工业出版社，2005：30-577.李建新，现代通信系统分析与仿真-MATLAB 通信

26、工具箱，西安：西安电子科技大学出版社，2000：13-458潘子宇，Matlab通信仿真设计指导书，南京工程学院，2011：24-369刘敏，MATLAB通信仿真与应用，北京：国防工业出版社：56-60附 录源程序代码如下源代码1：时变多径信道仿真（1）移动台在不同位置的多径信号clear allf=9e8; %发射信号频率v=0; %移动台速度，静止情况为0c=3e8; %电磁波速度，光速r0=1000; %移动台距离基站初始距离d=15000; %基站距离反射墙的距离t1=0:0.0000000000005:0.00000001; %时间E1=cos(2*pi*f*(1-v/c).*t1-

27、r0/c)./(r0+v.*t1); %直射径信号E2=cos(2*pi*f*(1+v/c)*t1+(r0-2*d)/c)./(2*d-r0-v*t1); %反射径信号subplot(2,3,1)plot(t1,E1,t1,E2,'-g',t1,E1-E2,'-r') %画出直射径、反射径和总的接收信号legend('直射径信号','反射径信号','移动台接收的合成信号')%axis(0 12 -0.5 0.5)subplot(2,3,4)plot(t1,E1-E2)f=9e8; %发射信号频率v=0; %移动台速

28、度，静止情况为0c=3e8; %电磁波速度，光速r0=3000; %移动台距离基站初始距离d=15000; %基站距离反射墙的距离t1=0:0.0000000000005:0.00000001; %时间E1=cos(2*pi*f*(1-v/c).*t1-r0/c)./(r0+v.*t1); %直射径信号E2=cos(2*pi*f*(1+v/c)*t1+(r0-2*d)/c)./(2*d-r0-v*t1); %反射径信号subplot(2,3,2)plot(t1,E1,t1,E2,'-g',t1,E1-E2,'-r') %画出直射径、反射径和总的接收信号%axis

29、(0 12 -0.5 0.5)subplot(2,3,5)plot(t1,E1-E2)f=9e8; %发射信号频率v=0; %移动台速度，静止情况为0c=3e8; %电磁波速度，光速r0=9000; %移动台距离基站初始距离d=15000; %基站距离反射墙的距离t1=0:0.0000000000005:0.00000001; %时间E1=cos(2*pi*f*(1-v/c).*t1-r0/c)./(r0+v.*t1); %直射径信号E2=cos(2*pi*f*(1+v/c)*t1+(r0-2*d)/c)./(2*d-r0-v*t1); %反射径信号subplot(2,3,3)plot(t1,

30、E1,t1,E2,'-g',t1,E1-E2,'-r') %画出直射径、反射径和总的接收信号%axis(0 12 -0.5 0.5)subplot(2,3,6)plot(t1,E1-E2)（2）仿真不同频率的信号的多径信号改变clear allf=3e8; %发射信号频率v=0; %移动台速度，静止情况为0c=3e8; %电磁波速度，光速r0=3000; %移动台距离基站初始距离d=15000; %基站距离反射墙的距离t1=0:0.0000000000005:0.00000001; %时间E1=cos(2*pi*f*(1-v/c).*t1-r0/c)./(r0+

31、v.*t1); %直射径信号E2=cos(2*pi*f*(1+v/c)*t1+(r0-2*d)/c)./(2*d-r0-v*t1); %反射径信号subplot(2,3,1)plot(t1,E1,t1,E2,'-g',t1,E1-E2,'-r') %画出直射径、反射径和总的接收信号legend('直射径信号','反射径信号','移动台接收的合成信号')%axis(0 12 -0.5 0.5)subplot(2,3,4)plot(t1,E1-E2)f=9e8; %发射信号频率v=0; %移动台速度，静止情况为0c=3

32、e8; %电磁波速度，光速r0=3000; %移动台距离基站初始距离d=15000; %基站距离反射墙的距离t1=0:0.0000000000005:0.00000001; %时间E1=cos(2*pi*f*(1-v/c).*t1-r0/c)./(r0+v.*t1); %直射径信号E2=cos(2*pi*f*(1+v/c)*t1+(r0-2*d)/c)./(2*d-r0-v*t1); %反射径信号subplot(2,3,2)plot(t1,E1,t1,E2,'-g',t1,E1-E2,'-r') %画出直射径、反射径和总的接收信号%axis(0 12 -0.5

33、0.5)subplot(2,3,5)plot(t1,E1-E2)f=27e8; %发射信号频率v=0; %移动台速度，静止情况为0c=3e8; %电磁波速度，光速r0=3000; %移动台距离基站初始距离d=15000; %基站距离反射墙的距离t1=0:0.0000000000005:0.00000001; %时间E1=cos(2*pi*f*(1-v/c).*t1-r0/c)./(r0+v.*t1); %直射径信号E2=cos(2*pi*f*(1+v/c)*t1+(r0-2*d)/c)./(2*d-r0-v*t1); %反射径信号subplot(2,3,3)plot(t1,E1,t1,E2,&

34、#39;-g',t1,E1-E2,'-r') %画出直射径、反射径和总的接收信号%axis(0 12 -0.5 0.5)subplot(2,3,6)plot(t1,E1-E2)源代码2：瑞利信道仿真function h=rayleigh(fd,t) %产生瑞利衰落信道fc=900*106; %选取载波频率v1=30*1000/3600; %移动速度v1=30km/hc=3*108; %定义光速fd=v1*fc/c; %多普勒频移ts=1/10000; %信道抽样时间间隔t=0:ts:1; %生成时间序列h1=rayleigh(fd,t); %产生信道数据v2=120*1

35、000/3600; %移动速度v2=120km/hfd=v2*fc/c; %多普勒频移h2=rayleigh(fd,t); %产生信道数据subplot(2,1,1),plot(20*log10(abs(h1(1:10000) title('v=30km/h时的信道曲线')xlabel('时间');ylabel('功率')subplot(2,1,2),plot(20*log10(abs(h2(1:10000)title('v=120km/h时的信道曲线')xlabel('时间');ylabel('功率

36、9;)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; %区别个条

37、路径的均匀分布随机相位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); %乘归一化功率系数得到传输函数源代码3：多径信道对单音频信号的影响仿真clc;close all;cl

38、ear all;%-参数定义-A=1;f=10;decay=0.8;dt=0.01;t=0:dt:10;L=20;fdelay=2*rand(1,L);wdelay0=rand(1,L)*2*pi;x=cos(2*pi*f*t);for i=1:L wdelay(i,:)=cos(2*pi*fdelay(i)*t); s(i,:)=decay*cos(2*pi*f*t+wdelay(i,:)+wdelay0(i);endy=sum(s)/sqrt(L);%-原始单频信号- figure(1);subplot(211);plot(t,x);xlabel('t/s');ylabel

39、('x(t)');title('单频信号时域图');axis(0 2 -1.5 1.5); subplot(212);f1=abs(fft(x,1024);plot(0:100/1024:100-100/1024,f1);xlabel('f/Hz');ylabel('H(f)');title('单频信号频域图');axis(0 30 0 500);%-经过多径信道的单频信号-figure(2);subplot(211);plot(t,y);xlabel('t/s');ylabel('y(t)

40、');title('经过20径后的时域图'); subplot(212);f2=abs(fft(y,1024);plot(0:100/1024:100-100/1024,f2);xlabel('f/Hz');ylabel('H(f)');title('经过20径后的频域图');axis(0 30 0 500);源代码4：多径信道对数字信号影响(1).作出信道的幅频响应以及相频响应clc;clear all;f=-2:0.001:2; %定义f的范围h1=0.5*exp(-j*2*pi*f*0);h2=0.707*exp(-j*2*pi*f*1);h3=0.5*exp(-j*2*pi*f*2);h=h1+h2+h3; %信道的响应subplot(2,1,1);plot(f,h);grid on;axis(-2 2 -0.5 2); %定义坐标范围subplot(2,1,2);plot(f,angle(h)/pi);grid on ;(2).作出信号输入前和经过信道后的波形以及其幅频响应TS=1时clc;clear