外部インタフェース/API

MAT-ファイルの例

この節では、MAT-ファイルの書き出し、読み込み、診断のCおよびFortranの例題を説明します。含まれるトピックスは以下の通りです。

• MAT-ファイルをCで作成
• MAT-ファイルをCで読み込む
• MAT-ファイルをFortranで作成
• MAT-ファイルをFortranで読み込む

MAT-ファイルをCで作成

この例題プログラムは、MATLABにロードされるMAT-ファイルを作成するライブラリルーチンの使用法を記述しています。プログラムは、読み込みまたは書き出しの失敗に対するMAT-ファンクション呼び出しの返り値のチェック方法も示しています。

/*
 * MAT-file creation program
 *
 * See the MATLAB API Guide for compiling information.
 *
 * Calling syntax:
 *
 *   matcreat
 *
 * Create a MAT-file which can be loaded into MATLAB.
 *
 * This program demonstrates the use of the following functions:
 *
 *  matClose
 *  matGetArray
 *  matOpen
 *  matPutArray
 *  matPutArrayAsGlobal
 *
 * Copyright 1984-2000 The MathWorks, Inc.
 */
/* $ Revision: 1.10 $ */
#include <stdio.h>
#include <string.h> /* For strcmp() */
#include <stdlib.h> /* For EXIT_FAILURE, EXIT_SUCCESS */
#include "mat.h"

#define BUFSIZE 256

int main() {
  MATFile *pmat;
  mxArray *pa1, *pa2, *pa3;
  double data[9] = { 1.0, 4.0, 7.0, 2.0, 5.0, 8.0, 3.0, 6.0, 9.0 };
  const char *file = "mattest.mat";
  char str[BUFSIZE];
  int status; 

  printf("Creating file %s...\n\n", file);
  pmat = matOpen(file, "w");
  if (pmat == NULL) {
    printf("Error creating file %s\n", file);
    printf("(Do you have write permission in this directory?)\n");
    return(EXIT_FAILURE);
  }

  pa1 = mxCreateDoubleMatrix(3,3,mxREAL);
  if (pa1 == NULL) {
      printf("%s : Out of memory on line %d\n",
              __FILE__, __LINE__); 
      printf("Unable to create mxArray.\n");
      return(EXIT_FAILURE);
  }
  mxSetName(pa1, "LocalDouble");

  pa2 = mxCreateDoubleMatrix(3,3,mxREAL);
  if (pa2 == NULL) {
      printf("%s : Out of memory on line %d\n",
              __FILE__, __LINE__);
      printf("Unable to create mxArray.\n");
      return(EXIT_FAILURE);
  }
  mxSetName(pa2, "GlobalDouble");
  memcpy((void *)(mxGetPr(pa2)), (void *)data, sizeof(data));
  
  pa3 = mxCreateString("MATLAB: the language of technical
                         computing");
  if (pa3 == NULL) {
      printf("%s :  Out of memory on line %d\n",
              __FILE__, __LINE__);
      printf("Unable to create string mxArray.\n");
      return(EXIT_FAILURE);
  }
  mxSetName(pa3, "LocalString");

  status = matPutArray(pmat, pa1);
  if (status != 0) {
      printf("%s :  Error using matPutArray on line %d\n",
              __FILE__, __LINE__);
      return(EXIT_FAILURE);
  }  
  
  status = matPutArrayAsGlobal(pmat, pa2);
  if (status != 0) {
      printf("Error using matPutArrayAsGlobal\n");
      return(EXIT_FAILURE);
  } 
  
  status = matPutArray(pmat, pa3);
  if (status != 0) {
      printf("%s :  Error using matPutArray on line %d\n",
              __FILE__, __LINE__);
      return(EXIT_FAILURE);
  } 
  
  /*
   * Ooops! we need to copy data before writing the array.  (Well,
   * ok, this was really intentional.) This demonstrates that
   * matPutArray will overwrite an existing array in a MAT-file.
   */
  memcpy((void *)(mxGetPr(pa1)), (void *)data, sizeof(data));
  status = matPutArray(pmat, pa1);
  if (status != 0) {
      printf("%s :  Error using matPutArray on line %d\n",
              __FILE__, __LINE__);
      return(EXIT_FAILURE);
  } 
  
  /* clean up */
  mxDestroyArray(pa1);
  mxDestroyArray(pa2);
  mxDestroyArray(pa3);

  if (matClose(pmat) != 0) {
    printf("Error closing file %s\n",file);
    return(EXIT_FAILURE);
  }

  /*
   * Re-open file and verify its contents with matGetArray
   */
  pmat = matOpen(file, "r");
  if (pmat == NULL) {
    printf("Error reopening file %s\n", file);
    return(EXIT_FAILURE);
  }

  /*
   * Read in each array we just wrote
   */
  pa1 = matGetArray(pmat, "LocalDouble");
  if (pa1 == NULL) {
    printf("Error reading existing matrix LocalDouble\n");
    return(EXIT_FAILURE);
  }
  if (mxGetNumberOfDimensions(pa1) != 2) {
    printf("Error saving matrix: result does not have two
            dimensions\n");
    return(EXIT_FAILURE);
  }

  pa2 = matGetArray(pmat, "GlobalDouble");
  if (pa2 == NULL) {
    printf("Error reading existing matrix GlobalDouble\n");
    return(EXIT_FAILURE);
  }
  if (!(mxIsFromGlobalWS(pa2))) {
    printf("Error saving global matrix: result is not global\n");
    return(EXIT_FAILURE);
  }

  pa3 = matGetArray(pmat, "LocalString");
  if (pa3 == NULL) {
    printf("Error reading existing matrix LocalDouble\n");
    return(EXIT_FAILURE);
  }
  
  status = mxGetString(pa3, str, sizeof(str));
  if(status != 0) {
      printf("Not enough space. String is truncated.");
      return(EXIT_FAILURE);
  }
  if (strcmp(str, "MATLAB: the language of technical
                    computing")) {
    printf("Error saving string: result has incorrect
            contents\n");
    return(EXIT_FAILURE);
  }

  /* clean up before exit */
  mxDestroyArray(pa1);
  mxDestroyArray(pa2);
  mxDestroyArray(pa3);

  if (matClose(pmat) != 0) {
    printf("Error closing file %s\n",file);
    return(EXIT_FAILURE);
  }
  printf("Done\n");
  return(EXIT_SUCCESS);
}

この例題プログラムの実行可能なバージョンを作成するには、ファイルをコンパイルし適切なライブラリをリンクしてください。種々のプラットフォームでのMAT-ファイルプログラムのコンパイルとリンクの方法の詳細は、「MAT-ファイルプログラムのコンパイルとリンク」に記述されています。

MAT-ファイルプログラムをコンパイルしリンクすると、作成したスタンドアロンアプリケーションを実行できます。このプログラムは、MATLABにロード可能なMAT-ファイルmattest.matを作成します。アプリケーションを実行するには、プラットフォームによって、アイコンをダブルクリックするか、システムプロンプトでmatcreatと入力してください。

matcreat
Creating file mattest.mat...

MAT-ファイルが作成されたことを確認するには、MATLABプロンプトでつぎを入力します。

whos -file mattest.mat
  Name               Size         Bytes  Class

  GlobalDouble       3x3             72  double array (global)
  LocalDouble        3x3             72  double array
  LocalString        1x43            86  char array

Grand total is 61 elements using 230 bytes

関連するファイル

MAT-ファイルをCで読み込む