多径时变信道模型仿真及性能分析分解

1、“*实践教学“*兰州理工大学计算机与通信学院2013年春季学期通信系统仿训练真课程设计题目：基于MATLAB的FIR滤波器语音信号去噪专业班级：姓 名：学号：指导教师：成绩：摘要本次课程设计做的是多径时变信道模型的仿真与性能分析，首先需要建立 信道模型，通过对输入信号和移动台的有些参数进行调整，使用MATLAB进行 仿真，得到时域和频域图，对比分析掌握多径信道的特点；其次，对瑞利衰落的 多径信道仿真，分析信道模型的特点；最后，观察单频和数字信号经过多径信道 后接收信号的情况。经过多次修改调试，最终完成了设计任务。关键词：多径时变信道；瑞利衰落；仿真；信道模型 TOC o 1-5 h z HYP

2、ERLINK l bookmark7 o Current Document 一多径信道的基本原理1 HYPERLINK l bookmark10 o Current Document 1.1移动通信1 HYPERLINK l bookmark13 o Current Document 1.2多径时变信道11.2.1信道模型的分类11.2.2时变信道的特点1 HYPERLINK l bookmark19 o Current Document 1.3瑞利信道衰落2 HYPERLINK l bookmark22 o Current Document 二实现框图3 HYPERLINK l bookma

3、rk25 o Current Document 2.1多径时变信道性能仿真实现框图3 HYPERLINK l bookmark28 o Current Document 2.2多径时变信道仿真实现4 HYPERLINK l bookmark31 o Current Document 三详细设计5 HYPERLINK l bookmark34 o Current Document 3.1瑞利信道的特性5 HYPERLINK l bookmark41 o Current Document 3.2多径时变信道的特性8 HYPERLINK l bookmark44 o Current Document

4、 3.3单频信号经过时变信道11 HYPERLINK l bookmark49 o Current Document 3.4数字信号经过多径时变信道13 HYPERLINK l bookmark52 o Current Document 总结15 HYPERLINK l bookmark55 o Current Document 参考文献16 HYPERLINK l bookmark64 o Current Document 附录17 HYPERLINK l bookmark67 o Current Document 致谢29前言在无线移动环境下进行高速可靠通信是具有挑战性的，电波通过物理媒体

5、 传播并与环境中的物体相互作用，因此，无线电波的传播是个复杂过程。在高频 (HF)频段范围内，电磁波经由天波传播时经常发生的问题是信号多径。电磁波的 多径传播主要是因为电磁波经电离层的多次折、反射，电离层的高度不同，电离 层不均匀性引起漫射现象等引起的。当信号的多径发生在发送信号经由传播路径 以不同的延迟到达接收机的时候，一般会引起数字通讯系统中的符号间干扰。而 且，由不同传播路径到达的各信号分量会相互削弱，导致信号能量衰减，造成信 噪比降低。移动无线信道是一个充满复杂干扰的信道。由环境中的各种障碍物所引起的 信号多径传播是其主要特点之一。另一个特点是多普勒效应。由于多径效应和移动台运动等影响

6、因素，使得移动信道对传输信号在时间、 频率和角度上造成了色散，即时间色散、频率色散、角度色散等等，因此多径信 道的特性对通信质量有着重要的影响，而多径信道的包络统计特性则是我们研究 的焦点。根据不同无线环境，接收信号包络一般服从几种典型分布，如瑞利分布、莱 斯分布等。在此专门针对服从瑞利分布的多径信道进行模拟仿真，进一步加深对 多径信道特性的了解。多径信道的基本原理1.1移动通信移动无线信道是一个充满复杂干扰的信道。由环境中的各种障碍物所引起的 信号多径传播是其主要特点之一。同一发射机发射的电磁波向各个方向辐射，不 同的波遇到不同的障碍物发生反射折射以及散射衍射等作用会使得波束到达接 收机的时

