當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

onexf

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1269.7
帖子: 173
在線: 155.3小時(shí)
蟲號: 454769

[交流] 【求助】幫忙作一個(gè)圖~

相關(guān)的mathematica 公式已經(jīng)給寫好，就是不清楚怎么做，請高手幫忙，給出源代碼和結(jié)果(最好用7.0的版本)~
C1 = 120;
C2 = 100;
C3= 60;

Aa = - (Cos[alpha] Cos[v] + Sin[alpha] Sin[v])/Tan [beta];
Dir = Cos[beta]^2 (Cos[alpha]^2 Cos[v]^2 + Sin[alpha]^2 Sin[v]^2)/(1 + Aa^2) + (Sin[beta] Aa)^2/(1 + Aa^2);

p = (Cos[alpha] Cos[beta])^4 + (Sin[alpha] Cos[beta])^4 + Sin[beta]^4;
You = 1/(1/3/(C1 + 2*C2) - (1 - 3*p)/3/(C1 - C2) + (1 - p)/2/C3);

rad = ((C1 + 2*C2)*(C1 - C2 - 2*C3)*Dir + 2*C2*C3)*You/(2*(C1 + 2*C2)*(C1 - C2)*C3);

[ Last edited by onexf on 2011-3-23 at 19:50 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

大開眼界

» 猜你喜歡

283求調(diào)劑（080500）已經(jīng)有3人回復(fù)
286求調(diào)劑已經(jīng)有7人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有4人回復(fù)
308求調(diào)劑已經(jīng)有6人回復(fù)
275求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
287求調(diào)劑已經(jīng)有9人回復(fù)
274求調(diào)劑已經(jīng)有23人回復(fù)
085602 化工專碩 338分求調(diào)劑已經(jīng)有6人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有16人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

求幫忙作圖，急�。。� 已經(jīng)有16人回復(fù)
液相色譜圖主峰本來很好，一個(gè)月后突然出現(xiàn)肩峰，求大俠幫忙，附圖。已經(jīng)有49人回復(fù)
請幫忙分析兩張色譜圖已經(jīng)有4人回復(fù)
實(shí)驗(yàn)得到一張質(zhì)譜圖，麻煩各位大哥幫忙解析一下這，含有結(jié)構(gòu)式。已經(jīng)有12人回復(fù)
已有核磁圖，但有一個(gè)峰不知是什么，求高手幫忙，謝謝！已經(jīng)有17人回復(fù)
幫忙用XCrySDen出個(gè)圖吧. 已經(jīng)有4人回復(fù)
大神來幫忙解釋一下，上下兩個(gè)圖是什么意思已經(jīng)有10人回復(fù)
幫忙看一下下面這個(gè)譜圖是不是我要的？為什么11ppm左右有一個(gè)峰？已經(jīng)有8人回復(fù)
求高手幫解一個(gè)圖譜，謝謝已經(jīng)有7人回復(fù)
幫忙分析一個(gè)HR-TEM圖！已經(jīng)有11人回復(fù)
以前沒做過，求幫忙看下DMA圖已經(jīng)有8人回復(fù)
幫忙分析一個(gè)XRD圖已經(jīng)有4人回復(fù)
急急急，哪位大神幫忙分析一下透鏡圖~~看不懂~汗一個(gè)先已經(jīng)有20人回復(fù)
編輯來信，讓提供一個(gè)圖作為TOC，這是什么意思？已經(jīng)有8人回復(fù)
做一個(gè)核磁圖譜，有提供原始數(shù)據(jù)幫忙做成圖！已經(jīng)有36人回復(fù)
請問XRD的局部放大的圖是怎么做出來的？已經(jīng)有4人回復(fù)
用CASTEP做出來的一個(gè)能帶圖，大家?guī)兔纯此膶?dǎo)電性已經(jīng)有21人回復(fù)
【請教】幫忙分析一張N2氣吸脫附的圖，不閉合已經(jīng)有9人回復(fù)
LaTex使用問題，幫忙畫兩個(gè)圖標(biāo) 已經(jīng)有3人回復(fù)
【求助】請高手幫忙分析一下這張固化DSC曲線圖，非常感謝已經(jīng)有18人回復(fù)
【求助】誰能幫忙分析一張多肽二級質(zhì)譜圖，必有重謝！已經(jīng)有17人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-03-23 16:00:38

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

