1.1 韦布尔幅度分布建模

瑞利分布会导致对杂波分布动态范围的过低估计,而韦布尔分布在更宽泛环境下可以提供比瑞利分布对杂波分布更精确描述。韦布尔分布是一种适用于多数环境下的杂波模型。

韦布尔分布的PDF为

式中,p为形状参数,反应函数的形状参数;q为尺度参数,反应函数的尺度参数。

在对杂波的分布参数进行估计时常用矩估计法,则韦布尔分布的n阶矩为

韦布尔分布的方差为

通常使用韦布尔分布的中值rm对参数进行稳健性估计,其中值为

则归一化变量R=r/rm的概率密度函数为

则韦布尔分布的归一化概率密度函数如图所示

由图可知,随着形状参数p的减小,大幅度样本出现的概率增加,韦布尔分布的拖尾变长。

1.2 高斯型多普勒谱建模

通常使用高斯模型对杂波的多普勒谱进行建模,其功率谱密度表达式为

式中,f为目标多普勒频率,f0为中心多普勒频率,Δf为3dB带宽。

经过傅里叶变换后,其自相关函数表达式为

则其自相关系数函数为

1.3 基于ZMNL的杂波仿真

无记忆非线性变换法实现相关韦布尔分布仿真原理框图如图所示

图中,x1x2为0均值、方差为1的标准高斯白噪声列向量,使用三角变换法将高斯白噪声转换为相关高斯分布随机变量y1y2,再经过非线性变换将相关高斯随机变量y1y2转换为相关韦布尔分布随机变量。详细步骤如下:

(1)根据韦布尔分布的功率谱密度函数S(f)计算自协方差系数函数r(t),并对r(t)离散化得到韦布尔分布的相关系数矩阵

(2)根据韦布尔分布与标准高斯分布相关系数的非线性关系计算标准高斯分布的自协方差系数函数λ(t),非线性关系为

式中,Γ()为伽马函数,2F1()为高斯超越几何函数,韦布尔分布相关系数r和高斯分布相关系数λ的关系如图所示

观察上图,韦布尔分布相关系数r和高斯分布相关系数λ的关系在0~1单调,可以通过设置一定精度,使用插值计算对应的高斯分布相关系数λ

通过插值后,计算高斯分布相关系数

(3)对λx进行 Cholesky分解,得到三角阵

由于计算误差,λx 可能不能进行Cholesky分解,可使用nearestSPD函数寻找与λx 最接近的正定对称矩阵。

(4)使用A分别与列向量x1x2进行矩阵乘法,得到所需的相关高斯随机变量y1y2。最后进行非线性变换,得到所需的相关韦布尔分布随机变量r

仿真参数:脉冲数:1024;韦布尔分布形状参数p=1.77;脉冲重复频率PRF=1000Hz;3dB带宽Δf=60Hz。

仿真结果:

仿真韦布尔分布归一化幅度概率密度函数如图所示

杂波信号时域信号如图所示。

杂波归一化功率谱密度如图所示

代码:MATLAB零记忆非线性方法(ZMNL)相关韦布尔分布

