通信系统仿真-论文

《通信系统仿真》专业：通信工程班级：06级姓名：学号：2009年12月20日目录21、前言………………………………………………………………………3-32、正文………………………………………………………………………4-83、参考文献…………………………………………………………………9-93摘要TLAB名字是由MATrix和LABoratory两个词的前三个字母组合而成的.它是MathWorks公司于1982年推出的一套高性能的数值计算和可视化数学软件.被誉为巨人肩上的工具.由于使用Matlab编程运算与人进行科学计算的思路和表达方式完全一致，所以不象学习其它高级语言--如Basic、Fortran和C等那样难于掌握，用Matlab编写程序犹如在演算纸上排列出公式与求解问题，所以又被称为演算纸式科学算法语言一般数值分析、矩阵运算、数字信号处理、建模和系统控制和优化等应用程序，并集应用程序和图形于一便于使用的集成环境中.在这个环境下，对所要求解的问题，用户只需简单地列出数学表达式，其结果便以数值或图形方式显示出来.引言MATLAB的含义是矩阵实验室（MATRIXLABORATORY），主要用于方便矩阵的存取，其基本元素是无须定义维数的矩阵.MATLAB自问世以来，就是以数值计算称雄.MATLAB进行数值计算的基本单位是复数数组（或称阵列），这使的MATLAB高度向量化.经过十几年的完善和扩充，现已发展成为线性代数课程的标准工具.由于它不需定义数组的维数，并给出矩阵函数、特殊矩阵专门的库函数，使之在求解诸如信号处理、建模、系统识别、控制、优化等领域的问题时，显得大为简捷、高效、方便，这是其它高级语言所不能比拟的.美国许多大学的实验室都安装有MATLAB供学习和研究之用.在那里，MATLAB是攻读学位的大学生硕士生博士生必须掌握的基本工具.MATLAB的应用范围很广，在医疗技术研究中，计算机技术研究中、嵌入式技术研究中、DSP技术研究中、电源技术研究中都有应用，本文主要从MATLAB在DSP中的应用来研究。4正文MATLAB的作用：基本的数据处理优化和解方程动态过程仿真：实时的和非实时的数据来源：Excel、数据库、A/D等等嵌入式控制：Pc/104和DSP神经元网络、小波分析、GA等等虚拟现实网络MATLAB辅助DSP设计的方法概述实现MATLAB辅助DSP设计的方法有两种，即CCSLink和EmbeddedTargetforTITMS320C6000DSP。前者可实现MATLAB与C2000/C5000/C6000DSP之间的数据交互和程序控制；后者可实现从Simulink模型自动生成DSP目标代码，但仅支持C6000系列的部分DSP器件。考虑到实验平台是以C5000DSP为基础的，并且第二种工具较容易操作等因素，下面主要研究利用CCSLink辅助DSP设计的方法。1、mbeddedTargetforTITMS320C6000DSPCodeComposerStudio(CCS)是TI提供的DSP的集成开发环境。对于C6000系列而言，利用工具EmbeddedTargetforTITMS320C6000DSP，用户可以从Simulink模型直接生成针对DSP目标板的可执行文件或者CCS工程文件，实现了MATLAB/Simulink与CCS集成开发环境(IDE)的无缝连接，从而把设计人员从繁杂的DSP编程中解脱出来，将主要的精力放在设计而不是DSP编程上，所以极大地缩短了系统的测试和开发周期，进而保证了所设计系统的高性能。2、MATLABLinkforCodeComposerStudio集成在MATLAB6.5或更高版本中的CCSLink工具把MATLAB、TIDSP及其集成开发环境CCS连接在一起，使得我们可以在MATLAB环境下就可以完成对CCS和DSP目标板的操作，包括与目标内存之间的数据交互，检测处理器的状态，控制DSP程序的运行等。它提供了MATLAB、CCS和目标DSP板的双向连接，开发者可以利用MATLAB中强大的可视化、数据处理和分析工具对DSP中的数据进行分析和再处理，以辅助DSP系统设计，这样就极大地降低了开发人员编写DSP代码的难度和工作量，提高了整个DSP系统的性能和可行性。我们利用CCSLink提供的相关函数能够实现MATLAB与目标DSP的存储器及寄存器间的信息交互，如同操作MATLAB变量一样来读、写TIDSP中的存储器或寄存器，即整个目标DSP对于MATLAB而言是透明的。CCSLink支持TI的C2000/C5000/C6000等多系列DSP，应用广泛。3、SLink的实现方式及工作原理（1）CCSLink的实现方式利用CCSLink即可实现MATLAB、CCS和目标DSP之间的信息交互。CCSLink共提供了3种实现方式：1）使用与CCSIDE的连接对象。利用此对象来创建CCSIDE和MATLAB的连接，从MATLAB命令窗中即可运行CCSIDE中的应用程序，查询目标DSP的状态信息，修改或读取目标DSP的存储器或寄存器中的数据，甚至可以调试DSP程序。其工作原理见3.2节。2）使用与RTDX(Real-TimeDataExchange)的连接对象。提供MATLAB和硬5件DSP之间的实时通信通道。利用此连接对象，可以打开、使能、关闭或禁止DSP的RTDX通道，利用此通道可以实时地向硬件目标DSP发送和取出数据而不必停止DSP上运行的程序。例如把原始数据发送给程序进行处理，并把处理结果取回到MATLAB空间中进行分析。RTDX连接对象实际上是CCS连接对象的一个子类，在创建CCS连接对象的同时创建RTDX连接对象，它们不能分别创建。3）使用嵌入式对象。在MATLAB环境中创建一个可以代表嵌入在目标C程序中的变量的对象。利用嵌入式对象可直接访问嵌入在目标DSP中的存储器和寄存器中的变量，即把目标C程序中的变量作为MATLAB的一个变量对待。（2）CCSLink的工作原理上面的三种方式在具体实现时略有不同，但是其原理相似，不失一般性，我们通过第一种方式——与CCSIDE的连接对象，来介绍CCSLink的工作机理。1）选择目标DSP并创建CCSIDE连接对象在MATLAB环境下使用函数ccsdsp并配置相应的属性名和属性值即可创建一个CCSIDE连接对象cc：cc=ccsdsp(‘boardnum’,x,‘procnum’,y,‘timeout’,z);其中，boardnum、procnum、timeout都是属性名。分别表示创建连接对象时所使用的目标板编号、DSP处理器编号及全局超时值。2）在MATLAB环境下加载CCS工程文件projfile=fullfile(matlabroot,‘toolbox’,‘tiddk’,…,‘ccsproject.pjt’);%取得工程文件projpath=fileparts(projfile);%取得工程文件路径open(cc,projfile);%加载工程文件，该函数实现把文件filename加载在CCSIDE中cd(cc,projpath);%改变CCS的工作路径3）在MATLAB环境下编译、连接并下载可执行文件创建了CCSIDE连接对象并加载工程文件后，即可通过函数build()编译和链接工程文件，以生成针对目标DSP的可执行文件（.out文件）；通过函数load()下载可执行文件到目标DSP中。具体方法如下：build(cc);%编译、链接工程文件，生成目标DSP可执行文件load(cc,‘ccsproject.out’);%下载可执行文件到目标DSP，可执行文件名与工程名相同4）在MATLAB环境下控制程序运行CCSLink提供了3个函数用于控制目标DSP的运行状态。通过这些函数，在MATLAB下对DSP的控制操作就像在CCSIDE中一样方便，实际上也确实是这样进行的。halt(cc);%中断处理器的运行restart(cc);%复位程序计数器PC，重新执行程序run(cc,‘runtohalt’,30);%执行程序到断点处，最后一个参数表示超时时间，此处为30s5）在MATLAB环境下对CCSIDE连接对象进行操作在MATLAB中可以方便的操作DSP存储器和寄存器中的数据，包括读操作和写操作。下面以读写目标DSP存储器和寄存器中的数据来说明。6date=read(cc,[256,0],‘double’,100);%从目标存储器目的地址100H(即十进制数256)处开始，读取100个双精度数并赋给MATLAB对象datereg=regread(cc,‘A0’,‘2scomp’);%读取DSP寄存器A0中的数据，并转化成双精度数分配给MATLAB对象regwrite(cc,[256,0],double([1.02.03.0]));%向目标存储器100H开始的地址处，写入三个双精度数据6）关闭CCSIDE连接对象当一切操作执行完成后，需要删除ccsdsp对象，以收回系统资源。clearcc;%关闭CCSIDE连接对象（4）CCSLink应用实例通过CCSLink工具，可以使用MATLAB来进行DSP的辅助设计。在此以FIR滤波器为例，介绍CCSLink在DSP设计中的应用。在该实例中，FIR滤波的核心内容由DSP实现，而算法仿真、系数产生、用户界面、结果验证等部分由MATLAB辅助完成。1）FIR滤波器设计设输入信号为x(n)，输出信号为y(n)，则FIR滤波器对应的差分方程可表示为：其中，h(k)(k=0,1,2,…,N-1)表示滤波器系数序列，是滤波器设计的关键。h(k)可以通过MATLAB工具箱提供的滤波器函数fir1()和fir2()来获取：首先通过建立MATLAB/Simulink模型，根据所需设计低通滤波器的性能要求来确定滤波器的阶数和截止频率等参数。本例所设计的FIR低通滤波器参数为：截止频率(Fc)2kHz，采样频率(Fs)8kHz，阶数(N)40。然后借助函数fir1(N,WN)即可返回对应滤波器的系数矩阵，其中N为阶数，WN为滤波器的归一化截止频率，即WN=2Fc/Fs，且0≤WN≤1。2)DSP滤波算法实现结合MATLAB产生的滤波器系数，就可以编写DSP滤波主程序了，关键代码如下：FIR_TASK：STM#K_FIR_BFFR,BK;设置FIR循环缓冲区的大小LD*INBUF_P+,-1,A;数据输入FIR_FILTER：STLA,*FIR_DATA_P+%;最新数据替代最旧数据RPTZA,(K_FIR_BFFR-1)MAC*FIR_DATA_P+0%,*FIR_COFF_P+0%,A;乘累加STHA,*OUTBUF_P;数据输出3)MATLAB辅助设计本设计中，算法仿真、系数产生、用户界面等部分由MATLAB辅助完成。算法仿真、系数产生方法在4.1中已详细介绍，本节重点介绍用户界面部分。本图形用户界面可以实现如下功能：修改参数，如滤波参数n和Wn，这种参数选取方法避免了反复修改DSP程序的弊端，可以满足不同层次使用者的使用；控制操作过程，如加载工程文件、下载可执行文件、运行程序等；绘制DSP输入数据的时域及频域信号；绘制滤波后DSP输出数据的时域及频域信号；使用MATLAB提供7的滤波函数对输入数据进行滤波计算，并把所得结果与DSP输出信号进行比较，以验证或改良DSP算法。用户界面的流程图如图1所示，主要由图形界面和回调函数两大部分组成。其中，图形界面由文件FIR_BPF.m实现；回调函数总共由6个文件实现，分别是input.m，input2.m，output.m，output2.m，lvbopin.m和lvbohou.m。当输入带有干扰的方波信号后，分别查看输入波形、MATLAB滤波输出信号、DSP滤波输出信号，结果如图2所示。比较MATLAB滤波输出信号和DSP滤波输出信号，据此可以采用修改滤波器系数或DSP程序等措施来改进DSP滤波算法。图1用户图形界面流程图8图2实验结果总结MATLAB中包括了被称作工具箱（TO