专利名称：一种基于数字信号处理器快速频谱细化的方法
技术领域：
本发明属于数字信号频域分析技术领域，在数字信号频域分析领域中频谱分析设备通常会需要使用数字信号处理器进行细化的频谱分析，本发明涉及快速、高效的频谱分析方法。
背景技术：
在数字信号频域分析领域中频谱分析设备通常会需要使用数字信号处理器进行细化的频谱分析，已有的频谱分析设备使用离散傅里叶变换分析或离散快速傅里叶变换分析。离散傅里叶变换分析的方式，虽然能够计算出任意点数任意频谱范围的频谱，但如果一次分析大量的数据，数字信号处理器工作会非常繁琐，运算使用的时间非常长。数字信号处理器进行离散快速傅里叶变换分析虽然运算快速，但只能得出2的N次方点数的频谱数据， 并且计算出的频谱是全景频谱，通常会超过频谱分析设备所关心的频率范围，很多频谱数据并不需要被计算出来。
发明内容
本发明目的是解决频谱分析设备中数字信号处理器进行大数据量的频谱分析的工作时间过长，运算复杂的问题，提供频域分析领域中一种基于数字信号处理器快速频谱细化的方法。能够实现本发明的硬件结构包括I)高速模拟与数字信号转换芯片用于采样需要分析的模拟信号，将模拟信号转换为数字信号处理器能够使用的数字信号，能够工作在频谱分析设备关心的频率范围的终止频率的两倍以上的采样率上。该芯片可以采用National Semiconductor公司12位170MSPS的高速模拟与数字信号转换芯片ADC12C170。2)、数字信号处理器用于接收高速模拟与数字信号转换芯片输出的数字信号，并对数字信号进行离散快速傅里叶变换分析和离散傅里叶变换分析，数字信号处理器需要有内部数据存储器，能够处理浮点数据，能够实现对外部扩展数据存储器进行的存储和读取操作，可以采用Analog Devices公司的SHARC系列32位浮点数字信号处理器芯片ADSP-21369。3)、外部扩展数据存储器用于存储数字信号处理器细化频谱分析过程中使用的窗系数、蝶形算子、采样波形、频谱细化波形等的一些大数据量的数据，可以采用hynix公司的16位同步动态随机存储器芯片H57V2562GTR-75C。本发明的硬件结构示意图如附图I所示，本发明的程序开发的编译环境可以为AnalogDevices 公司的 VisualDSP++5. O 软件。本发明方法的全部流程如附图2所示，其中步骤I到步骤4的详细过程如附图3所示，步骤5的详细过程如附图4所示，本发明技术方案的实现步骤如下步骤I、生成窗系数和蝶形算子数据
根据频谱分析设备作频谱分析的频率分辨率带宽值Frtw和高速模拟与数字信号转换芯片的数据采样率Fsam确定所需的采样数据量Lsam，使用窗函数的补偿系数为Pwin，所需的采样数据量如公式I所述Lsam = [ (Fsam+Frbw) X Pwin+0. 5](I)公式(I)中“ □”为数据取整的符号，如[2.6] = 2，使用编译器提供的窗函数生成数量为Lsam的窗系数，定义位于外部扩展数据存储器中的Lsam个浮点类型的数据区作为窗系数的存储器，命名为第一存储器(I)。确定使用频谱分析的点数L，本方法能够处理L=2aXN点的频谱分析，a和N均为正整数，确定N和a值，a —般可以取9至12，并确保数字信号处理器的内部存储器能够分配23个浮点数类型的数据区，根据L ^ Lsam推导出N的值，确保2aXN的值不小于所需的采样数据量Lsam,且2a远大于N。离散快速傅里叶变换的点数M为2的a次方，即M=2a，L=MXN=2aXN, (Μ >> N)，根据M和N确定以下两组的复数数据。
第一组数据为L个复数数据，为NX M蝶形算子(计算方法见公式(2 ))，在步骤5. 2使用，设η和m为正整数，将在步骤4、步骤5中继续使用，其中第(η — I) XM+m个数据为ぐ-t(Kn<NjKm〈M)，j 为虚数単位，f= — I響—ι,πι-1= e—i(T7)(n—⑵将公式(2)转换成能够让数字信号处理器实现的NXM蝶形算子计算公式(3)Re(W 1,m"1) = cos (2π(η - l)(m - 1)/L)= -sin (2ιτ(η - l)(m — 1)/L)(3)定义位于外部扩展数据存储器中LX2个浮点类型的数据区作为NXM蝶形算子数据存储器，命名为第二存储器(2)，从第二存储器(2)的第一个位置开始奇数位置存储复数数据的实部，偶数位置存储复数数据的虚部。第二存储器(2)的第I个位置存储Re(WL3 °)(n=l, m=l)，第二存储器(2)的
第2个位置存储im(Wt P)(n=l，m=l),第二存储器(2)的第3个位置存储Re^.+' ')
(n=l, m=2),第二存储器(2)的第4个位置存储丨m(v_ご')(n=l, m=2),......，第二存储
器(2)的第(n-l)X2XM+(m-l)X2+l个位置存储第二存储器⑵的第
(n-1) X2XM+On-1) X2+2个位置存储丨111(ノ<~1'塑ィ),……，第二存储器(2)的第2XL-I个
位置存储Re(wf_1‘ M_i)(n=N,m=M),第二存储器(2)的第2XL个位置存储Im(Vif_1 ‘ *卜1)(n=N，m=M)。第二组数据为NXN个复数数据，为NXN蝶形算子(计算方法见公式(4))，为步骤5. 3的N点离散傅里叶变换使用，设定其中第Oi1 - I) XN+n2个的数据为(Κη^Ν,Κη^Ν), j 为虚数単位，j2= — I:W^1"1,112"1 = e—D(4)将公式(4)转换成能够让数字信号处理器实现的NXN蝶形算子计算公式(5)
Re(W^= cos (2π(Πι - l)(n2 - 1)/N)
lni(Wf Mr1) = —sin (2π(η, - l)(n2 - 1)/N)定义位于数字信号处理器的内部数据存储器中NXNXZ个浮点类型的数据区作为N X N蝶形算子数据存储器，命名为第三存储器(3 )，从第三存储器(3 )的第一个位置开始奇数位置存储复数数据的实部，偶数位置存储复数数据的虚部。第三存储器(3)的第I个位置存储Re(W^ 0Xn1=I, n2=l)，第三存储器(3)的第2个位置存储Im(W^ cKn1=I, n2=l)，第三存储器(3)的第3个位置存储Re(W:' 1Kn1=I,n2=2)，第三存储器(3)的第4个位置存储Im(wJ 1Xn1=I, n2=2),……，第三存储器(3)的第(Ii1 — I) X2XN+(n2 — I) X2+1个位置存储Re(W^n:_1)，第三存储器(3)的第(Ii1 — I) X2XN+(n2 — I) X2+2个位置存储_1)，……，第三存储器(3)的第 2XNXN — I个位置存储ReOVifO1)U1=N, n2=N)，第三存储器(3)的第2XNXN个位置存储Im(W^nt)(Ii1=N, n2=N)。本步骤在确定了频谱分辨率带宽Frtw和数据采样率Fsam后只需操作一次，通过编译器提供的功能将数据分别存储在以上三个存储器上，供每次采样后的频谱分析使用。步骤2、生成与频谱分析设备关心的频谱范围有关的第四、第五存储器的数据设定频谱分析设备关心的频谱范围的起始频率为Fstart,终止频率为Fs_，起始频率Fstart在L点频谱上的位置为Istart (计算见公式(6))，终止频率Fsttjp在L点频谱上的位置为Is_ (计算方法如公式(7)):Istart = [(L + Fsam) XFstart-O. 5](6)Istop = [(L + Fsam) XFstop-O. 5](7)定义位于外部扩展数据存储器中M个整型类型的数据区的存储器，作为步骤5离散傅里叶变换的数量存储器，命名为第四存储器(4)，和位于外部扩展数据存储器中L个整型类型的数据区的存储器，作为步骤5离散傅里叶变换的参数存储器，命名为第五存储器(5 )，将第四存储器(4 )和第五存储器(5 )初始化为O值。设定整数I从Istart到Istop逐个数取值(Istart ( I ( Is_)，每次完成以下操作求取I整除M后得到的商数1_和余数Inwd，第四存储器(4)的第IMd+l个位置读取出来并加I后得到数值N·，再将Nnum存储入到第四存储器(4)的第Im()d+1个位置，将1_存储到第五存储器(5)的第I-XN+N·个位置。该步骤运算出的两个存储器的数值，即能够在步骤5运算步骤中实现只运算出从
Fstart到Fsttjp的频谱数据。如果频谱分析设备所关心的频谱范围不变，即起始频率Fstart的数值和终止频率Fstop的数值都没有变化，步骤2只在开始操作一次，不需要每次采样都操作，只有起始频率Fstart和终止频率Fsttjp有变化时，需要重新操作步骤2。步骤3、将模拟信号转换成数字信号，数字信号加窗后存储到外部扩展数据存储器定义位于外部扩展数据存储器中L个浮点类型的数据区作为加窗后的采样数据存储器，命名为第六存储器(6)，使用高速模拟与数字信号转换芯片将需要分析的模拟信号转换成数字信号，连续采样数量为Lsam的数字信号数据，将这些采样数据逐个乘以从第一存储器(I)的第一个位置开始读取的每个数据，从第一个位置开始逐个存储到第六存储器
(6)，共存储Lsam个位置，第六存储器(6)的第Lsam+1个位置到第L个位置共计L — Lsam个位置存储为数值O。步骤4、对数据抽取分组，在数字信号处理器的内部数据存储器中进行离散快速傅里叶变换定义位于数字信号处理器的内部数据存储器中的M个浮点数类型的数据区，用于存储M点离散快速傅里叶变换的输入数据，命名为第七存储器(7)，定义位于外部扩展数据存储器中LX2个浮点类型的数据区为本步骤的结果存储器，命名为第八存储器(8)，存储本步骤的N次M点离散快速傅里叶变换的复数結果。
·
将存储在第六存储器(6)中的L点浮点数据抽取成N组M点数据，依次作M点的离散快速傅里叶变换。正整数η的值为1，判断η的值如果η小于或等于N，从第六存储器(6)的数据中读取M点的数据，读取第六存储器(6)中的第n，N+n，2 X Ν+η，……，L 一 Ν+η个位置的数据，存储于第七存储器(7)，对这M点数据做M点的离散快速傅里叶变换，得出M点的复数結果，存储到第八存储器(8)的第(n— I) ΧΜΧ2+1到ηΧΜΧ2位置的共计MX2个位置，其中数据的实部存储在奇数位置，虚部存储在偶数位置，η的值加1，重新判断η的值，重复以上操作；如果η大于N，则结束本步骤。步骤5、生成频谱分析设备关心的频谱范围内的频谱细化波形数据定义位于数字信号处理器的内部数据存储器中ΝΧ2个浮点类型的数据区的存储器，作为步骤5. 3的运算输入存储器，命名为第九存储器(9)，定义位于外部扩展数据存储器中L个浮点类型的数据区为频谱数据存储器，命名为第十存储器(10)，正整数m的值为1，判断m的数值如果m小于或等于M，则完成以下步骤5. I到步骤5. 3操作；否则结束本步骤，如附图3所示步骤5. I.如果第四存储器(4)中的第m个位置读取的数据为0，跳出本轮循环，m加I进入下一轮循环，否则完成以下其余的操作，步骤5. 2.正整数η的值为1，判断η的值当η小于或等于N时，如果m或η有ー个数值为1，则从第八存储器(8)中的第(n-1) X2XM+(m — I) X2+1个位置开始的两个位置读取出两个浮点数直接存储到第九存储器(9)的第(η — I) X2+1个位置开始的两个位置，否则将第八存储器(8)中的第(η —I) X2XM+(m 一 I) X2+1个位置开始的两个位置读取出两个浮点数第一个数作为复数的实部和第二个数作为虚部与从第二存储器(2)同样位置读取到的复数相乘，将复数结果存储到第九存储器(9)的第(η — I) X 2+1个位置开始的两个位置，η的数值加1，重新判断η的值，重复以上操作，当η大于N吋，则进行步骤5. 3的操作。步骤5. 3.正整数η的值为1，判断η的值如果η小于或等于第四存储器(4)中的第m个位置读取的数值，则读取出第五存储器(5)的第(m — I) XN+n个位置的数据1_。从第九存储器(9)中依次读取出N个复数，第k次从第kX 2个位置读取出数为Xk( I < k < N)，
从第三存储器(3)读取出N个复数，第k次从第(kXN+Iqu。)X 2个位置读取出复数wf, (l^k^ N)，利用以下计算公式(8)求出离散傅里叶变换结果Dn，
K
k=1 (8)将Dn的模的平方|Dn|2 (实部数值和虚部数值的平方和)存储入第十存储器(10)的第Ι-ΧΜ+m位置，η的数值加1，重新判断η的值，重复步骤5. 3操作；如果η大于第四存储器(4)中的第m个位置读取的数值，m的值加1，重新判断m的值，重复完成步骤5. I到步骤5. 3的操作。复数求取模值的计算公式如(9)式 |Dj2= [Re (Dn)]2+[Im(Dn)]2(9)运算完成后，第十存储器(10)即得到基于L点离散傅里叶频谱分析的从Fstart到Fstop的线性频谱数据，对各个数据求取常用对数后乘以10，得到对数频谱数据，此数据即为频谱细化波形数据。本发明的优点和有益效果本发明为了解决频谱分析设备中数字信号处理器进行大数据量的频谱分析的工作时间过长，运算复杂的问题，所以提供了一种快速、简化的2的N次方的整倍数点数的频谱分析技术方案，能够省去大量的数据分析使用的时间。将L分解成L=MXN的形式，L点的离散傅里叶变换分解成作N个M点的离散快速傅里叶变换和M个N点的离散傅里叶变换，离散傅里叶变换的复数乘法计算量为L2，本方法的复数乘法计算量为匪log2M+L+MN2，能够按照频谱分析设备所关心的频率范围，只计算出部分频谱波形，大幅减少了数字信号处理器的运算量，复数乘法计算量最多为匪log2M+L+(MN2)/2，提高数字信号处理器的工作效率。假设作L=10240点的离散傅里叶变换，分解为M=2048，N=5，本方法比直接的离散傅里叶变换运算量减少了 99. 7%，比离散快速傅里叶变换(16384点，必须为2的N次方)运算量减少了 35.3%。本发明将需要在外部扩展数据存储器中进行的一次大数据量的运算分解成能够在数字信号处理器内部数据存储器中进行的多次小数据量运算，能够大幅度缩短数字信号处理器的运算时间。本发明的创新点是采用一种方法，使用数字信号处理器，将需要在外部扩展数据存储器中进行的一次大数据量的离散傅里叶变换细化频谱分析运算分解成能够在数字信号处理器内部数据存储器中进行的多次小数据量的离散快速傅里叶变换和离散傅里叶变换的运算。其特点在于使用这种方法能够减少数字信号处理器的细化频谱分析运算量，能够只分析得到频谱分析设备关心频率范围内的频谱，减少细化频谱分析的数据量，实现2的N次方的整倍数的点数的频谱分析，大大缩短数字信号处理器的运算处理时间，提高频谱分析设备工作效率。
图I本发明的硬件结构示意图；图2频谱细化流程图；图3步骤I至步骤4的流程图4步骤5的流程图；图5布莱特曼窗系数波形；图6采样到的数字信号；图7加窗后的数字信号；图8将加窗后的数字信号扩展到L点；图9本发明生成的频谱分析设备所关心频谱波形。
具体实施例方式实施例 如图I所示，结合一个实例说明本发明方法的硬件结构选择的数字信号处理器为Analog Devices公司的SHARC系列32位浮点数字信号处理器芯片ADSP-21369，外部扩展数据存储器为hynix公司的16位同步动态随机存储器芯片H57V2562GTR-75C，存储容量为256Mbit，高速模拟与数字信号转换芯片为National Semiconductor公司12位170MSPS的高速模拟与数字信号转换芯片ADC12C170。步骤I、设定频谱分析的频率分辨率带宽值Faw为30KHz (O. 03MHz)，高速模拟与数字信号转换芯片ADC12C170的数据采样率Fsam为166. 4MHz，使用布莱克曼窗函数(补偿系数Pwin为I. 62)，所需的采样数据量经公式(I)计算为=Lsam= [(166. 4 + 0. 03) X I. 62+0. 5]=8986个数据，使用布莱克曼窗函数生成数量为8986的布莱克曼窗系数，窗系数的波形如附图5所示，通过编译器提供的功能存储入第一存储器(I)。整数a—般可以取9至12，这里数值a取值ll，M=2n为2048，根据L彡8986推导出N的值，确定N取值5，L=2n X 5=2048 X 5=10240不小于所需的采样数据量8986,2048远大于5，进行10240的频谱分析。按照步骤I中所述的方法，求出5X2048蝶形算子的数据。定义位于外部扩展数据存储器中20480个浮点类型的数据区作为存储5 X 2048蝶形算子的第二存储器(2)，从第ニ存储器(2)的第一个位置开始奇数位置存储复数数据的实部，偶数位置存储复数数据的虚部。第二存储器(2)的第I 个位置存储Re(mr °)=cos(2X π XOX0/10240)=1 (n=l，m=l),第二存储器(2)的第2 个位置存储___)=—sin(2X π XOX0/10240)=()(η=1,m=l),第二存储器(2)的第3 个位置存储Re(じ:.')=005(2X i XOX 1/10240)= I (n=l,m=2),第二存储器(2)的第4 个位置存储 Im(Vこ')=-sin(2X π XOX I/10240)=0 (η=I,m=2),......，第二存储器(2)的第(η — I) X 4096+(m — I) X 2+1个位置存储Rei:.. )=cos(2 X πΧ (η— I) X (m — )/10240)，
第二存储器(2)的第(η — I) X 4096+(m 一 I) X 2+2个位置存储
Im(じ .)=—.sin(2X- X(η I)X(m— Ij/10240),......，第 ニ 存储器(2) 的第20479个位置存储 Re(wf"1, M^)=cos(2X π Χ4Χ2047/10240)=0 3 67(η=Ν, m=M),第二存储器(2)的第20480 个位置存储 Im(wf_1 M_1 )= — sin(2 X π X 4 X2047/10240)=0.9518(η=Ν，m=M)。按照步骤I中所述的方法，求出5X5蝶形算子的数据。定义位于数字信号处理器的内部数据存储器中50个浮点类型的数据区作为存储5X5蝶形算子的第三存储器(3)，从第三存储器(3)的第一个位置开始奇数位置存储复数数据的实部，偶数位置存储复数数据的虚部。第三存储器(3)的第I 个位置存储Re(Wご °)=cos(2X π XOX0/5)=1 Cn1=Ln2=I),第三存储器(3)的第2 个位置存储3)=—sin(2X π XOX0/5)=0Cn1=I,n2=l)，第三存储器(3)的第3 个位置存储Re(w^ VcosiSX π XOX 1/5)=1 U1=I, η2=2)，第三存储器(3)的第4 个位置存储Im(W^ sin(2X π XOX 1/5)=0U1=I,n2=2),......，第三存储器(3)的第(Ii1 — I) X2XN+(n2 — I) X 2+1个位置存储Re(W;;_1i.n:-i)=cos(2X π X(m — I)X(Ii1 —1)/5),第三存储器(3)的第(Ii1 — I) X2XN+(n2 — I) X 2+2个位置存储
^=-SinpX π X(η,-I)X(n,-1)/5),......，第三存储器(3)的第2XNXN— I 个位置存储x_1)=cos(2X π X4X4/5)=0.3090 Cn1=N, n2=N),第三存储器(3)的第2XNXN 个位置存储sin (2X π χ4X4/5)=— 0.9511(1^=1 η2=Ν)；步骤2、设定起始频率Fstart为5MHz，终止频率Fsttjp为65MHz，经公式(6) (7)计算
Istart 和 Istop 的值，Istart=[(10240 + 166. 4) X5+0. 5] = 308Istop=[ (10240^-166. 4) X65+0. 5] = 4000按照步骤2中所述的方法，求取第四存储器(4)和第五存储器(5)的数值I=Istart=308整除M=2048后得到的商数1_=0和余数Im()d=308，第四存储器(4)的第309个数值读取出来(为O)并加I后得到1，再将I存储入到第四存储器(4)的第309个数值，将O存储到第五存储器(5)的第1541个数值。1=309整除2048后得到的商数O和余数309，第四存储器(4)的第310个数值读取出来(为O)并加I后得到1，再将I存储入到第四存储器(4)的第310个数值，将O存储到第五存储器(5)的第1546个数值。 ......1=2356整除2048后得到的商数I和余数308，第四存储器(4)的第309个数值读取出来(为I)并加I后得到2，再将2存储入到第四存储器(4)的第309个数值，将I存储到第五存储器(5)的第1542个数值；步骤3、定义10240个浮点类型的数据区作为存储加窗后的采样数据的第六存储 器(6)，使用高速模拟与数字信号转换芯片连续采样数字信号数量为8986的数据，采样数据波形如附图6所示，将这些采样数据逐个乘以从第一存储器(I)读取到的数据，加窗后的采样的数字信号数据波形如附图7所示，从第一个位置开始逐个存储到第六存储器(6)，共存储8986个位置，第六存储器(6)的第8987个位置到第10240个位置共计1254个位置存储为数值O，第六存储器(6)数据波形如附图8所示；步骤4、正整数η的值为1，η小于5，从第六存储器(6)的数据中读取2048点的数据，读取第六存储器(6)中的第1，6，11，……，10236个位置的数据，存储于位于数字信号处理器内部数据存储器的第七存储器(7)，对这2048点数据做2048点的离散快速傅里叶变换，得出2048点的复数结果，存储到第八存储器(8)的第I到4096的位置，其中数据的实部存储在奇数位置，虚部存储在偶数位置，η的值加I。正整数η等于2，η小于5，读取第六存储器(6)中的第2，7，12，……，10237个位置的数据，存储于第七存储器(7)，对这2048点数据做2048点的离散快速傅里叶变换，得出2048点的复数结果，存储到第八存储器(8)的第4097到8192的位置，η的值加1，重新判断η的值，重复以上操作。......当正整数η等于5，读取第六存储器(6)中的第5，10，15，……，10240个位置的数据，存储于第七存储器(7)，对这2048点数据做2048点的离散快速傅里叶变换，得出2048点的复数结果，存储到第八存储器(8)的第16385到20480的位置，η的值加1，此时η大于5，结束本步骤；步骤5、定义4096个浮点类型的数据区的第九存储器(9)，定义10240个浮点类型的数据区的第十存储器(10)，正整数m的值为1，判断m的数值如果m小于或等于2048，则完成以下步骤5. I到步骤5. 3操作；否则结束本步骤步骤5. I.如果第四存储器(4)中的第m个位置读取的数据为0，跳出本轮循环，m加I进入下一轮循环，否则完成以下其余的操作。步骤5. 2.正整数η的值为I，判断η的值当η小于或等于5时，如果m或η有一个数值为1，则从第八存储器(8)中的第(n- I) X4096+(m-I) X2+1个位置开始的两个位置读取出两个浮点数直接存储到第九存储器(9)的第(n - I) X2+1个位置开始的两个位置，否则将第八存储器(8)中的第(η —I) X4096+(m 一 I) X 2+1个位置开始的两个位置读取出两个浮点数第一个数作为复数的实部和第二个数作为虚部与从第二存储器(2)同样位置读取到的复数相乘，将复数结果存储入到第九存储器(9)的第(n - I) X 2+1个位置开始的两个位置，η的数值加1，重新判断η的值，重复以上操作；
当η大于5时，则进行步骤5. 3的操作；步骤5. 3.正整数η的值为1，判断η的值如果η小于或等于第四存储器(4)中的第m个位置读取的数值，则读取出第五存储器(5)的第(m — l)X5+n个位置的数据1_。从第九存储器(9 )中依次读取出5个复数，第k次从第k X 2个位置读取出数为Xk( I < k < N)，从第三存储器(3)读取出5个复数，其中第k次从第(kX5+Iqu。)X2个位置读取出复数
W^ucClN)，利用公式(8)求出离散傅里叶变换结果Dn。将Dn的模的平方|Dn|2 (实部数值和虚部数值的平方和)存储入第十存储器(10)的第I_X2048+m位置，η的数值加1，重新判断η的值，重复以上操作；如果η大于第四存储器(4)中的第m个位置读取的数值，m的值加1，重新判断m的值，重复完成步骤5. I到步骤5. 3的操作。 运算完成后，第十存储器(10)即得到基于10240点离散傅里叶频谱分析的从5MHz到65MHz的线性频谱数据，对各个数据求取常用对数后乘以10，即得到对数频谱数据，频谱波形如附图9所示。
权利要求
1.一种基于数字信号处理器快速频谱细化的方法，能够实现本发明的硬件结构包括 1)高速模拟与数字信号转换芯片用于采样需要分析的模拟信号，将模拟信号转换为数字信号处理器能够使用的数字信号，能够工作在频谱分析设备关心的频率范围的终止频率的两倍以上的米样率上； 2)、数字信号处理器用于接收高速模拟与数字信号转换芯片输出的数字信号，并对数字信号进行离散快速傅里叶变换分析和离散傅里叶变换分析，数字信号处理器需要有内部数据存储器，能够处理浮点数据，能够实现对外部扩展数据存储器进行的存储和读取操作； 3)、外部扩展数据存储器用于存储数字信号处理器细化频谱分析过程中使用的窗系数、蝶形算子、采样波形、频谱细化波形的一些大数据量的数据； 其特征在于，该方法的步骤如下 步骤I、生成窗系数和蝶形算子数据 根据频谱分析设备作频谱分析的频率分辨率带宽值Frtw和高速模拟与数字信号转换芯片的数据采样率Fsam确定所需的采样数据量Lsam，使用窗函数的补偿系数为Pwin，所需的采样数据量如公式(I)所述； Lsam = [ (Fsam^-Frbw) X Pwin+0. 5](O 公式(I)中“ □”为数据取整的符号，使用编译器提供的窗函数生成数量为Lsam的窗系数，定义位于外部扩展数据存储器中Lsam个浮点类型的数据区作为窗系数的存储器，命名为第一存储器(I);确定使用频谱分析的点数L，本方法能够处理L=2aXN点的频谱分析，a和N均为正整数，确定N和a值，a—般取9至12，并确保数字信号处理器的内部存储器能够分配23个浮点数类型的数据区，根据L ≥ Lsam推导出N的值，确保2aXN的值不小于所需的采样数据量Lsam，且2a远大于N; 离散快速傅里叶变换的点数M为2的a次方，即M=2a，L=MXN=2aXN, M >> N，根据M和N确定以下两组的复数数据； 第一组数据为L个复数数据，为NXM蝶形算子，计算方法见公式(2)，在步骤5. 2中使用，设η和m为正整数，将在步骤4、步骤5中继续使用，其中第(η — I) XM+m个数据为K"' Kn<N, l〈m〈M，j 为虚数单位，j2= — I ；
全文摘要
一种基于数字信号处理器快速频谱细化的方法。该方法使用数字信号处理器，将大量数据的细化频谱分析分解成多次小量数据分析，得到频谱分析设备关心频率范围内的2的N次方的整倍数的点数的频谱分析，减少分析的数据量和运算量，缩短处理时间，提高频谱分析设备工作效率。本发明的方法包括生成窗系数和蝶形算子数据；生成与频谱范围有关的数据；模拟信号转换成数字信号，存储到外部扩展数据存储器；对数据抽取分组，在数字信号处理器内部数据存储器中进行离散快速傅里叶变换；生成频谱细化波形数据。本方法在10240点频谱分析时比离散傅里叶变换运算量减少了99.7%，比离散快速傅里叶变换(16384点，必须为2的N次方)运算量减少了35.3%。
文档编号G06F17/00GK102841880SQ20121032610
公开日2012年12月26日 申请日期2012年9月6日 优先权日2012年9月6日
发明者孟敬, 李晓光 申请人:天津市德力电子仪器有限公司