无线通信信道建模:时变与多径模型详解与比较(以及接收信号什么情况下能做N点FFT)

在无线通信的各个应用场景下,对无线信道建模往往是开展研究的第一步。

然而,对于无线信道建模的方式有多种多样。在什么条件下,该选择什么样的建模方式让人头疼。太复杂的建模难以推进研究,太简化的模型又无法表达信道特性。

本文针对 h ( τ ) h(\tau) h(τ) h ( t , τ ) h(t,\tau) h(t,τ) h ( n ) h(n) h(n) ∑ h i ( t ) \sum h_i(t) hi(t) ∑ h i ( τ ) \sum h_i(\tau) hi(τ) H H H(Töplitz矩阵)等多种常见的无线信道建模方式进行一一解析,分析其应用场景的全面性,从覆盖场景的全面化角度重新排序

此外,还针对经典的问题**“时域卷积=频域乘积”**进行讨论:

接收信号 y = h ∗ x + n y=h*x+n y=hx+n究竟在什么条件下,可以无失真的表示为 Y = H X + N Y=HX+N Y=HX+N

一、无线信道建模

以下信道模型按全面性降序排列

h ( t , τ ) h(t,\tau) h(t,τ)——时变信道冲激响应

h ( t , τ ) h(t,\tau) h(t,τ)是时变信道的冲激响应,表示信道在时间 t t t 和时延 τ \tau τ 上的响应。是最全面的模型
输出信号: y ( t ) = x ( t ) ∗ h ( t , τ ) = ∫ x ( t − τ ) h ( t , τ ) d τ y(t) = x(t) * h(t,\tau) = \int x(t - \tau) h(t,\tau) d\tau y(t)=x(t)h(t,τ)=x(tτ)h(t,τ)dτ

h ( t , τ ) h(t,\tau) h(t,τ)是关于时间 t t t 和时延 τ \tau τ 的二维函数,该怎么理解呢?

  • 在每个时刻 t 0 t_0 t0,信道响应 h ( t 0 , τ ) h(t_0,\tau) h(t0,τ)随着不同径的时延 τ \tau τ的变化都不同,捕获了多径效应
  • 对于任意一条径 τ 0 \tau_0 τ0,信道响应 h ( t , τ 0 ) h(t,\tau_0) h(t,τ0)随着时间 t t t的变化而变化,捕获了时变效应

∑ i h i ( t ) δ ( τ − τ i ) \sum_i h_i(t) \delta(\tau - \tau_i) ihi(t)δ(ττi)——时变多径信道

∑ i h i ( t ) δ ( τ − τ i ) \sum_i h_i(t) \delta(\tau - \tau_i) ihi(t)δ(ττi)是时变信道的离散多径模型,增益 h i ( t ) h_i(t) hi(t) 随时间变化,时延 τ i \tau_i τi 固定。
输出信号:
y ( t ) = ∑ i h i ( t ) x ( t − τ i ) y(t) = \sum_i h_i(t) x(t - \tau_i) y(t)=ihi(t)x(tτi)
∑ i h i ( t ) δ ( τ − τ i ) \sum_i h_i(t) \delta(\tau - \tau_i) ihi(t)δ(ττi) h ( t , τ ) h(t,\tau) h(t,τ)的一种特例,它表示信道的多径是有限且固定的,每条路径的增益随时间变化。

h ( τ ) h(\tau) h(τ)——时不变信道冲激响应

h ( τ ) h(\tau) h(τ)表示时不变信道的冲激响应,仅依赖时延 τ \tau τ,它表示,无论时间 t t t如何变化,信道在每条径上的响应都是不变的。

输出信号:
y ( t ) = x ( t ) ∗ h ( τ ) = ∫ x ( t − τ ) h ( τ ) d τ y(t) = x(t) * h(\tau) = \int x(t - \tau) h(\tau) d\tau y(t)=x(t)h(τ)=x(tτ)h(τ)dτ

h ( τ ) h(\tau) h(τ)是时变信道冲激响应模型 h ( t , τ ) h(t,\tau) h(t,τ)的一个特例,是从时变变成了时不变。

适用于固定或慢变信道。描述多径传播但不考虑时间变化。

∑ i h i δ ( τ − τ i ) \sum_i h_i \delta(\tau - \tau_i) ihiδ(ττi)——时不变多径信道

∑ i h i δ ( τ − τ i ) \sum_i h_i \delta(\tau - \tau_i) ihiδ(ττi)是时不变信道的离散多径模型,增益 h i h_i hi 固定,时延 τ i \tau_i τi 离散

输出信号:
y ( t ) = ∑ i h i x ( t − τ i ) y(t) = \sum_i h_i x(t - \tau_i) y(t)=ihix(tτi)

