c对mat文件的几个操作，打开mat文件 读取mat文件下的阵列，生成mat文件

#include "stdafx.h"
#include <stdlib.h>
#include "mat.h"
#include <iostream>

using namespace std;
int main()
{

    MATFile *pmat=NULL;
    MATFile *pmatFile = NULL;
    const char **dir;
    const char *file;
    //FILE *pp;
    const char *name;
    int         ndir;
    mxArray *pa;
    mxArray *pname;
    file = "E://alogrithmprogramming//C++//cmatdll//cmatdll//Y.mat";//双反斜杠防止转义  
    pmat = matOpen(file, "r");//打开文件，返回指向文件指针  
    //pp = matGetFp(pmat);
    if (pmat == NULL) {
        cout << "Error opening file" << endl;
        return 1;
    }

    dir = (const char **)matGetDir(pmat, &ndir); //ndir 表示mat文件中含有矩阵数目  
    if (dir == NULL) {
            cout << "Error reading directory of file" << endl;
            return(1);
    }

    else
    {
        cout << "opening file" << endl;
        cout << "Mat NUM: " << ndir << endl;
        for (int i = 0; i < ndir; i++)//mat文件下多个阵列名
            cout << dir[i] << endl;
    }

    pmat = matOpen(file, "r");//重新打开文件;  
    pa = matGetNextVariable(pmat, &name);//返回指向文件头文件信息的指针，指针类型为*mxArray;  
    cout << name << endl;//name 是矩阵的名字;  
    //pname=matGetVariable(pmat,name);// Read the array value for the specified variable name from a MAT-file.
    //cout << pname << endl;
    int ii = mxGetM(pa);//矩阵行数
    int jj = mxGetN(pa);//矩阵列数  
    cout << ii << "," << jj << endl;
    //double *pm = (double*)mxGetData(pa);//获取矩阵数值，返回指向矩阵第一数值的指针;
    //cout << *(pm++)<< endl;
    //mxSetData(pa, pm);

    // 生成.mat文件   
    pmatFile = matOpen("A.mat", "w");
    matPutVariable(pmatFile, "A", pa);
    matClose(pmatFile);

    mxFree(dir);
    matClose(pmat);
    system("pause");
    return 0;
}