[求助] 求助三元函数二阶偏微分C语言写法

2楼: Originally posted by huycwork at 2011-08-31 12:59:59:
你把算法直观地说一遍，我帮你写C语言要的不？

3楼: Originally posted by zhenzhen5516 at 2011-08-31 16:01:08:
非常感谢！为了您方便看，我把算法上传到附件里了，麻烦您下载一下！只要您写出一项我应该就能明白了~

4楼: Originally posted by huycwork at 2011-09-01 09:12:31:
呃，你太客气了……
我记得偏微分好像是常规微分在各个方向上微分之和，但你给出的符号都是啥意思呃？哪些是函数，哪些是变量，k是什么意思，各个函数跟哪些相关变量相关，那个下三角的意思又是什么呢？中括号括 ...

5楼: Originally posted by zhenzhen5516 at 2011-09-01 10:07:43:
哦，可能是我没说清楚。函数是u1=u1(x,y,z), K 是常量，u只是x y z的函数，下三角是二次微分算子，具体的我另上传了一个附件，可供参考~谢谢，等待您的回复

double func(double a, double b, double x){
    return a*x+b;
}

typedef struct{
    double a, b;
    double (*func)(double a, double b, double x);
} CFunc;
CFunc *factory(double a, double b){
    CFunc *f = malloc(sizeof(CFunc));
    f->a = a;
    f->b = b;
    f->func = func;
    return f;
}
double caller(CFunc *fc, double x){
    return fc->func(fc->a, fc->b, x);
}

6楼: Originally posted by huycwork at 2011-09-01 10:45:05:
我大概知道你的困难了，你看是不是这样：
混合偏导数的定义式df/dxdy = d(df/dy)/dx，df/dy部分返回一个函数，而C语言没有直接支持返回带参函数的机制。
如果是这样的话，你可以自己实现闭包函数即可。比如说返 ...

7楼: Originally posted by zhenzhen5516 at 2011-09-01 11:33:37:
呃，非常感谢您的详细回复。不过我没太看明白，其实我的问题就是df/dxdy = d(df/dy)/dx如何在程序里写进去，比如df2/dx2(表示函数F对x的二次微分）可以在程序里简单写成f(x+1,y,z)+f(x-1,y,z)-2f(x,y,z)，那么d( ...

Func: f(x, y, z);
Funcx: df/dx = Func(x+dx, y, z) - Func(x, y, z)   / dx;
Funcy: df/dy = Func(x, y+dy, z) - Func(x, y, z)   / dy;
Funcz: df/dz = Func(x, y, z+dz) - Func(x, y, z)   / dz;
Funcyx: Funcy/dx = Funcy(x+dx, y, z) - Funcy(x, y, z)  / dx;
Funczy: Funcz/dy
Funcxz: Funcx/dz

8楼: Originally posted by huycwork at 2011-09-01 13:53:10:
你没看明白我说的东西。你试试一些支持lambda表达式或者闭包的语言可能更能理解一些。
数学里面的函数完成的是一个映射，对应到C语言就是一个叫做函数的过程。两者的概念是一致的，即完成某种映射，但是C语言函 ...