7、间、幅度和相位均发生延迟与畸变，例如若发射一个窄脉冲经过无线信 道后将在接收端收到一连串幅度和相位均不同的脉冲串，如果在这个期间内连续 发射多个脉冲，将在接收端产生混叠发生误码的几率大大提高。另一个特点是多 普勒效应。移动通信中的终端基本处于移动的状态，这就导致了电磁波的多普勒 效应。由于到达接收机的杂散波的方向相位均不同所引起的多普勒效应也不尽相 同，更加恶化了接收信号。1.2多径时变信道1.2.1信道模型的分类按照调制信道模型，信道可以分为恒参信道和随参信道两类。部分无线信道 和各种有线信道，包括卫星链路（link）和某些视距传输链路，可以当做恒参信 道看待，因为它们的特性变化很小、很慢，

8、可以视作其参量恒定。恒参信道实际 上就是一个非时变线性网络。1.2.2时变信道的特点随参信道对信号传输的影响，依靠天波传播和地波传播的无线电信道、某些 视距传输信道和各种散射信道就是随参信道。随参信道的特性是“时变”的。例 如，在用天波传播时，电离层的高度和离子浓度随时间、季节、年份而在不断变 化，使信道特性随之变化。在移动通信中，由于移动台在运动，收发两点之间的 传输路径自然也在变化，从而使得信道参量也不断变化。一般来说，各种随参信道具有的共同特性是：（1）信号的传输衰减随时间而变；（2）信号的传输时延随时间而变；信号经过几条路径到达接收端，而且每条路径的长度(时延)和衰减 都随时间而变，即

9、存在多径传播(multipath propagation)现象。多径传播对信号的影响称为多径效应。他对信号传输质量的影响很大。1.3瑞利信道衰落在陆地移动通信中，移动台往往受到各种障碍物和其他移动体的影响，以致 到达移动台的信号是来自不同传播路径的信号之和。而描述这样一种信道的常用 信道模型便是瑞利衰落信道。和振幅恒定、单一频率的发射信号相比，接收信号波形的包络有了起伏，频 率也不再是单一频率，而有了扩展，成为窄带信号，信号包络因传播有了起伏的 现象称为衰落(fading)。多径传播使信号包络产生的起伏虽然比信号的周期缓慢， 但是仍然可能是在秒或秒以下的数量级，衰落的周期常能和数字信号的一个码

10、元 周期相比较，故通常将由多径效应引起的衰落称为快衰落。瑞利衰落信道(Rayleigh fading channel)是一种无线电信号传播环境的统计 模型。这种模型假设信号通过无线信道之后，其信号幅度是随机的，表现为“衰 落”特性，并且多径衰落的信号包络服从瑞利分布。由此，这种多径衰落也称为 瑞利衰落。这一信道模型能够描述由电离层和对流层反射的短波信道，以及建 筑物密集的城市环境。瑞利衰落只适用于从发射机到接收机不存在直射信号的情 况，否则应使用莱斯衰落信道作为信道模型。假设经反射(或散射)到达接收天线的信号为N个幅值和相位均随机的且统 计独立的信号之和。信号振幅为r,相位为，则其包络概率密度

11、函数为1-1所示：r _ r2e 2c 2P(r)= c 2(r z 0)(1-1)相位概率密度函数为1-2所示：P( )=1/2 兀(0 展 V 2兀)(1-2)二实现框图2.1多径时变信道性能仿真实现框图信号经过多径信道后，会产生码间干扰和衰落，其中衰落快慢取决于信道随 时间变化的快慢，而码间干扰的严重程度取决于码元间隔和多径间的时延差的相 对关系。多径效应总的来说有三点，即对单一正弦波产生频域弥散，对宽带信号 频率选择性衰落以及对数字信号产生时域弥散。时变信道仿真实现框图如2-1-1所示：图2-1多径时变信道仿真实现框图2.2多径时变信道仿真实现时变信道是指信道的参数随时间变化的信道，特

12、点是信号的传输衰减随时 间变化；信号的传输时延也是随时间而变的。时变信道对信号传输的影响是使输 入信号的频率弥散。多径信道是输入信号传输的传输路径不止有一条，接收端同时收到来自多 条传输路径的信号，这些信号可能是同相相加或反向相消。多径传播对信号的影 响称为多径效应，会对信号传输质量造成很大的影响。当输入为单频（振幅恒定，频率单一）信号时，经过多径时变信道的传输 后，接收信号的波形包络随时间随机起伏，输出不再是单频信号，而是一个窄带 信号，带宽大小随时变因素的快慢决定。经过多径时变信道传输以后，多径信道 的时延以及衰减均不相同，导致接收信号的幅度不同，频率也增多了。当输入为数字信号时，经过多径

