无线电定位原理与技术实验报告姓名:学号:班级:1105201****:**院系:电子与信息工程学院哈工大电子与信息工程学院电子工程系实验一 连续波雷达测速实验1、1 雷达测速原理雷达利用多普勒频率来提取目标的径向速度(即距离变化率),从而可以区分运动目标与固定目标及杂波。
多普勒效应描述了由于目标相对于辐射源的运动而引起发射信号的中心频率发生多普勒频移,目标的运动方向的不同决定了多普勒频移的正负。
00022d r vf v f f f f c v c --⎛⎫=-=≈ ⎪+⎝⎭(如果v c <<)雷达雷达λ>λλ'<图2- 多普勒效应1、2 连续波雷达测速实验仪器连续波雷达测速系统主要由三部分组成:微波发射与接收器件,差频放大与滤波电路,DSP 信号采集与处理电路。
其中微波发射与接收器件可以采用微波发射介质稳频振荡与微波接收混频器。
放大与滤波电路,在近距离时,测量直接由混频器输出的信号较大,由雷达方程可知,随着目标距离的增加,混频器输出会减小。
实验中采用三级放大电路,第一级射随阻抗匹配,第二三级可调增益放大。
其次由于背景噪声与扰动会引入杂波,对接收信号需要进行滤波。
DSP 信号采集与处理电路,采集多组回波数据,对数据进行分析得到相应的多普勒频率与速度值,由公式2rd v f λ=,算得速度r v 。
5402DSP 测速传感器混频器连续波发射机传感器输出信号放大滤波AD 串行接口PC 机FFT图2-2 连续波雷达测速实验仪器原理框图图2-3 测速雷达传感器1、3 实验要求本实验为演示实验,观察实验现象,并在PC 机使用Matlab 对实验数据进行分析。
实验要求:1、掌握雷达测速原理,2、了解连续波雷达测速实验仪器原理及使用,3、使用Matlab 对实验数据进行分析,得到回波多普勒频率与目标速度。
1、4 实验内容1、采集三组数据,每组数据 2048 点,采样频率为 2048Hz2、从每组数据中分别选取波形较好的 512 点,作出时域波形与频谱,并求出目标速度,其中,发射波频率为 10GHz 。
1、5实验结果分析首先将实验中TXT文件录取的数据读入到MATLAB中去,对其进行FFT变换并分析频域。
在对数据的处理过程中,我注意到数据的平均值不为零,这个说明存在着直流分量。
因而采用MATLAB强制地将零频出的值变为0,以免对实验的数据产生影响。
对于数据一,实验的截图如下所示。
上图为时域的信号,相对而言比较接近正弦波,采样点数N=434;下图为频域分析图,由于在[0,2π]之间存在两个峰值,因此通过 axis函数控制只显示了正半频域部分。
同时,为了精确求出幅值最大点,我在文件中利用 matlab 逻辑运算操作,执行以下命令由上可知峰值对应出的N值nd=10因为我在实验中采用的采样频率为2048Hz,发射波的频率为10GHz,波长为3米通过计算可知:多普勒频移 = 47、1889 Hz目标移动速度= 0、7078 m/s对于数据二,实验的截图如下所示。
上图为时域的信号,相对而言比较接近正弦波,采样点数N=396;下图为频域分析图,由于在[0,2π]之间存在两个峰值,因此通过 axis函数控制只显示了正半频域部分。
同时,为了精确求出幅值最大点,我在文件中利用 matlab 逻辑运算操作,执行以下命令由上可知峰值对应出的N值nd=12因为我在实验中采用的采样频率为2048Hz,发射波的频率为10GHz,波长为3米通过计算可知:多普勒频移 = 62、0606 Hz目标移动速度= 0、9309 m/s对于数据三,实验的截图如下所示。
上图为时域的信号,相对而言比较接近正弦波,采样点数N=398;下图为频域分析图,由于在[0,2π]之间存在两个峰值,因此通过 axis函数控制只显示了正半频域部分。
同时,为了精确求出幅值最大点,我在文件中利用 matlab 逻辑运算操作,执行以下命令由上可知峰值对应出的N值nd=11因为我在实验中采用的采样频率为2048Hz,发射波的频率为10GHz,波长为0、03米通过计算可知:多普勒频移 = 56、6030 Hz目标移动速度= 0、8490 m/s1、6实验源程序%实验一:连续波雷达测速实验clear all;t=load('C:\Users\dell\Desktop\实验\1110520115、1、txt');fs=2048;%抽样频率2048HzT=fft(t);T(1)=0;%由于所测数据含有直流成分,将其零频处强制置为0N=length(t);n=0:N-1;nd=sum(n(1:N/2)'、*(abs(T(1:N/2))==max(abs(T(1:N/2)))));%找出FFT 后最大值所在的位置,只算[0,pi]之间的正半部分fd=nd*fs/N;%求出多普勒频率v=0、5*fd*0、03;%求出目标速度figure(1);subplot(211);plot(n/fs,t);xlabel('时间/s');ylabel('电压幅值/U');grid on;title('多普勒差回波');subplot(212);plot(n(1:64)*fs/N,abs(T(1:64)));%利于观察,频域图像只显示一部分xlabel('频率/Hz');ylabel('电压/U');grid on;title('频率分析');axis([0 434 0 2300]);2实验二 线调频信号及匹配滤波仿真实验LFM 信号以其优越的频谱性能广泛应用于雷达与众多电子工程中,匹配滤波器在相参滤波分析中也得到广泛的应用。
2、1 线调频信号谱分析线调频(LFM)信号时域表达式:20()cos()2t kt S t Arect()t T ω=+ 式中:t rect()T就是矩形函数,k 就是调频斜率,并且与调制频偏ω∆的关系就是: 2f k T Tωπ∆∆== T 为时域波形宽度,简称时宽;f B ∆=2为调频范围。
简称频宽。
BT D =为时宽带宽积,就是线性调频信号一个很重要的参数。
LFM 信号的频谱近似为:20()[]}()240j S k ωωπω⎧--+⎪=⎨⎪⎩02others ωωω∆-≤ 近似程度取决于时宽带宽积D, D 越大,近似程度越高,即频谱越接近于矩形。
Time in u secFrequency in MHz图2-1 LFM 信号时域频域图(例)2、2 线调频信号匹配滤波 雷达发射LFM 脉冲信号,固定目标的回波时域表示:20()()cos(())2r r i r t t k t t S t Arect()t t T ω--=-+ 对应的匹配滤波器的传输函数近似(大时宽带宽积下)为: 20()()exp{[]}24H j k ωωπω-=- 02ωωω∆-≤ 匹配滤波器输出:0()()exp()d i d j t S S H j t ωωωω-=-= 02ωωω∆-≤ 代入相关参数,002,2,2B k B T f ωππωπ∆===匹配滤波器时域输出:02()1()()2d j t o o i f t t d S t S e d ωπωωπ+∞-∞-==⎰时宽带宽积:D BT =匹配滤波器的包络输出如下图3-2所示,所示,通常规定顶点下降到－4dB 处的宽度为输出脉冲的脉宽0T ,并且有01T B=,所以脉冲压缩比:0T BT D T ==Time in sec ⨯BA m p l i t u d eTime in sec ⨯BA m p l i t u d e ,d BTime in sec ⨯B A m p l i t u d e ,d B图2-2 LFM 通过匹配滤波器的时域图(例)对应的匹配滤波器的传输函数在大时宽带宽积下,如上图3-3所示,与辛格函数拟与很好,在主瓣与临近的几个旁瓣都没有偏差,但就是在小时宽带宽积下,仅在主瓣与辛格函数拟与无偏差,而在旁瓣偏差较大。
2、3 实验要求本实验为仿真实验。
实验要求:1、掌握线调频信号及其频谱特征,2、使用Matlab 对线调频信号及其频谱进行仿真3、掌握匹配滤波理论,4、使用Matlab 线调频信号进行匹配滤波仿真。
5、讨论时宽带宽积对线调频信号频谱与匹配滤波的影响。
2、4 实验内容用 matlab 编写源程序,实现上述实验要求。
2、5 实验结果与分析备注:1、在本实验当中,由于原先采集的数据与实际值相符合地不就是很好,因此在写实验报告的时候我决定放弃原先采集的数据,直接采用MATLAB来进行数据的模拟,实际上也确实模拟地不错。
2、根据信号与系统中所学的傅里叶变换的形式,可知在实验中进行匹配滤波时应该采用正交双通道处理。
在MATLAB中如果直接采用chirp函数产生余弦形式的信号进行匹配滤波,其效果与理想的sinc函数有一定的差距。
2、5、1线性调频信号仿真与分析线性调频信号仿真:我编写了一个chirp_m、m的子函数文件可以产色很难过exp形式的线性调频信号,一共有三组信号,各自的时域图像与频域图像如下所示(假设各自的抽样频率均为100MHz):第一组信号:实验信号x1,扫频周期为8us就是,fs=100MHz,扫频带宽8MHz,时域带宽积为64。
第二组信号:实验信号x2,扫频周期为8us就是,fs=100MHz,扫频带宽14MHz,时域带宽积为112。
第三组信号:实验信号x1,扫频周期为10us就是,fs=100MHz,扫频带宽20MHz,时域带宽积为200。
分析:通过对上面的三幅图的对比与分析可知,LFM信号在时域变化有规律随着时间疏密程度不同,在频域为一接近矩形的窗,并且近似程度与带宽B与扫描周期T无关,只与时宽带宽积D有关。
时宽带宽积越大,其近似程度也就越大。
不过从图中也可以瞧出,随着D的增大,其边沿处的上冲并不会消失,也就就是说Gibbs效应依然存在。
2、5、2匹配滤波仿真与分析匹配滤波仿真:分别对以上三组信号进行匹配滤波,并两两进行比较。
考虑到matlab中直接时域卷积conv计算比较慢,因此我采用的匹配滤波方法就是直接将信号延时作为回波信号,作傅立叶变换后并作共轭,与接收信号的傅立叶变换相乘后,再作傅立叶逆变换。
为了防止频域混叠,做fft变换时要对回波信号与原始信号补零,使之满足L>M+N-1的条件,之后再做ifft。
首先进行的就是第一组与第二组信号的匹配滤波,为了便于观察,在M文件中使用了axis函数控制直系那就是匹配滤波后部分时间段的图像,实验结果的截图如下:分析:通过对第一个图与第二图的对比可以发现,对应的匹配滤波器的传播函数在大时宽带宽积下(第二个图),与sinc函数相拟合地很好,相临近的两个旁瓣的误差并不就是很大,但就是对于小的时宽带宽积来说(第一个图),从第一个旁瓣就开始出现了一定的偏差。