MATLAB零记忆非线性方法(ZMNL)相关韦布尔分布相关推荐

  1. matlab 零速检测,一种基于车辆零速检测的惯性导航误差修正方法与流程

    本发明涉及车载导航与定位领域,尤其是涉及一种基于车辆零速检测的惯性导航误差修正方法. 背景技术: 惯性导航系统(inertialnavigationsystem,ins)能根据惯性传感器(陀螺仪.加速 ...

  2. matlab三角函数降次,Matlab实验-传递函数表示方法.ppt

    <Matlab实验-传递函数表示方法.ppt>由会员分享,可在线阅读,更多相关<Matlab实验-传递函数表示方法.ppt(57页珍藏版)>请在装配图网上搜索. 1. Matl ...

  3. 基于Matlab的结点电压法及相关定理验证的人机交互界面!

    基于Matlab的结点电压法及相关定理验证的人机交互界面! 前言 本篇主要针于无储能元件的电路分析 由结点电压法求得相关结点的电压 进而验证戴维宁定理.叠加定理等 提示:以下是本篇文章正文内容,下面案 ...

  4. 冯·诺伊曼奖得主Jorge Nocedal:增强学习中零阶优化方法及其应用(附演讲视频和PPT)

    2020 北京智源大会 本文属于2020北京智源大会嘉宾演讲的整理报道系列.北京智源大会是北京智源人工智能研究院主办的年度国际性人工智能高端学术交流活动,以国际性.权威性.专业性和前瞻性的" ...

  5. 【 MATLAB 】DFT的性质讨论(二)序列的循环移位及其 MATLAB 实现(频域方法)

    上篇博文:[ MATLAB ]DFT的性质讨论(二)序列的循环移位及其 MATLAB 实现(时域方法) 提到了对序列x(n)做循环移位后的DFT形式为: 上篇博文已经讨论过了第一种实现循环移位的方法, ...

  6. matlab中 三种方法计算 Ax b,在MATLAB中,方程Ax=B的解可以用哪个命令求得? matlab 求助 解方程组...

    matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=inv(A)*b - 采用求逆运算解方程组: (2)x=A\B - ...

  7. matlab拟合高次相,matlab 多元高次非线性函数拟合,回归,求教高手!

    问题描述: matlab 多元高次非线性函数拟合,回归,求教高手! 有这样几个数组 自变量n个: x=[...] y=[...] m=[.] w=[.] . 因变量一个: z=[.] 模型1(两个自变 ...

  8. matlab中怎么画五维图,使用MatLab绘制散点图的方法

    使用MatLab绘制散点图的方法: 1.首先从给出的excel表格中采用数据筛选方法,将所有某类型地区的相关数据全部筛选 出来.例如第1类区域的数据如下: 4043 1895 14 1 9.17 24 ...

  9. matlab零状态响应幅度频谱,matlab零状态、零输入响应

    <matlab零状态.零输入响应>由会员分享,可在线阅读,更多相关<matlab零状态.零输入响应(9页珍藏版)>请在人人文库网上搜索. 1.1.已知离散时间系统的差分方程为: ...

最新文章

  1. tinybert华为
  2. 从源码角度来读Handler
  3. 【Android 应用开发】BluetoothClass详解
  4. oracle自动选择索引
  5. python os模块方法_python os模块方法总结
  6. 条件队列java_Java并发系列(4)AbstractQueuedSynchronizer源码分析之条件队列
  7. Python可以这样学(第五季:Socket编程)-董付国-专题视频课程
  8. python编程入门视频-2020年5个经典python编程入门视频教程推荐学习
  9. rinetd 安装、配置方法 通过端口转发来访问内网服务
  10. Deploy a replica set
  11. 什么是表达能力?如何提高表达能力?
  12. 解析机器人视觉系统的神奇之处
  13. html5微课程制作,翟猛老师《微课开发及制作-基于H5课件制作模式》
  14. cad在线转换_超强在线转换工具,打开就能转换
  15. 读书笔记009:《伤寒论》- 足太阳膀胱经
  16. 三维计算机动画的特征是真实性,三维动画技术有哪些优势特征呢?
  17. 高德地图之周边信息查询
  18. 微信公众号推文内可以添加附件了吗?
  19. matlab推导二体问题,二体问题二体问题.ppt
  20. RNA 25. SCI文章中只有生信没有实验该怎么办?

热门文章

  1. 小松鼠是一只很可怜的动物(原创)
  2. c语言与java类型转换_C语言与JAVA理论区别
  3. cortex-A8的第一个裸跑程序——点亮LED灯
  4. 糖尿病常见并发症有哪些?
  5. 设置utools开机自启动
  6. 汇总js中的特殊符号
  7. 2008年网游业收入逾183亿元远超传统三大娱乐业
  8. win10下编译C代码:cmd 命令行 + notepad 记事本
  9. 移远BC26连接华为云
  10. QGIS 用于GIS中的土地退化评估和SDG监测