 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 576 | 回復(fù): 3 | ||
ilovexiaomu金蟲 (小有名氣)
|
[求助]
mathematica中如何實(shí)現(xiàn)c++程序轉(zhuǎn)換為mathematica語句 已有1人參與
|
|
各位蟲友： 我遇到了一個(gè)比較復(fù)雜的2重積分，我的積分函數(shù)表達(dá)式如下： 如圖片所示 在對(duì)上面的kx及ky求二重積分時(shí)無法用mathematica中的語句直接計(jì)算，因?yàn)榉e分函數(shù)相當(dāng)復(fù)雜，用直接的mathematica語句無法積分，在網(wǎng)上查閱了一個(gè)用c++編寫出來的版本，現(xiàn)在我想把他轉(zhuǎn)換成mathematica的語句。希望對(duì)此熟悉的朋友能幫忙轉(zhuǎn)換一下。 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 或者各位蟲友能有更好的解決這個(gè)積分的辦法，也可以麻煩告知，多謝 積分函數(shù).png |
回復(fù)此樓» 猜你喜歡
0856材料專碩353求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
293求調(diào)劑
已經(jīng)有4人回復(fù)
-
306求0703調(diào)劑一志愿華中師范
已經(jīng)有10人回復(fù)
-
【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】
已經(jīng)有7人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
材料與化工328分調(diào)劑
已經(jīng)有4人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有20人回復(fù)
-
一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271
已經(jīng)有6人回復(fù)
-
085404電子信息284分求調(diào)劑
已經(jīng)有3人回復(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ù)方程組的問題，跪求，快瘋了
已經(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問題？關(guān)于解微分方程組的和畫圖的一些問題。
已經(jīng)有19人回復(fù)
- 如何計(jì)算德拜函數(shù)？求助
已經(jīng)有10人回復(fù)
- 用Wofram Mathematica解方程的問題
已經(jīng)有9人回復(fù)
- 請(qǐng)教一個(gè)積分用留數(shù)定理計(jì)算的問題
已經(jīng)有15人回復(fù)
- mathematica請(qǐng)教計(jì)算結(jié)果問題
已經(jīng)有4人回復(fù)
- Delta函數(shù)如何積分
已經(jīng)有3人回復(fù)
- 【求助】【求助】mathematica中的一種數(shù)值積分問題
已經(jīng)有5人回復(fù)
- 【原創(chuàng)】在WORD中使用Mathematica畫的圖形的體會(huì)
已經(jīng)有13人回復(fù)
- 【求助】Mathematica怎么對(duì)不等式解的最小值集合作圖？
已經(jīng)有8人回復(fù)
- 論文投稿時(shí)，用Mathematica畫函數(shù)圖形如何更光滑
已經(jīng)有11人回復(fù)
xzczd
木蟲 (小有名氣)
- 應(yīng)助: 72 (初中生)
- 金幣: 1941.4
- 紅花: 11
- 帖子: 290
- 在線: 144.2小時(shí)
- 蟲號(hào): 3305534
- 注冊(cè): 2014-07-03
2樓2015-05-13 10:38:53
walk1997
金蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 4676.2
- 紅花: 22
- 帖子: 1066
- 在線: 798.1小時(shí)
- 蟲號(hào): 416039
- 注冊(cè): 2007-06-29
- 性別: GG
- 專業(yè): 粒子物理學(xué)和場(chǎng)論
3樓2015-05-13 16:23:42
ilovexiaomu
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1125.1
- 紅花: 1
- 帖子: 234
- 在線: 68.5小時(shí)
- 蟲號(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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +3 | 超級(jí)伊昂大王 2026-03-24 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +7 | A-哆啦Z夢(mèng) 2026-03-23 | 12/600 |
|
|
|
[考研] 335分 | 材料與化工專碩 | GPA 4.07 | 有科研經(jīng)歷 +4 | cccchenso 2026-03-23 | 4/200 |
|
|
|
[考研] 350求調(diào)劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 3/150 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|