 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (5009)
- >導(dǎo)師招生 (879)
- >考研 (668)
- >文獻(xiàn)求助 (330)
- >蟲(chóng)友互識(shí) (227)
- >考博 (162)
- >休閑灌水 (155)
- >基金申請(qǐng) (144)
- >碩博家園 (103)
- >招聘信息布告欄 (92)
- >博后之家 (75)
- >綠色求助(高懸賞) (46)
- >論文投稿 (32)
- >找工作 (28)
- >SciFinder/Reaxys (23)
- >外文書(shū)籍求助 (15)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 575 | 回復(fù): 3 | ||
ilovexiaomu金蟲(chóng) (小有名氣)
|
[求助]
mathematica中如何實(shí)現(xiàn)c++程序轉(zhuǎn)換為mathematica語(yǔ)句 已有1人參與
|
|
各位蟲(chóng)友： 我遇到了一個(gè)比較復(fù)雜的2重積分，我的積分函數(shù)表達(dá)式如下： 如圖片所示 在對(duì)上面的kx及ky求二重積分時(shí)無(wú)法用mathematica中的語(yǔ)句直接計(jì)算，因?yàn)榉e分函數(shù)相當(dāng)復(fù)雜，用直接的mathematica語(yǔ)句無(wú)法積分，在網(wǎng)上查閱了一個(gè)用c++編寫(xiě)出來(lái)的版本，現(xiàn)在我想把他轉(zhuǎn)換成mathematica的語(yǔ)句。希望對(duì)此熟悉的朋友能幫忙轉(zhuǎn)換一下。 c++語(yǔ)句如下： /* ============================================================= * * 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 或者各位蟲(chóng)友能有更好的解決這個(gè)積分的辦法，也可以麻煩告知，多謝 積分函數(shù).png |
回復(fù)此樓» 猜你喜歡
0854 考研調(diào)劑 招生了！AI 方向
已經(jīng)有14人回復(fù)
-
328求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑
已經(jīng)有5人回復(fù)
-
招08考數(shù)學(xué)
已經(jīng)有14人回復(fù)
-
一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑
已經(jīng)有4人回復(fù)
-
0854電子信息求調(diào)劑 324
已經(jīng)有3人回復(fù)
-
303求調(diào)劑
已經(jīng)有3人回復(fù)
-
336求調(diào)劑
已經(jīng)有3人回復(fù)
-
306求調(diào)劑
已經(jīng)有9人回復(fù)
-
工科0856求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- mathematica生成矩陣
已經(jīng)有17人回復(fù)
- mathematica如何給出不能積分的函數(shù)的積分結(jié)果
已經(jīng)有3人回復(fù)
- 如何在已知一階導(dǎo)數(shù)但不知原函數(shù)的情況下用mathematica求高階導(dǎo)數(shù)，求高手指導(dǎo)
已經(jīng)有6人回復(fù)
- 一個(gè)涉及廣義函數(shù)的積分
已經(jīng)有8人回復(fù)
- 用含有定積分的隱函數(shù)擬合數(shù)據(jù)，matlab或者origin
已經(jīng)有18人回復(fù)
- Mathematica關(guān)于繪制隱函數(shù)方程組的問(wèn)題，跪求，快瘋了
已經(jīng)有11人回復(fù)
- 數(shù)值積分方法求助
已經(jīng)有6人回復(fù)
- Mathematica里的compile怎么使用
已經(jīng)有27人回復(fù)
- 溫度的分段函數(shù)如何導(dǎo)入fluent
已經(jīng)有7人回復(fù)
- 請(qǐng)教mathematic問(wèn)題？關(guān)于解微分方程組的和畫(huà)圖的一些問(wèn)題。
已經(jīng)有19人回復(fù)
- 如何計(jì)算德拜函數(shù)？求助
已經(jīng)有10人回復(fù)
- 用Wofram Mathematica解方程的問(wèn)題
已經(jīng)有9人回復(fù)
- 請(qǐng)教一個(gè)積分用留數(shù)定理計(jì)算的問(wèn)題
已經(jīng)有15人回復(fù)
- mathematica請(qǐng)教計(jì)算結(jié)果問(wèn)題
已經(jīng)有4人回復(fù)
- Delta函數(shù)如何積分
已經(jīng)有3人回復(fù)
- 【求助】【求助】mathematica中的一種數(shù)值積分問(wèn)題
已經(jīng)有5人回復(fù)
- 【原創(chuàng)】在WORD中使用Mathematica畫(huà)的圖形的體會(huì)
已經(jīng)有13人回復(fù)
- 【求助】Mathematica怎么對(duì)不等式解的最小值集合作圖？
已經(jīng)有8人回復(fù)
- 論文投稿時(shí)，用Mathematica畫(huà)函數(shù)圖形如何更光滑
已經(jīng)有11人回復(fù)
xzczd
木蟲(chóng) (小有名氣)
- 應(yīng)助: 72 (初中生)
- 金幣: 1941.4
- 紅花: 11
- 帖子: 290
- 在線: 144.2小時(shí)
- 蟲(chóng)號(hào): 3305534
- 注冊(cè): 2014-07-03
2樓2015-05-13 10:38:53
walk1997
金蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 4676.2
- 紅花: 22
- 帖子: 1066
- 在線: 798.1小時(shí)
- 蟲(chóng)號(hào): 416039
- 注冊(cè): 2007-06-29
- 性別: GG
- 專業(yè): 粒子物理學(xué)和場(chǎng)論
3樓2015-05-13 16:23:42
ilovexiaomu
金蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1125.1
- 紅花: 1
- 帖子: 234
- 在線: 68.5小時(shí)
- 蟲(chóng)號(hào): 443234
- 注冊(cè): 2007-10-27
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
4樓2015-05-19 11:38:12
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ilovexiaomu 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +3 | 元夕元 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 333求調(diào)劑 +6 | 87639 2026-03-21 | 10/500 |
|
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +5 | vv迷 2026-03-22 | 5/250 |
|
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +6 | 皮皮皮123456 2026-03-22 | 6/300 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研]
 一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑
+3
|
Grand777 2026-03-21 | 3/150 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語(yǔ) 2026-03-20 | 3/150 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|