精通VC与MATLAB联合编程——综合举例二

【 tulaoshi.com - C语言心得技巧 】

精通VC与MATLAB联合编程<十>——综合举例二

作者：邓科

本文使用的工具：VC6.0+MATLAB6.5

　　在本章中实现的是在VC中调用MATLAB神经网络工具箱，文章中用到的程序是一篇硕士论文的一部分，因此没有提供源程序，目的是通过这个例子介绍如何在VC中调用利用com组件，实现VC和MATLAB数据的交换。首先在MATLAB中完成两个神经网络函数，一个用来检验，一个用来预测。

1、检验函数：

function [w1,b1,w2,b2,ET,a2]=SJHSJY(P,T,me,eg,lr,P2,T2,S1)[Pn,meanp,stdp,Tn,meant,stdt]=prestd(P,T);   %预处理[w1,b1,w2,b2]=initff(Pn,S1,''tansig'',Tn,''purelin'');  %数据初始化df=10;   tp=[df me eg lr];[w1,b1,w2,b2,tp,tr]=trainbpx(w1,b1,''tansig''); %样本训练        P2n=trastd(P2,meanp,stdp);   ......后处理   %必须为非负

2、预测函数：

function [a2]=SJHSYC(P,T,P2,me,eg,lr,S1)[Pn,meanp,stdp,Tn,meant,stdt]=prestd(P,T); [w1,b1,w2,b2]=initff(Pn,S1,''tansig'',Tn,''purelin'');........参数输出

3、在MATLAB中将该两个函数做成COM组件，具体方法可以参考浅析VC与MATLAB联合编程<三>。

4、新建基于对话框的VC工程，在VC中引入COM组件，具体方法可以参考浅析VC与MATLAB联合编程<三>,完成后可以在CLASSView中看到多了两个新类ISJHSJY和ISJHSYC，如图1所示，它们分别实现函数SJHSJY.m和SJHSYC.m的功能。


图1

仔细查看类ISJHSJY的参数，如图2所示


图2

在类ISJHSJY中 有调用函数格式说明：

 ISJHSJY __RPC_FAR * This,            /* [in] */ long nargout,            /* [out][in] */ VARIANT __RPC_FAR *w1,            /* [out][in] */ VARIANT __RPC_FAR *b1,            /* [out][in] */ VARIANT __RPC_FAR *w2,            /* [out][in] */ VARIANT __RPC_FAR *b2,            /* [out][in] */ VARIANT __RPC_FAR *ET,            /* [out][in] */ VARIANT __RPC_FAR *a2,            /* [in] */ VARIANT P,            /* [in] */ VARIANT T,            /* [in] */ VARIANT me,            /* [in] */ VARIANT eg,            /* [in] */ VARIANT lr,            /* [in] */ VARIANT P2,            /* [in] */ VARIANT T2,            /* [in] */ VARIANT S1);

5、关键代码：

SAFEARRAYBOUND rgsabound1[2];rgsabound1[0].lLbound=0;rgsabound1[1].lLbound=0;rgsabound1[0].cElements=4;    rgsabound1[1].cElements=6;VariantInit(&TTT);TTT.vt=VT_R8|VT_ARRAY;TTT.parray=SafeArrayCreate(VT_R8,2,rgsabound1);TTT.parray->pvData=HH2;

6、调用代码：

if(FAILED(CoInitialize(NULL)))//COM初始化{MessageBox("Error");}ISJHSJY *pclass=NULL;HRESULT hr=CoCreateInstance(CLSID_SJHSJY,NULL, CLSCTX_ALL, IID_ISJHSJY, (void **)&pclass);if(SUCCEEDED(hr)){long int u=6;HRESULT hr1=pclass->sjhsjy(u,&w1,&b1,&w2,&b2,&et,&a2,PP,TTT,me,eg,lr,P2,T2,S1);//函数调用   if(!SUCCEEDED(hr1)){MessageBox("请输入正确数据");return;}//函数调用失败   if(SUCCEEDED(hr1))   {   memcpy(W1,w1.parray->pvData,420*sizeof(double));//将Variant型结果数据转换成double型   memcpy(W2,w2.parray->pvData,240*sizeof(double));   memcpy(B1,b1.parray->pvData,60*sizeof(double));   memcpy(B2,b2.parray->pvData,4*sizeof(double));   memcpy(ET,et.parray->pvData,4*sizeof(double));   memcpy(A2,a2.parray->pvData,4*sizeof(double));      m_edit7.Format("%f rn %f rn %fr