∑ i h i δ ( τ − τ i ) \sum_i h_i \delta(\tau - \tau_i) ihiδ(ττi) h ( τ ) h(\tau) h(τ)的一种特例,它表示信道的多径是有限且固定的,每条路径的增益固定。

h ( n ) h(n) h(n)——离散时间信道冲激响应

h ( n ) h(n) h(n) h ( τ ) h(\tau) h(τ)的离散化版本,适用于数字通信系统,依赖离散时间 n n n
输出信号:
y ( n ) = ∑ k x ( n − k ) h ( k ) y(n) = \sum_k x(n - k) h(k) y(n)=kx(nk)h(k)
由于离散化的原因, h ( n ) h(n) h(n)不能表示连续时延 τ \tau τ的多径,受到采样间隔 T s T_s Ts的影响,它能表示的多径时延只能是采样间隔 T s T_s Ts的整数倍,无法表示小数倍。

H \mathbf{H} H——离散时间卷积矩阵

H \mathbf{H} H 是基于 h ( n ) h(n) h(n) 构建的Töplitz矩阵,表示离散时间信道的卷积操作,适用于线性时不变系统。

输出信号:
y = H x + n \bm{y} = \mathbf{H} \bm{x} + \bm{n} y=Hx+n
[ y ( 1 ) y ( 2 ) ⋅ ⋅ ⋅ y ( N ) ] = [ h ( 1 ) 0 ⋅ ⋅ ⋅ 0 h ( 2 ) h ( 1 ) ⋅ ⋅ ⋅ 0 ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ h ( N ) h ( N − 1 ) ⋅ ⋅ ⋅ h ( 1 ) ] [ x ( 1 ) x ( 2 ) ⋅ ⋅ ⋅ x ( N ) ] \begin{bmatrix} y(1) \\ y(2) \\ ···\\ y(N) \end{bmatrix} = \begin{bmatrix} h(1) & 0 & ··· & 0 \\ h(2) & h(1) & ··· & 0 \\ ··· & ···&···&··· \\ h(N) & h(N-1) & ··· & h(1) \end{bmatrix} \begin{bmatrix} x(1) \\ x(2) \\ ···\\ x(N) \end{bmatrix} y(1)y(2)⋅⋅⋅y(N) = h(1)h(2)⋅⋅⋅h(N)0h(1)⋅⋅⋅h(N1)⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅00⋅⋅⋅h(1) x(1)x(2)⋅⋅⋅x(N)

这么一看是不是很直观,把复杂的卷积运算换成矩阵乘法,方便处理。不过我这个写法还不对,这么写漏掉了前一个信号的拖尾效应。因为卷积是会拉长信号长度的,这么写成矩阵的形式就有偏差,严谨的写法如下:
在这里插入图片描述
在这里插入图片描述
其中, H a \mathbf{H}_a Ha为 信道矩阵,表示信道对第i个周期信号的影响, H b \mathbf{H}_b Hb也是 信道矩阵,表示信道对第i-1个周期信号 造成的拖尾影响。
则一组接收信号中既包含了当前周期目标信号的成分,也包含了前一个周期的信号 的成分,这是多径衰落信道中不可避免的多径干扰。

有点复杂,不过好处是做FFT的时候,信道矩阵可以对角化,此处不再赘述。

总结与比较

总结

在这里插入图片描述

比较

序号信道模型适用性模型说明局限性
1 h ( t , τ ) h(t,τ) h(t,τ)适合高速移动通信场景,如5G、UWB、复杂时变信道描述信号在时间t和时延τ上的响应,适用于时变信道模型复杂,计算成本高
2 ∑ i h i ( t ) δ ( τ − τ i ) \sum_i h_i(t) \delta(\tau - \tau_i) ihi(t)δ(ττi)适合高速移动通信场景,如城市移动通信、瑞利衰落描述时变信道中离散多径分量的时变增益假设时延固定,无法描述连续时延分布
3 h ( τ ) h(\tau) h(τ)适合固定或慢变信道,如室内Wi-Fi、有线信道描述时不变信道的时延响应,适用于多径传播无法描述时变特性
4 ∑ i h i δ ( τ − τ i ) \sum_i h_i \delta(\tau - \tau_i) ihiδ(ττi)适合城市环境固定无线通信,如GSM、CDMA描述时不变信道中的离散多径分量仅限于离散多径,无法描述连续时延
5 h ( n ) h(n) h(n)适合数字通信系统,如OFDM、数字基带处理离散时间域的信道冲激响应,基于采样受限于采样精度
6 H \mathbf{H} H (Töplitz矩阵)适合MIMO、信道均衡等矩阵运算场景基于h(n)构建的Töplitz矩阵,表示离散卷积仅适用于LTI系统,应用特定

