信号与系统论文

相关论文：Vc++下如何利用Matlab工具箱进行数字信号处理摘要： 本文详述了在Vc环境下如何利用Matlab工具箱进行数字信号处理，全文以Matlab工具箱中功率谱密度分析函数为例，介绍了通过Matlab自带的引擎、Matlab自身的编译器以及利用MathTools公司的Matcom进行对工具箱函数的调用。 关键词：Matlab M-文件 引擎 编译器 Matcom Vc++Matlab的信号处理工具箱是信号算法文件的集合，它处理的基本对象是信号与系统，信号处理工具箱位于目录、Toolbox\Signal下，利用工具箱中的文件可以实现信号的变换、滤波、谱估计、滤波器设计等。在其它的环境如Vc下如果能调用Matlab工具箱中的文件，会大大地加快一些算法的实现，同时其可靠性也很高。利用Matlab引擎 Matlab引擎采用客户和服务器计算方式，在运用中，Vc的C语言或C++语言的程序作为前端客户机，它向Matlab引擎传递命令和数据信息，并从Matlab引擎接收数据信息，它提供了下列几个函数： engOpen, engGetArray, engPutArray, engEvaString, engOutputBuffer ,engClose与客户机进行交互。下面例程是在Vc下建一个基于对话框的应用程序，在对话框中设置一个Button控件OnMatlabEngine.,在对话框 .cpp文件中加入”” 和“” 头文件，下面给出部分程序清单。Void CtestmatlabDlg::OnMatlabEngine(){Engine *ep;mxArray* T=NULL,*result=NULL,*mFs=NULL,*mnfft= NULL;double datax[1024];char buffer[1024];for(int j=0;j<1024;j++)//注：如通过采集卡采集数据可将采集的数据放在datax[]数组中,此循环就不需要{double samt=(double)(); datax[j]=sin(**samt**);}double *pPxx,*pFxx;if(!(ep=engOpen(" \0"))){//打开Matlab引擎，建立与本地Matlab的连接fprintf(stderr,"\n Can't start MATLAB engine\n");exit(-1);} double Fs[1]={1024};//因为Matlab所有参与运算的参数都是矩阵的形式，因而下列几行将参数转变double nfft[1]={1024};//成Matlab可接受的矩阵形式。T=mxCreateDoubleMatrix(1,1024,mxREAL);mnfft=mxCreateDoubleMatrix(1,1,mxREAL);mFs=mxCreateDoubleMatrix(1,1,mxREAL);mxSetName(T,"T");mxSetName(mnfft,"mnfft");mxSetName(mFs,"mFs");memcpy((char*)mxGetPr(T),(char*)datax, 1024*sizeof(double));memcpy((char*)mxGetPr(mnfft),(char*)nfft, sizeof(double));memcpy((char*)mxGetPr(mFs),(char*)Fs,1*sizeof(double));engPutArray(ep,T); //将转化的参数放入引擎中，此时可在Matlab command窗口下查看此参数engPutArray(ep,mnfft); engPutArray(ep,mFs);engEvalString(ep,"[pxx,fo]=psd(T,mnfft,mFs);"); //利用引擎执行工具箱中文件engOutputBuffer(ep,buffer,512); //如只想看显示图形，可将返回参数去掉，psd无返回参数缺省情况下会自动画图形result=engGetArray(ep,"pxx");//取出引擎中的数据放在所指的区域中供后续处理pPxx=mxGetPr(result); result=engGetArray(ep,"fo");pFxx=mxGetPr(result);engEvalString(ep,"plot(fo,10*log10(pxx));");//利用引擎画图engEvalString(ep,"title('功率谱分析');");engEvalString(ep,"xlabel('Hz');");engEvalString(ep,"ylable('db');");mxDestroyArray(T); //释放内存mxDestroyArray(mFs);mxDestroyArray(mnfft);mxDestroyArray(result);engEvalString(ep,"close;");engClose(ep);}上述程序在Vc下编译需要将 和两个动态库利用以下的命令：lib/def:<自己的Matlab的安装路径,下同>e:\ Matlab\extern\include\*.def /machine:ix86 /out:*.lib来生成程序所需的静态连接库和,将和所在的目录加入Vc++ project/link/object/library modules下即可。利用Matlab自身的编译器调用工具箱中的函数 Matlab的编译器可将Matlab的M文件转换为为C或C++的源代码以产生完全脱离Matlab运行环境的独立的运用程序，但Matlab本身的资料说明编译器如用来建立独立的运用程序，不能调用Matlab工具箱中的函数，这非常不利于搞一些特殊的算法。本人研究了一段时间发现，工具箱中的函数既然是M文件就一定可以用编译器来编译，以提供如Vc的调用函数，但是编译器只能编译一个独立的M文件，即这个 M文件不依赖于其他的M文件。如果M文件中又调用了其他的M文件，可将被调用的M文件拷贝到调用M文件的相应位置，作适当的改动就可以用于编译器编译。编译器不支持图形函数，所以M文件中如有图形函数需注释掉。当Matlab的编译器mcc加入适当的参数-e(mcc –e *.*)或-p(mcc –p *.*)就可生成将输入的M文件转换为适用于特定运用的C或C++源代码。这样如果要在Vc下编译通过,还需连入以下几个库, , , . . 以及Matlab C MATH库，建议采用前述的方法将动态连接改为静态连接。对于C/C++编译环境的设置，在Matlab command窗口下运行mex –setup 然后依提示操作，而对于C/C++连接环境的设置，运行mbuild –setup依提示操作即可。下面给出利用编译器将Matlab工具箱中文件生成可供Vc调用的函数。将文件拷贝一份至Matlab\bin目录下，改写相应调用的M文件如, 等。为生成的代码简洁，对于采集数据处理输入参数很明了的情况下可作大量的删减，最终使成为一个不依赖于其他M文件的独立的M文件，注意千万注释掉作图代码，最终改成如下形式，限于篇幅给出关键的几步： function [Pxx,f]=psd(Fs,nfft,noverlap,x)window=o,5*(1-cos(2*pi*(1:nfft)’/(nffft+1)));//hanning 窗dflag=’none’;window=window(;)………………………………….以上只要稍懂Matlab语言和信号处理知识就可完成这项工作。假设上述代码重新存为,在Matlab command 窗口下设置好环境参数运行mcc –e testwin,则可在Matlab\bin下生成 ,如运行mcc –p testwin 则生成. Vc下建立一个基于对话框的文件，然后在对话框里加一个Button控件OnButtonPsd 将上述生成的.c文件的头文件加入到工程的.cpp中，且将#ifdef_cplusplus extern “c”{#end ifc代码声明加入Vc的包含文件和生成的.C的包含文件之间将#ifdef_cplusplus}#end if加入.cpp文件未尾为了简洁且便于处理将生成的c函数稍改动，给出部分代码如下： void CTestpsdwinDlg::OnButtonPsd(){mxArray* x_rhs_;//指向采集数据存放区Fs=23510;//数据采集的频率 nfft=1024;//1024点的fftdouble datax[1024]//采集的数据x_rhs_mxCreateDoubleMatrix(1,1024,mxReal);memcpy(mxGetPr(x_rhs_),datax,1024*sizeof(double));noverlap=512;……………….……………….mccCopy(&Pxx,&Spec);mccCopy(&f,&frevgg_vector);for(int j=0;j<(int)(nfft/2+1);j++){datap[j]=mccGetRealVectorElement(&Pxx, (j+1));//功率谱密度存于datap[]数组dataf[j]=mccGetRealVectorElement(&f, (j+1));//相应频率存于数组dataf[]中}mccFreeMatrix(&Pxx);……………….SendMessageBox(WM_PAINT,0,0);//利用Vc下的图形函数画图Return;}如上生成的程序可读性不太好，而生成的c++代码则可读性较好，但千万注意只能用 Matlab的MATH库，不可用c++的MATH库，否则编译会出错，限于篇幅在此不述。3）利用Matcom调用工具箱中的函数Matcom编译M文件，先将M文件按照与Matcom的cpp库的对应关系翻译为cpp源代码，然后用对应版本的c编译器将cpp文件编译成相应的exe或dll文件，所以第一次运行要指定c编译器的路径，否则无法编译，指定好的编译信息就写在Matcom\bin\文件中，不过这一步按装matcom时，它自动寻找编译器并将其写入文件中，版中使用 OCX控件，因而它支持图形操作。我们依然用上述的文件，不要将图形函数注释掉，利用Mideva来生成可被Vc调用的信号处理程序。运行Mideva在主界面上直打开M文件，在菜单中选择compile to dll,输入testwin..在Matcom debug目录下可以找到这样的几个文件， ,等。 将上述和加入工程中，project/add to project/files并且在相应的文件中加入”” 加连接库：Tools\option\directory\ ， 选include选项,加入e:\matcom45\lib (包含) library选项，加入e:\matcom45\lib4) project\add to project\files 文件类型选项选（.lib）将e:\matcom45\lib\加入工程中编译运行。相应代码如下：void CtestmatcomDlg::OnpsdButton(){double datap[512],dataf[512];initM(MATCOM_VERSION);//初始化matcom库Mm Fs,nfft,noverlap;//创建矩阵Mm x=zeros(1,1024);Fs=1024;nfft=1024;noverlap=128;dMm(Pxx_o);dMm(f_o);//创建并命名矩阵datax[];//数据采集的数据存于此数组中for(int i=1;i<=1024;i++){(1,i)=datax[i+1];//给x阵赋值}testwin(Fs,nfft,noverlap,x,i_o,Pxx,f_o);//matcom生成的函数for(i=0;i<513;i++){//取出功率谱密度分析结果dataf[i]=(i+1,1);datap[i]=(i+1,1);}exitM();return;}可见利用Matcom进行M文件转换非常的容易，生成的代码可读性很好，以上的转换同时生成了可供Vc调用的动态连接库，其使用和一般的动态库一样使用。同时需指明Matcom不仅可转换独立的不依赖于其它M文件的M文件，同时可转换调用其它M文件的M文件嵌套。条件是这此M文件在同一个目录下面，如前所述的可直接用上述方法转换，生成了多个重载形式的psd函数结论： 利用Mtlab引擎调用工具箱中的函数可节省大量的系统资源，应用程序整体性能较好，但不可脱离Matlab 的环境运行。用Matlab编译器进行工具箱函数的调用，须转换相应的M文件使其成为独立的M文件，且不支持图形函数，转换的代码可读性不太好。用Matcom 进行转换非常方便，生成的代码可读性很好，支持图形函数，且代码执行的速度比不转换平均要快倍以上。以上程序在Vc++ 中调试通过,以上方法在工程实践中已得到很好的运用。仅供参考，请自借鉴希望对您有帮助

去找一本《信号与系统》，看看他的前言和绪论就基本上知道这本书是干什么的，这个课程是学习什么的了，而且有时候说明时还会带有例子的，易理解。

信号与系统整本书讲的都是信号的处理过程（模拟信号，数字信号……）包括时域分析与频域分析，对于时域模拟信号的处理主要方法是傅里叶变换分析，与拉普拉斯变换，前者现在应用特别广泛，计算机分析也是基于此的，后者主要是为了简化手工计算，其没有物理意义。对于数字信号的处理主要方法是离散傅里叶分析（现在计算机中的MATLAB处理时用的是FFT（快速散傅里叶变换）），对应于拉普拉斯变换，数字信号处理也会用Z变换来计算。信号与系统在日常生活中应用特别广泛，几乎每一个电子系统都会用到其知识。目前数字信号处理在实际应用中比较多，其方法也比较成熟。 这是我个人知道，因为我是读电信的，如果你感觉不够明白的话，你可以去问一下你身边读电子信息工程或电子科学技术的学生，这两个专业的学生必学信号与系统。

高速铁路信号是高速列车安全、高密度运行的基本保障。下面是我整理的高速铁路信号技术论文，希望你能从中得到感悟!

基于无线通信技术的高速铁路信号系统应用

摘 要

高速铁路信号系统是高速列车安全、高密度运行的基本保障。无线通信技术在铁路信号系统的应用，不但减少了高速铁路的信号系统成本，还较好的确保了高速铁路的安全。随着科学技术的进步，高速铁路不断的向着智能信息化转变，这就给无线通信技术领域提出了更加严格的要求，为了适应高速铁路的快速发展，各国都在潜心研究基于无线通信技术的新一代的铁路信号系统。本文介绍了国外无线通信系统在高速铁路信号系统中的发展情况，分析了运用无线通信技术的高速铁路信号系统的特点和问题，并探讨了无线通信技术在高速铁路信号系统中的应用。

【关键词】无线通信 高速铁路 信号系统

在整个高速铁路工程中，虽然信号系统的投资总额所占比率较小，但其起到的作用十分关键。由于轨道电路传输环境较差、传输信息的速率较低、设备更新维护费用高，所以基于轨道电路的列车控制系统已经不能满足高速铁路的快速发展要求。在80年代，国外开始研究基于无线通信的铁路信号系统TBS(Transmission Based Signalling)，希望通过无线通信技术的应用来提高铁路的管理职能、缩短列车间隔时间、节约能源、降低系统的成本。1995年在关于TBS的国际会议中，会议代表分析了无线通信技术在铁路信号系统应用的的可行性，并指出了无线通信技术可能给铁路信号系统带来的积极影响，表明了TBS将会成为未来铁路信号系统的发展方向。

1 国外TBS的发展情况

北美TBS的发展情况

1983年，美国铁道协会和加拿大铁道协会共同最早提出了基于无线通信的先进列车控制系统ATCS。ATCS主要是通过数字数据通信手段和先进的微处理器获取列车的精确位置和速度等信息，并对列车进行安全控制。ATCS的运用不仅避免了很多地面信号设备的安装，节省了系统成本，还消除信号盲区，增强了列车的安全系数。ATCS是由中央控制系统、无线数据通信网络、车载设备、路旁设备和线路维护人员移动终端五个子控制系统构成的。它的系统结构设计和功能模块的划分为以后基于无线通信的铁路信号系统奠定了基础。随着无线通信技术的发展，在ATCS之后北美又出现了很多基于无线通信的铁路信号系统，其中ARES可以提供非常可靠的检查和平衡手段，在很大程度上降低了人为操作失误造成的错误，使列车行驶更加安全。另外，PTS、PTC、AATC、ITCS等系统也是比较著名的。

欧洲TBS的发展情况

1992年国际铁盟下属的欧洲铁路研究机构提出了一套欧洲的铁路运输管理系统，包括车票发售、各国铁路互操作性等多个方面，ETCS就是其中非常重要的一部分。在欧共体委员会设立标准化欧洲铁路控制系统项目ETCS之前，欧洲各国铁路标准和模式不尽相同，轨距、信号设备、供电设备也不一样，因此各国只能使用自己的ATP、ATC系统。各国铁路制式上的差异使得欧洲铁路很难形成连续运输。在设立了标准化欧洲铁路控制系统项目ETCS后，各国的铁路开始逐渐按照统一标准进行规范，并逐渐取代各国不同的列车自动控制系统和防护系统。ETCS的目标就是要实现欧洲铁路的统一，提高各国铁路的互操作性，使铁路控制系统的功能和设备更加规范。

日本TBS的发展情况

在日本铁路信号系统的发展历程中，先后出现了ATS、现行ATC、数字式ATC、计算机和无线通信辅助信息控制系统等。其中现行ATC作为一种列车超速防护系统，以良好的自动制动功能保护了列车的安全。但在系统工作时，采用的最强的自动制动，影响了乘客的舒适程度。在1987年，日本开始基于无线通信的铁路信号系统的研究，为CARAT的出现奠定了坚实的基础。CARAT的使用能够使列车连续测定自身位置和行驶速度，使地面系统能够很好的了解列车运行情况，保证列车的运输安全。

2 TBS的特点和问题

在速度比较高的高速铁路上，距离比较近时，可以采用红外、蓝牙等无线通信技术实现对列车的控制;在距离比较远时，则可以通过全球定位控制系统、信标、计轴装置等来测定列车的速度和位置。车载计算机可以通过无线收发装置将列车的速度、位置信息发送给调度控制计算机，通过调度控制计算机的处理，再将列车允许的最大速度等信息通过无线通信发回给列车计算机。列车司机可以根据车载计算机的提醒进行相应的操作，如果列车司机没有及时作出反应，信息控制系统还可以自行将车速降低到允许范围以内。

TBS的特点

(1)在TBS中，主控中心可以根据列车的运行状态和操作状态通过车载计算机来调整列车的运行，加大了高速铁路信号系统的管理职能，保证了列车的安全，提高了铁路线路的通行能力。

(2)在无线通信信号系统控制下，列车和地面的可靠信息量增大，列车运行变得更加稳定，且避免了不必要的加速和制动，节约了能源，也让旅客乘车变得更加舒适。

(3)无线通信技术的运用，省掉了大量的地面信号装备，大大减少了设备的安装、维护、修整费用。

(4)无线通信信号系统的适应能力极强，通过软件上的调整就可以使列车的运行速度提高，且能够自动调整运行图，大大的提高了铁路运输管理能力。

(5)无线通信信号系统还可以通过车地间的双向信息通道实现列车的闭锁控。

TBS的问题

(1)高铁信号系统使用轨道电路只能使用较低的信息发送频率，传输环境恶劣，很难让电码的传送速率满足高速铁路的运行速度要求。

(2)TBS通过环线设备和应答器件接受数据信息，列车进行操作可能会有时间上的延迟，可能会给列车的运行造成不良的影响。

(3)轨道间的电缆电线作为车地之间的双向信息通道，虽然传输信息量大，抗干扰能力强，但设备费用较高，且防盗能力很差，一旦丢失，后果严重。

3 无线通信技术在高速铁路信号系统中的应用 微机联锁

无线通信技术在微机联锁方面运用的可行性还需进一步研究，但ATCS中提出，可以将检测到的道岔、信号机闭锁状态发送给主控中心，并利用道旁接口单元来接收主控中心的控制命令，以实现控制一组道岔、信号机动作的目的。另外道旁接口单元可以利用无线信道联系控制中心，通过电缆连接现场设备，从而检测并控制一些辅助的子系统。目前看来，无线通信技术用于微机联锁的现场设备可能会增加一些投资，且大型站场道岔众多，干扰较大，但还是具有较好的发展前景。

集中调度

在调度集中系统中，调度中心职要根据车站到发线占用情况和区段内闭塞分区大概了解列车运行的状况，并根据得到的信息排列进路。但利用TBS，控制系统就能够准确的了解列车运行的位置、速度，并根据沿线的信号系统情况发送列车控制命令，保证列车在最短的实践间隔内高速、安全、稳定的运行。无线通信技术赋予列车与控制中心的双线数据通信，给列车的运行带来了很大的方便，且实现了行车指挥自动化。

中继器

在高速铁路的实际运行中，我不可能在所有的高速铁路中都设这无线通信基站，这样不但增加了设备投资，还使无线通信铁路信号系统失去了存在的真正意义。有了中继器，基站就可以通过中继器接受和发送一些射频信号，从而使基站不仅可以管理基站区域范围内的站区，还能够将管理中继器管理的一些车辆和线路。

提高平交道口的通过效率

为了提高平交道口的防护能力和和通过效率，防止由于无线设备故障造成不必要的损失，主控中心按照时间间隔不断的查询道口的运行状态，并将查询信息及时反馈给接近道口的列车。另外主控中心通过接收的列车位置、速度信息，可以计算列车通过道口的时间，并根据实际情况设定列车的最大允许速度和列车运行线路参考。这样，列车通过平交道口就有了安全保障，而且还大大提高了道口的通过效率。

加强维修处防护

在高速铁路某路段需要进行维修时，维修部门可以通过移动终端将维修点输入到系统中，通过主控中心的传送，列车就可以很好的了解路段情况。在实际的运行中，列车可以根据了解到的维修点信息对列车进行操作，另外在列车接近维修点事，移动终端接受到地面系统的警报信号，以保证列车能够及时在维修段之前停车。

4 总结

随着高速铁路的不断发展，要确保列车的安全，先进的信号系统成了高速铁路运行的重中之重。在高速铁路信息系统中，无线通信的运用仍处于初期阶段，在具体的TBS规划时应充分考虑其与全路运输管理系统的接口，使无线通信技术更充分的运用在高速铁路的发展当中。

参考文献

[1]闵耀兴.我国铁路列车安全控制系统的现状[J].哈铁科技通讯，1997(04).

[2]姚丽娟.我国铁路信号系统的现状与发展[J].铁道通信信号，2003(04).

[3]步兵.基于通信的列车控制系统的可靠性分析方法[J].交通运输工程学报，2001(01).

[4]杨绚，陈德旺，陈荣高.速铁路列控系统主动安全控制的分析与思考[J].铁路计算机应用，2012(08).

作者简介

孙屹枫(1982-)，男，天津市人。中国民用航空大学大学本科毕业。研究方向：铁路信号。

作者单位

铁道第三勘察设计院集团有限公司电化电信处 天津市 300251

点击下页还有更多>>>高速铁路信号技术论文