Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

AaronSpark 发表于 2005-7-26 21:27

[转帖]*.mat文件的内容，如何分析？

发信人: FangQ (光芒之神~等Offer中), 信区: MathTools 标题: Re: *.mat文件的内容，如何分析？ 发信站: BBS 水木清华站 (Wed Apr 19 22:25:01 2000) 好想是matlab的C库里面的吧，下面两个函数完成mat文件的d的读入和保存 /* * loadmat - C language routine to load a matrix from a MAT-file. * * Here is an example that uses 'loadmat' to load a matrix from a MAT-file: * * FILE *fp; * char name; * int type, mrows, ncols, imagf; * double *xr, *xi; * fp = fopen("foo.mat","rb"); * loadmat(fp, &type, name, &mrows, &ncols, &imagf, &xr, &xi); * fclose(fp); * free(xr); * if (imagf) free(xi); * * The 'loadmat' routine returns 0 if the read is successful and 1 if * and end-of-file or read error is encountered.'loadmat' can be called * repeatedly until the file is exhausted and an EOF is encountered. * * Author J.N. Little 11-3-86 * * Modified by Steffen Torp, Servolaboratoriet 14-10-92. */ int loadmat(FILE *fp,int *type,char *pname,int *mrows,int *ncols,int *imagf, double **preal,double **pimag) /* FILE *fp; File pointer */ /* int *type; Type flag: see reference section of guide */ /* int *mrows; row dimension */ /* int *ncols; column dimension */ /* int *imagf; imaginary flag */ /* char *pname; pointer to matrix name */ /* double **preal;pointer to real data */ /* double **pimag;pointer to imag data */ { Fmatrix x; int mn, namlen, conv; /* * Get Fmatrix structure from file */ if (fread((char *)&x, sizeof(Fmatrix), 1, fp) != 1) { return(1); } /* Check if the mat-file is created on a platform similar to */ /* the actual platform. If yes, the flag 'conv' will be FALSE*/ /* Otherwise it will be TRUE */ if (SYSTEM == 0) { switch (x.type) { case 0L: conv = FALSE; break; case 1L: conv = FALSE; break; case 0xE8030000L: conv = TRUE;break; case 0xE9030000L: conv = TRUE;break; default: { printf("\nUnknown file type in loadmat!\n"); printf("The type is: %lx \n",x.type); exit(1); } } } else if (SYSTEM == 1000) { switch (x.type) { case 0L: conv = TRUE;break; case 0x1000000L: conv = TRUE;break; case 1000L: conv = FALSE; break; case 1001L: conv = FALSE; break; default: { printf("\nUnknown file type in loadmat!\n"); printf("The type is: %lx \n",x.type); exit(1); } } } if(conv==TRUE) { *type= (int)longconv(x.type); *mrows = (int)longconv(x.mrows); *ncols = (int)longconv(x.ncols); *imagf = (int)longconv(x.imagf); namlen = (int)longconv(x.namlen); mn = (*mrows)*(*ncols); } else { *type= (int)(x.type); *mrows = (int)(x.mrows); *ncols = (int)(x.ncols); *imagf = (int)(x.imagf); namlen = (int)(x.namlen); mn = (*mrows)*(*ncols); } /* * Get matrix name from file */ if (fread(pname, sizeof(char), namlen, fp) != namlen) { return(1); } /* * Get Real part of matrix from file */ if (mn) { *preal = (double *)malloc(mn*sizeof(double)); if (!(*preal)) { printf("\nError: Variable too big to load\n"); return(1); } if (fread(*preal, sizeof(double), mn, fp) != mn) { free(*preal); return(1); } if(conv==TRUE) matconv(mn,*preal); } /* * Get Imag part of matrix from file, if it exists */ if (x.imagf) { *pimag = (double *)malloc(mn*sizeof(double)); if (!(*pimag)) { printf("\nError: Variable too big to load\n"); free(*preal); return(1); } if (fread(*pimag, sizeof(double), mn, fp) != mn) { free(*pimag); free(*preal); return(1); } if(conv==TRUE) matconv(mn,*pimag); } return(0); } /* * savemat - C language routine to save a matrix in a MAT-file. * * Here is an example that uses 'savemat' to save two matrices to disk, * the second of which is complex: * * FILE *fp; * double xyz, ar, ai; * fp = fopen("foo.mat","wb"); * savemat(fp, 2000, "xyz", 2, 3, 0, xyz, (double *)0); * savemat(fp, 2000, "a", 5, 5, 1, ar, ai); * fclose(fp); * * Author J.N. Little 11-3-86 */ void savemat(FILE *fp,int type,char *pname,int mrows,int ncols,int imagf,dou ble *preal,double *pimag) /* FILE *fp; File pointer */ /* int type; Type flag: Normally 0 for PC, 1000 for Sun, Mac, and*/ /* Apollo, 2000 for VAX D-float, 3000 for VAX G-float */ /* Add 1 for text variables. */ /* See LOAD in reference section of guide for more info. */ /* int mrows; row dimension */ /* int ncols; column dimension */ /* int imagf; imaginary flag */ /* char *pname; pointer to matrix name */ /* double *preal; pointer to real data */ /* double *pimag; pointer to imag data */ { Fmatrix x; int mn; x.type = (long)type; x.mrows= (long)mrows; x.ncols= (long)ncols; x.imagf= (long)imagf; x.namlen = (long)strlen(pname) + 1L; mn = (int)(x.mrows * x.ncols); fwrite(&x, sizeof(Fmatrix), 1, fp); fwrite(pname, sizeof(char), (int)x.namlen, fp); fwrite(preal, sizeof(double), mn, fp); if (imagf) { fwrite(pimag, sizeof(double), mn, fp); } }

页: [1]

声振论坛's Archiver

[转帖]*.mat文件的内容，如何分析？