二、接收信号能否直接做N点FFT

对这样一个离散接收信号模型:
y ( n ) = x ( n ) ∗ h ( n ) = ∑ k = 1 L x ( n − k ) h ( k ) y(n)= x(n) * h(n)= \sum_{k=1}^{L} x(n - k) h(k) y(n)=x(n)h(n)=k=1Lx(nk)h(k)

信号长度N,信道长度L,输出长度N+L-1。即经过信道,信号被拉长了。

根据线性时不变系统的“时域卷积=频域乘积”的理论
(注意,一定要是线性时不变系统,所以信道响应必须是时不变的)
应该有:
Y ( k ) = x ( k ) ⋅ h ( k ) Y(k)= x(k) · h(k) Y(k)=x(k)h(k)

信号 x ( n ) x(n) x(n)与信道响应 h ( n ) h(n) h(n)的卷积*是线性卷积,因此输出 y ( n ) y(n) y(n) 的长度会变长至N+L-1
理论上,接收信号要做N+L-1点的FFT才能保证频谱不失真,那什么情况下能做N点FFT,强行做N点FFT又有什么影响?

周期信号

做N点FFT就是默认将线性卷积变为N点的循环卷积,这就涉及到一个重要知识

循环卷积假设信号在长度 M M M 的窗口内是周期性的。
我们知道,如果 M ≥ N + L − 1 M≥N+L-1 MN+L1,循环卷积输出与线性卷积一致,无失真。

但是如果 M < N + L − 1 M<N+L-1 MN+L1呢?此时需要分情况讨论

如果 x ( n ) x(n) x(n)是周期信号, M = L M=L M=L点循环卷积在一个周期内等价于线性卷积,无失真。(OFDM通过添加CP实现伪周期性,从而应用循环卷积)
如果 x ( n ) x(n) x(n)是非周期信号, M < N + L − 1 M<N+L-1 MN+L1点循环卷积就会导致时间域混叠,出现失真

因此周期信号,以及添加了循环前缀CP的OFDM信号满足循环卷积的条件,可以直接做N点的FFT,并且频谱不失真。

非周期信号

经过刚刚的分析,我们知道非周期信号不满足循环卷积的条件,必须对N+L-1点的线性卷积进行FFT才能保证不失真。

如果强行做N点FFT,就需要从接收信号中丢弃L-1个样本,这会导致频谱出现怎样的失真?

