XPPAUT使用C语言产生的动态链接库

<gcc.exe 所在目录>\gcc -m32 -shared -o <目标文件名> <源文件名> 如果 gcc 找不到源文件，那么用源文件的全路径名。 3.指定.so 文件的查找路径 可用添加环境变量 LD_LABRARY_PATH 的值指定，DOS 命令行语句为 set LD_LABRARY_PATH=% LD_LABRARY_PATH %;<动态链接库所在的路径（注意 在路径最后加一个反斜杠"\"）> 4.ode 文件中嵌入调用动态链接库的操作语句（export 语句） 嵌入语句为： export {<实参数 1>，<实参数 2>,...}{<输出项 1>，输出项 2，...} 5.启动 xpppaut 后计算前加载动态链接库并指定调用的动态链接函数 加载方式为：依次单击 File->Edit->Load DLL 在弹出的对话框中输入编译的动态链接 库名，单击 ok，最后在 xppaut 主界面的命令行输入需要调用的函数名。 6.启动 xppaut 计算 下面给出一个简单的计算实例说明 xppaut 使用 C 生成的动态链接库的操作。
@ xlo=-5,ylo=-5,xhi=12,yhi=12 ====================================================================== 动态链接库源代码可写成 ====================================================================== /*这里文件名取为 mf.c*/ void myfun(double *in,double *out,int nin,int nout,double *var,double *con) { double x,y,alpha,beta; x=var[1]; y=var[2]; alpha=in[2]; beta=in[3]; out[0]=alpha*x+beta*y; out[1]=-beta*x+alpha*y; } ====================================================================== 再在 xpp.bat 文件开头加上下面 3 行代码并保存 xpp.bat 文件(也可以通过其他方式编译动态 链接库和设置动态链接库搜索路径)： ========================================= cd /d %~dp0 D:\CFree\mingw\bin\gcc -m32 -shared -o .\mf.so .\mf.c set LD_LIBRARY_PATH=% LD_LIBRARY_PATH %;. ========================================== 这里我们把 mf.c、usedll.ode、xpp.bat 放在相同目录，第一行表示跳转到当前目录，.\mf.so 和.\mf.c 斜杠前面的点和第 3 行结尾的点号都表示当前目录(DOS 命令行中点号是当前目录 的简写)。在我自己的电脑中 gcc.exe 所在目录为 D:\CFree\mingw\bin，这个目录要根据使用 者的 gcc.exe 的实际目录修改。 接下来进行如下操作 （1）把 usedll.ode 文件拖放操作到 xpp.bat 启动 xppaut （2）在 xppaut 主界面依次单击 File->Edit->Load DLL 在弹出的对话框中输入 mf.so 单击 ok 按钮 （3）在 xppaut 主界面的命令行输入需要调用的函数名 myfun (4) 依次单击 Initialconds->Go 启动计算 计算后可以得到如下的图示计算结果.
dx x y dt 对于求解动力系统 （ 1, 2 ， x (0) 10, y (0) 10 ） dy x y dt
我们可以书写普通的 ode 文件 ====================================================================== dx/dt=alpha*x+beta*y dy/dt=-beta*x+alpha*x par alpha=-1,beta=2 init x=10,y=10 @ total=50,dt=0.001,xp=x,yp=y @ xlo=-5,ylo=-5,xhi=12,yhi=12 ====================================================================== 现在我们用动态链接函数计算方程右边的表达式 x y 和 x y ，那么对应的 ode 文件可以写成： ====================================================================== # usedll.ode dx/dt=fx dy/dt=fy fx=0 fy=0 par alpha=-1,beta=2 init x=10,y=10 export {x,y,alpha,beta}{fx,fy} @ total=50,dt=0.001,xin[2],...分别对应 export 语句的<实参数 1>，<实参数 2>,<实参数 3>... out[0],out[1],out[2],...分别对应 export 语句的<输出项数 1>，<输出项 2>,<输出项 3>... nin 表示实参个数 nout 表示输出项的个数 var[0]表示时间变量 t，var[1],var[2],var[3],...依次对应 ode 文件顺序生成的存储变量（这里的 存储变量是指由方程默认定义的变量及由临时变量定义语句定义的临时变量） 如某个 ode 文件内容包含如下代码 ================================ dx[1..10]/dt=xp[j] dy[1..10]/dt=yp[j] xp[1..10]=0 yp[1..10]=0 par a=1,b=2,alpha=0.5 ... ================================ 那么 var[1],var[2],...var[10] 依次对应 x1,x2,...x10 var[11],var[12],...var[20] 依次对应 y1,y2,...y10 var[21],var[22],...var[30] 依次对应 xp1,xp2,...xp10 var[31],var[32],...var[40] 依次对应 yp1,yp2,...yp10 con[2],con[3],...依次顺序对应 ode 文件中的参数,上例中 con[2],con[3],con[4]依次对应 a,b, alpha（但实际操作来看库函数的 con 指针指向的对象不太明确，从而一般不要用这个指针） 按照函数的定义（按址传递） ，这里的"对应"应理解为共用的内存地址。 2.编译动态链接库为.so 文件 用 gcc 的编译代码为（前提是电脑上必须有 gcc 编译器,如果没有需要下载一个） ：
XPPAUT 允许读取 C 语言产生的（32 位）动态链接库，作用是： 1、简化 ode 文件，优化 ode 文件结构，使得 ode 文件可读性更强 2、动态链接函数实现 ode 文件代码不易或不能实现的操作 3、加快运算效率（这对耗时久的计算很有必要） 使用动态链接库的基本操作步骤为： 1.用 C 语言编写动态链接库源代码 可被调用的动态链接库源代码一般具有如下格式 =========================================== 文件包含语句（如果需要的话） 定义用户函数语句（如不需要可以不定义） 定义可能被 xppaut 调用的函数语句(可定义多个函数) ============================================= 其中：可能被 xppaut 调用的函数必须具有如下格式： void <函数名>（double *in,double *out,int nin,int nout, double *var, double *con）{函数体}