• C++调用MATLAB函数


    MATLAB具有着很高的计算性能,一些算法用MATLAB很容易实现,而用C++很难实现,如果想在C++中调用MATLAB编写的函数,可以将该函数编译成库文件,之后在C++中对其进行调用。

    本文详细讲解如何将MATLAB函数编译成库,并在C++中进行调用。方法可行,本人亲测。

    第一步、在MATLAB中的命令行窗口中输入mex -setup,会出现如下界面。

    第二步、点击界面中的“mex -setup C++”,出现如下提示。

    第三步、之后就是进行对MATLAB中的function的编译了。这里,我的函数名为phasecong,实现的是图像相位一致性的计算。如下图所示。

    这个函数是一位大牛学者写的,原本是有多个输入和多个输出的,但是为了编译方便和后续调用方便,我把其他几个参数删掉了,只保留了一个输入矩阵和一个输出矩阵。

    这样做的好处是可以减少编译和调用出错的概率。编译方法见第四步。

    第四步、下面就是对function进行编译了,编译的命令为:mcc -W cpplib:pc -T link:lib phasecong

    这条命令中pc表示的是想要生成的库的名称,可以自己定,最后面的phasecong表示的是要编译的函数的名,注意函数名后面不要加.m后缀。其他的参数都不用改。

    如下图所示。

    第五步、输入完上面的命令后就是等待MATLAB编译完成了,我当时编译好像用了30秒左右的时间。

    编译完成后就会出现编译成功的提示,如下图所示。

    需要注意的是,在编译时,当前文件夹一定要为function所在的文件夹。

    第六步、经过上面的编译后,在该文件夹下就会生成.h .cpp .lib .dll这四个文件,因为我设置的输出的库名称为pc,因此我这里生成pc.h pc.cpp pc.lib pc.dll这四个文件。

    第七步、接下来就是将生成的这四个文件加入到我们的C++工程中了。

    将这四个文件放在C++工程的文件夹中,最好在该文件夹中新建两个文件夹,分别命名为lib和include,因为后面我们还需要引入其他的库文件和头文件,这样做可以让文件更整齐。然后将pc.lib放入lib文件夹,将pc.h放入include文件夹。然后在VS中右键项目,点击属性,在属性中选择“配置属性”->“VC++目录”->“包含目录”,将刚才建立的include文件夹包含进去,如下图所示。其他的属性都不用改。

    这样我们就将头文件包含了进去,之后进行加载lib文件,这里我们使用静态加载方式,在工程文件中加入这样一条语句:

    #pragma comm#pragma comment(lib, ".\lib\pc.lib" )ent(lib, ".\lib\pc.lib" )

    对于pc.dll和pc.cpp文件,值需要将其放在工程目录下就可以,不用进行配置。
    这样就完成了对这四个文件的配置。

    第八步、配置好生成的这四个文件后就可以运行C++的程序了,但是点击运行后应该会报错(通常一定会报错),不要怕,这是因为我们前面生成的这四个文件是需要依赖其他的库和头文件才能工作的。我们根据报错的信息可以知道缺少的头文件,例如下图这个错误:

    通过这个错误我们就可以知道,我们的项目中缺少mat.h这个头文件,这时我们就要找到这个头文件并将其加入到include文件夹中。到哪里去找呢,在这个文件夹中“D:Program FilesMATLABR2014bexterninclude”。重复这一操作,直到不再提示这种错误。这一操作需要重复很多次,像我这个项目需要加入下图这么多的头文件。

    之后还需要加入额外的lib库文件,具体加入哪些lib文件,需要根据上面的这些头文件而定,我们到文件夹“D:Program FilesMATLABR2014bexternlibwin64microsoft”中(因为我的MATLAB是64位的,所以在win64文件夹下,如果大家的MATLAB是32位的,那么估计是在win32文件夹下,大家自己找一下),寻找上面这些头文件同名的lib文件,不过并不是每个头文件都能找到同名的lib文件,不过不要紧,我们只要把能找到的lib文件复制到工程目录下的lib文件夹中即可,除了这些lib文件之外,还需要加上“libmat.lib”和“libmx.lib”这两个文件,这两个文件好像包括着必要的类型的定义之类的吧。搞定之后我的lib文件夹中有下图这些个文件。

    这些lib的加载方式与之前我们说的pc.lib的加载方式相同,如下图所示。

    进行到此,我们所有的配置就都搞定了,如果没有意外的话,C++的项目编译就可以通过了,

    但是虽然编译通过了,在运行的过程中还是有可能提示缺少某个库文件,我们只需要采用同样的方法把那些库文件加到项目里即可,比如我这个程序还提示我缺少mclmcr.dll这个文件,我就把这个文件复制到工程目录下就可以了,所需的dll文件通常都在“D:Program FilesMATLABR2014binwin64”这个文件夹可以找到。

    至此我们就完成了MATLAB函数的编译和C++调用的整个过程。

    最后是几点提示:

    1、我这里使用的是MATLAB2014a和vs2013。

    2、如果使用的MATLAB是64位的版本,那么C++的程序也要在x64下配置!!!!!!!

          如果使用的MATLAB是32位的版本,那么C++的程序也要在win32下配置!!!!!!

          这点非常重要!!!!!!!!

    3、MATLAB函数的输入和输出参数通常都是矩阵,在C++中矩阵变量使用mwArray这个类型,这个类的使用方法可以看我下面贴的代码,也可以在网上找教程。

    4、在C++程序中,我们调用这个函数的名称与MATLAB中定义的函数名称相同,这个函数是在MATLAB编译后生成的头文件中进行的声明,例如我这里是在pc.h中。

    5、在C++程序中,在定义输入数组和输出数组时,最好使用动态数组,因为传统的数组是在栈中分配内存,这样会造成内存溢出。

    最后为了更直观的展示在C++中如何对这些个库文件进行的调用以及如何使用这个函数,我在这里贴上了C++程序的源代码。

     1 // pcproj.cpp : 定义控制台应用程序的入口点。
     2 //
     3  
     4 #include "stdafx.h"
     5 #include "pc.h"
     6 #include <iostream>
     7 #include <fstream>
     8  
     9 using namespace std; 
    10  
    11 #pragma comment(lib, ".\lib\pc.lib" )
    12 #pragma comment(lib, ".\lib\mclmcrrt.lib")
    13 #pragma comment(lib, ".\lib\libmx.lib")
    14 #pragma comment(lib, ".\lib\libmat.lib")
    15 #pragma comment(lib, ".\lib\mclmcr.lib")
    16  
    17 #define WIDTH 512
    18 #define HEIGHT 512
    19  
    20 int _tmain(int argc, _TCHAR* argv[])
    21 {
    22     //mclInitializeApplication(NULL, 0);
    23     //mclmcrInitialize();
    24  
    25     UINT8 *img_befor = new UINT8[WIDTH*HEIGHT];
    26     DOUBLE *img_after = new DOUBLE[WIDTH*HEIGHT];
    27  
    28     int data;
    29     double ddata;
    30     FILE *fp;
    31  
    32     errno_t err = fopen_s(&fp, "pixel.txt", "r");
    33     if (err!=0){
    34         cout << "open pic fail!" << endl;
    35         return 0;
    36     }
    37  
    38     for (int i = 0; i < WIDTH; i++){
    39         for (int j = 0; j < HEIGHT; j++){
    40             fscanf_s(fp, "%d", &data);
    41             img_befor[i*WIDTH + j] = data;
    42         }
    43     }
    44  
    45     fclose(fp);
    46     fp = NULL;
    47  
    48     pcInitialize();
    49  
    50     mwArray img_input_array(HEIGHT, WIDTH, mxUINT8_CLASS, mxREAL);
    51     mwArray img_output_array(HEIGHT, WIDTH, mxDOUBLE_CLASS, mxREAL);
    52     int nargout = 1;
    53  
    54     img_input_array.SetData(img_befor, HEIGHT*WIDTH);
    55  
    56     phasecong(nargout,img_output_array,img_input_array);
    57     
    58     img_output_array.GetData(img_after, HEIGHT*WIDTH);
    59  
    60     pcTerminate();
    61  
    62     //将结果写到文件
    63     err = fopen_s(&fp, "pixel_after.txt", "w");
    64     if (err != 0){
    65         cout << "open pic_after fail!" << endl;
    66     }
    67     else{
    68         //写文件
    69         for (int i = 0; i < WIDTH; i++){
    70             for (int j = 0; j < HEIGHT; j++){
    71                 ddata = img_after[i*WIDTH + j];
    72                 fprintf(fp, "%.4f ", ddata);
    73             }
    74             fprintf(fp, "
    ");
    75         }
    76     }
    77     fclose(fp);
    78     fp = NULL;
    79  
    80  
    81     return 0;
    82 }

  • 相关阅读:
    JSTL标签用法
    (转)javascript正则实战
    Spring中几种数据源的配置
    (转)常用正则表达式
    EL表达式详解
    JAVA环境变量配置(转)
    S2SH整合步骤
    正则表达式30分钟教程(转)
    在Flash Player 10.2中使用原生鼠标指针
    time
  • 原文地址:https://www.cnblogs.com/ybqjymy/p/14535459.html
Copyright © 2020-2023  润新知