13、时变信道，针对数字信号体现在码间干扰上。 由于各径时延不同，通过个路径的衰减也是不同的，信号经过多条路径后到达接 收端形成码间干扰。三详细设计本次课程设计中多径信道的信道模型建立，性能分析，瑞利衰落的多径信 道模型等仿真分析都使用了 MATLAB平台。MATLAB的数据分析和处理功能 十分强大，运用它来进行语音信号的分析、处理和可视化相当便捷。而且编程易 学、直观，代码非常符合人们的思维习惯。MATLAB几乎可以在各种机型和操 作系统上运行，所以在可移植性和可扩充性上MATLAB远优越于其他的高级编 程语言。3.1瑞利信道的特性瑞利分布分析幅度与相位的分布特性：包络r服从瑞利分布，。在02n内

14、服从均匀分布。瑞利分布的概率分布密 度如图3-1-2所示：多径衰落信道基本模型离散多径衰落信道模型为3-1所示：N (t)J(t) = r (t)%(t -t )k=i * k(3-1)r (t)t其中，k 复路径衰落，服从瑞利分布；k是多径时延。多径衰落信道模型框图如图3-1-1所示：图3-1-1多径衰落信道模型框图产生服从瑞利分布的路径衰落r(t)利用窄带高斯过程的特性，其振幅服从瑞利分布，其振幅可以使用3-2所示的公式求出:(t)=w (t)2+n (t)2(3-2)上式中气、ns，分别为窄带高斯过程的同相和正交支路的基带信号。模型的适用瑞利衰落模型适用于描述建筑物密集的城镇中心地带的无

15、线信道。密集的建 筑和其他物体使得无线设备的发射机和接收机之间没有直射路径，而且使得无线 信号被衰减、反射、折射、衍射。通过电离层和对流层反射的无线电信道可用瑞 利衰落来描述，因为大气中存在的各种粒子能够将无线信号大量散射。瑞利衰落 属于小尺度的衰落效应，它总是叠加于如阴影、衰减等大尺度衰落效应上。信道衰落的快慢与发射端和接收端的相对运动速度的大小有关。相对运对导 致接收信号的多普勒频移。瑞利分布的概率密度曲线及瑞利信道的包络如图3-1-2所示：BO Figure 1File Edit View Tnseirt T o o IsDesktop WindcuA/ HelpD CJ d 怀七孩巧0

16、 X威j槌n目丁回图3-1-2瑞利分布的概率密度曲线及瑞利信道的包络移动台速度变化时对信道的瑞利衰落会产生一定的影响，图3-1-3所示为V=30千米/小时以及V=120千米/小时时的信道曲线图:图3-1-3不同速度下的瑞利信道曲线从图3-1-3中可以看出，速度越大对信道瑞利衰落影响就会越大。3.2多径时变信道的特性移动无线信道是一个充满复杂干扰的信道。由环境中的各种障碍物所引起 的信号多径传播是其主要特点之一。多径导致频率选择性，在同一位置，由于反射径信号的存在，发射不同频 率的信号时，在接收机处接收到信号有的频率是被增强了，有的频率是被削弱了。 频率选择性由此产生。把那些受到的影响基本一致的

17、频率范围叫做相干带宽。时变信道是指信道的参数随时间变化的信道，特点是信号的传输衰减随时 间变化；信号的传输时延也是随时间而变的；信号经过多条路径到达接收端，而 且每条路将的长度（时延）和衰减都是随时间而变化的。时变信道对信号传输的 影响是使输入信号的频率弥散。多径信道是输入信号传输的传输路径不止有一条，接收端同时收到来自多 条传输路径的信号，这些信号可能是同相相加或反向相消。多径传播对信号的影 响称为多径效应，其会对信号传输质量造成很大的影响。多普勒频移反映了信道的时变性，多普勒频移越大，信道的相干时间就越 短，合成信号包络变化越快。相干时间只是说信道在这段时间内特性基本不变。 至于这段时间内

