CN1971306A - 调频信号频谱反转的系统及方法 - Google Patents

调频信号频谱反转的系统及方法 Download PDF

Info

Publication number
CN1971306A
CN1971306A CN 200510086970 CN200510086970A CN1971306A CN 1971306 A CN1971306 A CN 1971306A CN 200510086970 CN200510086970 CN 200510086970 CN 200510086970 A CN200510086970 A CN 200510086970A CN 1971306 A CN1971306 A CN 1971306A
Authority
CN
China
Prior art keywords
signal
frequency
spectrum
cos
πkt
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN 200510086970
Other languages
English (en)
Other versions
CN100487483C (zh
Inventor
王一丁
杨虹
赵实
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Graduate School of CAS
University of Chinese Academy of Sciences
Original Assignee
University of Chinese Academy of Sciences
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Chinese Academy of Sciences filed Critical University of Chinese Academy of Sciences
Priority to CN 200510086970 priority Critical patent/CN100487483C/zh
Publication of CN1971306A publication Critical patent/CN1971306A/zh
Application granted granted Critical
Publication of CN100487483C publication Critical patent/CN100487483C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明涉及数据采集技术领域,特别是一种用于调频信号调频斜率反转的调频信号频谱反转的系统及方法。调频信号频谱反转的系统使用了欠采样的系统设计方法。当一路调频信号输入该系统后,通过正交解调方法获得两路位于奇次频带中的带通信号,然后进行带通信号的欠采样获得两路基带信号,最后通过正交调制方法获得一路输出调频信号,输出调频信号具有与输入调频信号相反的调频斜率。

Description

调频信号频谱反转的系统及方法
技术领域
本发明涉及数据采集技术领域,特别是一种用于调频信号调频斜率反转的调频信号频谱反转的系统及方法。
背景技术
编码辐射定标技术是合成孔径雷达(SAR)辐射定标技术的发展方向之一,其关键技术在于有源编码转发器的研制。在有源编码转发器的设计中需要完成输入调频信号调频斜率的反转。
目前,欠采样方法被广泛应用于高速数据采集,无线电信号分析等高速,超高速数字处理领域。虽然欠采样之后的数字频谱结构复杂,但是欠采样能够保留完整的信号特征,且具有频谱搬移的能力。也就是说欠采样方法具有信号频谱结构反转的能力。
若某一带通信号S(t)的上截止频率为fh,下截止频率为fl,考虑采样后的频谱无混叠,根据奈奎斯特采样定理,最小的采样频率为fs=2fh。然而根据欠采样定理,允许无误地重建带通信号S(t)的最小采样频率为:
F s = 2 ( f h - f l ) ( 1 + M N ) = 2 B ( 1 + M N ) ≈ 2 B - - - ( 1 )
其中,B=fh-fl,M=[fh/(fh-fl)]-N,N为不超过fh/(fh-fl)的最大正整数。
当信号S(t)的带宽B<<fh时,欠采样定理确定的最小采样频率fs将远小于奈奎斯特采样定理确定的最小采样频率fs。除了采样频率的降低外,欠采样还能实现信号频谱结构的反转,下面给出证明。
利用满足式(1)的欠采样频率Fs对信号S(t)进行等间隔采样,其频谱将位于[kB,(k+1)B],k=0,1,2,...的不同频带上,其中基带频谱位于[0,B]。
设位于[0,B]基带中偏离其中心频率f00分别为±Δf的两个余弦信号为:
s 1 0 ( t ) = A cos [ 2 π ( f 00 - Δf ) t ] - - - ( 2 )
S 2 0 ( t ) = A cos [ 2 π ( f 00 + Δf ) t ] (3)
式中 f 00 = B 2 ,其采样信号为:
S 1 0 ( n ) = A cos [ 2 π ( B 2 - Δf ) n f s ] (4)
S 2 0 ( n ) = A cos [ 2 π ( B 2 + Δf ) n f s ] (5)
取fs=mB,将其带入式(4)和式(5)其中fs满足式(1)的最小要求。
S 1 0 ( n ) = A cos ( 1 m nπ - 2 nπΔf mB ) - - - ( 6 )
S 2 0 ( n ) = A cos ( 1 m nπ + 2 nπΔf mB ) - - - ( 7 )
同理可得位于第k个频带[kB,(k+1)B],k=1,2,...上偏离其中心频率 f 0 k = ( 2 k + 1 ) B 2 分别为±Δf的两个采样信号为:
S 1 k ( n ) = A cos ( 2 m knπ + 1 m nπ - 2 nπΔf mB ) - - - ( 8 )
S 2 k ( n ) = A cos ( 2 m knπ + 1 m nπ + 2 nπΔf mB ) - - - ( 9 )
取m=2的情况,容易得出当k为偶数时:
S 1 0 ( n ) = S 1 k ( n ) , S 2 0 ( n ) = S 2 k ( n ) - - - ( 10 )
而当k为奇数时:
S 1 k ( n ) = A cos ( 3 2 nπ - nπΔf B ) - - - ( 11 )
S 2 k ( n ) = A cos ( 3 2 nπ + nπΔf B ) - - - ( 12 )
根据三角公式,以上二式可以化简为
S 1 k ( n ) = A cos ( 1 2 nπ + nπΔf B ) - - - ( 13 )
S 2 k ( n ) = A cos ( 1 2 nπ - nπΔf B ) - - - ( 14 )
参照式(6)和式(7)则:
S 1 0 ( n ) = S 2 k ( n ) , S 2 0 ( n ) = S 1 k ( n ) - - - ( 15 )
对位于[0,B]基带中偏离其中心频率f00分别为±Δf的两个正弦信号,结果如下:
S 1 0 ( n ) = - S 2 k ( n ) , S 2 0 ( n ) = - S 1 k ( n ) - - - ( 16 )
根据式(15)和式(16)的结果,我们可以得到下面结论:
·在k为奇数时,基带[0,B]中余弦信号的频谱与其在[kB,(k+1)B],k=1,3,...中的频谱是关于中心频率反转的。
·在k为奇数时,基带[0,B]中正弦信号的频谱与其在[kB,(k+1)B],k=1,3,...中的频谱是关于中心频率反转的,同时幅度乘以-1。
·在k为偶数时,基带[0,B]中的频谱与[kB,(k+1)B],k=2,4,...中的频谱是相同的。
根据上述理论分析,选择合适的欠采样频率,可以在基带中实现带通信号的频谱结构反转,然后利用低通滤波器提取出基带信号,最后经过混频得到所需要的中心频率,即完成了信号频谱结构的反转。
发明内容
为了实现该功能,本发明提出了基于欠采样方法的调频信号频谱反转系统的设计方案和原理框图。该方案通过正交解调、欠采样和正交调制方法有效实现了转发调频信号频谱结构的反转,使得转发调频信号具有与输入调频信号相反的调频斜率。
本发明的目的在于设计一种简单、实用的输入调频信号频谱反转的系统。
本发明的技术要点在于首先通过正交解调方法获得位于奇次频带的输入调频信号,然后使用欠采样方法得到位于基带的调频信号。基带调频信号的频谱结构与奇次频带调频信号的频谱结构相反,提取基带调频信号并经过正交调制获得输出调频信号,该输出调频信号与输入调频信号具有相反的调频斜率。
上面通过理论推导和仿真实验证明了欠采样方法能够实现输入调频信号频谱的反转。下面研究欠采样方法在有源编码转发器设计中的具体应用。由于合成孔径雷达发射的是线性调频信号,所以地面有源编码转发器需要通过频谱结构的反转实现输入线性调频信号调频斜率的反转。
为了有效实现输入线性调频信号调频斜率的反转,我们提出了基于欠采样方法的设计方案。为了证明本欠采样方案的正确性,我们采用加拿大Radarsat-1合成孔径雷达的参数进行仿真实验,具体参数如下:
发射信号载波频率:fc=5.3GHz。发射线性调频信号:脉冲宽度为T=20μs,调频带宽为B=30MHz,调频斜率为 K = B T 。有源编码转发器的输入调频信号为:S(t)=cos(2πfct+πKt2)。有源编码转发器中的本振频率:fos=5.03GHz。有源编码转发器中欠采样前的奇次跟踪滤波器的中心频率为: f 0 = ( 2 k + 1 ) 2 B = 285 MHz , k = 9
一种调频信号频谱反转的方法,通过欠采样方法实现输入调频信号基带频谱的搬移,获得偶次频带频谱和奇次频带频谱,其中输入调频信号偶次频带中的频谱与其基带中的频谱关系一致,输入调频信号奇次频带中的频谱与其基带中的频谱关于中心频率反转。
所述的调频信号频谱反转的方法,采用正交解调方法获得两路带通信号,通过欠采样方法,在基带中分别实现两路信号频谱结构的反转。
所述的调频信号频谱反转的方法,采用正交调制方法获得输出调频信号,通过正交调制方法能够获得与输入调频信号相同的中心频率和带宽,但调频斜率相反的输出调频信号。
本发明的优点在于:
1.欠采样方法降低了输入调频信号的采样速率,同时也降低了硬件实现的难度。
2.欠采样方法使宽带高速信号的低速处理成为可能。
3.正交解调方法能够保留完整的输入调频信号频谱。
4.正交调制方法能够保证输出调频信号与输入调频信号具有相同的中心频率。
5.基于欠采样的调频信号频谱反转系统简单、实用。
6.调频信号调频斜率反转系统的稳定性好、精度高。
附图说明
图1是调频信号频谱反转的欠采样系统图。
图2是输基于欠采样方法的有源编码转发器系统图。
图3是输入线性调频信号频谱图。
图4输入线性调频信号的频谱图。
图5是输出线性调频信号频谱图。
图6是输出调频信号的频谱图。
具体实施方式
图1的调频信号频谱反转的欠采样系统,由两条通道组成,每条通道包括:调频信号输入的乘法器、奇次跟踪滤波器、欠采样电路、低通滤波器、D/A转换器、加法器,其中,调频信号输入的乘法器、奇次跟踪滤波器、欠采样电路、低通滤波器、D/A转换器依次顺序连接,D/A转换器通过一D/A转换器连接到加法器,加法器连接到带通滤波器,再由带通滤波器输出调频信号。
图2是具体的设计方案。是图1的具体实施电路,首先接收天线收到的射频信号S(t)(见观测点1)与本振信号fos=5.03GHz进行正交解调获得两路位于奇次频带的中频带通信号:
Si(t)=cos(2πf0kt+πKt2)    (17)
Sq(t)=sin(2πf0kt+πKt2)
式(17)中,f0k=kB是中频带通信号的起始频率,因为k=9,所以式(17)可以进一步简化为:
Si(t)=cos(2πBt+πKt2)    (18)
Sq(t)=sin(2πBt+πKt2)
然后利用采样频率fs=2B进行奇次频带中频带通信号的欠采样,在基带获得反转的调频信号波形为:
Sis(t)=cos(2πBt-πKt2)    (19)
Sqs(t)=-sin(2πBt-πKt2)
上面的两路正交信号与中心频率为f0k的中频参考信号cos(2πf0kt)正交调制获得输出中频带通信号(见观测点2):
So(t)=Sis(t)cos(2πf0kt)+Sqs(t)sin(2πf0kt)
=cos(2πBt-πKt2)cos(2πf0kt)-sin(2πBt-πKt2)sin(2πf0kt)  (20)
=cos(2πf0kt+2πBt-πKt2)
最后该中频带通信号与本振信号fos=5.03GHz上变频得到转发的射频信号:
St(t)=cos(2πfct+2πBt-πKt2)    (21)
此时转发射频调频信号与输入有源编码转发器的调频信号的调频斜率正好相反,这样就实现了输入调频信号调频斜率的反转。
图3是输入线性调频信号的波形,该信号的调频斜率为正,信号频率由慢到快。
图4是输入线性调频信号的频谱结构。
图5是输出线性调频信号,比较图3和图5可以看到该输出信号频率由快到慢,其调频斜率为负与输入调频斜率相反。
图6是输出调频信号的频谱结构,比较图6和图4可以看到它们是关于中心频率反转的。

Claims (5)

1.一种调频信号频谱反转的方法,通过欠采样方法实现输入调频信号基带频谱的搬移,获得偶次频带频谱和奇次频带频谱,其中输入调频信号偶次频带中的频谱与其基带中的频谱关系一致,输入调频信号奇次频带中的频谱与其基带中的频谱关于中心频率反转。
2.根据权利要求1所述的调频信号频谱反转的方法,其特征在于,采用正交解调方法获得两路带通信号,通过欠采样方法,在基带中分别实现两路信号频谱结构的反转。
3.根据权利要求1所述的调频信号频谱反转的方法,其特征在于,采用正交调制方法获得输出调频信号,通过正交调制方法能够获得与输入调频信号相同的中心频率和带宽,但调频斜率相反的输出调频信号。
4.根据权利要求2所述的调频信号频谱反转的方法,正交解调的过程是,首先接收天线收到的射频信号S(t)与本振信号fos进行正交解调获得两路位于奇次频带的中频带通信号:
Si(t)=cos(2πf0kt+πKt2)                           (17)
Sq(t)=sin(2πf0kt+πKt2)
式(17)中,f0k=kB是中频带通信号的起始频率,因为k=9,所以式(17)可以进一步简化为:
Si(t)=cos(2πBt+πKt2)                             (18)
Sq(t)=sin(2πBt+πKt2)
然后利用采样频率fs=2B进行奇次频带中频带通信号的欠采样,在基带获得反转的调频信号波形为:
Sis(t)=cos(2πBt-πKt2)                            (19)
Sqs(t)=sin(2πBt-πKt2)
上面的两路正交信号与中心频率为f0k的中频参考信号cos(2πf0kt)正交调制获得输出中频带通信号:
So(t)=Sis(t)cos(2πf0kt)+Sqs(t)sin(2πf0kt)
=cos(2πBt-πKt2)cos(2πf0kt)-sin(2πBt-πKt2)sin(2πf0kt) (20)
=cos(2πf0kt+2πBt-πKt2)
最后该中频带通信号与本振信号fos上变频得到转发的射频信号:
St(t)=cos(2πfct+2πBt-πKt2)                     (21)
此时转发射频调频信号与输入有源编码转发器的调频信号的调频斜率正好相反,这样就实现了输入调频信号调频斜率的反转。
5.一种调频信号频谱反转的系统,由两条通道组成,每条通道包括:调频信号输入的乘法器、奇次跟踪滤波器、欠采样电路、低通滤波器、D/A转换器、加法器,其中,调频信号输入的乘法器、奇次跟踪滤波器、欠采样电路、低通滤波器、D/A转换器依次顺序连接,D/A转换器通过一D/A转换器连接到加法器,加法器连接到带通滤波器,再由带通滤波器输出调频信号。
CN 200510086970 2005-11-24 2005-11-24 调频信号频谱反转的系统及方法 Expired - Fee Related CN100487483C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200510086970 CN100487483C (zh) 2005-11-24 2005-11-24 调频信号频谱反转的系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200510086970 CN100487483C (zh) 2005-11-24 2005-11-24 调频信号频谱反转的系统及方法

Publications (2)

Publication Number Publication Date
CN1971306A true CN1971306A (zh) 2007-05-30
CN100487483C CN100487483C (zh) 2009-05-13

Family

ID=38112219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200510086970 Expired - Fee Related CN100487483C (zh) 2005-11-24 2005-11-24 调频信号频谱反转的系统及方法

Country Status (1)

Country Link
CN (1) CN100487483C (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101793952A (zh) * 2010-03-22 2010-08-04 南京航空航天大学 一种基于欠采样的电子侦察宽带数字接收信号处理方法
CN102571658A (zh) * 2012-01-16 2012-07-11 中国船舶重工集团公司第七二三研究所 一种高中频欠采样条件的相位编码快速解调方法
CN107831476A (zh) * 2017-10-13 2018-03-23 中国科学院电子学研究所 雷达信号降采样率的装置及方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101793952A (zh) * 2010-03-22 2010-08-04 南京航空航天大学 一种基于欠采样的电子侦察宽带数字接收信号处理方法
CN101793952B (zh) * 2010-03-22 2012-09-05 南京航空航天大学 一种基于欠采样的电子侦察宽带数字接收信号处理方法
CN102571658A (zh) * 2012-01-16 2012-07-11 中国船舶重工集团公司第七二三研究所 一种高中频欠采样条件的相位编码快速解调方法
CN107831476A (zh) * 2017-10-13 2018-03-23 中国科学院电子学研究所 雷达信号降采样率的装置及方法

Also Published As

Publication number Publication date
CN100487483C (zh) 2009-05-13

Similar Documents

Publication Publication Date Title
CN102325058B (zh) 一种变频系统群时延测试方法
CN103728593B (zh) 一种实现地波超视距雷达同时多频发射/接收的方法
CN107064883A (zh) 一种基于压缩感知的线性调频雷达信号处理方法
CN100575984C (zh) 一种使用混沌信号的合成孔径雷达系统
CN105242243A (zh) 基于去斜处理和两次延时的宽带接收数字波束形成方法
CN102707263A (zh) 一种多频多基地高频地波雷达系统及其操作方法
CN202166734U (zh) 一种宽带调频连续波雷达收发子系统
CN108008369A (zh) 一种非线性调频信号欠采样处理方法
CN101825698B (zh) 微波宽频段多极化单抛物面天线散射测量系统
CN103399301B (zh) 一种宽带sar信号的接收装置及接收方法
CN103209147A (zh) 一种多频段毫米波接收机及方法
CN102751998A (zh) 一种基于软件无线电接收机的数字中频模块
CN100487483C (zh) 调频信号频谱反转的系统及方法
CN101672914B (zh) 一种圆锥扫描高分辨率微波散射计的信号处理方法
CN201555948U (zh) 子带合成实现超大带宽信号产生与接收系统
CN103529429B (zh) 用于线性调频信号的频域中的脉冲压缩方法
CN105245262A (zh) 一种俯仰向数字波束合成的实时处理方法和装置
CN102694609B (zh) 一种卫星无线电测定业务通道零值的标定方法
CN104181521A (zh) 一种发射多频载波的高距离分辨率雷达
CN102497225B (zh) 一种适用于反向天线阵的求取数字化共轭相位值的装置
Zongbo et al. Design and application of DRFM system based on digital channelized receiver
CN206773180U (zh) 宽带雷达目标回波信号模拟设备
CN105242242A (zh) 一种基于参数拟合的超大带宽信号预失真补偿方法
Li et al. A novel UHF radar system design for river dynamics monitoring
CN102104564B (zh) 一种通用超宽带接收信道化处理方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090513

Termination date: 20111124