DSP试验大作业精

1、inDSP实验课程序设计报告学院:学号:姓名:指导教师:DSP实验课大作业设计实验目的在DSP上实现线性调频信号的脉冲压缩、动目标显示（MTI ）和动目标检测（MTD，并将结果与MATLAB上的结果进行误差仿真。二实验内容2.1 MATLAB 仿真设定带宽、脉宽、采样率、脉冲重复频率，用 MATLAB产生16个脉冲的LFM，每个脉冲有4个目标（静止，低速，高速），依次做2.1.1脉压2.1.2相邻2脉冲做MTI，产生15个脉冲2.1.3 16个脉冲到齐后，做MTD，输出16个多普勒通道2.2 DSP实现将MATLAB产生的信号，在visual dsp中做脉压，MTI、MTD，并将结果与MATL

2、AB 作比较。三实验原理3.1线性调频线性调频信号是指频率随时间而线性改变（增加或减少的信号。本实验中用的是线 性调频矩形脉冲信号，其表达式可写为：5(0=严 mT其中信号的复包络为Ir/2叶卜r/2T为脉冲宽度，M =型T为调频斜率，B为调频带宽，也称频偏。3.2脉冲压缩原理脉冲压缩处理将发射的宽脉冲信号压缩成窄脉冲信号，它既可以发射宽脉冲以提高平均功率和雷达的检测能力，又能保持窄脉冲的距离分辨率。本实验中是对线性调频脉冲信号的压缩处理，是通过在发射的长脉冲上加上频率调 制，并通过匹配滤波器接收机来压缩接受到的信号，其过程如下图所示：图1 LFM脉冲压缩处理框图（频率域）3.3 MTI （动

3、目标显示）原理MTI是指利用用杂波抑制滤波器来抑制各种杂波，提高雷达信号的信杂比, 以利于运动目标检测的技术。实验中使用的是单延迟对消器，其原理如下图所示:图2单延迟线对消.讹）-川7）3.4 MTD （动目标检测）原理MTD是一种利用多普勒滤波器组来抑制各种杂波，以提高雷达在杂波背景下 检测运动目标能力的技术。四实验步骤带宽以下是该实验中设定的几个参数Ban dWidth=2.0e6TimeWidth=42.0e-6脉宽Fs=2.0e6采样率P RT=240e-6脉冲重复周期TargetDistance= 3000 8025 8025 12000 目标距离TargetVelocity= 50

4、 0 -120 225目标速度假设接收到的回波数是 16个，噪声为高斯随机噪声。4.1在MATLAB产生线性调频信号。4.2根据目标距离得出目标回波在时域的延迟量，根据目标速度得出多普勒相移，从而在MATLAB产生4个目标的16个回波串（接收到的回波含噪声）。4.3由匹配滤波理论产生对应于目标回波的滤波系数（脉压系数），在时域中做线性卷积，实现匹配滤波（时域脉压）；在频域中做回波信号和脉压系数的FFT，点乘后作逆FFT，实现频域脉压。两者进行比较，讨论其差别。4.4对16个去暂态点后的脉冲串按接收顺序进行排列，用次相消器（一种滤波方式）实现 MTI。4.5 做16通道的FFT，实现 MTD4.

5、6在DSP中对MATLAB生的回波数据和脉压系数进行处理，实现频域脉压。导入DSP的回波数据为时域数据，而脉FFT压系数为频域数据。将导入 DSP的时域回波数据进行一次变换到频域，然后将其与频域脉压系数进行点积，得到频域脉压结果。对该结果再做一次IFFT，将频域转换成时域。在这一步中需要调用库函数fft_flp 32.asm。该子程序可实现8192点复数的FFT功能。由于C语言中无法实现复数运算，因此，对8192个复数按照实部虚部交替的顺序进行重排列， 用长度为16384的数组来存放时域回波数据。频域相乘后， 做乘积结果的IFFT，得到脉压结果。做IFFT，仍需调用库函数fft_flp32.a

6、sm，此时要通过 FFT子程序实现IFFT的功能，要对频域的数据进行处理，才能达到这一目的。4.7对脉压后的数据按照脉冲号重排，相邻序列的数据相减（滑动对消），实现 MTI。4.8调用子程序fft 16.asm，做16通道FFT，实现MTD入口参数为16通道的脉压数据。五实验结果及讨论5.1 脉压、MTI、MTD结果分析：5.1.1脉压图2脉压2的情况下依然不能分辨。对于幅度目标1的回波由于雷达闭锁 的影响，幅度较低；目标 2和目标3的回波叠加在一起，且 由于相位的影响，两个回波同相叠加而增强或者反相抵消而 减弱，所以在16个回波中目标2和目标3的回波的和幅度有 时增强，有时减弱；目标 4回波

7、的幅度是一般情况。结果如图所示(修改程序，使横轴表示目标距离，以便分 析)，试验初试是设置的目标距离参数为TargetDistance=3000 8025 8025 ，12000，现在图中只看到三个目标分别是(300, 41.94 )、(8025, 122.7 )、( 12000，84.74)，因为 目标2和目标3的距离参数相同，即使距离参数不同，在小Mr*于雷达的距离分辨力sqrt(Sig Po wer(k= 1,1,0.5,1,5.1.2 MTIMTI18016014012010080604020015Y: 9Z: 163.4Y: 1510'-_.*._X： 412輕_X: 161

8、500图2(a MTI三维图MTD图2(b MTI正面图在实验中做的是单延迟对消 MTI，即相邻两个脉冲做差，最后 只剩下15个脉冲，对于静止目标，两次回波相差不大，互相 抵消，回波2此时只剩下目标3的回波；由于 MTI滤波器幅 频响应特点的影响，目标 1和目标3的回波幅度减小大于目 标4的回波幅度的减小，所以出现上图的幅度现象；另，由 于随机噪声的影响，各次脉冲的幅度有一定的起伏。5.1.3 MTD.-I140012001000X: 107丫:一 0- Z: 1330800 -一600<:40Z: 640.1 一一一X200Z6420,-2-4-6-8图3(a MTD三维图图 3(c

9、MTD Y-Z View对图3中的幅度分析,MTD中的目标的幅度是由不同脉冲相干积累得到，同时由于FFT存在频谱泄漏，主瓣展宽等，所以结果中目标的幅度近似为脉压的16倍。对图3进行多普勒和时延分析，可以看图目标1的距离和速度都已经测出，目标 2和目标3由于速度不同，落在了不同的FFT通道，从而区分出来。在该实验参数的设置下目标4的速度-174.16m/s，这显然是和实际不符的，即速度出现了模糊，究其原 因为在该实验参数下的测速范围是-199199m/s，超出这个范围将会产生 速度模糊（多普勒模糊）。本实验雷达的多普勒分辨率为"八 介如£川对应的径向速度为:其中TotalNu

10、mber=Sam pleNumber* Pu lseNumber=7680 ，在PulseNumber不变的情况下，更改 PRT可以改变SampleNumber，进而改变雷达的测速范围，下面将更改PRT 为 150e-6，此时4800，测速范围为-318318m/s，可满足试验要求。结果如下图:1400X: 160Y: 6Z: 10721200X: 1071000Y: 0Z: 1351800600X: 40_ =Y: 1400Z: 573.8X: 41200Z: 375.420='-2-4图4解多普勒模糊此时目标4的速度V=228.8，存在定的误差，但较原参.-6-8数已有很大的改善。

11、以上方法在仿真实验中可以解决存在的问题，然而在实际 中由于目标的速度是未知的，所以上述方法显然是不可行 的。实际中为了解决多普勒模糊，采用两种脉冲重复频率， 即和来解决多普勒模糊，且选择 <1和4，使得它们相互成质数，一种选择是使 $广必4珂山1）几，在所希望的一个脉冲重复间隔）周期内，两个脉冲重复频率只在某一个位置上重合，此时：如果厶> 切，则匚 ，如果 d厶，贝，且真实多普勒为坳+亿十芒，最后，如果厶5 ,则A S心, 此时，可能发生盲多普勒，可用第三种脉冲重复频率加以解 决。5.2 DSP结果与Matlab结果对比，误差分析。图5频域脉压的相对误差图6 MTI数据的相对误差图

12、7 MTD数据的相对误差从上面的图中可以看出，MATLA审DSP的相对误差均在10-4之内，误差在可 以接受的范围内。这样的误差可能是由于算法和不同字长的舍入造成的误差。六实验思考题及心得6.1 MTI结果中看不到静止目标答：MTI中利用相邻脉冲差值来检测目标，当目标静止时，两次回波的幅度和相位变化不大（杂波的影响），相差后几乎为零，所以MTI结果中看不到静止的目标。6.3 MTD速度/多普勒通道的含义答：不同的多普勒通道对应着不同的速度，目标速度的不同，导致多普勒频移不同，落到不同的多普勒通道。6.2速度模糊如果目标的多普勒频移不超过单个滤波器带宽的一半，即小亍十？,多普勒滤波组就可以分辨出

13、目标的多普勒频移，否则就会产生多普勒模糊。也就是 说，雷达发射脉冲的重复频率PRF不能低于目标最大多普勒频移的2倍，否则雷 达无法分辨目标的多普勒信息练习五、编写C程序，建立磁盘文件，练习读写数据实验目的在VisualDSP4.5中学习并掌握对磁盘文件的读写，并对数据进行简单的运 算。实验内容2.1编写C程序，产生1000个随机浮点数（随机类型和参数自行确定），建立一个磁盘文件，把所有随机数依序、按照相应格式写入文件。2.2关闭文件，再打开此文件，从文件中把数据读入，计算其均值、方差，并 写入另一个文件中。三实验原理3.1随机数的产生利用rand（产生所需要的数据;运用格式化文件输入、输出函数

14、:3.2磁盘文件的读写int fprin tf(FILE *stream, const char *format, /*args*/ .;int fsca nf(FILE *stream, const char *format, /* args */.;3.3均值、方差叫工t ；1 、 rD广U站-叫y；N J I四实验步骤4.1熟悉实验原理4.2编写C程序4.3仿真、调试程序，直至程序满足实验要求五实验结果及讨论5.1产生的随机数据和从外部文件中调入所产生的数据Al ''J r' - 1K'lbir. -iri宀* p/-hk L “ . 1*"诵讯

15、rji汕心E 平aI!和一益"ss. - +G1 " 1 r ' J：*1 *iH 八：u-L _ H Lllil r- r *h£ . ' "-;mrGhqVJHi:中1 b I b dl 1 I1 1 > ： 4ii fZ »i r ,'F/ j-r*!： T，:r ' 1n I r rfTarvi.r：和JII z irEri 1"W八："-' " gi _，山1Hi * L吨!i； T：1 lW".F!：Jl! -<?",. . i 订

16、 J汕I 1«!： 4L-.，-.朮如k-.p11l、S .Vie- c El b h- _h ”h u p aj 卜 u. 4臣雷於血于 无-T 一 a f - r J t L - b n I- L- ' h九咗门注血产rg=L I . ' L3I FT. ur :-i- ir-I J* 區 6心1叶I iit->KaFUJI4Cl»*i'+lSpHF"-IT *1MlEr- in-d .LbiC'-'.-a cE-i c n.-iwI a1d<ahhl-. M.J*：-HI.lE!il-T-b-：.I / b

17、 ri < lid I n n- - T ; L < ' 4- fl ' i r b- L - - - B r r ht7hL” F' J hin"u5.e F+ 1-T4Fprrp1 c Jr 'Jf 77c 吋“ilWlZ l 'hr 吗J iiTEr 71.>/ L P H E|.b LITtr5.2均值和方差图1数据对比图2均值和方差六实验思考题及心得6.1此练习中，计算机和DSP#自起什么作用？如果DSP脱离计算机，此程序还能运行吗?答：在此练习中，计算机提供文件存储的磁盘空间，DSP向计算机磁盘写入或读出数据，并进

18、行运算；在此过程中 DSP需要借助计算机的磁盘完成功能，因此,DSP脱离计算机后该程序将无法正常运行。练习六、调试器和EZ-KIT板的性能比较实验目的对比程序在仿真环境下和EZ KIT板上的执行性能，并分析原因。二实验内容2.1把练习二改为执行5次、10次，在软仿真环境下运行，并用记录执行的时间;实验原理运用VisualDSP中的linear profiling窗口观察并记录程序在仿真环境和EZKIT板上执行的时间。四实验步骤4.1修改练习二的程序，在程序中加入控制次数的for循环，以控制程序执行5次、10次。4.2编译成功后，打开 VisualDSP的linear profiling窗口观察

19、并记录程序在仿真环境下的执行时间，通过tool/linear profiling打开lin ear p rofili ng 窗口。，按照默认project），其余默，并将上述 窗口中观4.3 生成 LDF文件，点击 toolsexpert linker'creat LDF选项，依次选择1：文件保存的路径（与原工程一致）、(sin gle p rocesstype选择assembly ； 2：选择单处理器4.3在断电情况下连接EZ KIT板，创建LDF文件导入 DSP中，运行程序，在认；3， finishEZ KIT sessi on lin ear p rofili ng 察并记录程序

20、执行的时间。五实验结果及讨论结果:硬件平台5次64点DFT耗时0.94ms10次64点DFT耗时1.88ms软件平台5次64点DFT耗时16s10次64点DFT耗时33s讨论：软件仿真比硬件仿真慢，原因主要有两个方面：第一， 软件仿真是在计算机上模拟程序在DSP上的运行，这就需要一定的系统开销，同时计算机内部的数据传输叫DSP中复杂，这也造成了系统开销；第二，DSP内部存在专用的硬件算法如乘累加，在计算机上是不存在这样的结构的。练习七、用ADSP_TS201 EZ-KIT板实现定时器实验目的利用ADSP-TS201 EZ-KIT板上的硬件资源，通过对板上 A B两个DSP 的FLAGrn示灯进

21、行编程，来完成定时器的操作。二实验内容编译讲义上的程序到没有错误，观察flag丁的亮灭，并分析。三实验原理利用ADSP-TS201 EZ-KIT板上的定时器timer1，产生中断，在中断服务子程序中实现对FLAG的控制四实验步骤4.1连接硬件首先给EZ-KIT板加电，可以看到板上的 FLAG丁在闪，板上的power灯（绿色）持续亮一段时间，RESETT（红色）持续亮一段时间又熄灭。接下来板 子上的USB_monitor灯一直亮。4.2 打开 VisualDSP+，选择 Session/Select Session/ADSP_TS201SEZ_KIT Lite via Debug Age nt，

22、选择 ADS P- TS201 EZ-KIT session。4.3打开工程，编译链接该工程，没有错误后运行程序。观察并记录结果。4.4通过改变中断服务程序xr4的值，可以看到灯闪亮的速度在变化。五实验思考题及心得5.1灯的闪灭是通过什么来控制的?答：灯的闪灭是通过中断服务函数控制的，每次中断到来时，利用异或取反, 使FLAG丁的状态改变，从而控制灯的闪灭。5.2代表秒的灯闪一次，定时器发生了多少次中断?答：代表秒的灯闪一次，定时器发生了 10次中断。练习九、中断响应、标志查询一实验目的1.1模拟产生外部中断和可编程I/O （ Flag03）；1.2学会使用中断服务子程序及查询方式;1.3熟练

23、使用断点调试。实验内容及步骤实验一1、添加外部中断，以IRQO为例。如图1所示。IrntM-nlpt TmnLiUNOKCancelName1 Mm1 MaxJ D4IAGO110000ItDOOitnoiHIJ dFt图1外部中断IRQ2、编写C语言程序，在主程序中以查询标志位 en的方式进行加法计算，en为1则开始进行加法运算（data1n+data2n），在中断服务子程序中加入断点观察每次进入中断服务子程序的指令周期间隔。数组data1和data2，每个长度20 （这个可以自己设定），用in elude函数将外部“ .dat ”文件包含进 来。在中断服务子程序中设置标志位 en为1。加法

24、完成后程序停止，此时在VisualDS P Mermoy中察看计算结果是否正确。根据如下代码提示部分完成程序的编写:void mai n(voidin itial(; /初始化寄存器函数/ irq0_isr（外部IRQ0中断服务子程序，在子程序中断中置en标志位in terru pt(SIGIRQ0,irq0_isr;while(1while(标志位=1resultn=input1n+input2n;/ 实现加法运算;en=0;/设置标志位为0;n+;if(n= N+1 /数组长度break;/跳出最外层while循环;void in itial(配置系统配置寄存器SYSCON外部总线接口寄存

25、器_builtin_sysreg_write(_SYSCON,0x279e7;_builtin_sysreg_write _INTCTL,0x0; /设置外部中断触发方式为边沿触发,关闭定时实验二设置可编程I/O的输入，在程序中查询FLAG灯的变化状态，加断点观察中断 产生的时间间隔。1设置FLAG输入，本实验以FLAG1为例，如图2所示:OftCM!SSSiSaiiManif 1 M'k fl 忖* I Di*W 1EU 1IOODO 100图2 FLAG1输入2、编写如下程序void ma in (voidinitial(; /初始化寄存器函数while(1查询flagflag仁built in _sysreg_read(_SQSTAT; /SQSTAT 程序控制状态寄存器flag仁flagi &0x