莱斯信道中QPSK信号Matlab仿真_第1页
莱斯信道中QPSK信号Matlab仿真_第2页
莱斯信道中QPSK信号Matlab仿真_第3页
莱斯信道中QPSK信号Matlab仿真_第4页
莱斯信道中QPSK信号Matlab仿真_第5页
已阅读5页,还剩2页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、1 移动无线信道的定义及分类各类信号从发射端发送出去以后,在到达接收端之前经历的所有路径统称为信道。如果其中传输的是无线电信号,电磁波所经历的路径称之为无线信道。与其他通信信道相比,无线信道是最为复杂的一种。无线传播环境是影响无线通信系统的基本因素。发射机与接受机之间的无线传播路径,因从经历简单的视距传播,到遭遇各种复杂的地物(如建筑物、山脉和树林等)所引起的反射、绕射和散射传播等而显得非常复杂。另外,移动台相对于发射台移动的方向和速度,甚至收发双方附近的移动物体也对接受信号有很大的影响。因此,这使得无线信道具有极度的随机性。移动通信信号在空间传播中所经历的衰落大体可以分为2类, 即大尺度衰落

2、和小尺度衰落,如图1所示。大尺度衰落是因为发射机与接收机之间的距离和两者之间障碍物(如山丘、森林、建筑物等)的遮蔽影响而造成的信号强度的衰减,它反映了移动信号在较大区域中的平均能量的减少或称为路径损失。而小尺度衰落是指当移动台在一个较小的范围运动时,引起接收信号的幅度、相位和到达角等的快速变化。图1 移动无线信道的分类信号在传播的过程中,受各种环境的影响会产生反射、衍射和散射,这样就使得到达接收机的信号是许多路径信号的叠加,因而这些多径信号的叠加在没有视距传播情况下的包络服从瑞利分布。当多径信号中包含一条视距传播路径时,多径信号就服从莱斯分布。根据信号多径附加时延的大小,小尺度衰落又可以分为平

3、坦衰落和频率选择性衰落。另外,由于移动台的移动性而导致接收到的信号产生多普勒频移(频率色散),根据多普勒扩展的大小,信道又可以分为快衰落信道和慢衰落信道。2 瑞利分布和莱斯分布的特性在实际情况中对数字通信系统来说,调制符号的周期比由多径传播引起的时延扩展要大,因此在一个符号周期内的所有频率分量都会经历相同的衰减和相移。信道对于所有频率分量来说是平坦的, 因而定义这类信道为平坦衰落信道。理论分析和实测试验结果表明:平坦衰落的幅度在大多数情况下,符合瑞利分布(ray leighdistribution)或莱斯分布(ricedistribution)。由于移动通信信道的复杂性,其仿真一般是以平坦衰落

4、信道建模为基础的,然后在此基础上,再对频率选择性信道等进行建模和仿真。用Matlab中能生成伪随机序列的randn语句就可以得到期望的莱斯衰落序列。瑞利衰落序列可以由K=0得到。图2是一个当K=7dB时典型的莱斯衰落信号包络,衰落幅度用分贝表示。图2 当K =7dB时莱斯衰落信号的包络3 MATLAB仿真接下来,利用Matlab对莱斯分布的累积分布函数(CDF)进行近似估计。莱斯分布的累积分布函数是通过迭代法得到的,在每一步的迭代中利用Matlab中的find和length函数来得到符合要求的衰落序列,并使用上面产生莱斯分布的M文件rice_fading.m得到K=7dB时的莱斯分布的累积分布