假设取 y lin [ n ] y_{\text{lin}}[n] ylin[n] 的前 L L L 点,定义截断信号:
y trunc [ n ] = { y lin [ n ] , 0 ≤ n ≤ N − 1 0 , otherwise y_{\text{trunc}}[n] = \begin{cases} y_{\text{lin}}[n], & 0 \leq n \leq N-1 \\ 0, & \text{otherwise} \end{cases} ytrunc[n]={ylin[n],0,0nN1otherwise

这等价于对 y lin [ n ] y_{\text{lin}}[n] ylin[n] 施加矩形窗 w [ n ] w[n] w[n](长度 N N N):
y trunc [ n ] = y lin [ n ] ⋅ w [ n ] , w [ n ] = { 1 , 0 ≤ n ≤ N − 1 0 , otherwise y_{\text{trunc}}[n] = y_{\text{lin}}[n] \cdot w[n], \quad w[n] = \begin{cases} 1, & 0 \leq n \leq N-1 \\ 0, & \text{otherwise} \end{cases} ytrunc[n]=ylin[n]w[n],w[n]={1,0,0nN1otherwise

在频域中,窗函数 w [ n ] w[n] w[n] 的 DTFT W ( f ) W(f) W(f) Y lin ( f ) Y_{\text{lin}}(f) Ylin(f) 卷积:
Y trunc ( f ) = Y lin ( f ) ∗ W ( f ) Y_{\text{trunc}}(f) = Y_{\text{lin}}(f) * W(f) Ytrunc(f)=Ylin(f)W(f)
W ( f ) W(f) W(f) 是矩形窗的频谱(类似 sinc 函数),具有主瓣和旁瓣,导致频谱展宽和泄漏。

频谱具体来说会出现如下失真:

  • 频谱展宽:截断相当于施加矩形窗, Y trunc ( f ) Y_{\text{trunc}}(f) Ytrunc(f) 的频谱被 W ( f ) W(f) W(f) 的主瓣展宽,主瓣宽度约为 2 / ( N T s ) 2/(N T_s) 2/(NTs)。这使得 Y [ k ] Y[k] Y[k] N N N 点 FFT 的结果)在每个频率点包含来自邻近频率的贡献,降低频率分辨率。
  • 旁瓣泄漏:矩形窗的旁瓣导致频谱能量泄漏到其他频率分量,尤其在 y lin [ n ] y_{\text{lin}}[n] ylin[n] 的尾部( n = N , … , N + L − 2 n = N,\dots, N + L - 2 n=N,,N+L2)包含显著能量时,泄漏更明显。 幅度失真:由于丢失了尾部 L − 1 L - 1 L1 个样本, Y [ k ] Y[k] Y[k] 的幅度可能偏离真实频谱 Y lin ( f ) Y_{\text{lin}}(f) Ylin(f),特别是在低频和高频区域。
  • 相位失真:截断可能改变 y lin [ n ] y_{\text{lin}}[n] ylin[n] 的时间结构,导致 Y [ k ] Y[k] Y[k] 的相位信息不准确,影响相位敏感的应用(如相干检测)。

频谱的失真从时域上来等效观察,相当于时间域混叠。

如果直接对 x [ n ] x[n] x[n] h [ n ] h[n] h[n] 进行 N N N 点循环卷积,等价于将尾部的L-1个样本折叠到开头。

结论

  1. 如果信道响应h(n)长度为L=1,无论什么信号都可以直接做N点FFT
  2. 如果信道响应h(n)长度L>1,周期信号可以做N点FFT,也可以做N+L-1点FFT,后者频谱分辨率更高。
  3. 如果信道响应h(n)长度L>1,非周期信号可以做N+L-1点FFT保证不失真,也可以通过添加CP 将非周期信号的线性卷积转化为循环卷积,做N点FFT。但是如果直接做N点FFT,就会使得频谱展宽,旁瓣泄露。
  4. 核心差异在于:非周期信号的线性卷积扩展了输出长度,而周期信号的周期性和 OFDM 的 CP 机制限制了输出长度为 L L L,从而简化 FFT 点数需求。
【评估目标跟踪方法】9个高度敏捷目标在编队中的轨迹和测量研究(Matlab代码实现)内容概要:本文围绕“评估目标跟踪方法”,重研究9个高度敏捷目标在编队飞行中的轨迹生成测量过程,并提供完整的Matlab代码实现。文中详细模拟了目标的动态行为、运动约束及编队结构,通过仿真获取目标的状态信息观测数据,用于验证和比较不同目标跟踪算法的性能。研究内容涵盖轨迹建模、噪声处理、传感器测量模拟以及数据可视化等关键技术环节,旨在为雷达、无人机编队、自动驾驶等领域的目标跟踪系统提供可复现的测试基准。; 适合人群:具备一定Matlab编程基础,从事控制工程、自动化、航空航天、智能交通或人工智能等相关领域的研究生、科研人员及工程技术人员。; 使用场景及目标:①用于目标跟踪算法(如卡尔曼滤波、粒子滤波、GM-CPHD等)的性能评估对比实验;②作为无人机编队、空中交通监控等应用场景下的轨迹仿真传感器数据分析的教学研究平台;③支持对高度机动目标在复杂编队下的可观测性跟踪精度进行深入分析。; 阅读建议:建议读者结合提供的Matlab代码进行实践操作,重关注轨迹生成逻辑测量模型构建部分,可通过修改目标数量、运动参数或噪声水平来拓展实验场景,进一步提升对目标跟踪系统设计评估的理解。
本软件实现了一种基于时域有限差分法结合时间反转算法的微波成像技术,旨在应用于乳腺癌的早期筛查。其核心流程分为三个主要步骤:数据采集、信号处理三维可视化。 首先,用户需分别执行“WithTumor.m”“WithoutTumor.m”两个脚本。这两个程序将在模拟生成的三维生物组织环境中进行电磁仿真,分别采集包含肿瘤模型不包含肿瘤模型的场景下的原始场数据。所获取的数据将自动存储为“withtumor.mat”“withouttumor.mat”两个数据文件。 随后,运行主算法脚本“TR.m”。该程序将加载上述两组数据,并实施时间反转算法。算法的具体过程是:提取两组仿真信号之间的差异成分,通过一组专门设计的数字滤波器对差异信号进行增强净化处理,随后在数值模拟的同一组织环境中进行时间反向的电磁波传播计算。 在算法迭代计算过程中,系统会按预设的周期(每n次迭代)自动生成并显示三维模拟空间内特定二维切面的电场强度分布图。通过对比观察这些动态更新的二维场分布图像,用户有望直观地识别出由肿瘤组织引起的异常电磁散射特征,从而实现病灶的视觉定位。 关于软件的具体配置要求、参数设置方法以及更深入的技术细节,请参阅软件包内附的说明文档。 资源来源于网络分享,仅用于学习交流使用,请勿用于商业,如有侵权请联系我删除!
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值