引言
Visual C++自诞生以来,一直是Windows环境下最主要的应用开发系统,利用Visual C++开发系统可以完成各种应用程序的开发,从底层软件直到上层直接面向用户的软件都可以用Visual C++来开发,而且Visual C++强大的调试功能也为大型复杂软件的开发提供了有效的排错手段[1],但是其数值计算能力和二、三维图形显示方面却不理想。Matlab是由MathWorks公司于1984年推出的一套数值计算软件,其功能十分强大,可以实现数值分析、优化、统计、偏微分方程数值解、自动控制、信号处理、图像处理等若干个领域的计算和图形显示功能[2]。但是,Matlab在实现人机对话、数据的传输方面却不是很方便,将两者结合起
来,取长补短,发挥两种语言的优势,是开发人员研究的一个方向。
1 VC与Matlab混合编程的方法:
Visual C++和Matlab混合编程有多种方法[3-7],其中最主要的有以下三种方法:
1.1 在VC++中利用Matlab Engine函数
Matlab Engine函数库是MathWorks公司提供的一组函数库,它提供了一种在用户程序进程中与独立的Matlab进程通讯的方法,它包括了和 Matlab进行交互所必需的全部功能,用户不用去关心Matlab Engine是如何实现的,只需利用这些函数,开发者就可以方便的操纵Matlab完成所需的功能。具体应用中往往在Visual C++中设计程序框架,以编译的程序作为前端客户机,通过调用Matlab Engine在后台与Matlab服务器建立连接,实现动态通讯。
1.2 利用Matlab丰富的数学函数库
Matlab 中包含内容丰富的数学库函数,同时它还提供了C语言和C++语言的数学函数接口,从Matab5.1版本开始,MathWorks公司推出了一系列的 Matlab自带编译器来解决Matlab与C++接口问题,LCC编译器可以将Matlab的C/C++数学库编译成VC++能识别的代码嵌入VC++ 环境,用户可以方便的在VC++的IDE(集成开发环境)中调用这些代码。
1.3 利用Matcom
Matcom 是MathTools公司推出的一个能将M文件转化成相同功能C++代码的工具,是为Matlab中的M文件进行高效解释和调试的集成开发环境, Matcom编译M文件,先将M文件按照与Matcom的cpp库的对应关系,翻译为cpp源代码,然后用对应版本的C编译器将.cpp文件编译成相应 的.exe文件。
本文采用的是第一种方法,利用VC++强大的可视化功能编辑应用程序对话框,方便对采集卡的控制,采集数据,在VC++环境中调用Matlab Engine函数,对采集到的数据进行处理,并显示处理后的图形。
2 本文所用采集卡简介[8]:
本 文采用UA302型USB采集卡,这种采集卡具有16bit的分辨率,满量程时精度优于0.02%,最高采样频率可达100KHz,输入通道为16或32 模入通道,可以采用定时器触发和软件触发两种触发方式。UA302采集卡可以使用各种Windows编程工具编程,专用的动态链接库UA300.DLL提 供了简洁高效的采集和控制函数,支持UA302采集卡的各种功能,用户可简单方便的调用这些函数完成各种采集工作。本文要用到的采集卡函数:
OpenUA300 打开UA302设备
CloseUA300 关闭UA302设备
minit 单或多通道多点采集初始化(第一种方式)
readdata 单或多通道多点采集(第一种方式)
3 混合编程实例:
3.1 用VC++6.0编写如图1所示的基于对话框的应用程序界面
对话框上的Button控件用来控制采集卡采集信号并对其进行处理,以灰色显示的按钮在程序运行时不可用,Picture控件用来显示采集到的信号。
图1 应用程序界面
3.2本文用到的MatlabEngine函数:
engopen 打开Matlab引擎
engcolse关闭Matlab引擎
mxCreateDoubleMatrix创建双精度数据类型的Matlab矩阵
engPutVariable向Matlab引擎传输一个矩阵
engOutputBuffer创建字符缓冲区以获取Matlab的文本输出
engEvalString执行Matlab命令
3.3程序的具体实现:
本程序使用了多线程编程技术,除了主线程外采用了两个辅线程,一个用来采集信号数据,另一个用来显示采集到的信号。
采集数据线程程序:
UINT AdoptDataThreadProc(LPVOID pParam)//采集数据线程函数,该线程不能停否则会有漏点
{
//初始化端口,用于单次采样
minit(husb,0,1,1);
readdata(husb,singleaddata,6000000/SampFrequency,1024);
return 1;
}
绘制波形线程程序:
UINT DrawThreadProc(LPVOID pParam) //该线程负责绘制波形
{
do
{
CDC *pDC=pWnd->GetWindowDC();
CBitmap *m_pBitmap;
CDC *m_pdcmem;
m_pdcmem=new CDC;
m_pBitmap=new CBitmap;
m_pdcmem->CreateCompatibleDC(pDC);
m_pBitmap->CreateCompatibleBitmap(pDC,rect1.right,rect1.bottom);
CBitmap* poldbitmap=m_pdcmem->SelectObject(m_pBitmap);
NCCDraw(m_pdcmem,1);
Sleep(100);
pDC->BitBlt(0,0,rect1.Width(),rect1.Height(),m_pdcmem,0,0,SRCCOPY); m_pdcmem->SelectObject(poldbitmap);
Sleep(1);//可有可无
delete m_pBitmap;
delete m_pdcmem;
pWnd->ReleaseDC(pDC);
}while(1);
return 1;
}
这两个线程在单击“采集”按钮时被启用,程序如下:
void CVcMatlabDlg::OnBtnAdoptData() //首先打开断端口,然后再开线程用来采集
{
husb=OpenUA300();
AdoptDataThread=AfxBeginThread(AdoptDataThreadProc,0,THREAD_PRIORITY_NORMAL); DrawThread=AfxBeginThread(DrawThreadProc,0,THREAD_PRIORITY_NORMAL);
GetDlgItem(IDC_BTN_ADOPT)->EnableWindow(FALSE);
GetDlgItem(IDC_BTN_STOP)->EnableWindow(TRUE);
}
用上面的程序采集频率为1000Hz的正弦信号后的图形如图2:
采集到的数据在单击“FFT”按钮时被处理,Matlab Engine就是在这个函数里被调用的,程序如下:
图2 采集到的正弦信号
void CVcMatlabDlg::OnBtnFft()
{
// TODO: Add your control notification handler code here
Engine *ep;
mxArray *X=NULL;
mxArray *Num=NULL;
mxArray *Frequency=NULL;
char buffer[1024];
if (!(ep=engOpen("\0")))
//打开Matlab引擎,建立与本地Matlab的连接
{
AfxMessageBox("Can't start MATLAB engine!");
exit(-1);
}
//参与matlab运算的参数都是矩阵形式,因而要将参数转换为Matlab可接收的矩阵形式
X=mxCreateDoubleMatrix(1,1024,mxREAL);
Num=mxCreateDoubleMatrix(1,1,mxREAL);
Frequency=mxCreateDoubleMatrix(1,1,mxREAL);
double var1=(double)SampNum;//强制转换SamNum和SampFrequency从short到double
double var2=(double)SampFrequency;
memcpy((double*)mxGetPr(Num),(double*)&var1,sizeof(double));
memcpy((double*)mxGetPr(Frequency),(double*)&var2,sizeof(double));
double var3[1024];
for (int i=0;i<1024;i++)
{
var3[i]=(double)(singleaddata[i]/1);
} memcpy((char*)mxGetPr(X),(char*)var3,1024*sizeof(double));
//将转换的参数放入引擎中,可在Matlab commond窗口下察看参数值
engPutVariable(ep, "X", X);
engPutVariable(ep, "Num",Num);
engPutVariable(ep, "Frequency", Frequency);
//下面开始执行MATLAB命令
engEvalString(ep,"pxx=fft(X,Num);");//engEvalString()函数是在vc中执行matlab函数
engOutputBuffer(ep,buffer,1024);
engEvalString(ep,"ff1=-Frequency/2:1/Num*Frequency:Frequency/2-Frequency/Num;");
engEvalString(ep,"subplot(2,1,1);");
engEvalString(ep,"plot(ff1,fftshift(abs(pxx)));");
engEvalString(ep,"title('信号频谱图')");
engEvalString(ep,"xlabel('频率/Hz')");
engEvalString(ep,"ylabel('幅值谱')");
engEvalString(ep,"p=pxx.*conj(pxx)/1024");
engEvalString(ep,"ff2=Frequency*(0:511)/1024;");
engEvalString(ep,"subplot(2,1,2);");
engEvalString(ep,"plot(ff2,p(1:512));");//利用引擎画图
engEvalString(ep,"title('信号功率谱图')");
engEvalString(ep,"xlabel('频率/Hz')");
engEvalString(ep,"ylabel('功率谱密度')");
mxDestroyArray(X);//释放内存
mxDestroyArray(Num);
mxDestroyArray(Frequency);
engClose(ep);//关闭引擎,图片随之关闭
GetDlgItem(IDC_BTN_FFT)->EnableWindow(FALSE);
}
编译运行程序后,对图2所示的正弦信号进行处理,得到的图形如图3:
图3 分析后的频域波形
4 结束语
在Visual C ++中利用多线程编程可以同时采集和显示采集卡采集到的信号,防止采集过程中的掉点;调用Matlab Engine函数可以发挥两者的优势,利用Visual C++强大可视化功能,方便的对采集卡进行控制,利用Matlab Engine强大的数据处理能力,对采集到的数据进行数字信号处理,并可以方便的显示处理后的结果。
参考文献:
[1]David J.Kruglinski著.潘爱民,王国印译. Visual C++技术内幕[M]. 北京:清华大学出版社,1998
[2]董长虹主编. Matlab信号处理及应用[M]. 北京:国防工业出版社,2005.1
[3] 何晓涛,于春田. VC调用Matlab的方法[J].河北科技大学学报,2003,24(1):35-39
[4] 王安红,孙志毅. 一种VC++与Matlab混合编程的实现方法[J]. 计算机应用与软件,2003,20(6):12-13,77
[5]高崇明. VC++6.0与Matlab混合编程技术的原理与实现[J]. 无线电工程,2000,30(2):53-56
[6]肖永韧. VC与Matlab混合编程之DLL实现方法[J]. 计算机工程与应用,2001(13):174-176
[7]4亓波. 实现VC++6.0与Matlab的混合编程[J]. 电脑编程技巧与维护,2000(12):62-64
[8] UA302/H型A/D采集器使用说明。
本文转自博客园知识天地的博客,原文链接:多线程VC++和Matlab混编在信号采集和处理中的应用,如需转载请自行联系原博主。