CN106019249B - 一种基于三角级数拟合的通道均衡滤波器设计方法 - Google Patents
一种基于三角级数拟合的通道均衡滤波器设计方法 Download PDFInfo
- Publication number
- CN106019249B CN106019249B CN201610573013.7A CN201610573013A CN106019249B CN 106019249 B CN106019249 B CN 106019249B CN 201610573013 A CN201610573013 A CN 201610573013A CN 106019249 B CN106019249 B CN 106019249B
- Authority
- CN
- China
- Prior art keywords
- sequence
- channel
- frequency spectrum
- signal
- fitted
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于三角级数拟合的通道均衡滤波器设计方法,通过对各单元特显点回波数据采集,在频域与基准理想线性调频信号频谱相除,然后填充频谱边缘,再通过三角级数拟合获得均衡滤波器的频谱特性,然后经过IFFT获得均衡滤波器的时域单位冲击响应。相对于传统频域均衡方法,此方法能够有效降低频谱边缘的起伏,使得均衡滤波器逼近理想幅相特性,大大提高宽带数字阵列雷达空域及时域性能。
Description
技术领域
本发明属于宽带数字阵列雷达信号处理领域,具体涉及一种基于三角级数拟合的宽带雷达通道均衡滤波器设计方法,可以广泛用于宽带数字阵列、宽带SAR/ISAR等产品。
背景技术
宽带数字阵列雷达相对传统模拟相控阵雷达来讲,具有动态范围大、精度高、多功能、多波束、低副瓣、抗干扰能力强等优势,但是发射及接收阵列通道的模拟器件存在幅相特性不一致性,即通道失配,使得宽带数字阵列雷达副瓣抬高、测角精度变差、自适应抗干扰能力下降、雷达成像分辨力下降等问题潜在优势难以发挥。
针对此问题,通常采用自适应通道均衡技术补偿通道失配,即设计通道均衡滤波器,补偿系统工作带宽内的幅度相位波动。
目前设计方法有时域算法与频域算法。时域算法是基于维纳滤波理论与最小二乘理论方法;频域算法首先将参考通道与工作通道变换到频域,然后求取通道比值,通过最小二乘拟合求取通道均衡滤波器。
传统时域算法受采样率、均衡器阶数、信噪比的影响非常大;传统频域算法总体上性能优于时域算法,但是通道频率响应两端的边缘频带区域失配比中心频带区域严重,使得性能相对理论有所下降,同时传统算法较为依赖带宽延迟积,对系统参数设计较为苛刻。
本发明总体上属于频域算法,在充分考虑发射及接收通道幅相特性的基础上,首先录取特显点目标回波数据,通过三角级数拟合方法优化均衡滤波器频谱的边缘频带区,消除了边缘频带的振荡效应,降低了系统带宽延迟积设计要求,使得实际均衡滤波器性能逼近理论均衡滤波器特性,大大提高宽带数字阵列雷达空域及时域性能指标。
发明内容
要解决的技术问题
为了避免传统通道均衡滤波器设计过程中存在的频谱边缘起伏剧烈、对带宽时延积要求苛刻等问题,本发明提出一种基于三角级数拟合的通道均衡滤波器设计方法,以满足宽带数字阵列体制成像质量的要求。
技术方案
实现本发明的技术思路是:首先,在雷达的远场放置一个较大RCS的三角形角反,采集信噪比较大的基带回波,完成系统收发通道幅相特性采集;其次将基带时域回波与理想线性调频信号变换到频域,获得频域商,求取频域商的包络与相位,再对频谱边缘填充,然后分别对包络与相位进行三角级数拟合,最后获的均衡滤波器的幅频序列,然后进行IFFT,即可获得均衡滤波器的时域响应。
一种基于三角级数拟合的通道均衡滤波器设计方法,其特征在于步骤如下:
步骤1:信号与场景设置
1a)设计线性调频信号:信号调频带宽B和雷达系统工作带宽一致;信号时宽T大于100,线性调频信号为:
其中n=0,1,2...N-1,N为采样点的数量,fs为采样率,要求fs≥B;
1b)设置场景:在弱杂波区或无杂波区放置一个点目标,目标距离大于距离盲区c为光速,通过选择RCS合适的点目标,使得回波信杂噪比不小于20dB;
步骤2:回波采集
采集雷达各通道基带回波数据sm(n),sm(n)要包含目标回波区域,m为通道号,m=0,1,2…M-1,n为采样点,n=0,1,2,…N-1;
步骤3:计算频谱:对基带回波数据sm(n)做N点FFT,获得频谱Sm(n);对线性调频信号sLFM(n)做N点FFT,获得频谱SLFM(n);
步骤4:频域相除:求第一个通道频域商:将频域商表示为:H1(n)=|H1(n)|·exp(∠(H1(n)));其中,|H1(n)|为幅度序列,∠(H1(n))为相位序列;
步骤5:频谱填充:
频谱支撑区为:
步骤6:傅立叶级数拟合
6a)对幅度序列|H1(n)|进行三角级数拟合;
即获得新的幅度序列ai为余弦级数系数,bi为正弦级数系数,x取值为0,1,2,...N-1,w为角频率,L为级数;
6b)求新序列的端点差,即如果小于等于第一门限时,拟合停止,然后转入6c),否则继续拟合,对新序列进行FFTSHIFT,即转入6a);
6c)获得拟合后的幅度序列
6d)对相位序列∠H进行三角级数拟合;
即获得新的相位序列ai为余弦级数系数,bi为正弦级数系数,x取值为0,1,2,...N-1,w为角频率,L为级数;
6e)求新序列∠H的端点差,即|∠H(N-1)-∠H(0)|,如果|∠H(N-1)-∠H(0)|小于等于第二门限时,拟合停止,转入6f),否则继续拟合,对新序列进行FFTSHIFT,即转入6d);
6f)获得拟合后的相位序列∠H;
6g)获得新的复数序列:对此序列进行IFFT,获得均衡滤波器时域序列h(n);
6h)对h(n)进行FFTSHIFT,求取能量95%以上的中间区域,作为均衡滤波器系数;
步骤7:转入步骤4~步骤6,求取其它通道的均衡滤波器系数。
所述的点目标为三角形角反。
步骤6b)中所述的第一门限为δ为10-2。
步骤6e)中所述的第二门限为δ为10-2。
有益效果
本发明提出的一种基于三角级数拟合的均衡滤波器设计方法,由于充分利用了三角级数的特点,对实测数据通过多次反复拟合获得较为理想的均衡滤波器频谱特性,同时结合特显点实测回波获得的收发通道带内幅相波动特性,使得均衡后的通道特性逼近理想通道特性,提高了宽带数字阵列雷达空域与时域性能。
附图说明
图1是本发明的实现流程图;
图2是频谱填充示意图;
具体实施方式
现结合实施例、附图对本发明作进一步描述:
图1是本发明的实现流程图,本发明的实现步骤如下:
1)信号与场景设置
1a)设计线性调频信号:信号调频带宽B和系统工作带宽一致;信号时宽T尽量宽,时宽带宽积至少大于100,距离盲区为c为光速,线性调频信号为:
其中n=0,1,2...N-1,fs为采样率,要求fs≥B。
1b)设置场景:在弱杂波区或无杂波区放置一个强点目标,比如三角形角反,目标距离大于距离盲区r,回波信杂噪比不小于20dB;
2)回波采集
采集各通道基带回波数据sm(n),sm(n)要包含目标回波区域,m为通道号,m=0,1,2…M-1,n为采样点,n=0,1,2,…N-1;
3)计算频谱:对基带回波sm(n)做N点FFT,获得频谱Sm(n);对参考信号sLFM(n)求其频谱SLFM(n);
4)频域相除:求第一个通道频域商:将频域商表示为:H1(n)=|H1(n)|·exp(∠(H1(n)));
5)频谱填充:
频谱支撑区为:
支撑区数据长度为M,支撑区平均值为:对频谱边缘进行支撑区平均值填充,即令
6)傅立叶级数拟合
6a)对幅度序列|H1(n)|进行三角级数拟合;
即获得新的幅度序列ai为余弦级数系数,bi为正弦级数系数,x取值为0,1,2,...N-1,w为角频率,L为级数;
6b)求新序列的端点差,即如果
时,拟合停止,δ可以根据需要选择,通常为10-2或更小,然后转入6c),否则继续拟合,对新序列进行FFTSHIFT,转入6a);
6c)获得拟合后的幅度序列
6d)对相位序列∠H进行三角级数拟合;
即获得新的相位序列ai为余弦级数系数,bi为正弦级数系数,x取值为0,1,2,...N-1,w为角频率,L为级数;
6e)求新序列∠H的端点差,即|∠H(N-1)-∠H(0)|,如果
时,δ可以根据需要选择,通常为10-2或更小,拟合停止,转入6f),否则继续拟合,对新序列进行FFTSHIFT,转入6d);
6f)获得拟合后的相位序列∠H;
6g)获得新的复数序列:对此序列进行IFFT,获得均衡滤波器时域序列h(n);
6h)对h(n)进行FFTSHIFT,求取能量95%以上的中间区域,作为均衡滤波器系数;
步骤7:转入步骤4~步骤6,求取其它通道的均衡滤波器系数。
本发明提出了一种基于三角级数拟合的通道均衡滤波器设计方法,通过对各单元特显点回波数据采集,在频域与基准理想线性调频信号频谱相除,然后填充频谱边缘,再通过三角级数拟合获得均衡滤波器的频谱特性,然后经过IFFT获得均衡滤波器的时域单位冲击响应。相对于传统频域均衡方法,此方法能够有效降低频谱边缘的起伏,使得均衡滤波器逼近理想幅相特性,大大提高宽带数字阵列雷达空域及时域性能。
Claims (2)
1.一种基于三角级数拟合的通道均衡滤波器设计方法,其特征在于步骤如下:
步骤1:信号与场景设置
1a)设计线性调频信号:信号调频带宽B和雷达系统工作带宽一致;信号时宽T大于100,线性调频信号为:
其中n=0,1,2...N-1,N为采样点的数量,fs为采样率,要求fs≥B;
1b)设置场景:在弱杂波区或无杂波区放置一个点目标,目标距离大于距离盲区c为光速,通过选择RCS合适的点目标,使得回波信杂噪比不小于20dB;
步骤2:回波采集
采集雷达各通道基带回波数据sm(n),sm(n)要包含目标回波区域,m为通道号,m=0,1,2…M-1,n为采样点,n=0,1,2,…N-1;
步骤3:计算频谱:对基带回波数据sm(n)做N点FFT,获得频谱Sm(n);对线性调频信号sLFM(n)做N点FFT,获得频谱SLFM(n);
步骤4:频域相除:求第一个通道频域商:将频域商表示为:H1(n)=|H1(n)|·exp(∠(H1(n)));其中,|H1(n)|为幅度序列,∠(H1(n))为相位序列;
步骤5:频谱填充:
频谱支撑区为:
M为支撑区数据长度,支撑区平均值为:对频谱边缘进行支撑区平均值填充,即令
步骤6:傅立叶级数拟合
6a)对幅度序列|H1(n)|进行三角级数拟合;
即获得新的幅度序列ai为余弦级数系数,bi为正弦级数系数,x取值为0,1,2,...N-1,w为角频率,L为级数;
6b)求新序列的端点差,即如果小于等于第一门限时,拟合停止,然后转入6c),否则继续拟合,对新序列进行FFTSHIFT,即转入6a);
6c)获得拟合后的幅度序列
6d)对相位序列∠H进行三角级数拟合;
即获得新的相位序列ai为余弦级数系数,bi为正弦级数系数,x取值为0,1,2,...N-1,w为角频率,L为级数;
6e)求新序列∠H的端点差,即|∠H(N-1)-∠H(0)|,如果|∠H(N-1)-∠H(0)|小于等于第二门限时,拟合停止,转入6f),否则继续拟合,对新序列进行FFTSHIFT,即转入6d);
6f)获得拟合后的相位序列∠H;
6g)获得新的复数序列:对此序列进行IFFT,获得均衡滤波器时域序列h(n);
6h)对h(n)进行FFTSHIFT,求取能量95%以上的中间区域,作为均衡滤波器系数;
步骤7:转入步骤4~步骤6,求取其它通道的均衡滤波器系数;
步骤6b)中所述的第一门限为δ为10-2;
步骤6e)中所述的第二门限为δ为10-2。
2.根据权利要求1所述的一种基于三角级数拟合的通道均衡滤波器设计方法,其特征在于所述的点目标为三角形角反。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610573013.7A CN106019249B (zh) | 2016-07-20 | 2016-07-20 | 一种基于三角级数拟合的通道均衡滤波器设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610573013.7A CN106019249B (zh) | 2016-07-20 | 2016-07-20 | 一种基于三角级数拟合的通道均衡滤波器设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106019249A CN106019249A (zh) | 2016-10-12 |
CN106019249B true CN106019249B (zh) | 2018-08-31 |
Family
ID=57116659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610573013.7A Active CN106019249B (zh) | 2016-07-20 | 2016-07-20 | 一种基于三角级数拟合的通道均衡滤波器设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106019249B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112147590B (zh) * | 2020-08-25 | 2023-04-07 | 西安电子科技大学 | 一种基于响应估计频域拟合的通道均衡方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1494285A (zh) * | 2002-11-01 | 2004-05-05 | 上海奇普科技有限公司 | 一种自适应均衡器 |
JP4593959B2 (ja) * | 2004-04-09 | 2010-12-08 | ソニー株式会社 | 適応等化装置及び方法 |
JP5096179B2 (ja) * | 2008-01-28 | 2012-12-12 | ルネサスエレクトロニクス株式会社 | 適応等化装置 |
CN101527697B (zh) * | 2009-03-16 | 2012-08-29 | 上海高清数字科技产业有限公司 | 时域自适应均衡器及其均衡方法 |
-
2016
- 2016-07-20 CN CN201610573013.7A patent/CN106019249B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106019249A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105259552B (zh) | 一种基于非线性调频信号的合成孔径雷达成像方法和装置 | |
CN108387900B (zh) | 一种直升机载旋转式合成孔径雷达的振动误差补偿方法 | |
US20130002488A1 (en) | Wideband beam forming device; wideband beam steering device and corresponding methods | |
US8760340B2 (en) | Processing radar return signals to detect targets | |
CN104777467B (zh) | 基于频率扫描天线的目标检测方法 | |
CN109581362B (zh) | 合成孔径雷达高度计在可变脉冲簇模式下的信号处理方法 | |
CN111090080A (zh) | 基于空时编码阵列的超宽带雷达单通道数字波束形成方法 | |
CN110879391B (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
Yu et al. | Ground moving target motion parameter estimation using Radon modified Lv's distribution | |
JP4444150B2 (ja) | フィルタ装置 | |
JP4723880B2 (ja) | 電波誘導装置 | |
CN110045338A (zh) | 基于希尔伯特变换的幅相误差估计与校正方法 | |
CN106019249B (zh) | 一种基于三角级数拟合的通道均衡滤波器设计方法 | |
CN109143235B (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
JP6323156B2 (ja) | 画像レーダ装置 | |
CN110531354B (zh) | 一种频控扫描雷达色散信号的二维成像方法 | |
CN110308445A (zh) | 一种基于车载数字阵调频连续波雷达的成像方法 | |
CN111458677A (zh) | 一种双通道单脉冲比幅测角方法、装置 | |
CN109343056A (zh) | 一种针对非线性调频sar的rd成像方法和装置 | |
CN109633641B (zh) | 一种太赫兹频段的旋翼叶片逆合成孔径雷达成像算法 | |
Wang et al. | A novel signal processing algorithm for staggered SAR with low oversampling factors | |
CN106093907A (zh) | 一种基于特显点的宽带lfmcw体制雷达通道均衡方法 | |
Tran et al. | An experimental study of radar tomographic imaging in a multi-bistatic scenario | |
Jin et al. | Optimal mismatched filter design by combining convex optimization with circular algorithm | |
Izadi et al. | Weighted Capon beamformer combined with coded excitation in ultrasound imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |