matlab调用 fortran,matlab和fortran混合编程

matlab2016b+vs2010+ivf2013+f90

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

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

找到：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无关。

1 #include "fintrf.h" !必须有的头文件,里面有mxGetPr, mxGetM, mxGetN,mxCreateDoubleMatrix等函数的申明 http://yuehongliang2000.lofter.com/post/1d001431_ab72071

2

3

4 subroutine mexFunction(OutSum,OutVar,InSum,InVar)!函数接口名称必须为mexFunction,5

6

7 !OutSum:输出参数个数8

9

10 !OutVar:输出参数数组指针11

12

13 !InSum:输入参数个数14

15

16 !InVar:输入参数数组指针17

18

19 !参数顺序不能随意更改20 Integer InSum,OutSum21

22

23 mwPointer InVar(*),OutVar(*) !mwPointer专门用于表示指针变量,这个不能随意用Integer代替24

25

26 mwPointer mxGetPr, mxGetM, mxGetN,mxCreateDoubleMatrix !这个对返回指针函数的再次申明,27 integer, parameter :: fp = selected_real_kind(15,300)28 Integer , parameter :: myINT = SELECTED_INT_KIND( 8)29

30

31 Real(fp),Allocatable::x(:,:),y(:,:)32 real(fp)::z33

34

35 Integer m,n36

37

38

39

40

41 If(InSum/=2)Then42

43

44 call mexErrMsgIdAndTxt('MATLAB:InputTooBig','输入参数个数必须为2个')45

46

47 Return48

49

50 EndIf51

52

53

54

55

56 m=mxGetM(InVar(1))!获取第1个输入参数的行数57

58

59 n=mxGetN(InVar(1))!获取第1个输入参数的列数60

61

62 Allocate(x(m,n),y(m,n))63

64

65 Call mxCopyPtrToReal8(mxGetPr(InVar(1)),x,m*n)!将第1个参数数组赋值给x变量66

67

68 Call mxCopyPtrToReal8(mxGetPr(InVar(2)),z,1)!将第2个整数变量赋值给z69

70

71

72

73 Call Add(x,y,z,m,n)!调用内部函数74

75

76

77

78

79 OutVar(1)=mxCreateDoubleMatrix(m,n,0)!给返回参数分配内存80

81

82 Call mxCopyReal8ToPtr(y,mxGetPr(OutVar(1)),m*n)!将返回参数赋值给分配的内存83

84

85

86

87

88 DeAllocate(x,y)!释放临时分配的内存89

90

91 Return92

93

94 End SubRoutine95

96

97

98

99 ! SubRoutine Add(x,y,z,m,n)100 !101 ! Integer,Intent(In)::m,n102 ! integer, parameter :: fp = selected_real_kind(15,300)103 !104 ! Real(fp),Intent(In)::x(m,n),z105 !106 ! Real(fp),Intent(Out)::y(m,n)107 !108 ! Integer i,j109 !110 ! Do i=1,m111 !112 ! Do j=1,n113 !114 ! y(i,j)=x(i,j)+z115 !116 ! EndDo117 !118 ! EndDo119 !120 ! Return121 !122 !End SubRoutine123

124 SubRoutine Add(x,y,z,m,n)125

126

127 Integer,Intent(In)::m,n128 integer, parameter :: fp = selected_real_kind(15,300)129

130

131 Real(fp),Intent(In)::x(m,n),z132

133

134 Real(fp),Intent(Out)::y(m,n)135

136

137 Integer i,j138

139

140 Do i=1,m141

142

143 Do j=1,n144

145

146 y(i,j)=x(i,j)+z147

148

149 EndDo150

151

152 EndDo153

154

155 Return156

157

158 End SubRoutine

将matlab当前路径改为Add.f90同目录路径，然后 ： mex Add.f90   ,之后会生成Add.mex64

在matlab命令行窗口输入：

x=[1 2 3];

a=2;

y=Add(x,a)

---

我把Add.f90改为multiple.f90的话：

效果一样。

同样的功能的m文档，会优先调用.mexw64，且速度很快。所以可以写一个同名的m文档，只有函数名，之后全是注释文档，用以查阅。