《基于Matlab测速仪设计(多普勒效应回波信号处理)》由会员分享，可在线阅读，更多相关《基于Matlab测速仪设计(多普勒效应回波信号处理)（14页珍藏版）》请在金锄头文库上搜索。
1、1 “数字信号处理课程设计”任务书 题目题目 5基于 Matlab 测速仪设计(多普勒效应回波信号处理) 主要主要 内容内容 利用离散付里叶变换进行频谱分析，设计数字滤波器组，并对滤波器组加窗 降低副瓣，通过频域方法测量速度。 设计设计 要求要求 1.产生输入回波信号及加性噪声信号。 2.计算回波频率偏移于目标速度的对应关系。 3.按技术指标要求对输入信号进行傅立叶变换，对输入信号进行频域分析。 4.写出所设计的窄带滤波器组的系统函数并显示出频率响应曲线。 5.设计适当的加权窗函数，写出设计过程，画出加权后滤波器频率响应曲线 6.正确显示信号频谱和对应的滤波器响应间的关系。 7.采用门限等方法
2、确定并显示目标速度。 主要主要 仪器仪器 设备设备 1、计算机 1 台，安装 MATLAB 软件 主要主要 参考参考 文献文献 美数字信号处理使用MATLABM.西安：西安交通大学出版社，2002. 课程设计进度计划（起止时间、工作内容）课程设计进度计划（起止时间、工作内容） 本课程设计共安排 7 个题目，这是其中题目之一。整个课程设计共 45 学时，分 1.5 周安排， 具体进度如下： 8 学时 复习题目相关知识，掌握实现的原理； 20 学时 用 MATLAB 语言实现题目要求； 9 学时 进一步完善功能，现场检查、答辩； 8 学时 完成课程设计报告。 课程设计开始日期课程设计开始日期 20
3、10.12.6 课程设计完成日期课程设计完成日期 2010.12.15 课程设计实验室名称课程设计实验室名称信号与信息处理实验室地地 点点实验楼 3-503 2 资料下载地址资料下载地址http:/59.64.74.111/实践环节/数字信号处理课程设计 目录目录 一、一、概述概述44 1、目的4 二、二、相关知识简介相关知识简介44 1、背景简介4 2、多普勒频移4 3、 程序中的相关计算5 4、 测速仪的技术规格6 5、 设计内容6 6、 测量方法6 三、三、设计过程设计过程77 1、设计步骤7 2、 设计流程图 3、 程序7 四、四、实验结果及分析实验结果及分析1212 1、实验结果12
4、 五、五、总结总结1515 参考文献资料参考文献资料1616 3 一、一、 概述概述 1 1、 目的：目的： 1、熟悉 DFT 的基本性质，利用离散付里叶变换对信号和系统进行频域分析； 2、熟悉并加深对采样定理的理解； 3、设计数字滤波器组，实现对存在加性干扰的时域离散信号的滤波； 4、对滤波器组进行加权设计，降低副瓣，通过频域方法测量速度。 二、二、相关知识简介相关知识简介 1 1、背景简介：背景简介： 车辆监测用测速雷达是依据多普勒效应原理实现的。由雷达发出一束微波，遇被测目标 （车辆）时微波被反射回来，再由雷达接收反射波。如果目标相对于雷达有距离上的运动， 反射波的频率将与发射频率发生差
5、异，而这种差异的大小正比于目标与雷达的相对运动速 度。于是，雷达通过检测反射波频率和发射波频率的差，就可以计算出被测车辆的移动速 度。 2 2、多普勒频移：多普勒频移： 多普勒频移，当运动在波源后面时，会产生相反的效应。波长变得较长，频率变得较低 （红移 red shift） 。波源的速度越高，所产生的效应越大。根据光波红（蓝）移的程度， 可以计算出波源循着观测方向运动的速度。把声波视为有规律间隔发射的脉冲，可以想象 若你每走一步，便发射了一个脉冲，那么在你每走一步时，面前的声源发出的脉冲相对于 你的传播距离比你站立不动时近了一步，而在你后面的声源则比原来不动时远了一步。或 者说，在你之前的脉
6、冲频率比平常变高，而在你之后的脉冲频率比平常变低了。 所谓多普勒效应就是当发射源与接收体之间存在相对运动时，接收体接收的发射源发射信 息的频率与发射源发射信息频率不相同，这种现象称为多普勒效应，接收频率与发射频率 之差称为多普勒频移。声音的传播也存在多普勒效应，当声源与接收体之间有相对运动时， 接收体接收的声波频率 f与声源频率 f 存在多普勒频移 f(doppler shift)即 f=f-f 当接收体与声源相互靠近时，接收频率 f大于发射频率 f 即：f0 当接收体与声源相互远离时，接收频率 f小于发射频率 即： f0 可以证明若接收体与声源相互靠近或相互远离的速度为 v，声速为 c，则接
7、收体接收声波的 多普勒频率为： f= f(c+-v1)/(c-+v2) 括号中分子和分母的加、减运算分别为“接近”和“远离”之意。 发射信号表示为： 回波信号： 4 回波与发射信号之间有固定相位差： 当目标与雷达之间有相对运动时，则距离 R 随时间变化： 由于通常雷达和目标间的相对运动速度 v，远小于电磁波速度 c，故时延： 相位差： 频差（多普勒频移）： 多普勒测速系统是利用多普勒效应测定飞行器径向速度的无线电跟踪测量系统。由地面向 飞行器或由飞行器向地面发射固定频率的等幅电磁波，因飞行器与测控站存在相对运动， 接收信号的频率与发射信号的频率互不相同。其频率之差就是多普勒频移。多普勒频移正
8、比于测控站航天器方向上的速度分量，所以测出多普勒频移的大小，就可获得目标对测 控站的相对径向速度。 3 3、 程序中的相关计算程序中的相关计算 测速范围： 精度： 70km/h - 150km/h 5km/h 换算： 19.4m/s - 41.6m/s 1.387m/s fd: 606Hz - 1300Hz 43.3Hz wp=2*pi*f/Fs=0.233*pi ws=2*pi*f/Fs=0.249*pi Bt=2*pi*f/Fs=0.016*pi=ws-wp Fs=45*fd 4 4、 5 带通滤波器的设计原理带通滤波器的设计原理 一个理想的滤波器应该有一个完全平坦的通带，例如在通带内没有
9、增益或者衰 减，并且在通带之外所有频率都被完全衰减掉，另外，通带外的转换在极小的频率范围 完成。实际上，并不存在理想的带通滤波器。滤波器并不能够将期望频率范围外的所有 频率完全衰减掉，尤其是在所要的通带外还有一个被衰减但是没有被隔离的范围。这通 常称为滤波器的滚降现象，并且使用每十倍频的衰减幅度dB 来表示。通常，滤波器 的设计尽量保证滚降范围越窄越好，这样滤波器的性能就与设计更加接近。然而，随着 滚降范围越来越小，通带就变得不再平坦 开始出现“波纹” 。这种现象在通带的边 缘处尤其明显，这种效应称为吉布斯现象。 除了电子学和信号处理领域之外，带通滤波器应用的一个例子是在大气科学领 域，很常见
10、的例子是使用带通滤波器过滤最近3 到 10 天时间范围内的天气数据，这 样在数据域中就只保留了作为扰动的气旋。 在频带较低的剪切频率 f1 和较高的剪切频率 f2 之间是共振频率，这里滤波器的 增益最大，滤波器的带宽就是f2 和 f1 之间的差值。 布莱克曼窗计算公式：布莱克曼窗计算公式： 增加一个二次谐波余弦分量，可进一步降低旁瓣，但主瓣宽度进一步增加，增 加 N 可减少过渡带。 频谱的幅度函数为： +0.04 5 5、设计内容：设计内容： 基于 Matlab 设计车辆测速仪的信号处理部分，要求测速范围：18255Km/h，测速精度： 优于 1Km/h，响应时间：小于 30ms。 6 6、测
11、量方法：测量方法： a)测量方法大致有两种方法:时域法和频域法。 b)因环境影响，接收到的多普勒信号的“背景”十分复杂，信噪比大大降低，采 用传统的时域处理方法对被淹没在干扰和噪声中的多普勒信号检出或识别往 6 往是困难的，使得测频精度明显下降。 c)采用频域谱分析方法，选择合适的采样频率及适当的窗口，可以大大提高测 频精度和可靠性。 三、三、设计过程：设计过程： 1 1、 设计步骤：设计步骤： 1.按技术指标要求确定抽样频率，分析车辆反射波的信号形式、载波、频偏范围，仿真 车辆反射波的抽样信号（线性调频信号） 。 2.按照对实际背景噪声分析，仿真产生加性随机噪声输入信号。 3.分析输入信号的
12、频谱，根据载波频率和车辆速度范围要求，分析多普勒频率与车速的 对应关系，确定滤波器系统的处理带宽。 4.选择合适的滤波器，按照技术指标及测速精度要求和测量范围要求，分析分辨率与滤 波器之间的关系，设计和实现滤波器组，显示出窄带滤波器组的频率响应图形。 5.选用合适的窗函数加权，以降低副瓣，实用中，为避免副瓣淹没附近较小的车辆回波 的主信号，引起漏测，必须压低副瓣电平。加权处理可以抑制副瓣电平。加权处理可 以在时域或频域进行。时域加权会降低发射机的利用率，且特定形状的波形调制困难。 一般在频域进行，频域加权处理不必增加实时处理的运算量。实际采用的加权函数包 括汉明、汉宁、余弦平方、余弦立方等。从
13、信噪比损失、主瓣宽度和旁瓣电平方面分 析选择较优的加权函数，显示出加权后滤波器响应曲线。 6.采用门限等方法判明信号所出滤波器位置，确定目标速度。 7.利用 Matlab 将测速结果仿真显示出来。 2 2、设计流程图设计流程图 7 3 3、 程序：程序： clc;clear all; fs=5213;N=1200; n=0:N-1;t=n/fs; v1=75; v2=(v1*1000)./3600; fd=(2*v2)./(6.4*0.01); y=sin(2*pi*fd*t)/800+randn(1,N)/N; B=2*pi*100/fs;%一个过渡带 wc1=2*pi*584.55/fs;
14、wc2=2*pi*627.95/fs; figure(1) plot(t,y); 8 M=ceil(11*pi/B); a=(M-1)/2; k=0:1:M-1; m=k-a; ep=0.0000001; hd1=sin(wc1*m+ep)./(pi*m+ep); hd2=sin(wc2*m+ep)./(pi*m+ep); hd=hd2-hd1; Ha=abs(fft(hd,1024); f=0:1023*fs/1024; figure(2) plot(f,20*log10(Ha/max(Ha);grid;title(a) 带通滤波器幅频特性) axis(0,fs/3,-60,20); xlabel(f/Hz);ylabel(幅度) %*
