matlab和fortran混合编程

提示：扫一扫查出行
【扫一扫了解最新限行尾号】
复制提示

matlab2016b+vs2010+ivf2013+f90

其实默认是f77语法，但通过配置可以改变为自由格式。

默认只能f77代码，怎样修改：

https://ww2.mathworks.cn/matlabcentral/answers/333875-i-am-trying-to-create-a-mex-file-for-fortran90-code-i-get-an-error-because-of-and-long-sentence

找到：mex_FORTRAN_win64，删掉fixed就可以用f90了。

大体过程：

1.在matlab中配置fortran编译环境： mex -setup FORTRAN
2.写好源代码，自定义子函数仍可以像往常一样调用别的子函数
3.将matlab当前路径调整到当前源代码文档所在路径
4.在matlab命令行输入编译命令：mex csuherror.f90
5.像使用matlab内置函数一样使用此函数，error=csuherror(X),其中X为csuh的9个参数（包括C），输出的为计算误差
6.若需要调试，需要这样编译：mex -g csuherror.f90 ，之后通过vs附加matlab到进程，并打开源代码加断点，在matlab运行时就可以调试了，具体请参考帮助文档。

下面是一个例子：

先写fortran代码：传入一个数组和一个数字，输出一个数组

Add.f90 ,注意此文档名为将来matlab调用的函数名，调用函数名时与其内的子函数名称Add无关。

#include "fintrf.h" !必须有的头文件,里面有mxGetPr, mxGetM, mxGetN,mxCreateDoubleMatrix等函数的申明 http://yuehongliang2000.lofter.com/post/1d001431_ab72071 subroutine mexFunction(OutSum,OutVar,InSum,InVar)!函数接口名称必须为mexFunction, !OutSum:输出参数个数 !OutVar:输出参数数组指针 !InSum:输入参数个数 !InVar:输入参数数组指针 !参数顺序不能随意更改
Integer InSum,OutSum mwPointer InVar(*),OutVar(*) !mwPointer专门用于表示指针变量,这个不能随意用Integer代替 mwPointer mxGetPr, mxGetM, mxGetN,mxCreateDoubleMatrix !这个对返回指针函数的再次申明,
integer, parameter :: fp = selected_real_kind(,)
Integer , parameter :: myINT = SELECTED_INT_KIND( ) Real(fp),Allocatable::x(:,:),y(:,:)
real(fp)::z Integer m,n If(InSum/=)Then call mexErrMsgIdAndTxt('MATLAB:InputTooBig','输入参数个数必须为2个') Return EndIf m=mxGetM(InVar())!获取第1个输入参数的行数 n=mxGetN(InVar())!获取第1个输入参数的列数 Allocate(x(m,n),y(m,n)) Call mxCopyPtrToReal8(mxGetPr(InVar()),x,m*n)!将第1个参数数组赋值给x变量 Call mxCopyPtrToReal8(mxGetPr(InVar()),z,)!将第2个整数变量赋值给z Call Add(x,y,z,m,n)!调用内部函数 OutVar()=mxCreateDoubleMatrix(m,n,)!给返回参数分配内存 Call mxCopyReal8ToPtr(y,mxGetPr(OutVar()),m*n)!将返回参数赋值给分配的内存 DeAllocate(x,y)!释放临时分配的内存 Return End SubRoutine ! SubRoutine Add(x,y,z,m,n)
!
! Integer,Intent(In)::m,n
! integer, parameter :: fp = selected_real_kind(,)
!
! Real(fp),Intent(In)::x(m,n),z
!
! Real(fp),Intent(Out)::y(m,n)
!
! Integer i,j
!
! Do i=,m
!
! Do j=,n
!
! y(i,j)=x(i,j)+z
!
! EndDo
!
! EndDo
!
! Return
!
!End SubRoutine SubRoutine Add(x,y,z,m,n) Integer,Intent(In)::m,n
integer, parameter :: fp = selected_real_kind(,) Real(fp),Intent(In)::x(m,n),z Real(fp),Intent(Out)::y(m,n) Integer i,j Do i=,m Do j=,n y(i,j)=x(i,j)+z EndDo EndDo Return End SubRoutine