雷达系统中杂波信号的建模与仿真.doc
《雷达系统中杂波信号的建模与仿真.doc》由会员分享,可在线阅读,更多相关《雷达系统中杂波信号的建模与仿真.doc(51页珍藏版)》请在咨信网上搜索。
1、(完整版)雷达系统中杂波信号的建模与仿真1. 雷达系统中杂波信号的建模与仿真目的雷达的基本工作原理是利用目标对雷达波的散射特性探测和识别目标。然而目标存在于周围的自然环境中,环境对雷达电磁波也会产生散射,从而对目标信号的检测产生干扰,这些干扰就称为雷达杂波。对雷达杂波的研究并通过相应的信号处理技术可以最大限度的压制杂波干扰,发挥雷达的工作性能.雷达研制阶段的外场测试不仅耗费大量的人力、物力和财力,而且容易受大气状况影响,延长了研制周期。随着现代数字电子技术和仿真技术的发展,计算机仿真技术被广泛应用于包括雷达系统设计在内的科研生产的各个领域,在一定程度上可以替代外场测试,降低雷达研制的成本和周期
2、。长期以来,由于对杂波建模与仿真的应用己发展了多种杂波类型和多种建模与仿真方法。然而却缺少一个集合了各种典型杂波产生的成熟的软件包,雷达系统的研究人员在需要用到某一种杂波时,不得不亲自动手,从建立模型到计算机仿真,重复劳动,造成了大量的时间和人力的浪费.因此,建立一个雷达杂波库,就可以使得科研人员在用到杂波时无需重新编制程序,而直接从库中调用杂波生成模块,用来产生杂波数据或是用来构成雷达系统仿真模型,在节省时间和提高仿真效率上的效益是十分可观的。从七十年代至今已经公布了很多杂波模型,其中有几类是公认的比较合适的模型。而且,杂波建模与仿真技术的发展己有三十多年的历史,己经有了一些比较成熟的理论和
3、行之有效的方法,这就使得建立雷达杂波库具有可行性。为了能够反映雷达信号处理机的真实性能,同时为改进信号处理方案提供理论依据,雷达杂波仿真模块输出的杂波模拟信号应该能够逼真的反映对象环境的散射环境。模拟杂波的一些重要散射特性影响着雷达对目标的检测和踉踪性能,比如模拟杂波的功率谱特性与雷达的动目标显示滤波器性能有关;模拟杂波的幅度起伏特性与雷达的恒虚警率检测处理性能有关。因此,杂波模拟方案的设计是雷达仿真设计中极其重要的内容,杂波模型的精确性、通用性和灵活性是衡量杂波产生模块的重要指标。2.Simulink简介Simulink是MATLAB最重要的组件之一,它提供一个动态系统建模、仿真和综合分析的
4、集成环境。在该环境中,无需大量书写程序,而只需要通过简单直观的鼠标操作,就可构造出复杂的系统.Simulink具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点,并基于以上优点Simulink已被广泛应用于控制理论和数字信号处理的复杂仿真和设计.同时有大量的第三方软件和硬件可应用于或被要求应用于Simulink。 Simulink是MATLAB中的一种可视化仿真工具, 是一种基于MATLAB的框图设计环境,是实现动态系统建模、仿真和分析的一个软件包,被广泛应用于线性系统、非线性系统、数字控制及数字信号处理的建模和仿真中3。Simulink可以用连续采样时间、离散采样时间或两种
5、混合的采样时间进行建模,它也支持多速率系统,也就是系统中的不同部分具有不同的采样速率.为了创建动态系统模型,Simulink提供了一个建立模型方块图的图形用户接口(GUI) ,这个创建过程只需单击和拖动鼠标操作就能完成,它提供了一种更快捷、直接明了的方式,而且用户可以立即看到系统的仿真结果。S函数是系统函数(System Function)的简称,是指采用非图形化的方式(即计算机语言,区别于Simulink的系统模块)描述的一个功能块。用户可以采用MATLAB代码,C,C+,FORTRAM或Ada等语言编写S-函数.S-函数提供了函安息代码与Simulink之间的接口,使得用户编写的代码既能像
6、Simulink模型库中的模块那样具有统一的仿真接口,同时能够实现各种灵活的控制和计算功能.S函数是一种特定的语法构成,用来描述并实现连续系统、离散系统以及复合系统等动态系统;S函数能够接收来自Simulink求解器的相关信息,并对求解器发出的命令做出适当的响应,这种交互作用非常类似于Simulink系统模块与求解器的交互作用。一个结构体系完整的S-函数包含了描述动态系统所需的全部能力,所有其它的使用情况都是这个结构体系的特例。往往S函数模块是整个Simulink动态系统的核心。3. 主要内容1首先综述常规脉冲雷达杂波产生的机理,具体阐述了在指定的杂波功率谱下,幅度服从Rayleigh分布、L
7、ogNormal分布、Weibull分布和K分布的杂波建模与仿真的方法。指出杂波散射现象可以理解为与地、海面随机形态相关的一种随机过程,因此通常用杂波幅度分布统计模型和杂波相关模型来描述。(1)Rayleigh分布适用于描述气象杂波、箔条干扰、低分辨率雷达的地杂波。当在一个杂波单元内含有大量相互独立的、没有明显贡献的散射源时,雷达杂波包络服从Rayleigh分布。(2)LogNormal分布使用于低入射角,复杂地形的杂波数据或者平坦区高分辨率的海杂波数据。(3) Weibull分布的动态范围介于上述两种分布之间,能在更广阔的范围内精确表示实际的杂波分布.通常,在高分辨率雷达,低入射角的情况下一
8、般海情的海浪杂波能够用Weibull分布精确地描述,地杂波也能用Weibull分布描述。(4)K分布适用于描述高分辨率雷达的非均匀杂波,多见于对海杂波、地杂波的描述。K分布是一种符合分布模型,它可由一个均值是慢变化的Rayleigh分布来表示,其中这个慢变化的均值服从分布。除了上述具有特定概率密度函数的非相关雷达杂波仿真外,在雷达信号处理的有些场合还需要知道雷达杂波功率谱分布,常见的有Gaussian型、Cauchy型、AllPole型等。2针对以上论述的十二种杂波分布给出了在MATLAB平台下用m语言实现的算法实现流程。3将这十二种杂波模型生成Simulink自定义模块,并添加到Simuli
9、nk仿真模型库中供仿真调用。4.杂波建模与模拟方法 杂波可以说是雷达在所处环境中接收到的不感兴趣的回波。就像目标回波一样,杂波也是极为复杂的。为了有效地克服杂波对信号检测的影响,需要知道杂波的幅度特性以及频谱特性。除独立的建筑物、岩石等可以认为是固定目标外,大多数地物、海浪杂波都是极为复杂的,它可能既包含有固定的部分又包含有运动的部分,而每一部分反射回来的回波,其振幅和相位都是随机的.通常采用一些比较接近而又合理的数学模型来表示杂波幅度的概率分布特性,这就是雷达杂波模型。目前描述杂波模型主要有三种方式:(1)描述杂波散射单元机理的机理模型;(2)描述杂波后向散射系数的概率密度函数的分布模型;(
10、3)描述由实验数据拟合与频率、极化、俯角、环境参数等物理量的依赖关系的关系模型。4。1 雷达杂波幅度分布模型到目前为止,人们已经提出了许多杂波模型,有关描述杂波后向散射系数的概率密度函数的分布模型,比较公认的幅度概率密度函数分布模型为Rayleigh分布、LogNormal、Weibull分布、K分布等。(1) Rayleigh(瑞利)分布在雷达可分辨范围内,当散射体的数目很多的时候,根据散射体反射信号振幅和相位的随机特性,它们合成的回波包络振幅是服从瑞利分布的6。以x表示杂波回波的包络振幅,以2表示它的功率,则x的概率密度函数为: (41)相对应的概率密度函数分布曲线如图4.1所示.图4.1
11、 瑞利分布概率密度函数分布曲线图瑞利分布与每个散射体的振幅分布无关,只要求散射体的数目足够多,并且所有散射体中没有一个起主导作用。需要说明的是,瑞利分布只能代表同一个距离单元上杂波从这次扫描到下次扫描的变化规律,它不能用来表示同一个扫描过程中杂波回波的振幅分布,因为杂波的强度一般都是随着距离的增大而减弱的。对于低分辨力雷达,当高仰角和平稳环境时,瑞利分布的杂波模型可以得到较为精确的结果。但是,随着对雷达杂波分布特性分析的逐步深入,人们发现,对于海浪杂波和地物杂波,瑞利分布模型并不能给出令人满意的结果。特别是随着距离分辨力的提高,杂波分布出现了比瑞利分布更长的“尾巴,即出现高振幅的概率相当大。因
12、而,如果继续采用瑞利分布模型,将出现较高的虚警概率.海浪杂波的分布不仅是脉冲宽度的函数,而且也与雷达极化方式、工作频率、天线视角以及海情、风向和风速等因素有关,地物杂波也受类似因素的影响。对于高分辨力雷达,在低仰角或恶劣海情下,海浪杂波己不服从瑞利分布,而通常能用韦布尔分布来描述。类似地,地物杂波通常能用LogNormal分布来描述7。(2) LogNormal(对数一正态)分布设x代表杂波回波的包络分布,则x的LogNormal分布是: (4-2)其中代表lnx的标准差,xw是x的中值。相对应的概率密度函数分布曲线如图4.2所示:图4。2 LogNormal分布概率密度函数分布曲线图LogN
13、ormal分布的严重缺点是在最影响虚警和灵敏度的区域里,吻合程度反而较差。对数一正态分布和瑞利分布之间的主要差别在于前者的“尾巴”较长,也就是说,大幅度的概率要比后者大一些。(3) Weibull(韦布尔)分布一般来说,对于大多数试验和理论所确定的杂波幅度分布,瑞利分布模型和对数一正态分布模型仅适用于它们中的有限分布。瑞利分布模型一般地倾向于低估实际杂波分布的动态范围,而对数一正态分布倾向于高估实际杂波分布的动态范围8.韦布尔杂波分布模型比瑞利分布模型、对数一正态杂波分布模型常常能在更广的环境内精确的表示实际的杂波分布。适当地调整韦布尔分布的参数,能够使它成为瑞利分布或接近于对数一正态分布。通
14、常,在使用高分辨力雷达,低入射角的情况下,海浪杂波能够用韦布尔分布模型精确地描述,地物杂波也能够用韦布尔分布模型描述。设x代表杂波回波的包络振幅,则x的韦布尔分布为: (43)其中:xm为尺度参数,是分布的中值;a为分布的形状(斜度)函数.相对应的概率密度函数分布曲线如图4。3所示:图4。3 Weibull分布概率密度函数分布曲线图如果把式(2-3 )形状参数a固定为2,并把改写成22,则式(2-3 )变为: (44)这就是瑞利分布。所以,瑞利分布是韦布尔分布的特例。如果a=1,并把改写成22,则韦布尔分布变为: (4-5)这就是指数分布。从信号检测的观点来说,对数一正态分布代表着最恶劣的杂波
15、环境:瑞利分布代表最简单的杂波环境;韦布尔分布是中间模型9。在许多情况下,它是一种比较合适的杂波模型,因此,它比瑞利分布能适应更宽的杂波范围。(4)K分布:它的表达式如下: (46)其中x为幅度,a为量化参数,v为形状参数,它的取值范围为一1v=maxN,Nc (412)式(4-12)中,N是所要产生的杂波过程的长度;Nc是杂波过程的相关长度.利用IFFT变换就可以得到相关随机序列x(kt),其过程如图4。5所示。图4.5相关高斯序列产生原理方框图图4。5中, (413)当且仅当相位因子n各分量的概率分布是高斯分布时,x(kt)的各分量才是高斯分布。这样便产生了符合特定功率谱要求的相关高斯分布
16、序列.ZMNL法需要找到高斯序列与所需序列相关系数之间的非线性关系g(.),且它随不同的分布而不同,故不能对协方差矩阵和概率密度函数进行独立控制12。4.3。2 相关对数正态分布杂波的建模利用已产生的相关高斯随机序列,可得到相关对数正态分布的随机序列,分析如下:对于随机变量yN(lnu,2) ,做线性变换x=exp(y),得到具有两个参数的对数正态分布随机变量X的PDF为: (414)式(214)中,u为比例因子,根据具体的环境而不同。由于在时域幅度上的变换会对功率谱产生影响,所以图5。4中的H(f)不再是单纯的高斯谱采样,而是要经过非线性补偿。设Si,j为y(nt)的自相关函数,对应高斯谱P
17、i,j为经过补偿的功率谱序列P(nt)对应的自相关函数,设其对应的随机序列为w(nt),则Si,j与Pi,j有对应关系: (4-15)其中,形状参数决定了非线性变换造成的谐失真程度,Pi,j经过FFT得到经线性补偿的功率谱序列P(nt).令图4。5中的 (416)则可得到非线性补偿功率谱高斯分布的随机序列y(nt),在经非线性变换x=exp(y),可得高斯谱对数正态分布的随机序列x(nt),此过程如图4。6所示:图4。6相关对数止态分布序列产生方框图4。3。3 相关韦布尔分布杂波的建模韦布尔分布的概率密度函数为: (417)式(2-17)中:q为比例因子,p为形状因子.若p=2,就成为瑞利分布
18、,若P=1就是指数概率密度函数。Si,j, Pi,j的关系为: (4-18)式(218)中:(.)为gamma函数,2F1(。)为高斯超几何分布函数。高斯超几何分布函数的定义如下: (4-19)由于无法由上式(4-19)解出Pij的表达式,因此无法象对数一正态分布一样直接计算Pij。实际的应用中可以把两者的关系存入一个表中,通过查表法来计算。此模拟过程如图4.7所示图4.7相关韦布尔分布序列产生方框图4。3。4 相关K分布杂波的建模在相关K分布模型中,杂波回波的幅度被描述为两个因子的乘积,第一部分是斑点分量(即快变化分量),它是由大量散射体的反射进行相参叠加而成的,符合瑞利分布,即第二部分是基
19、本幅度调制分量(即慢变化分量),它反映了与大面积结构有关的散射束在空间变化的平均电平,具有长相关时间,服从x分布13。相关K分布的概率密度函数为: (4-20)其中,x0是杂波幅度,(.)为Gamma函数,Kv(.)为第二类修正的v阶Bessel函数v一1是形状参数,a是标度参数。形状参数v规定了与平均值有关的较大的矩,并通过混合模型规定了在平均值中的不均匀量;对于大多数杂波,形状参数v的取值范围是0。1v。对于较小的取值(即v0。1)杂波有较长的拖尾,并意味着有尖峰杂波;而当v-时,杂波的分布接近于瑞利分布。K分布随机序列的n阶矩为: ,则杂波功率为: 。K分布的混合模型包含了杂波起伏的两个
20、部分,它们具有不同的相关时间。第一部分是基本幅度调制分量(即慢变化分量),它是由与散射体结构有关的散射束在空间变化的平均电平所致。它用取平方根的伽马分布表示,并且具有相关性;该相关性依赖于风速等环境条件,有秒量级长的相关时间,且不受频率捷变的影响。这种基本调制分量可通过对数据进行1/4秒的积分而分离出来。第二部分是斑点分量(即快变化分量),它是由任何距离单元中杂波的多路径散射性质产生的,它是通过散射体内部运动或通过频率捷变而去相关.斑点分量可以从相当短的时间序列的杂波数据中分离出来,这个杂波数据序列时间在单一距离选通波门的清况下约为200毫秒。所以杂波功率z2可以表示成两个独立随机变量的乘积:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 雷达 系统 中杂波 信号 建模 仿真
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。