matlab 实现线性调频信号以及分析处理
里面有关于实现matlab的算法以及分析处理山国科技记文在线分布的时频平面作直线积分投影的变换,统称对信号作变换在分布的时频平面里惯用轴的截距和斜率为参数表小直线。因此,当需要沿作直线积分时,可将积分路径(直线)的参数(u,a)替换成()日两对参数之间的关系为:m=-cot,w=! sina。若求信号的变换,并以参数表示积分路径,则有:D.a=PQ线w, (t, wB u-u du∫r(,n)ma(w-mn-m)nh∫m(,w[一(m+motcw lt, wo +mt dt/sinaWo=u/sina上式表明,若是参数为和的信号,则积分值最大;而当参数偏离与或时,积分值迅速减小,即对‘定的信号,其变换会在对应的参数处呈现尖峰。我们自然会想到:多分量的信号的特性在平面里更加突出。即表现为各个尖峰,因而更有利于区别交叉项和噪声。利用变换一定能够获得更好的性能。作为时频分析方法之一,分数阶傅里叶变换ˉ与分布()变换()分别有着一定的数学关系,借助它们的联系,可进一步说明分数阶傅里叶变换的物理意义。信号的分布函数的定义为t+=xtde作为能量型时频表示满足许多期望的数学性质,这里给出其边缘特性X tt wdvXw=wtwat对WD旋转C角度,即对分布实施变换,其结果是RWIW=∫f山国技记文在线而信号的阶分薮阶傅里叶变换X。t的就是将信号的旋转c角度,即对于分数阶傅里叶变换只有旋转不变性,所以有X u= wtP可以看出,对时间轴与频率轴的积分分别是信号在时刻的瞬时功率和信号在频率的谱密度,而信号的对与时间成c角度的轴的积分投影对应着角度为a的分数阶傅里叶变换的幅度平方,这进步从能量的角度说明分数阶傅里叶变换作为广义傅里叶变换的含义。正弦信号在时频平面是一条平行于时间轴的直线,即它的频率不随时间变化,可视为旋转角度为°的完全时间域表示;冲击朕数在时频平面是一条平行于频率轴的直线可视为旋转角度为°的完全频率域表示;信号在时频平面是一条斜率为调频率的直线,当该信号的某一角度的分数阶傅里叶变换与其调频率一致时,在无限长度的理想情况下,表现为幅度为无穷大的冲击,在信号长度有限的情况下,其分数阶傅里卟变换呈现极大值这就是信号在分数阶傅里叶变换域的特点。离散 Chirp fourier变换是最近提出的一种有效的线性调频信号检测技术,它 Fourier变换的一种推广形式,可同时匹配 chirp信号的中心频率和调频率。本文利用修正离散Chirp- Fourie交换( MDCFT)实现干扰信号的检测和参数估计,从而实现对干扰的自适应抑制。分析和仿真表明,该方法可对FM干扰有着极好的抑制效果;同时,由于 Chirp- Fourie变换是维的线性变换,可借助快速傅里叶变换(FFT〕实现,与基于WVD的算法相比,不仅避免了交叉项十扰,而且降低了计算的复杂度,其实现更为简使3.基于Mat1ab的上机仿真过程及结果分析3.1对单分量信号的仿真及结果分析():输入解析信号为x()=eb的分布:40,图单分量信号的分布山国科技论文在线在上述解析信号中加入噪声后,用分布分析其性能图加入噪声的单分量信号的分布由图可以看出实际结果与前面的理论推导致。在实际应用中,信号长度总是有限长的,此时分布呈背鳍状。由图可以得到变换对噪声不太敏感,时频变换后信噪比较高。但当干扰的幅度大到一定程度时,变换的结果会严重变差,甚至分析不出结果。():前两个图是输入解析信号为x(t)=em的变换,后两个图是在这个解析信号中加入噪声以后用变换对其进行的分析:400C501m01501020100150图单分量信号的变换由理论分析可知,当旋转角度与线性调频信号的斜率相這应时,变换将出现一个峰值。这个分析在图中得到了证实。():图前两个图是输入解析信号为x()=e的分数阶傅里叶变换,后两个图是在山国科技论文在线这个解析信号中加入噪声以后用分数阶傅甲叶变换对其进行的分析:分数阶傅甲叶变换变换与变换的紧密联系在图和图的仿真中也可以得到证实HOD50图单分量信号的分数阶傅里叶变换():图的前两个图是输入中心频率是,调频率是的单分量线性调频信号后的Chirp- Fourier变换,后两个图是在这个信号中加入噪声以后用 Chirp-Fourier变换对其进行的分析。通过这个仿真,还将证明一个重要性质: Chirp- Fourier变换可同时匹配线性调频信号的中心频率和调频率的82a图单分量信号的 Chirp fourier变换比较结论:从以上几个仿真图形可以看出,对单分量的信号而言,上述几个变换山国科技论文在线都有非常好的时频聚集性,特别是分布与理论结果完仝一致。在抗噪声方面,对比几个图可知,变换和 Chirp- Fourier变换要比分布和分数阶傅里叶变换吏好。而对于分数阶傅里叶变换和分布,分数阶傅里叶变换的抗噪声性能要好3.2对多分量信号的仿真及结果分析个多分量的线性调频信号的D15020心Dm图多分量信号的一个多分量的线性调频信号的变换50.540多分量信号的变换山国科技论文在线个多分量的线性调频信号的分数阶傅甲叶变换:图多分量信号的分数阶傅里叶变换个多分量的线性调频信号(含两个分量,中心频率和调频率分别为k=)的 Chirp- Fourier变换50299,Q图多分量信号的 Chirp-fourier变换比较结论:从以上四个图可以看出,对于多分量信号,分布由于存在交叉项,时频面模糊不清,而其他三种变换则可以检测到两个信号。从图中还可以看到,Chirp- Fourier变换的效果是最好的。而且我们从图中还可以清楚地看到线性调频信号的中心频率和调频率。4LFM信号的应用线性词频)信号广泛地应用于雷达、声纳和通信等信息系统中。在这类系统中,信号的检测与参数估计是个重要的研究课题,受到特别的关注。下面给出一个基于FRT的MTD雷达信号处理过程的防真实例。假设有一个运动目标,回波信号为Stjn∫t-jwt+nt,其中nt为杂波信号,信号参数为nt是均值为零,方差为的高斯白噪声,信噪比为,观测时间为,采样频率为采样点数为N采用分数阶域的扫描上算法对该冋波信号作计算机仿真,仿真结果如图所从图中可以清楚看到一个LFM信号的存在,而闬目标的峰值非常突出,受杂波的影响相对较小。因此采用FRT的MTD雷达的抗干扰能力较强。另外由于日标的特征非常明显,可以通过适当提高杂波门限的方法来减小虚警概率山国科技论文在线图基于ⅣRFT的MTD雷达信号处理过程的防真5结束语非平稳信号是现代信号处理的主要研究对象之一,对其有很多种理论分析方法。本文介绍的分布,变换,分数阶傅里叶变换,变换是其中比较常用和重要的几种。本文对这几种变换做了初步的介绍,进而对它们进行了一些比较这有助于进一步了解各种变换的性能和作信号分析时选择合适的变换。时频分布之所以受到很多研究人员和信号处理领域的工程人员的重视,是因为它有很多传统傅立叶变换所不具备的性质。由时频分析的定义可知时频表示能给出信号在时域和频域的信息。经过儿年的发展,时频分析理论趋于成熟,并遂渐在实际应用中崭露头角,近年来已在实际的非平稳信号处理中获得了十分广泛的应用。如:信号检测与分类,吋频域滤波,信号综合,系统辩识和谱估计等。在的期刊和国际会议上发表的与采用时频工具处理非平稳干扰有关的论文及研究报告共有余篇,其中以美国大学教授的成果最为显著。时频分析是一个前景很广阔的研究方向,虽然取得了一定的成就,但理论体系尚不十分完备,需要进一步的发展。参考文献[1ˉ张贤达,保铮《非平稳信号分析与处理》[M1998年9月第1版国防工业出版社[2ˉ沈民奋,孙丽莎《现代随机信号与系统分析》M年月第版科学出版社[3丁凤芹,曹家麟《基丁分数阶傅里叶变换的多分量 Chirp信号的检测与参数估计》《语音技术》2004年第1期[4_孙泓波,郭欣,顾红,苏上民,刘国岁《修正 Chirp- Flourier变换及其在SAR运动目标检测中的应用》《电子学报》2003年第1期山国技记文在线[5董永强,陶然,思永,王越《基丁分数阶傅里叶变换的SAR运动目标检测与成像》《兵工学报》1999年第2期L6_陶然,齐林,王越《分数阶 Fourier变奂的原理与应用》LM」2004年8月第1版清华大学出版社[7董永强,陶然,周思永,王越《含未知参数的多分量 chirp信号的分数阶傅里叶分析》《北京理工大学学报》1999年第5期[8ˉ陈辉,王永良《利用离散 Chirp- Flourier变换技术估计调频信号参数》《空军雷达学院学报》2001年第1期[9ˉ齐林,穆晓敏,朱春华《系统中基于 Chirp- Fourier变换的扫频干扰抑制算》《电讯技术》年第期[10]李勇,徐震等《 MATLAB辅助现代工程数字信号处理》[M2002年10月鷥1版西安电子科技人学出版社「111胡昌华,周淘,夏启兵,张伟《基于 MATLAB的系统分析与设计—时频分析》「M12001年7月第1[2]干小宁,许家栋《离散调频-傅里叶变换及其作雷达成像中的应用》《系统工稈与电子技术》2002年第3期
- 2020-12-02下载
- 积分：1
RS纠错编码原理及其实现方法.pdf
RS纠错编码原理及其实现方法。Zhengzhou Oriole Xinda Electronic Information Cc., Ltd前言随着越来越多的系统采用数字技术来实现,纠错编码技术也得到了越来越广泛的应用。RS码既可以纠正随机错误,又可以纠正突发错误,具有很强的纠错能力,在通信系统中应用广泛。近些年来,随着软件无线电技术的发展,RS编码、译码一般都在通用的硬件平台上实现。通常采用基于FPGA的ⅦHDL编码硬件实现,或者在DSP、单片机上用C和汇编编程软件实现。RS纠错编码涉及的领域很广,特别是设计到很多数学知识。这对那些对数学不太感冒的工程技术人员来书是个不小的挑战。尽管讲RS编码的书籍很多但是那些书都是采用循序渐进,逐步引人的方式从汉明码到循环码,从循环码到BCH码,BCH码再引入悶S码。对亍工程技术人员他们需要的是简明扼要的讲解,和详细的实现方法。本人写这篇文章的宗旨就是尽量最简单的语言最简短的篇幅来讲RS纠错编码原理,把重点来放在实现方法上。为了便于读者仿真,本文采样MLAB程序实现,程序尽量符合硬件C语言写法,读者经过简单修改即可应用到工程中去。本文读者对象本文是为那些初识瑙编码的学生、工程技术人员而写,并不适合做理论研究,如果你是纠错编码方面的学者、专家,那么本文并不适合你。由于作者水平有限,错误在所难免,恳请读者批评指正。不得更改陈文礼2008-01于郑州Zhengzhou Oriole Xinda Electronic Information Cc., Ltd必备的一些代数知识1、在纠错编码代数中,把以二进制数字表示的一个数据系列看成一个多项式。例如二进制数字序列1010111,可以表示成:M(x)=ax+a5x0+a5不5+a+4 TasK +ax+a,x+ank式中的x表示代码的位置,或某个二进制数位的位置,X前面的系数表示码的值。若a;是一位二进制代码,则取值是0或1。dM()称为信息代码多项式多项式次数称系数不为0的x的最高次数为多项式/(x)的次数,记为Of(x)2、域域在R编码理论中起着至关重要的作用。简单点说域GF(2)有2设2个符号[0,n,a2…22且具有以下性质域中的每个元素都可以用a",a,a2,om的和来表示。a←la为本原多项式p(x)的根。运算规则有:在纠错编码运算过程中,加减、乘和除的运算是在伽罗华域中进行。现以GF(2)域中运算为例:加法例:a+a=0010+0110101(模2加法相当于0005与011或减法运算与加法相同乘法例:a·a0=a(8+10)modl5除法例:cs/a0=a-2=a-2+5=a不理解没关系,下面的例子也许对你有帮助。例:mF=4,p(x)=x4+x+1求GF(2")的所有元素因为a为p(x)的根得到a4+a+1=0或a4=a+1(根据运算规则)Zhengzhou Oriole Xinda Electronic Information Cc., Ltd由此可以得到域的所有元素元素二进制对应十进制对应码值000000101000a+100l⊥0110a(a+1)=a+a(mod p(a))12a(a+a=a+a(mod p(a)1011a(a+l(modula))+a+1)10C(a+1=a+a(mod p(a )a(a23+a)a+I(mod p(a)1110a(a+a+D=aa+a(modp(a)tatI(mod p(a))11a(a3+a2+a+1)=a34a2+1(modp(a)1001a(a+a+1=a+l(mod p(a)a(a+1=l(mod(a))由此可以看岀本原多项式是求解域的全部元素的关键。读者也许会有这样的疑问我们如何得到p(x)呢?本原多城式p(x)的特性是2+得到的余式等于0O(X由于作者也是工程技术人员,具体怎么得到p(x),也没有深究过。Zhengzhou Oriole Xinda Electronic Information Cc., Ltd作者在设计RS编码时候都是根据 MATLAB指令rsgeηpoly来得到p(x)。其格式为 rsgenpoly(n,k)参数n为码长一般n=2"-1,k为信息码元个数。例如m4,码长n=15,信息码元长度为9GF(2)的本原多项式可以根据指令>>rsgenpoly(15, 9)得到ans= GF(2 4)array. Primitive polynomial =D 4+D+1 (19 decimal)有读者来信问:我要做一个(158的RS编码,在 MATLAB中输入命令 rsgenpoly(158,128),结果MAB报错Error using =- rsgenpolyN must equal 2m-1 for some integer m这里做一下解释我们S编码时普先要根据码长选取mλ选择原则是2若码长为6那么我们可以选择n=8, rsgenpey命令的第少个参数必须为2"-1,第二个参数司以随便选择只要小于2”-1就形了在此给出m∈(2,16)的所有本原多项式(m=2)P[m+1]={1,1,1}/米1+x+x3*/P[m+1]-{1,1,0,1}/米1+x+x4*/P[m11]={1,1,0,0,1}/米1+x2+x5*/P|m+1={1,0,1,0,0,1};Zhengzhou Oriole Xinda Electronic Information Cc., Ltd(m=6)/米1+x+x6*/P[m+1]={1,1,0,0,0,0,1}7)/来1+x3+x7*P[m+1]={1,0,0,1,0,0,0,1}(m=8)/米14x2+x31x4+x8*/P[m+1]-{1,0,1,1,1,0,0,0,1/*1+x4+x9半P[m1]={1,0,0,0,1,0,0,0,(m=10)/1+x3+x10*/P|m+1={1,0,0,1,0,0,0,0,/*1+x2+x11P[m+1]={1,0,0,0,0,0,0,1}(m=12)/*1+x+x4+x6+x12P[m+1]-{1,1,0,0,、1,0,0,(m=13)/*1+x+x^3+x4+x^13*/P[m+1]={1,1,0,1,1,0,0,00,0,1};(m=14)/*1+x+x6+x10+x14来P[m+1]={1,1,0,0,0,0,1,0,0,0,1,0,0,0,1}(m=15)/米14x+x15*/P[m+1]={1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,1};(m=16)/*1+x+x3+x12+x16*/P[m+1]={1,1,0,1,0,0,0,0,0,0,0,0,1,0,0,0,1};Zhengzhou Oriole Xinda Electronic Information Cc., Ltd二、线性分组码的一些基本概念1、线性分组码一般用(n,)或(n,k,d)表示n为码长,k为信息码元的数目,n-k为监督码元的数目。d表示码元距离。定义:两个码组上对应位置上数字不同的个数称为码组的距离。发送的码字C=(1,C2C3,…C接收的矢量r=(,2,信道错误图样:e=c+r例如c=(1,1,0,0,0)(1,0,001)e=(1+1,1+0,0+0,0+0,0+1)(0,1,0,0,1)从而可以看出从左端起第2位和第5位是错误的2、校验矩阵概念码长为n,信息数为k,监督数为r。这样的一组码形式为:m:m2,P,P2Pm表示第个信息码,P表示第j个校验码各个校验码可从下列线性方程组求得hm+h2m2+…+n+1B1+012+0h2m1+2m2+…+h2m+0p1p20hmn+h,2m2+…+hm+O+0+…+1p,=0式中h;是常数校验方程组可写成校验矩阵100h21h2…,h2k010h000该矩阵具有r行和n列故式(1-1)可以写成c=0或c=08Zhengzhou Oriole Xinda Electronic Information Cc., LtdH矩阵称为[n,k,r码的校验矩阵。发送矢量为C接收矢量为F若rH≠0则说明接收到的码有错误。设错误图样为e则可写成以下关系式r=c+e为了纠错必须知道那些位上存在错误。这可由校正子(又称伴随式)s来确定s=rH=cH +eh=eh译码器的主要任务就是如何从中得到最像e的错误图样e从而译出c=r-e设第讠个是错误的因此e=(00..0第个有错误s=rH=(00…0、100000)00计算出的矢量示出i是出错误的位置。3、生成矩阵概念生成矩阵G,它是一个k行,n列的矩阵若已知信息组m,通过生存矩阵可求得相应的码字。c=mxG(m是k个信息元组成的信息组)这个应该比较容易理解,在此就不做过多解释。、RS码的一些重要性质1、RS码生成多项式:码长n=2”-1,监督元数目r=n-k=2t,能纠正t个错误。Zhengzhou Oriole Xinda Electronic Information Cc., Ltd定义:在(n,k,d)的RS码中,存在唯一的n-k次多项式g(x),使得每一个码多项式c(x)都是g(x)的倍式。g(x)称为n,k,d]RS码的生成多项式一般情况下g(x)=(x-a)(x-a2)…(x-a2)2、定理:在GF(2m)中,每个非0元素(1,a,a2…a22)均满足x2=1,反之x21-1=0的根必在GF(2")中。所以x-1=(x-a)(x-a)x3、RS码的校验多项式由于生成多项式g(x)是x-1的因式g(rh(g(x)为n-k次多项式,则h(x)为k次多项式,k3x+g)hx+…+x+4)由右式可以看出x"1,x2,x的系数均等于0即gg0010h1+g1bo=0g0h+g1h11+…+8nkh2(2k)=0∴.+n-kk-10n-kk式中g0+81h1+…+8nkh1(n=k)(表示X的系数10
- 2020-12-08下载
- 积分：1