18、是增强信号还是消弱信号则没有体现。（1）在其它参数不变的情况下画出移动台距离基站初始距离r0=3000时的接收信号情况，如图3-2-1所示:图3-2-1移动台距离基站初始距离为r0=3000的接收信号蓝色线：直射径的信号；绿色线：反射径的信号；红色线：移动台接收到的 第1径和第2径的合成信号，从图3-2-1可以看出，即使移动台是静止的（V=0）， 由于反射径的存在，使得接收到的合成信号最大值要小于直射径的信号。（2）移动台在不同位置的多径信号移动台距离基站初始距离分别设置为r0=1000，r0=9000，r0=14000其它条 件保持不变时，接收到的信号的情况如图3-2-2所示：1移动蹦 距离

19、基站初始B距离为r0=10001甘1 111.511 1110.50.50 .0-0.5-0.5,1 1 | r-11移动削基站初O距离为r0=9000-1-1.5|-100.20.40.60.8100.20.40.60.81-100.20.40.60.81图 3-2-2r0取不同值时接收信号的情况从图3-2-2中可得出结论：使移动台静止，由于反射径的存在，使接收信号 要比没有反射径时的信号弱，衰落由此产生。（3）不同频率的信号经过多径信道f=3e8, f=9e8, f=27e8，接收信号的情况如图3-2-3所示:420-2-2发射信号频率亍=98 .-400.20.40.60.81-400.

20、20.40.60.81x 108图 3-2-3f取不同值时接收信号的情况若f逐步变大，有些频率被削弱，f充分大时，看出合成信号被削弱了，那 些受到影响基本保持一致的频率范围称为相干带宽，在同一位置，由于反射镜信 号的存在，发射不同信号的频率，在接收机收到的信号有的被加强，有的信号被 减弱，这就是所谓的频率选择性衰落。（4）仿真移动台不同速度的信号的多径信号改变移动台速度不同的信号经过多径信道v=0，v=300其它条件不变时接收信 号的情况如图3-2-4所示：图3-2-4 V值不同时接收信号的情况移动台有速度时，发现即使同一频率，同一位置，在不同的时间点，合成信 号的强度也是不一样的，有的地方信

21、号衰减，有的地方信号增强。当速度由0 增加到300时，直射径信号减弱，反射径信号增强，合成信号减弱。3.3单频信号经过时变信道多径传播对信号的影响称为多径效应，其会对信号传输质量造成很大的影 响。当输入为单频信号时，经过多径时变信道的传输后，接收信号的波形包络随 时间随机起伏，输出不再是单频信号，而是一个窄带信号，带宽大小随上事变因 素的快慢决定。经过多径时变信道传输以后，多径信道的时延以及衰减均不相同， 导致接收信号的幅度不同，频率也增多了。一个幅度为1,频率为10Hz的单频信号经过20条路径传输得到的波形 及频谱，并且这20条路径的衰减相同，但时延的大小随时间变化，每径的时延 变化规律为正

22、弦型，变化的频率是从02Hz随机均匀抽取的。输入的单频正弦信号时域及频域图如3-3-1所示：图3-3-1单频信号的时域及频域图从图3-3-1中可看出，原信号为单一正弦信号，时域标准正弦，频域单一冲 激。信号的幅值为1，频率为10Hz。信号经多径传播后在接收端所得信号时域及频域图图3-3-2经过20径后信号的时域及频域图由图3-3-2可知，过多径传播后，不同路径时延不同，不同时延的信号叠加， 导致时域图形不再是单一正弦，出现了其它的频率。（由于采样频率的关系，看 起来还是光滑曲线）频域出现了毛刺，即频域扩散。3.4数字信号经过多径时变信道多径时变信道对数字信号产生时域弥散。数字信号经过多径时变信

23、道，发生 码间干扰上。由于各径时延不同，通过个路径的衰减也是不同，信号经过多条路 径到达接收端形成码间干扰。一条三径传输的信道s （t） = Z3 u b （t 一 t ,），其参数如下:i = 1u = 0.5, u = 0.707, u = 05；t = 0,t = 1,t = 2123123画出信道的幅频响应和相频响应如3-4-1所示:图3-4-1三径信道的幅频、相频响应图画出Ts=1时输出信号波形及输出信号的幅度谱如图3-4-2所示:图3-4-2数字信号进入三径信道输出信号的波形及幅度谱画出Ts=4时输出信号波形及输出信号的幅度谱如图3-4-3所示:图3-4-3改变Ts后画出信号波形及

