| 查看: 514 | 回复: 3 | |||
ilovexiaomu金虫 (小有名气)
|
[求助]
mathematica中如何实现c++程序转换为mathematica语句 已有1人参与
|
|
各位虫友： 我遇到了一个比较复杂的2重积分，我的积分函数表达式如下： 如图片所示 在对上面的kx及ky求二重积分时无法用mathematica中的语句直接计算，因为积分函数相当复杂，用直接的mathematica语句无法积分，在网上查阅了一个用c++编写出来的版本，现在我想把他转换成mathematica的语句。希望对此熟悉的朋友能帮忙转换一下。 c++语句如下： /* ============================================================= * * trapzd2d. Computes the nth stage of refinement of an extended * * 2d trapezian rule. func is a pointer to a function with two * * real variables, it is to be integrated * * between ax<=x<=ay and ay<=y<=by. Works analogous to the one- * * dimensional trapzd (Numerical Recipes, p.137). * * ============================================================= */ typedef double real; /* can be put e.g. to float if single precission is sufficient */ #define FUNC(x,y) ((*func)(x,y)) real trapzd2d(real (*func)(real, real), real ax, real bx, real ay, real by, int n) { real x, y, tnm, sum, delx, dely; static real s; int it, j, k; if (n==1) { /* 1. stage: Function is evaluated only at the corners. */ return (s=0.25*(bx-ax)*(by-ay)*(FUNC(ax,ay)+FUNC(bx,ay)+FUNC(ax,by)+FUNC(bx,by))); } else { for (it=1, j=1; j<n-1; j++) it <<= 1; tnm = it; /* delx, dely are the spacings of the points to be added */ delx = (bx-ax)/tnm; dely = (by-ay)/tnm; x = ax + 0.5*delx; /* Evaluates function at the borders of the integration area */ for (sum=0.0, j=1; j<=it; j++, x+=delx) { sum += 0.5*(FUNC(x,ay) + FUNC(x,by)); } y = ay + 0.5*dely; for (j=1; j<=it; j++, y+=dely) { sum += 0.5*(FUNC(ax,y) + FUNC(bx,y)); } /* Interiour of the area */ x = ax + 0.5*delx; for (j=1; j<=it; j++, x+= delx) { y = ay + 0.5*dely; for (k=1; k<=it; k++, y+=dely) { sum += FUNC(x,y); } } x = ax + delx; for (j=1; j<it; j++, x+=delx) { y = ay + 0.5*dely; for (k=1; k<=it; k++, y+=dely) { sum += FUNC(x,y); } } x = ax + 0.5*delx; for (j=1; j<=it; j++, x+=delx) { y = ay + dely; for (k=1; k<it; k++, y+=dely) { sum += FUNC(x,y); } } /* Replace s by its refined value */ s = 0.25*(s + (bx-ax)*(by-ay)*sum/(tnm*tnm)); return s; } } #undef FUNC 或者各位虫友能有更好的解决这个积分的办法，也可以麻烦告知，多谢  积分函数.png |
回复此楼» 猜你喜欢
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有7人回复
-
需要合成515-64-0，50g，能接单的留言
已经有4人回复
-
自荐读博
已经有4人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
-
临港实验室与上科大联培博士招生1名
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- mathematica生成矩阵
已经有17人回复
- mathematica如何给出不能积分的函数的积分结果
已经有3人回复
- 如何在已知一阶导数但不知原函数的情况下用mathematica求高阶导数，求高手指导
已经有6人回复
- 一个涉及广义函数的积分
已经有8人回复
- 用含有定积分的隐函数拟合数据，matlab或者origin
已经有18人回复
- Mathematica关于绘制隐函数方程组的问题，跪求，快疯了
已经有11人回复
- 数值积分方法求助
已经有6人回复
- Mathematica里的compile怎么使用
已经有27人回复
- 温度的分段函数如何导入fluent
已经有7人回复
- 请教mathematic问题？关于解微分方程组的和画图的一些问题。
已经有19人回复
- 如何计算德拜函数？求助
已经有10人回复
- 用Wofram Mathematica解方程的问题
已经有9人回复
- 请教一个积分用留数定理计算的问题
已经有15人回复
- mathematica请教计算结果问题
已经有4人回复
- Delta函数如何积分
已经有3人回复
- 【求助】【求助】mathematica中的一种数值积分问题
已经有5人回复
- 【原创】在WORD中使用Mathematica画的图形的体会
已经有13人回复
- 【求助】Mathematica怎么对不等式解的最小值集合作图？
已经有8人回复
- 论文投稿时，用Mathematica画函数图形如何更光滑
已经有11人回复
2楼2015-05-13 10:38:53
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
3楼2015-05-13 16:23:42
ilovexiaomu
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1125.1
- 红花: 1
- 帖子: 234
- 在线: 68.5小时
- 虫号: 443234
- 注册: 2007-10-27
- 专业: 凝聚态物性I:结构、力学和
4楼2015-05-19 11:38:12