调用Matlab对矩阵做奇异值分解的程序源代码.docxVIP

??

?

??

调用Matlab对矩阵做奇异值分解的程序源代码

?

??

?

?

?

?

?

?

?

???

?

?

?

?

?

为计算针对多目标节点集的相对连通系数需要对矩阵进行奇异值分解，为此可直接调用Matlab里的对应矩阵计算模块。从VC++调用Matlab进行计算有两种方法：一种是通过Matlab引擎直接调用Matlab，另一种是将Matlab函数编译为dll在VC中使用。

这里采用第二种方法，通过Matlab引擎调用矩阵计算功能。其实现代码如下。

voidConnectivitys：：svd_Matlab（void*destMatrixP）

{

inti，j；

double*vP=NULL，*sP=NULL；

mxArray*Xin，*Vout，*Sout；

Engine*ep；

//显式加载Matlab动态链接库，定义其中的函数

HINSTANCEhInst_libeng=LoadLibrary（“libeng.dll”）；

HINSTANCEhInst_libmat=LoadLibrary（“libmat.dll”）；

HINSTANCEhInst_libmx=LoadLibrary（“libmx.dll”）；

typedefEngine*（*_EngOpen）（constchar*startcmd）；

_EngOpenengOpen=（_EngOpen）GetProcAddress（hInst_libeng，“engOpen”）；

typedefint（*_EngClose）（Engine*ep）；

_EngCloseengClose=（_EngClose）GetProcAddress（hInst_libeng，“engClose”）；

typedefmxArray*（*_EngGetArray）（Engine*ep，constchar*name）；

_EngGetArrayengGetArray=（_EngGetArray）GetProcAddress（hInst_libeng，“engGetArray”）；

typedefint（*_EngPutArray）（Engine*ep，constmxArray*ap）；

_EngPutArrayengPutArray=（_EngPutArray）GetProcAddress（hInst_libeng，“engPutArray”）；

typedefint（*_EngEvalString）（Engine*ep，constchar*string）；

_EngEvalStringengEvalString=（_EngEvalString）GetProcAddress（hInst_libeng，“engEvalString”）；

typedefmxArray*（*_MxCreateDoubleMatrix）（intm，intn，mxComplexityflag）；

_MxCreateDoubleMatrixmxCreateDoubleMatrix=（_MxCreateDoubleMatrix）

GetProcAddress（hInst_libmx，“mxCreateDoubleMatrix”）；

typedefvoid（*_MxDestroyArray）（mxArray*pa）；

_MxDestroyArraymxDestroyArray=（_MxDestroyArray）

GetProcAddress（hInst_libmx，“mxDestroyArray”）；

typedefvoid（*_MxSetName）（mxArray*pa，constchar*s）；

_MxSetNamemxSetName=（_MxSetName）GetProcAddress（hInst_libmx，“mxSetName”）；

typedefdouble*（*_MxGetPr）（constmxArray*pa）；

_MxGetPrmxGetPr=（_MxGetPr）GetProcAddress（hInst_libmx，“mxGetPr”）；

//调用MATLAB引擎

if（！（ep=engOpen（NULL）））{

AfxMessageBox（_T（“Can’tstartMATLABengine”））；

exit（-1）；

}

//将平均偏差矩阵载入Matlab

Xin=mxCreateDoubleMatrix（des