24、输出信号幅度谱由图3-4-2及3-4-3可以看出，随着Ts的增大，输出信号的波形与输入信 号的波形越接近。因为信道幅度特性不是很理想，会造成输入信号失真，针对模 拟信号体现在波形失真上;针对数字信号体现在码间干扰上。由于各径时延不同， 通过个路径的衰减也是不同的，信号经过多条路径后到达接收端形成码间干扰。总结本次课程设计中多径信道的信道模型建立，性能分析，瑞利衰落的多径信道 模型等仿真分析都使用了 MATLAB平台。从仿真图形可以形象的反应出多径时 变信道的特点及性能。根据不同无线环境，接收信号包络一般服从几种典型分布， 如瑞利分布、莱斯分布等。在这次课程设计中专门针对服从瑞利分布的多径信道

25、进行模型仿真，进一步加深对多径信道特性的了解。信号经过多径信道后，会产生码间干扰和衰落，其中衰落快慢取决于信道 随时间变化的快慢，而码间干扰的严重程度取决于码元间隔和多径间的时延差的 相对关系。时变信道是指信道参数随时间变化的信道，特点是信号的传输衰减随时间 变化；信号的传输时延也是随时间而变的。时变信道对信号传输的影响是使输入 信号的频率弥散。多普勒频移反映了信道的时变性，多普勒频移越大，信道的相 干时间就越短，合成信号包络变化越快。多径效应总的来说有三点，即对单一正弦波产生频域弥散，对宽带信号频率 选择性衰落以及对数字信号产生时域弥散。单一频率的信号经多径时变信道传输 后，由于不同路径时延

26、不同，不同时延的信号叠加，导致时域图形不再是单一正 弦，从频谱图中可以看出，它的频带展宽出现了其它频率，即就是频域弥散。宽 带信号经过多径时变信道，在同一位置，由于反射径信号的存在，发射不同频率 的信号时，接收机处收到信号有的频率被增强了，有的频率被削弱了。频率选择 性由此产生。数字信号经过多径时变信道，发生码间干扰上。由于各径时延不同， 通过个路径的衰减也是不同，信号经过多条路径到达接收端形成码间干扰。通过理论与仿真的结合，我更加深入的了解了多径时变信道的性质和特点， 带给我很大的收获。参考文献覃团发、姚海涛、覃远年、陈海强.移动通信.重庆大学出版社23:48Gordon.Stuber .移

27、动通信原理.电子工业出版社 51:56宋荣方.矢量多径信道的衰落相关特性.南京邮电学院学18:22樊昌信.通信原理M .第5版.北京：国防工业出版社，2001. 44:48Bernard Sklar.数字通信一基础与应用.徐平平、宋铁成,译.北京:电子工 业出版社，2002. 77:80梁斌、朱洪波.移动通信Rician信道中的多普勒影响分析66:69附录fc=900*10.人6;wc=2*pi*fc;v1=30*1000/3600;c=300*10人6;wm=wc*(v1/c);fm=wm/(2*pi);N =128*100;%Carrier frequency%接收端速度%Maximum

