| 查看: 688 | 回复: 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 |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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 (幼儿园)
- 金币: 4677.7
- 红花: 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
40