Vc++下如何利用Matlab工具箱进行数字信号处理_计算机论文.docVIP

Vc++下如何利用Matlab工具箱进行数字信号处理_计算机论文 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文件中加入”engine.h” 和“math.h” 头文件，下面给出部分程序清单。 Void CtestmatlabDlg::OnMatlabEngine(){ Engine *ep mxArray* T=NULL,*result=NULL,*mFs=NULL,*mnfft= NULL double datax[1024] char buffer[1024] for(int j=0j1024j++)//注：如通过采集卡采集数据可将采集的数据放在datax[]数组中,此循环就不需要 { double samt=(double)(1.0/1024) datax[j]=sin(2.0*63.0*samt*3.1415926+1.15*3.1415926) } double *pPxx,*pFxx if(!(ep=engOpen( \0))){//打开Matlab引擎，建立与本地Matlab的连接 fprintf(stderr,\n Cant 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