28、shift%Doppler shift% generate Doppler power spectrumdeltaf = 2*fm/(N-1);T = 1/deltaf;sf0 = 1.5/(pi*fm);for n = 1:(N-2)/2sf(n) = 1.5/(pi*fm*sqrt(1-(n*deltaf/fm)A2);end classicf = fliplr(sf),sf0,sf;% generate two normally distributed random variablesgaussN_re1 = randn(1,(N-2)/2);gaussN_im1 = randn(1,(

29、N-2)/2);gaussN_pos1 = gaussN_re1 + i*gaussN_im1;gaussN_neg1 = conj(gaussN_pos1);gaussN1 = fliplr(gaussN_neg1),0,gaussN_pos1;gaussN_re2 = randn(1,(N-2)/2);gaussN_im2 = randn(1,(N-2)/2);gaussN_pos2 = gaussN_re2 + i*gaussN_im2;gaussN_neg2 = conj(gaussN_pos2);gaussN2 = fliplr(gaussN_neg2),0,gaussN_pos2;

30、% generating flat Rayleigh fading channelx = ifft(sqrt(classicf).*gaussN1);y = ifft(sqrt(classicf).*gaussN2);rayleigh_amp = sqrt(abs(x).A2+abs(y).A2);rayleigh_db = 20*log10(rayleigh_amp);figure(1);subplot(211)plot(rayleigh_db);title(瑞利分布的包络) r = sqrt(0.5*(gaussN_re1.A2 + gaussN_re2.A2);step = 0.1; r

31、ange = 0:step:3;h = hist(r, range);fr_approx = h/(step*sum(h);fr = (range/0.5).*exp(-range.A2);subplot(212)plot(range, fr_approx,ko, range, fr,k);title(瑞利分布的概率密度曲线)grid;function h=rayleigh(fd,t)fc=900*10A6;v1=30*1000/3600;c=3*10A8;fd=v1*fc/c;ts=1/10000;t=0:ts:1;h1=rayleigh(fd,t);v2=120*1000/3600;fd=

32、v2*fc/c;h2=rayleigh(fd,t);%产生瑞利衰落信道%选取载波频率%移动速度v1=30km/h%定义光速%多普勒频移%信道抽样时间间隔%生成时间序列%产生信道数据%移动速度v2=120km/h%多普勒频移%产生信道数据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(功率) functi

33、on 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;for n=1:M

34、%第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);clear allf=9e8;v=0;c=3e8;r0=3000;d=15000;%移动台速度，静止情况为0

35、%电磁波速度，光速%移动台距离基站初始距离 %基站距离反射墙的距离t1=0:0.00000000005:0.00000002;% 时间%直射径信号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); % 反射径信号figureplot(t1,E1,t1,E2,-g,t1,E1-E2,-r)%画出直射径、反射径和总的接收信号legend(直射径信号,反射径信号,移动台接收的合成信号)figureplot(t1,E1-E2) clear allf=9e8;v=

36、0;c=3e8;r0=1000;d=15000;%发射信号频率%移动台速度，静止情况为0%电磁波速度，光速 %移动台距离基站初始距离 %基站距离反射墙的距离%时间%直射径信号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,1)plot(t1,E1,t1,E2,-g,t1,E1-E2,-r)%画出直射径、反射径和总的接收信号legend(直射径信号,

37、反射径信号,移动台接收的合成信号)%axis(0 12 -0.5 0.5)subplot(2,3,4)plot(t1,E1-E2)f=9e8;%发射信号频率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,2)plot(t1,E1,t1,E2,-g,t1,E1-E2,-r)%axis(0 1

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

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

40、-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;%移动台速度，静止情况为0c=3e8;%电磁波速度，光速r0=3000;%移动台距离基站初始距离d=15000;%

41、基站距离反射墙的距离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 0.5)subplot(2,3,5)； plot(t1,E1-E2)f=27e8;%发射信号频率v=0;%移动台速度，静止情况为0c

42、=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,-g,t1,E1-E2,-r)%axis(0 12 -0.5 0.5)subplot(2,3,6)plot(

43、t1,E1-E2)clear allf=90;%发射信号频率v=0;%移动台速度，静止情况为0c=3e8;%移动台距离基站初始距离d=150;%基站距离反射墙的距离t1=0:0.0005:0.1;%时间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,2,1)plot(t1,E1,t1,E2,-g,t1,E1-E2,-r)%画出直射径、反射径和总的接收信号legend(直射径信号,反射径信号,移动台接收的合

44、成信号)subplot(2,2,2)plot(t1,E1-E2)clear allf=90;%发射信号频率v=300;%移动台速度，静止情况为0c=3e8;%电磁波速度，光速r0=30;%移动台距离基站初始距离d=150;%基站距离反射墙的距离t1=0:0.0005:0.1;%时间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,2,3)plot(t1,E1,t1,E2,-g,t1,E1-E2,-r)%画出

45、直射径、反射径和总的接收信号legend(直射径信号,反射径信号,移动台接收的合成信号)subplot(2,2,4)plot(t1,E1-E2)clc;close all;clear 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:Lwdelay(i,:)=cos(2*pi*fdelay(i)*t);s(i,:)=decay*cos(2*pi*f*t+wdelay(i,:)+wdelay0(i);endy=sum(s)/

46、sqrt(L);figure(1);subplot(211);plot(t,x);xlabel(t/s);ylabel(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);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径后信号频域图);axi