Matlab以MEX方式“调用”C源代码
- 格式:doc
- 大小:59.00 KB
- 文档页数:12
C语言与matlab接口MATLAB作为世界顶尖的数学应用软件,以其强大的工程计算、算法研究、工程绘图、应用程序开发、数据分析和动态仿真等功能,在航空航天、机械制造和工程建筑等领域发挥着越来越重要的作用。
而C语言功能丰富,使用灵活方便,目标程序效率高。
既有高级语言的优点,又有低级语言的特点。
因此,C语言是目前应用最广的编程语言。
虽然MATLAB是一个完整的、功能齐全的编程环境,但在某些情况下,与外部环境的数据和程序的交互是非常必须而且有益的。
众所周知,MATLAB是用M语言编程,不能在M文件中直接调用C语言程序。
可以通过MATLAB提供的应用编程接口(API)来实现与外部的接口,在MATLAB 环境中实现调用C语言或Fortran程序、输入或输出数据以及和其他软件程序间建立客户/服务器关系等功能。
MATLAB中调用C语言程序,必须通过MEX文件来实现。
一、C语言的MEX文件的结构C语言的MEX文件是一种动态连接子程序,可以象调用M文件一样调用它。
MEX文件主要有以下几方面的应用:(1)在MATLAB中,M文件的计算速度特别是循环迭代的速度远比C语言慢,因此可以把要求大量循环迭代的部分用C语言编写为MEX文件,提高计算速度。
(2)已经开发的C语言程序,则不必将其转化为M文件而重复劳动,通过添加入口程序mexFunction,可以由MA TLAB调用。
(3)直接控制硬件,如A/D采集卡,D/A输出卡等,以用于数据采集或控制应用。
C语言的MEX文件的源程序由两个非常明显部分组成:(1)计算程序,即在MEX文件中完成计算功能的程序代码,计算程序可以是普通的C语言程序,按照C语言规则编写即可(2)入口程序,将计算程序与MATLAB连接的入口函数mexFunction。
入口程序mexFunction相对要复杂一些,函数中有四个参数nlhs、plhs、nrhs和prhs。
这里nlhs是输出数据的个数,plhs是指向mxArray(MATLAB中所有数据都由mxArray定义)的输出数据的指针,nrhs是输入数据的个数,prhs是指向mxArray 的输入数据的指针。
Matlab调用c程序方法:第一步:要先在matlab中安装c程序编译器,步骤如下:键入命令:mex -setup;选择c/c++编译器;选择c/c++编译器版本;确认。
第二步:键入:mex *.c实例介绍:【转】如果我有一个用C语言写的函数,实现了一个功能,如一个简单的函数:double add(double x, double y) {return x + y;}现在我想要在Matlab中使用它,比如输入:>> a = add(1.1, 2.2)3.3000要得出以上的结果,那应该怎样做呢?解决方法之一是要通过使用MEX文件,MEX文件使得调用C函数和调用Matlab的内置函数一样方便。
MEX文件是由原C代码加上MEX文件专用的接口函数后编译而成的。
可以这样理解,MEX文件实现了一种接口,它把在Matlab中调用函数时输入的自变量通过特定的接口调入了C函数,得出的结果再通过该接口调回Matlab。
该特定接口的操作,包含在mexFunction这个函数中,由使用者具体设定。
所以现在我们要写一个包含add和mexFunction的C文件,Matlab调用函数,把函数中的自变量(如上例中的1.1和2.2)传给mexFunction的一个参数,mexFunction把该值传给add,把得出的结果传回给mexFunction的另一个参数,Matlab通过该参数来给出在Matlab语句中调用函数时的输出值(如上例中的a)。
比如该C文件已写好,名为add.c。
那么在Matlab中,输入:>> mex add.c就能把add.c编译为MEX文件(编译器的设置使用指令mex -setup),在Windows 中,MEX文件类型为mexw32,即现在我们得出add.mexw32文件。
现在,我们就可以像调用M函数那样调用MEX文件,如上面说到的例子。
所以,通过MEX文件,使用C函数就和使用M函数是一样的了。
关于MEX函数的说明.MEX⽂件是⼀种可在matlab环境中调⽤的C(或fortran)语⾔衍⽣程序,MEX⽂件的后缀名按32位/64位分别为 .mexw32/.mexw64。
MEX ⽂件是由C或Fortran语⾔编写的源代码,经matlab处理⽽⽣成的。
它是可以被matlab解释器⾃动装载并执⾏的程序,类似windows下的dll⽂件。
MEX⽂件实现了⼀种其它语⾔与MATLAB的接⼝,通过MEX⽂件可以在MATLAB中像调⽤内嵌函数⼀样调⽤使⽤和等语⾔编写的函数,实现了代码重⽤,同时也能提⾼MATLAB环境中数据处理的。
MEX⽂件是⼀种计算机⽂件,它提供MATLAB或Octave之间的接⼝以及⽤C,C ++或Fortran编写的函数。
它代表“MATLAB可执⾏⽂件”。
编译时,MEX⽂件是动态加载的,允许从MATLAB或Octave中调⽤外部函数,就像它们是内置函数⼀样。
为了⽀持MEX⽂件的开发,MATLAB和Octave都提供了外部接⼝功能,便于在MEX⽂件和⼯作区之间传输数据。
除了MEX⽂件之外,Octave还有⾃⼰的格式,使⽤⾃⼰的原⽣API,性能更好。
#include "mex.h"void timestwo_alt(double *y, double x){*y = 2.0*x;}void mexFunction ( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[] ){double *M;int m,n; //指针指向第⼀个参数的数据地址M = mxGetPr(prhs[0]);m = mxGetM(prhs[0]);n = mxGetN(prhs[0]);plhs[0] = mxCreateDoubleMatrix(m,n, mxINT32_CLASS,mxREAL); //⽣成m x n 的实矩阵,分配内存空间double *A;A = mxGetPr(plhs[0]);timestwo_alt(A,*M); //调⽤并直接赋值到指针指向的输出变量}编译链接C语⾔的MEX⽂件源程序,在MATLAB的控制窗⼝中输⼊:mex timestwoalt.c⽣成⼀个名为timestwoalt.mexw32的MEX⽂件;运⾏:在MATLAB的控制窗⼝中输⼊timestwoalt(10)结果为ans =20.0000 + 0.0000iMEX⽂件的编写和编译需要两个基本条件:⼀是必须按照MATLAB要求配置接⼝组件和相关⼯具,⼆是要有C语⾔或Fortran语⾔的。
Matlab以MEX方式调用C源代码(2009-04-23 20:10:15)转载标签:杂谈如果我有一个用C语言写的函数,实现了一个功能,如一个简单的函数:double add(double x, double y){return x + y;}现在我想要在Matlab中使用它,比如输入:>> a = add(1.1, 2.2)3.3000要得出以上的结果,那应该怎样做呢?解决方法之一是要通过使用MEX文件,MEX文件使得调用C函数和调用Matlab的内置函数一样方便。
MEX文件是由原C代码加上MEX文件专用的接口函数后编译而成的。
可以这样理解,MEX文件实现了一种接口,它把在Matlab中调用函数时输入的自变量通过特定的接口调入了C函数,得出的结果再通过该接口调回Matlab。
该特定接口的操作,包含在mexFunction这个函数中,由使用者具体设定。
所以现在我们要写一个包含add和mexFunction的C文件,Matlab调用函数,把函数中的自变量(如上例中的1.1和2.2)传给mexFunction的一个参数,mexFunction把该值传给add,把得出的结果传回给mexFunction的另一个参数,Matlab 通过该参数来给出在Matlab语句中调用函数时的输出值(如上例中的a)。
值得注意的是,mex文件是与平台有关的,以我的理解,mex文件就是另类的动态链接库。
在matlab6.5中使用mex -v 选项,你可以看到最后mex阶段有类似如下的信息:--> "del _lib94902.obj"--> "del "test.exp""--> "del "test.lib""也就是说,虽然在matlab6.5生成的是dll文件,但是中间确实有过lib文件生成。
1.Maltab 编译Mex文件(不可调试)1.1hello word在matlab命令行中输入:mex -setup选择合适编译器,完成设置。
新建一个文件,输入:#include <mex.h>void mexFunction(int nlhs ,mxArray *plhs[] ,int nrhs, const mxArray *prhs[]){mexPrintf("Hello word!");}把上面程序保存为一个C文件,如hello.c在command window下运行>> mex hello.c>>结束之后在current folder 窗口里会有一个hello.mexw32文件。
在command window下运行>> hello()Hello word!1.2如果是多维数组并不能用Arr[i][j]的格式寻址,在参数在C里只是一维数组,一维数组按列排列,C语言是按行排列的,如[a b c d;e f g h]C语言在内存中是按a b c d e f g h排列,而matlab是按 a e b f c g d h排列的。
如果要变成C语言可Arr[i][j]索引的格式,可用for循环转化打印函数:mexPrintf()相当于matlab的disp函数,mexPrintf函数里需要打印的部分必须用""而不能用''子函数的问题:子函数有两种方式,一种采用return值得方式返回值另一种方式为利用参数返回值,把需要return的值直接赋给形参。
多个.c文件的编译:多个.c文件中只能存在一个mexFunction函数可用mex ('a.c','b.c','d.c')的格式编译2.在VS中调试mex(1)在vs中建立dll工程,建立project.cpp文件,输入#include "project.h"#include <mex.h>void mexFunction(int nlhs ,mxArray *plhs[] ,int nrhs, const mxArray *prhs[]){mexPrintf("Hello word!");}建立.h文件,输入#include "matrix.h"#include "mex.h"#define TESTMATLAB_EXPORTS#ifdef TESTMATLAB_EXPORTS#define MEX_FUNCTION_API __declspec(dllexport)#else#define MEX_FUNCTION_API __declspec(dllimport)#endifMEX_FUNCTION_API void mexFunction(int nlhs, mxArray* plhs[], int nrhs, mxArray* prhs[]);建立project.def文件,输入LIBRARYEXPORTSmexFunction(2)添加include和lib路径若matlab为32位,liblujing为:D:\Program Files\MATLAB\R2014a\extern\lib\win3264位则为:D:\Program Files\MATLAB\R2014a\extern\lib\win64\microsoftinclude路径为:D:\Program Files\MATLAB\R2014a\extern\include(3)添加链接器附加依赖项libmx.lib libeng.lib libmat.lib libmex.lib,将连接器输出的文件名称后缀,project.dll改为project.mexw32,64位改为project.mexw64(4)设置C/C++-代码生成-运行库为多线程调试(/Mtd)(5)若为64位matlab,还需将vs的编译平台设置为x64(6)生成解决方案,会在debug目录生成project.mexw64.(7)打开matlab,将工作目录改到如上debug目录,包含project.mexw64文件(8)开始调试,在vs中cpp文件设置端点,选择工具-附件到进程-matlab-添加,再在matlab中运行该文件,即可跳入端点处.(9)调试完成后,将project.mexw64放入matlab的目录中,可以直接调用。
使用MATLAB的MCC命令生成CC程序使用MATLAB的MCC命令生成C/C++程序首先,配置MATLAB用如下命令:mex –setupmbuild -setup例一:将m文件转化成库文件使用1、建立一个名为ceshidll.m的M函数文件,该函数的功能是输入两组数完成两组数据的插值拟合,并将结果用图形表示:ceshidll.m文件内容如下:function ceshidll(x,y)a=min(x):0.1:max(x);b = interp1(x,y,a,''spline'');%一维插值函数plot(x,y,''*'',a,b);2、在MATLAB Command中用如下命令编译函数ceshidll.m:>> mcc-t -W libhg:dlltest-T link:lib -h libmmfile.mlib ceshidll.m参数说明:-t 将m文件编译为C\C++代码-W libhg:dlltest 生成使用C函数图形库的文件,生成的文件名为dlltest-T link:lib 生成库文件-h 辅助选项,可以将任何被调用的辅助函数都包含到编译的文件中libmmfile.mlib连接任何需要的共享函数库ceshidll.m 被编译的文件名编译完成后在MATLAB当前目录下会生成以下文件:ceshidll.c、ceshidll.h 、dlltest.c 、dlltest.exports、dlltest.h、dlltest.mlib、dlltest.exp、dlltest.lib、dlltest.dll。
其中dlltest.h 、dlltest.lib和dlltest.dll文件是我们需要的。
使用方法:#include "matlab.h"#include "dlltest.h"#pragma comment(lib,"dlltest")关键代码:UpdateData(TRUE); //刷新输入数据double X[100],Y[100];CString AA,BB,a;int i=1;mxArray *A=NULL; //初始化矩阵mxArray *B=NULL;AA=m_edit1; //字符串赋值BB=m_edit2;.....//将字符转化为数字mlfEnterNewContext(0, 0);//自动管理内存dlltestInitialize();mlfCeshidll(A,B); //调用dll文件中函数mxDestroyArray(A); //释放矩阵内存mxDestroyArray(B);mlfRestorePreviousContext(0, 0);例二:将m文件转换成对应的C\C++文件1、在MATLAB中编写如下函数: function [x]=gjfcz(A,b)%A=[-1.5 1 2; 4 2 3 ; -3 2 8]%b=[3;5;6]x=A\b保存名为gjfcz.m,该函数的功能为求解线形方程组,可参考第四章的内容。
matlab调⽤c语⾔步骤以及传递图像数组时遇到的bug@⽬录前⾔该⽂章为介绍了matlab调⽤c语⾔函数的步骤以及遇到的⼀个坑,开始⽬的是⽤c语⾔去处理图像,matlab显⽰图像的处理结果,⽅便调参,然后再把C代码移植到单⽚机上,但是最后发现很不⽅便,所以最后还是放弃了,但是也学习了⼀点东西,所以记录⼀下.matlab版本为matlab2020a.⼀、MEXmex是matlab提供的⼀个接⼝,简单来说,MEX-file是⼀种预编译的,⽤其他语⾔(C/C++,Fortran)编写的函数库,可以直接被Matlab调⽤。
⼆、使⽤步骤1.安装c/c++编译器在附加功能中安装MinGW-w642.新建.c⽂件在⽂件头部include ''mex.h"3.⾃定义函数在.c⽂件中写⼊⾃⼰想写的c语⾔函数4.在.c⽂件中定义MexFunction函数MexFunction是c语⾔和MATLAB的接⼝函数,这⼀步可将c语⾔和matlab对接上,该函数是在.c⽂件中。
整个c程序由⼀个接⼝⼦过程 mexFunction构成,前⾯提到过,Matlab的mex函数有⼀定的接⼝规范,就是指:void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])nlhs:输出参数数⽬plhs:指向输出参数的指针nrhs:输⼊参数数⽬mxGetScalar(prhs[0]) :把通过prhs[0]传递进来的mxArray类型的指针指向的数据(标量)赋给C程序⾥的变量;mxGetPr(prhs[0]) :从指向mxArray类型数据的prhs[0]获得了指向double类型的指针mxGetM(prhs[0]):获得矩阵的⾏数mxGetN(prhs[0]):获得矩阵的列数mxCreateDoubleMatrix(int m, int n, mxComplexity ComplexFlag) :实现内存的申请,m:待申请矩阵的⾏数; n:待申请矩阵的列数以我为例,我在matlab输⼊的是110*180的⼆维图像数组,输出也为⼆维图像数组.这⼀步最重要的是要检查参数传递的正确性,这⾥我遇到⼀个⼤坑,花了我半天的时间才找到,这后边单独说。
Matlab中mcc, mbuild和mex命令详解MATLAB编译及接口Matlab的mcc命令将m文件转换成C/C++文件和相应的MEX包裹文件(需要Matlab编译器)和exe文件,并在特定条件下可以自动调用mbuild或者mex。
Mex命令将入口函数为mexfunction()的C/C++文件编译成MEX(DLL)文件(需要C/C++编译器)。
Mbuild命令将入口函数为main()的C/C++文件生成可独立运行的C应用程序(需要C/C++编译器)。
对于64位的Windows平台,MATLAB没有现成的C编译器。
需要免费下载Microsoft Windows SDK编译器。
Microsoft Windows SDK for Windows 7 and .NET Framework 4下载地址:(需要联网安装):microsoft/en-us/download/details.aspx?id=8279安装前一定要卸载Microsoft Visual C++ 2010 x86 Redistributable - 10.0.30319以上版本的C++运行库,装完SDK7.1后需要安装VC-Compiler-KB2519277补丁::microsoft/zh-cn/download/details.aspx?id=4422 在使用mcc, mbuild和mex命令之前,首先要配置MATLAB使用Windows SDK进行C/C++语言编译,配置命令如下:mex -setup17mbuild -setup配置命令运行结果:>> mex -setupMEX 配置为使用'Microsoft Visual C++ 2012 (C)' 以进行C 语言编译。
警告: MATLAB C 和Fortran API 已更改,现可支持包含2^32-1 个以上元素的MATLAB 变量。
不久以后,您需要更新代码以利用新的API。
Matlab以MEX方式“调用”C源代码阅读() 评论() 发表时间:2008年09月26日 00:46本文地址:/blog/3311493-1222361217如果我有一个用C语言写的函数,实现了一个功能,如一个简单的函数:double add(double x, double y){return x + y;}现在我想要在Matlab中使用它,比如输入:>> a = add(1.1, 2.2)3.3000要得出以上的结果,那应该怎样做呢?解决方法之一是要通过使用MEX文件,MEX文件使得调用C函数和调用Matlab的内置函数一样方便。
MEX文件是由原C代码加上MEX文件专用的接口函数后编译而成的。
可以这样理解,MEX文件实现了一种接口,它把在Matlab中调用函数时输入的自变量通过特定的接口调入了C函数,得出的结果再通过该接口调回Matlab。
该特定接口的操作,包含在mexFunction这个函数中,由使用者具体设定。
所以现在我们要写一个包含add和mexFunction的C文件,Matlab调用函数,把函数中的自变量(如上例中的1.1和2.2)传给 mexFunction的一个参数,mexFunction把该值传给add,把得出的结果传回给mexFunction的另一个参数,Matlab通过该参数来给出在Matlab 语句中调用函数时的输出值(如上例中的a)。
值得注意的是,mex文件是与平台有关的,以我的理解,mex文件就是另类的动态链接库。
在matlab6.5中使用mex -v 选项,你可以看到最后mex阶段有类似如下的信息:--> "del _lib94902.obj"--> "del "test.exp""--> "del "test.lib""也就是说,虽然在matlab6.5生成的是dll文件,但是中间确实有过lib文件生成。
比如该C文件已写好,名为add.c。
那么在Matlab中,输入:>> mex add.c就能把add.c编译为MEX文件(编译器的设置使用指令mex-setup),在Windows中,MEX 文件类型为mexw32,即现在我们得出add.mexw32文件。
现在,我们就可以像调用M函数那样调用MEX文件,如上面说到的例子。
所以,通过MEX文件,使用C函数就和使用M函数是一样的了。
我们现在来说mexFunction怎样写。
mexFunction的定义为:void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) {/*....................................*/}可以看到,mexFunction是没返回值的,它不是通过返回值把结果传回Matlab的,而是通过对参数plhs的赋值。
mexFunction的四个参数皆是说明Matlab调用MEX文件时的具体信息,如这样调用函数时:>> b = 1.1; c = 2.2;>> a = add(b, c)mexFunction四个参数的意思为:nlhs = 1,说明调用语句左手面(lhs-left hand side)有一个变量,即a。
nrhs = 2,说明调用语句右手面(rhs-right hand side)有两个自变量,即b和c。
plhs是一个数组,其内容为指针,该指针指向数据类型mxArray。
因为现在左手面只有一个变量,即该数组只有一个指针,plhs[0]指向的结果会赋值给a。
prhs和plhs类似,因为右手面有两个自变量,即该数组有两个指针,prhs[0]指向了b,prhs[1]指向了c。
要注意prhs是const的指针数组,即不能改变其指向内容。
因为Matlab最基本的单元为array,无论是什么类型也好,如有double array、 cell array、struct array……所以a,b,c都是array,b = 1.1便是一个1x1的double array。
而在C 语言中,Matlab的array使用mxArray类型来表示。
所以就不难明白为什么plhs和prhs 都是指向mxArray类型的指针数组。
完整的add.c如下:#include "mex.h" // 使用MEX文件必须包含的头文件// 执行具体工作的C函数double add(double x, double y){return x + y;}// MEX文件接口函数void mexFunction(int nlhs,mxArray *plhs[], int nrhs,const mxArray *prhs[]){double *a;double b, c;plhs[0] = mxCreateDoubleMatrix(1, 1, mxREAL);a = mxGetPr(plhs[0]);b = *(mxGetPr(prhs[0]));c = *(mxGetPr(prhs[1]));*a = add(b, c);}mexFunction的内容是什么意思呢?我们知道,如果这样调用函数时:>> output = add(1.1, 2.2);在未涉及具体的计算时,output的值是未知的,是未赋值的。
所以在具体的程序中,我们建立一个1x1的实double矩阵(使用 mxCreateDoubleMatrix函数,其返回指向刚建立的mxArray的指针),然后令plhs[0]指向它。
接着令指针a指向plhs [0]所指向的mxArray 的第一个元素(使用mxGetPr函数,返回指向mxArray的首元素的指针)。
同样地,我们把prhs[0]和prhs [1]所指向的元素(即1.1和2.2)取出来赋给b和c。
于是我们可以把b 和c作自变量传给函数add,得出给果赋给指针a所指向的mxArray中的元素。
因为a是指向plhs[0]所指向的mxArray的元素,所以最后作输出时,plhs[0]所指向的mxArray赋值给output,则 output便是已计算好的结果了。
上面说的一大堆指向这指向那,什么mxArray,初学者肯定都会被弄到头晕眼花了。
很抱歉,要搞清楚这些乱糟糟的关系,只有多看多练。
实际上mexFunction是没有这么简单的,我们要对用户的输入自变量的个数和类型进行测试,以确保输入正确。
如在add函数的例子中,用户输入char array便是一种错误了。
从上面的讲述中我们总结出,MEX文件实现了一种接口,把C语言中的计算结果适当地返回给Matlab罢了。
当我们已经有用C编写的大型程序时,大可不必在 Matlab里重写,只写个接口,做成MEX文件就成了。
另外,在Matlab程序中的部份计算瓶颈(如循环),可通过MEX文件用C语言实现,以提高计算速度。
以上是对mex文件的初步认识,下面详细介绍如何用c语言编写mex文件:1 为什么要用C语言编写MEX文件 MATLAB是矩阵语言,是为向量和矩阵操作设计的,一般来说,如果运算可以用向量或矩阵实现,其运算速度是非常快的。
但若运算中涉及到大量的循环处理,MATLAB的速度的令人难以忍受的。
解决方法之一为,当必须使用for循环时,把它写为MEX文件,这样不必在每次运行循环中的语句时MATLAB都对它们进行解释。
2 编译器的安装与配置要使用MATLAB编译器,用户计算机上应用事先安装与MATLAB适配的以下任何一种ANSIC/C++编译器:5.0、6.0版的MicroSoft Visual C++(MSVC)5.0、5.2、5.3、5.4、5.5版的Borland C++LCC(由MATLAB自带,只能用来产生MEX文件)下面是安装与配置MATLAB编译器应用程序MEX的设置的步骤:(1)在MATLAB命令窗口中运行mex–setup,出现下列提示:Please choose your compiler for building external interface (MEX) files:Would you like mex to locate installed compilers [y]/n?(2)选择y,MATLAB将自动搜索计算机上已安装的外部编译器的类型、版本及所在路径,并列出来让用户选择:Select a compiler:[1] Borland C++Builder version 6.0 in C:\Program Files\Borland[2] Digital Visual Fortran version 6.0 in C:\Program Files\Microsoft Visual Studio[3] Lcc C version 2.4 in D:\MATLAB6P5P1\sys\lcc[4] Microsoft Visual C/C++ version 6.0 in C:\Program Files\Microsoft Visual Studio [0] NoneCompiler:(3)选择其中一种(在这里选择了3),MATLAB让用户进行确认:Please verify your choices:Compiler: Lcc C 2.4Location: D:\MATLAB6P5P1\sys\lccAre these correct?([y]/n):(4)选择y,结束MATLAB编译器的配置。
3 一个简单的MEX文件例子【例1】用m文件建立一个1000×1000的Hilbert矩阵。
ticm=1000;n=1000;a=zeros(m,n);for i=1:1000for j=1:1000a(i,j)=1/(i+j);endendtoc在matlab中新建一个Matlab_1.cpp 文件并输入以下程序:#include "mex.h"//计算过程void hilb(double *y,int n){int i,j;for(i=0;i<n;i++)for(j=0;j<n;j++)*(y+j+i*n)=1/((double)i+(double)j+1);}//接口过程void mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[]) {double x,*y;int n;if (nrhs!=1)mexErrMsgTxt("One inputs required.");if (nlhs != 1)mexErrMsgTxt("One output required.");if (!mxIsDouble(prhs[0])||mxGetN(prhs[0])*mxGetM(prhs[0])!=1)mexErrMsgTxt("Input must be scalars.");x=mxGetScalar(prhs[0]);plhs[0]=mxCreateDoubleMatrix(x,x,mxREAL);n=mxGetM(plhs[0]);y=mxGetPr(plhs[0]);hilb(y,n);}该程序是一个C语言程序,它也实现了建立Hilbert矩阵的功能。