摘要:MATLAB对于矩阵处理是非常高效的,而C++对于矩阵操作是非常麻烦的,因而可以采用C++与MATLAB混合编程求解矩阵问题。
主要思路就是,在MATLAB中编写函数脚本并使用C++编译为dll文件(在C++中可以调用编译的函数),然后对VS项目进行文件配置,编写C++代码调用MATLAB中定义的函数。
问题描述:对于一个多项式

需要求解c0到c5的值,由相关条件已知c0=c1=0,且
、
、
、

可得如下矩阵式:

对比类似AX=B,可求X=A\B。
1、写出MATLAB代码如下

运行结果:

2、使用C++编译器将ComputeC2_C5.m文件编译成C++可以调用的.dll文件,步骤:
2.1、在MATLAB控制栏输入"mex -setup"

2.2、点击"mex -setup C++",使用C++进行编译

2.3、输入"deploytool",弹出编译(Compiler)窗口,选择Library Compiler(库编译)。

2.4、进入MATLAB Compiler窗口,按照如下步骤,1)选择C++ Shared Library,2)点击"+"号按钮添加将要编译的.m文件,3)填写Library的信息,默认自动填写脚本文件名称,4)点击Package"√"按钮,进行编译打包处理。

2.5、选择一个文件夹存放打包生成的文件,默认和脚本文件名相同的.prj文件。

2.6、打包处理过程(大概需要几十秒)

2.7、生成后自动存放在与脚本同名的文件夹下,如下图所示,第一个文件夹中的exe文件用于程序发布时用,第二个文件夹中包含相当于release下得.h,.lib和 .dll文件,第三个文件包含测试用地.h,.lib和 .dll文件等文件,是测试需要的;至此MATLAB生成dll文件完成。

3、新建C++项目,并将打包生成的文件导入C++项目。
3.1、将for_testing文件夹中的.h、.lib、.dll文件复制到项目中并和.h、.cpp、.cu(CUDA项目文件)放在一起


3.2、并将.h(ComputeC2_C5.h)文件添加到项目中



3.3、配置电脑系统环境变量和VS工程属性
1)在系统环境变量的Path中添加"Matlab安装目录\bin\win64;"

2)在项目属性->配置属性->VC++ 目录->包含目录,添加"Matlab安装目录\extern\include"

3)在项目属性->配置属性->VC++ 目录->库目录,添加"Matlab安装目录\extern\include\lib\win64\microsoft"

4)在项目属性->配置属性->链接器->输入->附加依赖项,添加
" libmat.lib
libmex.lib
mclmcr.lib
mclmcrrt.lib
libemlrt.lib
libeng.lib
libfixedpoint.lib
libcovrt.lib "

4、编写代码

代码(.cu文件是CUDA项目的文件,和.cpp文件类似,一样可以使用):
#pragma once
#pragma comment(lib,"mclmcrrt.lib") //MATLAB中的库文件
#pragma comment(lib,"ComputeC2_C5.lib") //必须要链接ComputeC2_C5.lib库才能正常编译运行#include"mclcppclass.h" //使用mwArray类的头文件
#include "ComputeC2_C5.h" //MATLAB生成的头文件
void ComputeC(doubleamplitude)
{if (!ComputeC2_C5Initialize())//初始化; 注意这个初始化中 !ComputeC2_C5Initialize() 字段中的 ComputeC2_C5与matlab编译成的.lib同名
{
cout<< "Could not initialize addfunction!";
exit(0);
}
mwArray X(4, 1, mxDOUBLE_CLASS, mxREAL); //申请4行1列的矩阵数组
mwArray A(1, 1, mxDOUBLE_CLASS, mxREAL); //申请1行1列的矩阵数组
double *C2ToC5 = new double[4]; //保存X矩阵中的元素
A.SetData(&litude, 1); //设置参数
ComputeC2_C5(1, X, A); //调用函数
X.GetData(C2ToC5,4); //将函数返回值保存在C2ToC5数组中
cout <
{
cout<< "c" << i + 2 << "=" << C2ToC5[i] <
}
ComputeC2_C5Terminate();//注意结束中 !ComputeC2_C5Terminate() 字段中的 ComputeC2_C5也与matlab编译成的.lib同名
}intmain()
{
ComputeC(10);
system("pause");return 0;
}
5、结果(上边为VS项目运行结果,下边为MATLAB运行结果,结果一致):

6、常见问题:
6.1、MATLAB和VS项目需要位数相同,同为32位或64位。
6.2、编译通过,运行提示"mclmcr.lib",因为系统环境变量没有配置
6.3、编译通过,运行不出错但是窗口一闪而过,则是没有设置初始化和结束语句。即

