MATLAB外部接口编程(张威)第4章.pptVIP

图4-5 设置可执行文件的对话框 4.4 MEX 文 件 示 例 　　 例4-3 创建稀疏矩阵。 　　首先向读者演示一个创建稀疏矩阵的例子，该例子其实是第2章创建稀疏矩阵C语言例子的翻版，这里使用Fortran语言将其重新写了一遍，通过这个例子重点了解一下mxCopy函数的具体使用方法，下面是该例子的源代码： 001 C fcreatesparse.f 002 C 创建稀疏矩阵 003 C 入口函数 004 subroutine mexFunction(nlhs, plhs, nrhs, prhs) 005 C------------------------------------------ 006 C 参数声明 007 integer plhs(*), prhs(*) 008 integer nlhs, nrhs 009 C 指针 010 integer ir, jc , pr 011 C 实际数据 012 real*8 pdata(6) 013 data pdata /3, 1, 10, -7, 2, -2/ 014 integer irdata(6) 015 data irdata / 0, 3, 0, 1, 2, 3 / 016 integer jcdata(5) 017 data jcdata / 0, 2, 2, 5, 6 / 018 C------------------------------------------ 019 integer mxGetIr,mxGetJc,mxGetPr 020 integer mxCreateSparse 021 integer mxCopyReal8ToPtr 022 integer mxCopyInteger4ToPtr 023 C 代码行 024 plhs(1) = mxCreateSparse(4 , 4 , 6 , 0) 025 ir = mxGetIr(plhs(1)) 026 jc = mxGetJc(plhs(1)) 027 pr = mxGetPr(plhs(1)) 028 C 内存复制 029 call mxCopyReal8ToPtr(pdata,pr,6) 030 call mxCopyInteger4ToPtr(irdata,ir,6) 031 call mxCopyInteger4ToPtr(jcdata,jc,5) 032 C 函数文件尾部 033 return 034 end 　　本例子的核心部分从024行开始，首先使用mxCreateSparse函数创建了稀疏矩阵的mxArray数据对象，然后分别使用mxGetIr、mxGetJc、mxGetPr三个函数获取稀疏矩阵的三个重要的指针ir、jc和pr，这里是实数的稀疏矩阵，所以没有pi指针。接着，使用mxCopy函数完成数据内存的复制。这里需要用户注意的是整数类型的数据拷贝，在MATLAB的外部接口函数库中，具有三个不同的函数用于整数类型数据的拷贝，分别是： ● mxCopyInteger1ToPtr。 ● mxCopyInteger2ToPtr。 ● mxCopyInteger4ToPtr。 同样，还有三个函数用于将指针数据拷贝给相应的整数类型数据，这里就不一一列举了。需要读者注意的是，针对具体的整数数据类型须选取不同的拷贝函数。下面编译运行该文件： mex fcreatesparse.f A = fcreatesparse A = (1,1) 3 (4,1) 1 (1,3) 10 　 (2,3) -7 　 (3,3) 2 (4,4) -2 full(A) ans = 3 0 10 0 0 0 -7 0 0 0 2 0 　 1 0 0 -2 　　 例4-4 复数数据。 Fortran语言的数据类型是比较丰富的，其中，其独有的复数数据类型是该语言的一个优势，使其成为了早期计算机数值编程领域首选的编程语言。在MATLAB中，同样有相应的数据类型，不过Fortran语言的复数和MATLAB的复数之间进行交互的时候，只能够通过mxCopy函数完成数据的传递。例4-4演示了Fortran语言MEX文件处理复数的方法。下面是程序的源代码： 00