目录摘要 11、设计原理 2 1.1设计目的 2 1.2仿真原理 2 1.2.1瑞利分布简介 2 1.2.2多径衰落信道基本模型 2 1.2.3产生服从瑞利分布的路径衰落r(t 3 1.2.4产生多径延时 4 1.3仿真框架 4 2、设计任务 4 2.1设计任务要求 4 2.2 MATLAB 仿真程序要求 4 3、DSB 调制解调分析的MATLAB 实现 53.1 DSB调制解调的 MATLAB 实现 5 3.2瑞利衰落信道的MATLAB实现 6 4、模拟仿真及结果分析 7 4.1模拟仿真 7 4.1.1多普勒滤波器的频响 7 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 11 页1 / 11 4.1.2多普勒滤波器的统计特性 7 4.1.3信道的时域输入 /输出波形 8 4.2仿真结果分析 84.2.1时域输入 /输出波形分析 84.2.2频域波形分析 8 4.2.3多普勒滤波器的统计特性分析 95、小结与体会 96、参考文献 9MATLAB 通信仿真设计摘要主要运用 MATLAB 进行编程，实现采用对输入信号进行抑制载波的双边带调幅；而后将调幅波输入信道，研究多径信道的特性对通信质量的影响；最后将信道内输出的条幅波进行同步解调，解调出与输入信号波形相类似的波形，精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 11 页2 / 11 观测两者差别。同时输出多普勒滤波器的统计特性图及信号时域和频域的输入、输出波形。关键字： 双边带调幅瑞利衰落相干解调MATLAB1、设计原理1.1设计目的由于多径和移动台运动等影响因素，使得移动信道对传输信号在时间、频率和角度上造成了色散，如时间色散、频率色散、角度色散等等，因此多径信道的特性对通信质量有着至关重要的影响，而多径信道的包络统计特性成为我们研究的焦点。根据不同无线环境，接收信号包络一般服从几种典型分布，如精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 11 页3 / 11 瑞利分布、莱斯分布和Nakagami-m 分布。在设计中，专门针对服从瑞利分布的多径信道进行模拟仿真，进一步加深对多径信道特性的了解。1.2仿真原理1.2.1瑞利分布简介1）环境条件：通常在离基站较远、反射物较多的地区，发射机和接收机之间没有直射波路径，存在大量反射波；到达接收天线的方向角随机且在02 ）均匀分布；各反射波的幅度和相位都统计独立。 其中复路径衰落，服从瑞利分布；是多径时延。多径衰落信道模型框图如图 2 所示：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 11 页4 / 11 图 2 多径衰落信道模型框图1.2.3产生服从瑞利分布的路径衰落r(t利用窄带高斯过程的特性，其振幅服从瑞利分布，即 2）上式中，、分别为窄带高斯过程的同相和正交支路的基带信号。首先产生独立的复高斯噪声的样本，并经过FFT 后形成频域的样本，然后与 S。如下图 3所示: 图3 瑞利衰落的产生示意图其中， Average power (dB 1 0 0 2 310 -1.0 3 710 -9.0 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页，共 11 页5 / 11 4 1 090 -10.0 5 1 730 -15.0 6 2 510 -20.0 1.3仿真框架根据多径衰落信道模型 见图 2），利用瑞利分布的路径衰落见图 3）和多径延时参数见表 1），我们可以得到多径信道的仿真框图，如图4 所示：图4 多径信道的仿真框图2、设计任务2.1设计任务要求1）查找资料，了解瑞利衰落信道模型的分类，结合某种模型，掌握瑞利分布的多径信道仿真原理，用MATLAB 仿真实现瑞利分布的多径信道的仿真；2）根据已学的知识，实现一种基带信号的模拟调制并做出仿真；3）结合 1）2）步，观察已调信号通过瑞利信道后的时域波形图和频谱图；4）对仿真结果做适当分析。2.2MATLAB 仿真程序要求1）参数设计准确、合理；2）关键语句加注释；。%DSB 调制SignalInput=sin(t/50。%+cos(t/65。 % 调制信号精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 11 页6 / 11 c=cos(0.2*pi*t。% 载波信号y_in=SignalInput.*c。 % 调制delay=0 31 71 109 173 251。%10nspower=0 -1 -9 -10 -15 -20。%dBy_in=zeros(1,delay(6 y_in。% 为时移补零y_out=zeros(1,LengthOfSignal。% 存放经信道未解调的信号+1-delay(i:(delay(6+LengthOfSignal-delay(i*10(power(i/20。end 。% S(t*cos(w*t=m(t*cos(w*t*cos(w*t=0.5*m(t*(1+cos(2*w*t% 用一个低通滤波器将上式中的第一项和第二项分离，无失真的恢复出原始的调制信号。% 这种调制方法又称为同步解调或相干解调% 同步解调y_out_end=y_out.*c。% 同步解调或相干解调% 低通滤波wp=0.1*pi。ws=0.12*pi。Rp=1 。As=15 。N,wn=buttord(wp/pi,ws/pi,Rp,As。b,a=butter(N,wn。y_out_end =filter(b,a,y_out_end。% 滤波y_out_end =2* y_out_end。% 恢复幅度% 原信号的频谱K=fft(SignalInput。%DSB 调制后信号的频谱L=fft(y_in。%y_out 的频谱 。% 最终解调的频谱N=fft(y_out_end。% 输出figure(1。subplot(4,2,1。plot(SignalInput(delay(6+1:LengthOfSignal。axis(0,3000,-2,2。title( 原始输入信号 。subplot(4,2,2。plot(abs(fftshift(K。 axis(4900,5100,0,6000。title( 原始输入信号的频谱 。subplot(4,2,3。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页，共 11 页7 / 11 plot(y_in(delay(6+1:LengthOfSignal。 axis(0,3000,-2,2。%去除时延造成的空白信号title( 进入瑞利信道前，DSB调制后的信号 。subplot(4,2,4。plot(abs(fftshift(L。 axis(3500,6500,0,3000。title( 进入瑞利信道前，DSB调制后的信号的频谱 。subplot(4,2,5。plot(y_out(delay(6+1:LengthOfSignal。axis(0,3000,-0.08,0.08。% 去除时延造成的空白信号title( 经瑞利信道后，DSB解调前的信号 。subplot(4,2,6。plot(abs(fftshift(M。axis(3500,6500,0,100。title( 经瑞利信道后，DSB解调前的信号的频谱 。subplot(4,2,7。plot(y_out_end(delay(6+1:LengthOfSignal。axis(0,3000,-0.08,0.08。%去除时延造成的空白信号title( 最终解调后的信号 。subplot(4,2,8。plot(abs(fftshift(N。axis(4900,5100,0,200。title( 最终解调后的信号的频谱 。figure(2。subplot(3,1,1。hist(r,256。%绘制直方图title( 瑞利信道的幅度分布 subplot(3,1,2。hist(angle(r0。title( 瑞利信道的相位分布 。subplot(3,1,3。plot(Sf1。title( 多普勒滤波器的频响特性 。3.2瑞利衰落信道的 MATLAB 实现%Rayl.m 参考【 1】f=1:2*fm-1。%通频带长度y=1.5./(1-(f-fm/fm.2.(1/2/pi/fm。% 多普勒功率谱 ( 基带 图3Sf=zeros(1,LengthOfSignal。Sf1=y。Sf(fc-fm+1:fc+fm-1=y。%(把基带映射到载波频率x1=randn(1,LengthOfSignal。x2=randn(1,LengthOfSignal。nc=ifft(fft(x1+1i*x2.*sqrt(Sf。%同相分量% 首先产生独立的复高斯噪声的样本，并经过FFT 后形成频域的样本，% 然后与 Sx3=randn(1,LengthOfSignal。x4=randn(1,LengthOfSignal。ns=ifft(fft(x3+1i*x4.*sqrt(Sf。%正交分量r0=(real(nc+1i*real(ns。 %瑞利信号r=abs(r0。%瑞利信号幅值nc 、ns 分别为窄带高斯过程的同相和正交支路的基带信号）4、模拟仿真及结果分析4.1模拟仿真4.1.1多普勒滤波器的频响图5多普勒滤波器的频响4.1.2多普勒滤波器的统计特性图6 多普勒滤波器的统计特性4.1.3信道的时域输入 /输出波形精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共 11 页9 / 11 图7信道的时域 /频域输入 /输出波形4.2仿真结果分析4.2.1时域输入 /输出波形分析次实验主要是通过MATLAB 仿真瑞利衰落信道的传输过程，通过双边带调幅的调制与解调实现信号的传输。正如右图所示：图中第一、第二个波形是在进入瑞利衰落前，第三、第四个波形是在进入瑞利衰落后，有明显的噪声的存在。由第一个图输入，第四个图输出，信号的传递在存在干扰的情况下基本实现。第一个波形到第二个波形是实现了抑制载波的双边带调幅；第三个波形到第四个波形是运用同步解调或相干解调实现对载波信号的解调功能。图 8 信道的时域输入 /输出波形4.2.2频域波形分析精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页，共 11 页10 / 11 图 9 信道的频域变化分析图 9 中的第二、第四幅图可以发现输出信号的频谱图上一段频率内出现了多余的小频率，但总体频率没有很大变化。说明信道对输入波形存在影响，但本题中，输出的波形还是可以基本反映输入情况的。4.2.3多普勒滤波器的统计特性分析如图 6 中显示出了瑞利信道的幅度分布和相位分布情况，而在简介中提到包络 r 服从瑞利分布， 在 02 内服从均匀分布。瑞利分布的概率分布密度如图 1，经过分析对照可以发现，输出的多普勒滤波器的统计特性完全符合要求。5、小结与体会这次实习是自由定题，当听到这个时让我产生一种不知所措感，因为，我从没有真正讲学习的东西运用到稍微大一点的实际中。不过，越是这样才能让我们锻炼的更多，在用程序编出想要的波形时，我们不仅能体会成功的快乐，而且对通信的传输原理有了更加清晰的认识，这次设计要求我们学会综合运用我们所学的通信知识及用Matlab 实现，可以上网查资料，也可以图书馆里找。设计一开始的时候我充满了恐惧感，面对好多的书，好多的代码，好混乱的脑中知识总之，一万个不想做，但是想到如果不实践是不会有进步的，只得硬着头皮来，然而随着一点一点进行，发现一切都不像想象的那么艰巨，慢慢的我完成了我的课程设计瑞利衰落信道的仿真。其实做什么事都是那样，不开始，永远害怕打破那层膜，只用勇敢的迈出第一步，一切问题都会在面临的时候解决。这次设计最深刻的感受就是：只要开始，一切都可以变得简单。6、参考文献1胡宴如，耿苏燕主编，高频电子线路，北京：高等教育出版社，2009。2沈卫康，宋宇飞，宋红梅，数字信号处理，北京：清华大学出版社，2018。3樊昌信，曹丽娜编著，通信原理，北京：国防工业出版社，2018。4潘子宇， Matlab 通信仿真设计指导书，南京工程学院，2018。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页，共 11 页