登录
首页 » Others » 基于MATLAB复调制ZOOM-FFT算法的分析和实现

基于MATLAB复调制ZOOM-FFT算法的分析和实现

于 2020-12-08 发布
0 672
下载积分: 1 下载次数: 7

代码说明:

基于MATLAB复调制ZOOM-FFT算法的分析和实现2006年第4期舰船电子工程121滤波;使用函数来实现傅立叶变换次复数乘法。设数字滤波器的阶数为K,滤波器系数离线生成,则滤波需要DNK次复数乘法,则总4 Matlab仿真和验证的运算量为为验证上述算法及分析过程的正确性,在MatZFFTNloN+2N+DN·K(3)中产生一个正弦组合信号3随着细化倍数的增加,基带FFT和ZFFT的运算量x(t)=30cos(2m110t)+30cos(2x11145t)都会大幅度增加;zFF只有当细化频带较窄(此时+25cos(2x112.3t)+48cos(2m113.8t)无需数字滤波)或长序列的情况下,与基带FT相+50co(2x114.5t)比才具有运算量上的优势。分别利用基带FT和ZT对其进行谱分析ZFT算法存在自身的局限性,其存在的问题仿真条件:f=2048H,F点数N=1024,细化倍数D=50。基带FFT的频率分辨率4f=2H,历如下:(1)需要存放中间数据的内存空间巨大限制ZF的频率分辨率△f=0.04H。仿真结果如图了最大细化倍数2和图3所示。(2)采用具有线性相位的FIR数字滤波器实igure(n现抗混叠滤波,由于有限阶滤波器的吉布斯效应( Gibbs effect),滤波器截止频率处的频谱不可避免020040060080010001200会出现局部失真。(3)细化倍数越高,重釆样的选抽比越高,则细化带宽越窄。当需要细化的带宽较大时,必须进5行多次细化,这势必会增加计算量。Figure(4)频率成分调整较复杂。将FT和谱分析105110115130得到的频率成分调整到所选频带的频率成分式较Frequency(Hz复杂的过程,特别是为了避免低通抗混滤波器的边图3FF幅值频谱缘误差造成的频率混叠为了比较频率细化的效果,对图中谱线作了归化处理。图2中fgme(a)为原始信号,fgme(c)6小结为基带FYT处理后的幅值谱线,fgre(d)为移频后ZFT算法的关键在于利用傅立叶变换的移频基带FFT处理后的幅值谱线。由此图可以看出,基特性将感兴趣的高频段频率移至频谱原点,降低采带FFT的几个谱峰叠加为一个谱峰,各频率成分不可分辨。图3中fge(g)为重新采样后F处理样率重新釆样,从而获取较高的频率分辨率。它对后的幅值谱线,gure(h)为频率调整到实际频率处于获得某些特殊频段而不是整个带宽的信号细微的幅值谱线。此图中,因频率分辨率降低了D倍谱结构十分有用。该算法在实际工程技术中有较zF的幅值谱线中5条谱线清晰可见,说明ZF广泛的应用效果明显。参考文獻5ZF运算量和局限性讨论[1]胡广书.数字信号处理-理论、算法与实现[M]北京:清华大学出版社,1997当采用时域抽取FFT算法时,N点DT的复数[2] Vinay K ingle, John g proakis.数字信号处理及其乘法次数为l2N,复数加法次数为NN。为MATLAB实现[M].北京:电子工业出版社,1998[3]赵霞,熊小伏,郭珂.用细化频谱技术分析断路器简单起见,仅比较两种算法的复数乘法次数。操动机构振动信号[J.电力系统自动化,2003,(12):37设频率分辨率4f=fN,细化倍数D=△/404」f。要获得4/的分辨率,基带FFT的运算量为[4]丁康,谢明,张彼德等.基于复解析带通滤波器的FrTdN)lo复调制细化谱分析原理和方法[J.振动工程学报,2001,62(D14(1):30~35采用ZF算法,在复调制时只计算重采样的[5]宗孔德.多抽样率信号处理[M].北京:清华大学点,需N次复数乘法。同样,调制系数的计算也需N出版社,19基于 MATLAB复调制Z00M-FT算法的分析和实现旧WANFANG DATA文献链接作者:王力,张冰,徐伟, Wang li, Zhang bing, Xu Wei作者单位:王力,张冰, Wang Li, Zhang bing(江苏科技大学,镇江,212003),徐伟, Xu Wei(船舶系统工程部,北京,100036)刊名:舰船电子工程英文刊名SHIP ELECTRONIC ENGINEERING年,卷(期)2006,26(4)被引用次数:次参考文献(5条)1.宗孔德多抽样率信号处理19962.丁康;谢明;张彼德基于复解析带通滤波器的复调制细化谱分析原理和方法[期刊论文]振动工程学报2001(013.赵霞;熊小伏;郭珂用细化频谱技术分析断路器操动机构振动信号[期刊论文]电力系统自动化2003(12)4.陈怀琛数字信号处理教程- MATLAB释义与实现19985.胡广书数字信号处理一理论、算法与实现1997本文读者也读过(6条江波.唐普英基于复调制的ZooⅷFFT算法在局部频谱细化中的研究与实现[期刊论文]-大众科技2010(7)2.丁康.谢明.张彼德.赵玲.张晓飞. Ding Kang. Xie ming. Zhang bide. Zhao ling. ZHANG Xiaofei基于复解析带通滤波器的复调制细化谱分析原理和方法[期刊论文]-振动工程学报2001,14(1)3.罗利春. LUo Lic- hun zoom-FFT的改进、频谱反演与时-频局部化特性[期刊论文]-电子学报2006,34(1)4.戴振华.纪海林.徐运涛.DAⅠZhen-hua. JI Hai-1in.ⅫUYun-taoZ00MFFT算法在数字音频分析仪中的实现[期刊论文]-兵工自动化2007,26(10)5.黄镔.许婧.高峰.束洪春Z0OM-FFT在水电机组振动信号分析中的应用[期刊论文]-昆明理工大学学报(理工版)2002,27(5)6.王卫江改进的自适应Zoom-FFT算法研究[期刊论文]一电子技术应用2006,32(7)证文献(10条1.程兆刚.唐力伟.张淑琴.曹洪娜基于复调制Z0OM-FFT算法下阻尼比识别的研究[期刊论文]计算机与数字工程2012(1)2.刘树强.罗天.王宁.潘栋基于 Labview的异步电机转子断条检测[期刊论文]电子设计工程2011(3)3.王文森.邱宏安高精度超声流量检测系统设计[期刊论文]电声技术2011(2)4刘树强.罗天.谭兴文基于 Labview的笼型异步电动机转子断条故障在线检测系统[期刊论文]西南大学学报:自然科学版2011(9)5.王乐.苏小敏.杜林.李春化复白噪声中复正弦波频率估计方法硏究[期刊论文]火控雷达技术2011(36.周红霞.江佩勤.伍洲基于嵌入式系统的ZFFT移频轨道检测算法[期刊论文]通信技术2010(37.焦玮琦.陈特放基于局部频谱细化的轨道移频信号高精度检测[期刊论文]机车电传动2009(28.史瑞根.姚金杰基于 Labview的数字变频FFT设计[期刊论文]现代电子技术2009(7)9武中奇.杨世武丌FT算法在铁路移频信号分析中的应用及其DSP实现[期刊论文]铁道通信信号2008(7)10.时献江.张春喜.邵俊鹏异步电机断条故障诊断的细化包络方法[期刊论文]电机与控制学报2008(2)本文链接http://d.g.wanfangdata.com.cn/periodicaljcdzgc200604033.aspx

下载说明:请别用迅雷下载,失败请重下,重下不扣分!

发表评论

0 个回复

  • buckboost双向变换器Simulink仿真模型
    buckboost双向变换器Simulink仿真模型,完成DC-DC的双向转换
    2020-11-28下载
    积分:1
  • 语音增强报告(谱减法和维纳滤波
    语音增强报告(谱减法和维纳滤波) matlab编程实现
    2020-11-30下载
    积分:1
  • 反向传播算法推导—全连接神经网络
    反向传播算法是人工神经网络训练时采用的一种通用方法,在现代深度学习中得到了大 规模的应用。全连接神经网络(多层感知器模型,MLP),卷积神经网络(CNN),循环神 经网络(RNN)中都有它的实现版本。算法从多元复合函数求导的链式法则导出,递推的 计算神经网络每一层参数的梯度值。算法名称中的“误差”是指损失函数对神经网络每一层 临时输出值的梯度。反向传播算法从神经网络的输出层开始,利用递推公式根据后一层的误 差计算本层的误差,通过误差计算本层参数的梯度值,然后将差项传播到前一层(w, x,)+b这个神经元接受的输入信号为向量(),向量()为输入向量的组合权重,为徧置项,是标量。神经儿对输入冋量进行加权求和,并加上偏置项最后经过激活函数变换产生输出为表述简洁,我们把公式写成向量和矩阵形式。对每个神经元,它接受的来自前一层神经元的输入为向量,本节点的权重向量为,偏置项为,该神经元的输出值为先计算输入向量与权重向量的内积,加上偏置项,再送入一个函数进行变换,得到输出这个函数称为激活函数,典型的是函数。为什么需要激活函数以及什么样的函数可以充当激活函数,在之前的公众号文章“理解神经网终的激活函数”中已经进行了介绍。神绎网络一般有多个层。第一层为输入层,对应输入向量,神绎元的数量等于特征向量的维数,这个层不对数据进行处理,只是将输入向量送入下一层中进行计算。中间为隐含层,可能有多个。最后是输出层,神经元的数量等于要分类的类别数,输出层的输岀值被用来做分类预测。下面我们来看一个简单神经网络的例了,如下图所示这个网络有层。第一层是输入层,对应的输入向量为,有个神经元,写成分量形式为(),它不对数据做任何处理,直接原样送入下一层。中间层有个神经元,接受的输入数据为向量,输出向量为,写成分量形式为。第三个层为输出层,接受的输入数据为向量,输出向量为,写成分量形式为()。第一层到第层的权重矩阵为(,第二层到第三层的权重矩阵为()。权重矩阵的每一行为一个权重向量,是层所有神经元到本层某一个神经儿的连接权重,这里的上标表小层数如果激活函数选用函数,则第二层神经元的输出值为+(-(+0)+(1+(0)(-(()第三层神经元的输出值为如果把代入上面二式中,可以将输出向量表示成输出向量的函数。通过调整权重矩阵和偏置项可以实现不同的函数映射,因此神经网终就是一个复合函数需要解决的·个核心问题是·旦神经网络的结构(即神经元层数,每层神经元数量)桷定之后,怎样得到权重矩阵和偏置项。这些参数是通过训练得到的,这是本文推导的核心任务个简单的例子首先以前面的层神经网络为例,推导损失函数对神经网络所有参数梯度的计算方法假设训练样本集中有个样本()。其中为输入向量,为标签向量。现在要确定神经网络的映射函数:什么样的函数能很好的解释这批训练栟本?答案是神经网络的预测输出要尽可能的接近样本的标签值,即在训练集上最小化预测误差,如果使用均方误差,则优化的目标为:∑‖()-其中()和都是向量,求和项内部是向量的范数平方,即各个分量的平方和。上面的误差也称为欧氏距离损失函数,除此之外还可以使用其他损失函数,如交叉熵、对比损失等。优化目标函数的自变量是各层的权重矩阵和梯度向量,一般情况下无法保证目标函数是凸函数,因此这不是一个凸优化问题,有陷入局部极小值和鞍点的风险(对于这些概念和问题之前的公众号文章“理解梯度下降法”,“理解凸优化”中己经做了详细介绍)这是神经网络之前一直被诟病的一个问题。可以使用梯度下降法进行求解,使用梯度下降法需要计算出损失函数对所有权重矩阵、偏置向量的梯度值,接下来的关键是这些梯度值的计算。在这里我们先将问题简化,只考虑对单个样本的损失函数()-‖后面如果不加说明,都使用这种单样木的损失函数。如果计算出了对单个样木损失函数的棁度值,对这些梯度值计算均值即可得到整个目标函数的梯度值。和(要被代入到网络的后一层中,是复合函数的内层变量,我们先考虑外层的和。权重矩阵是一个x的矩阵,它的两个行分别为向量(和是个维的列向量,它的两个元素为()和()。网络的输入是向量,第一层映射之后的输出是向量首先计算损失函数对权重矩阵每个元素的偏导数,将欧氏距离损尖函数展开,有((+))(())6(如果,即对权重矩阵第行的元素求导,上式分了中的后半部分对来说是常数。根据链式法则有S()+()O如果,即对矩阵第二行的元素求导,类似的有:可以统一写成可以发现,第一个下标决定了权重矩阵的第行和偏置向量的第个分量,第二个下标决定了向量的第个分量。这可以看成是一个列向量与一个行向量相乘的结果,写成矩阵形式为上式中乘法⊙为向量对应元素相乘,第二个乘法是矩阵乘法。是个维列向量,+也是一个维列向量,两个向量执行⊙运算的结果还是个维列向量。是一个元素的列向量,其转置为维行向量,前面这个:维列向量与的乘积为的矩阵,这正好与矩阵的尺寸相等。在上面的公式中,权重的偏导数在求和项中由部分组成,分别是网络输出值与真实标签值的误差激活区数的导数+(),本层的输入值。神经网络的输出值、激活函数的导数值本层的输入值都可以在正向传播吋得到,因此可以晑效的计算出来。对所有训练样本的偏导数计算均值,可以得到总的偏导数对偏置项的偏导数为:如果上式分子中的后半部分对来说是常数,有:()⊥()如果类似的有这可以统写成:写成矩阵形式为偏置项的导数由两部分组成,分别是神经网络预测值与真实值之间的误差,激活函数的导数值,与权重矩阵的偏导数相比唯一的区别是少了。接下来计算对和的偏导数,由于是复合函数的内层,情况更为复杂。()是个的短阵,它的个行向量为(),(,(,(。偏置项()是维向量,个分量分别是(),(,(),(。首先计算损失函数对的元素的偏导数:而上式分子中的两部分都有,因此都与有关。为了表述简活,我们令:根据链式法则有:其巾((和和都是标量和()是两个()向量的内积,的每一个分量都是()的函数。接下来计算和这里的一是个向量,衣示的每个分量分别对求导。当时有:后面个分量相对于求导变量(都是常数。类似的当时有:()0)(()和时的结果以此类推。综合起来有:同理有:()十如果令合并得到()()[()-)。()。()写成矩阵形式为()最后计算偏置项的偏导数()类似的我们得到:合并后得到()写成矩阵形式为:(0)至此,我得到了这个简单网络对所有参数的偏导数,接下来我们将这种做法推广到更般的情况。从上面的结果可以看岀一个规律,输出层的权重矩阵和偏置向量梯度计算公式中共用了()-)()对」隐含层也有类似的结果完整的算法现在考虑一般的情况。假设有个训练样本(),其中为输入向量,为标签向量。训练的目标是最小化样木标签值与神经网络预测值之闩的误差,如果使用均方误差,则优化的目标为:其中为神经网络所有参数的集合,包括各层的权重和偏置。这个最优化问题是·个不带约束条件的问题,可以用梯度下降法求解。上面的误差函数定义在整个训练样本集上,梯度下降法每一次迭代利用了所有训练样本,称为批量棁度卜降法。如果样木数量很大,每次迭代都用所有样木进计算成木太高。为了解决这个问题,可以采用单样本梯度下降法,我们将上面的损失函数写成对单个样本的损失函数之和:定义对单个样本()的损失函数为)=-()如果采用单个样本进行迭代,梯度下降法第次迭代时参数的更新公式为:nV如果要用所有样本进行迭代,根据单个样本的损失函数梯度计算总损失梯度即可,即所有样本梯度的均值用梯度下降法求解需要初始化优化变量的值。一般初始化为一个随机数,如用正态分布(a)产生这些随机数,其中G是一个很小的正数到日前为止还有一个关键问题没有解决:日标函数是一个多层的复合函数,因为神经网络中每一层都有权重矩阵和偏置向量,且每一层的输出将会作为下一层的输入。因此,直接计算损失函数对所有权重和偏置的梚度很复杂,需要使用复合函数的求导公式进行递推计算几个重要的结论在进行推导之前,我们首先来看下面几种复合函数的求导。又如下线性映射函数:其中是维向量,是×的矩阵,是维向量。问题:假设有函数,如果把看成常数,看成的函数,如何根据函数对的梯度值Ⅴ计算函数对的梯度值Ⅴ?根据链式法则,由于只和有关,和其他的≠无关,因此有:c∑(对于的所有元素有:写成矩阵形式为:问题:如果将看成常数,将看成的函数,如何根据V计算Ⅴ?由于任意的和所有的都有关系,根据链式法则有写成矩阵形式为这是一个对称的结果,在计算函数映射时用矩阵乘以向量得到,在求梯度时用矩阵的转置乘以的梯度得到的梯度。问题:如果有向量到向量的映射:
    2020-12-09下载
    积分:1
  • 200SMART的MODBUS主站通讯变频器
    【实例简介】
    2021-07-06 00:31:05下载
    积分:1
  • RainFlow雨流计数法计算材料疲劳强度
    用雨流计数法计算材料疲劳强度的一种方法,这里提供了具体的算法
    2020-07-03下载
    积分:1
  • 华为风格手册 V2.0-CH.chm
    【实例简介】本文的制订旨在给所有从事华为产品资料写作工作的人员提供统一的写作规范,使大家在实际工作中有共同标准可以遵循。
    2021-10-29 00:33:47下载
    积分:1
  • 四参数Logistic拟合y=(a-d)/[1+(x/c)^b]+d
    在开发医学分析的软件中经常用到四参数Logistic拟合算法,例如化学发光、定量分析等。找了很久,咨询了很多人都敝帚自珍。我终于找到了ImageJ这个好东东,不过是Java的代码,有一天,我终于把它改成Delphi的了。再一天我把它传到网上去,大家分享。
    2020-11-27下载
    积分:1
  • 基于神经网络的手写数字识别MATLAB源码
    基于BP神经网络的手写数字识别MATLAB源码,包括GUI界面
    2020-12-02下载
    积分:1
  • STM32+ESP8266机智云MCU方案APP控制
    【实例简介】STM32+ESP8266机智云MCU方案,手机APP端控制LED亮灭,显示DHT11温湿度,STM32内部温度,NTC温度
    2021-11-12 00:36:11下载
    积分:1
  • x264库(64位and32位)
    x264编译库,Windows环境下包括64位和32位资源,可使用
    2020-12-03下载
    积分:1
  • 696516资源总数
  • 106611会员总数
  • 19今日下载