5、函数(CDF)的近似估计,如图3所示。然后通过Matlab中的hist函数得到的瑞利分布PDF的估计值与式子分析求得的PDF进行比较,结果如图4所示,所得的估计值与式子分析求得的PDF非常的接近。程序如下:function r = rice_fading(Kdb, N,Mi)K = 10(Kdb/10); const = 1/(2*(K+1);x = randn(1,N); y = randn(1,N);r = sqrt(const*(x + sqrt(2*K).2 + y.2);rt = zeros(1,Mi*length(r); ki = 1;for i=1:length(r)rt(ki:

6、i*Mi) = r(i); ki = ki+Mi;endr = rt;Kdb=7; N=100000; Mi=1;r=rice_fading(Kdb,N,Mi);RdB = 20*log10(r);Rt = min(RdB):max(RdB);for m = 1:length(Rt)fade = find(RdB < Rt(m);Nm = length(fade); AF(m) = Nm/N;endsemilogy(Rt,AF,'k-o'); set(gcf,'paperunits','centimeters');set(gcf,'

7、papersize',5 5); 设置图像大小为5cm*5cmgrid;运行结果为:图3 K=7dB时莱斯分布的CDF程序如下:N = 100000;x = randn(1, N); y = randn(1, N);r = sqrt(0.5*(x.2 + y.2);step = 0.1; range = 0:step:3;h = hist(r, range);fr_approx = h/(step*sum(h);fr = (range/0.5).*exp(-range.2);plot(range, fr_approx,'ko', range, fr,'k'

8、;);set(gcf,'paperunits','centimeters');set(gcf,'papersize',5 5); 设置图像大小为5cm*5cmgrid;图4瑞利分布的PDFQPSK信号传输过程的框图如图5所示。图5 QPSK 信号在莱斯衰落信道中传输过程的框图在仿真过程中主要用到了Matlab工具箱中这样一些函数:产生同分布随机变量矩阵的randint函数,基带数字调制器dmodce函数,基带数字解调器ddcmodce函数,以及计算误比特数和误比特率的biterr函数。最后得到当K=5dB与K=-inf时的误比特率(BER)的仿真

9、曲线,如图6所示。对于莱斯衰落信道,由于多径信号中含有视距传播分量,从图6可以看出随着K值的增大,接收信号中直视波的能量增大,使得误码率的性能大为改善。分析可知当K值逐渐增大时,所得误码率的曲线会越来越接近理想曲线。程序如下:function r = ray_doppler(fm, M, dt, N)T = N*dt-dt; t = 0:dt:T;c = sqrt(2/M); w = 2*pi*fm; x = 0; y = 0;for n = 1:Malpha = (2*pi*n-pi+(2*pi*rand-pi)/(4*M);ph1 = 2*pi*rand - pi;ph2 = 2*pi*r

10、and - pi;x = x + c*cos(w*t*cos(alpha) + ph1);y = y + c*cos(w*t*sin(alpha) + ph2);endr = sqrt(x.2 + y.2)/sqrt(2);function BER = psk_rice(logEbNo,M,Kdb)k = log2(M); EbNolin = 10.(logEbNo/10)Fd = 1; Fs = 1; Nit = 100000; Ns = 1000;Tstop = 100;for it_snr = 1:length(EbNolin)nstd = sqrt(1/(2*k*EbNolin(it_

11、snr);bit_err = 0;for it = 1:Nita = randint(Ns,1,M);s = dmodce(a,Fd,Fs,'psk',M);r = rice_fading(Kdb,Ns,1); v_r = s.*r'v = v_r + nstd*(randn(Ns,1) + .j*randn(Ns,1);z = ddemodce(v, Fd, Fs,'psk',M);errors = biterr(a,z);bit_err = bit_err + errors;if bit_err >= TstopbreakendendBER(it_snr) = bit_err/(it*Ns*k);Endsemilogy(logEbNo, BER,'k');xlabel('Eb/No(dB)');ylabel('BER');set(gcf,'paperunits','centimeters');set(gcf,'papersize',5 5); 设置图像大小为5cm*5cmgrid;图6 莱斯衰落信道中QPSK信号传输性能的仿真4 结 论对移动无线通信系统而言, 因为传播环境的复杂性和

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论