C++读写.mat文件的方法

最近在编写C++程序的过程中，需要使用matlab生成的.mat文件，于是查找了基本使用方法。 

我的计算机环境是win7 64位系统，VS2010，matlab R2010b。 

一、工程配置： 
1.VC++目录->包含目录 添加: 
MATLAB\R2010b\extern\include 
MATLAB\R2010b\extern\include\win64 
 ->库目录 添加： 
  MATLAB\R2010b\extern\lib\win64\microsoft 
  MATLAB\R2010b\extern\lib\win32\microsoft 
2.C/C++->常规->附加包含目录 添加: 
  MATLAB\R2010b\extern\include 
  MATLAB\R2010b\extern\include\win64 
3.链接器->输入->附加依赖库 添加: 
libmat.lib 
libmx.lib 
libmex.lib 
libeng.lib 

由于所安装的matlab为64位，要调用其函数，需要将工程转换为X64 
4.顶菜单->生成->配置管理器->平台：X64 
5.链接器->高级->目标计算机: 
MachineX64(/MACHINE:X64) 

6.计算机环境变量->path 添加： 
E:\DevTools\MATLAB\R2010b\extern\lib\win64\microsoft; 
E:\DevTools\MATLAB\R2010b\bin\win64; 

二、所用到的基本读取写入mat文件代码示例： 

首先需要包含头文件： 

[cpp] view plaincopy 

    MATFile *pmatFile = NULL;   
    mxArray *pMxArray = NULL;   

    // 读取.mat文件（例：mat文件名为"initUrban.mat"，其中包含"initA"）   
    double *initA;   

    pmatFile = matOpen("initUrban.mat","r");   
    pMxArray = matGetVariable(pmatFile, "initA");   
    initA = (double*) mxGetData(pMxArray);   
    M = mxGetM(pMxArray);   
    N = mxGetN(pMxArray);   
    Matrix<double> A(M,N);   
    for (int i=0; i<M; i++)   
        for (int j=0; j<N; j++)   
            A[i][j] = initA[M*j+i];   

    matClose(pmatFile);   
    mxFree(initA);   

    // 生成.mat文件   
    double *outA = new double[M*N];   
        for (int i=0; i<M; i++)   
            for (int j=0; j<N; j++)   
                outA[M*j+i] = A[i][j];   
    pmatFile = matOpen("A.mat","w");   
    mxSetData(pMxArray, outA);   
    matPutVariable(pmatFile, "A", pMxArray);   
    matClose(pmatFile);   

</div>

三、程序说明 

1.使用matOpen函数打开mat文件 

MATFile *matOpen(const char *filename,const char *mode) 

mode： 
r:以只读方式打开 
u:更新模式，可读可写，但是如果要打开的数据文件不存在，也不会创建新的文件 
w:以写入方式打开，只能写如入，如果要打开的数据文件不存在，则创建一个新的文件 

2.使用matGetVariable函数读取mat文件中的变量 

mxArray * matGetVariable(MATFile * pMF, const char * name); 

读取名为name的变量，返回一个数据阵列指针 

3.使用mxGetData函数获取数据阵列中的数据 

void *mxGetData(const mxArray *pa ); 

返回时需要使用强制类型转换。 

4.使用mxGetM和mxGetN函数获取数据阵列矩阵的维数 

size_t mxGetM(const mxArray *pa); 

size_t mxGetN(const mxArray *pa); 

5.使用mxSetData函数将变量存入数据阵列 

void mxSetData( mxArray *pa, void  *newdata ); 

6.使用matPutVariable函数将数据阵列存入到mat文件中 

int matPutVariable(MATFile * pMF, const char * name, const mxArray * pA); 

存入成功返回0，存入出错返回非0 

7.由于malab中矩阵的存储是按列存储，与c中不同，所以从mat文件中获取的矩阵需要重新排列。生成mat文件时，同样要注意。