leimiao_hit

木蟲之王 (文學(xué)泰斗)

★ ★
onexf(金幣+1):謝謝參與
ben_ladeng(金幣+1): 謝謝參與 2011-09-05 20:56:44

引用回帖:

2樓: Originally posted by mshwangg at 2011-03-23 23:45:40:
源程序如下:

C1 = 120;
C2 = 100;
C3 = 60;
Aa = -(Cos[alpha] Cos[v] + Sin[alpha] Sin[v])/Tan[beta];
Dir = Cos[
beta]^2 (Cos[alpha]^2 Cos[v]^2 + Sin[alpha]^2 Sin[v]^2)/(1 +
Aa^ ...

速度太慢了啊

回復(fù)此樓

13樓2011-09-05 11:08:19

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 17 個(gè)回答

mshwangg

至尊木蟲 (正式寫手)

★ ★ ★ ★ ★ ★
onexf(金幣+1):謝謝參與
cenwanglai(金幣+5): 謝謝，呵呵，很有耐心。等lz給你bb哈！ 2011-03-23 23:52:26
onexf(金幣+30): 多謝，麻煩繼續(xù)給力！ 2011-03-24 19:53:57
ben_ladeng: 2011-09-05 20:55:54

引用回帖:

Originally posted by onexf at 2011-03-23 16:00:38:

相關(guān)的mathematica 公式已經(jīng)給寫好，就是不清楚怎么做，請高手幫忙，給出源代碼和結(jié)果(最好用7.0的版本)~
C1 = 120;
C2 = 100;
C3= 60;

Aa = - (Cos[alph ...

源程序如下:

C1 = 120;
C2 = 100;
C3 = 60;
Aa = -(Cos[alpha] Cos[v] + Sin[alpha] Sin[v])/Tan[beta];
Dir = Cos[
   beta]^2 (Cos[alpha]^2 Cos[v]^2 + Sin[alpha]^2 Sin[v]^2)/(1 +
   Aa^2) + (Sin[beta] Aa)^2/(1 + Aa^2);
p = (Cos[alpha] Cos[beta])^4 + (Sin[alpha] Cos[beta])^4 + Sin[beta]^4;
You = 1/(1/3/(C1 + 2*C2) - (1 - 3*p)/3/(C1 - C2) + (1 - p)/2/C3);
rad = ((C1 + 2*C2)*(C1 - C2 - 2*C3)*Dir + 2*C2*C3)*
You/(2*(C1 + 2*C2)*(C1 - C2)*C3);
Plot3D[NMaximize[{rad, 0 < v < 2 \[Pi]}, v][[1]], {alpha, 0,
  2 \[Pi]}, {beta, -\[Pi]/2, \[Pi]/2}]
Plot3D[NMinimize[{rad, 0 < v < 2 \[Pi]}, v][[1]], {alpha, 0,
  2 \[Pi]}, {beta, -\[Pi]/2, \[Pi]/2}]
Quit[]
計(jì)算結(jié)果圖:

就是算得很慢

贊一下(3人)

回復(fù)此樓

2樓2011-03-23 23:45:40

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

circlewu

禁蟲 (職業(yè)作家)

★
onexf(金幣+1):謝謝參與

本帖內(nèi)容被屏蔽

4樓2011-03-24 08:21:38

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

onexf

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1269.7
帖子: 173
在線: 155.3小時(shí)
蟲號: 454769

引用回帖:

Originally posted by mshwangg at 2011-03-23 23:45:40:
源程序如下:

C1 = 120;
C2 = 100;
C3 = 60;
Aa = -(Cos[alpha] Cos[v] + Sin[alpha] Sin[v])/Tan[beta];
Dir = Cos[
beta]^2 (Cos[alpha]^2 Cos[v]^2 + Sin[alpha]^2 Sin[v]^2)/(1 +
Aa ...

謝謝你，我用SphericalPlot3D[rad,alpha,beta]做出來的一個(gè)3D球
最大值這個(gè)還可以（如下）：

最小會(huì)值這個(gè)怎么會(huì)有些奇點(diǎn)出現(xiàn)呢？（如下）

另我看別人的文獻(xiàn)中還可以將這種圖型切開（球的里面和外面都不一樣，還有立體感，如下），非常直觀漂亮 (和文章作者交流過，他說用mathmatica做的，但源碼已經(jīng)丟失了)。請教高手怎么實(shí)現(xiàn)？

[ Last edited by onexf on 2011-3-24 at 19:56 ]

贊一下

回復(fù)此樓

8樓2011-03-24 19:52:58

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡單回復(fù)

zhouliangyu3樓

2011-03-24 00:00

onexf(金幣+1):謝謝參與

引用回帖:

Originally posted by mshwangg at 2011-03-23 23:45:40: 源程序如下: C1 = 120; C2 = 100; C3 = 60; Aa = -(Cos[alpha] Cos[v] + Sin[alpha] Sin[v])/Tan[beta]; Dir = Cos[ beta]^2 (Cos[alpha]^2 Cos[v]^2 + Sin[alpha]^2 Sin[v]^2)/(1 + Aa ...

牛

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當(dāng)當(dāng)
[考研] 求調(diào)劑推薦材料 304 +12	荷包蛋hyj 2026-03-26	12/600	2026-03-27 14:46 by laoshidan
[考研] 321求調(diào)劑 +4	璞玉~~ 2026-03-25	4/200	2026-03-27 11:41 by sanrepian
[考研] 材料學(xué)碩，求調(diào)劑 6+5	糖葫蘆888ll 2026-03-22	10/500	2026-03-27 08:18 by hypershenger
[考研] 325求調(diào)劑 +5	李嘉圖·S·路 2026-03-23	5/250	2026-03-27 00:42 by wxiongid
[考研] 349求調(diào)劑 +5	杰斯塔里斯 2026-03-21	5/250	2026-03-27 00:31 by wxiongid
[考研] 349求調(diào)劑 +4	李木子啊哈哈 2026-03-25	4/200	2026-03-26 22:49 by fmesaito
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8	AZMK 2026-03-25	10/500	2026-03-26 22:37 by 學(xué)員8dgXkO
[考研] 0703化學(xué)求調(diào)劑 +3	丹青奶蓋 2026-03-26	5/250	2026-03-26 20:11 by macy2011
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17	RichLi_ 2026-03-25	17/850	2026-03-26 19:44 by plmuchong
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4	JKSOIID 2026-03-26	4/200	2026-03-26 15:58 by 不吃魚的貓
[考研] 尋找調(diào)劑 +5	倔強(qiáng)芒? 2026-03-21	8/400	2026-03-26 13:25 by 0906ljy
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 各位老師您好：本人初試372分 +5	jj涌77 2026-03-25	6/300	2026-03-25 14:15 by mapenggao
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4	為學(xué)666 2026-03-24	4/200	2026-03-25 11:27 by BruceLiu320
[考研] 0703化學(xué)求調(diào)劑 +6	奶油草莓. 2026-03-22	7/350	2026-03-25 10:00 by shangxh
[考研] 070300化學(xué)求調(diào)劑 +9	苑豆豆 2026-03-20	9/450	2026-03-24 17:15 by licg0208
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky
[論文投稿] 急發(fā)核心期刊論文 +3	賢達(dá)問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