版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

onexf

金蟲 (小有名氣)

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

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

相關(guān)的mathematica 公式已經(jīng)給寫好，就是不清楚怎么做，請(qǐ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ù)此樓

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

大開眼界

» 猜你喜歡

求調(diào)劑已經(jīng)有7人回復(fù)
272求調(diào)劑已經(jīng)有7人回復(fù)
352分化工與材料已經(jīng)有3人回復(fù)
0703化學(xué) 已經(jīng)有6人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有3人回復(fù)
求調(diào)劑已經(jīng)有3人回復(fù)
266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑已經(jīng)有4人回復(fù)
一志愿南師大0703化學(xué) 275求調(diào)劑已經(jīng)有4人回復(fù)
085701環(huán)境工程求調(diào)劑已經(jīng)有3人回復(fù)
266求調(diào)劑已經(jīng)有10人回復(fù)

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

求幫忙作圖，急�。�！已經(jīng)有16人回復(fù)
液相色譜圖主峰本來很好，一個(gè)月后突然出現(xiàn)肩峰，求大俠幫忙，附圖。已經(jīng)有49人回復(fù)
請(qǐng)幫忙分析兩張色譜圖已經(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ù)
請(qǐng)問XRD的局部放大的圖是怎么做出來的？已經(jīng)有4人回復(fù)
用CASTEP做出來的一個(gè)能帶圖，大家?guī)兔纯此膶?dǎo)電性已經(jīng)有21人回復(fù)
【請(qǐng)教】幫忙分析一張N2氣吸脫附的圖，不閉合已經(jīng)有9人回復(fù)
LaTex使用問題，幫忙畫兩個(gè)圖標(biāo) 已經(jīng)有3人回復(fù)
【求助】請(qǐng)高手幫忙分析一下這張固化DSC曲線圖，非常感謝已經(jīng)有18人回復(fù)
【求助】誰能幫忙分析一張多肽二級(jí)質(zhì)譜圖，必有重謝！已經(jīng)有17人回復(fù)

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

查看全部散金貼

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

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

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)給寫好，就是不清楚怎么做，請(qǐ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的回帖

zhouyadan2

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 330.6
帖子: 136
在線: 20.4小時(shí)
蟲號(hào): 1217449

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

哇，好牛�。�

回復(fù)此樓

6樓2011-03-24 08:47:20

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

onexf

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1269.7
帖子: 173
在線: 155.3小時(shí)
蟲號(hào): 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)丟失了)。請(qǐ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的回帖

mshwangg

至尊木蟲 (正式寫手)

★ ★
ben_ladeng(金幣+2): 可以在辛苦點(diǎn) 2011-09-05 20:56:25

引用回帖:

Originally posted by onexf at 2011-03-24 19:52:58:
謝謝你，我用SphericalPlot3D[rad,alpha,beta]做出來的一個(gè)3D球
最大值這個(gè)還可以（如下）：

最小會(huì)值這個(gè)怎么會(huì)有些奇點(diǎn)出現(xiàn)呢？（如下）
[eimg]de/bc/45 ...

最小值的奇點(diǎn)問題恐怕需要你自己去解決了，我看過在你給定的取值范圍內(nèi)最小值（極小值）有兩個(gè)，但數(shù)值是一樣的，你再考慮下計(jì)算過程和參數(shù)是否正確合理吧。
關(guān)于圖切開的問題，要是簡單的切開（或者說畫一半）倒是簡單，就是c和d圖那種，只需要將三維球坐標(biāo)范圍設(shè)定到0到Pi就可以的，剩下的就是改變顏色的事了，看看幫助應(yīng)該能搞定。但是a和b那種剖面的沒搞出來。原本想畫個(gè)平面然后和a圖合到一起的，然而用球坐標(biāo)我不會(huì)畫平面

那位蟲子教教我？最近比較忙，腦筋有點(diǎn)不夠用了

贊一下(1人)

回復(fù)此樓

9樓2011-03-29 22:13:09

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

rachel9610

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1455.9
帖子: 58
在線: 60.7小時(shí)
蟲號(hào): 587705

★ ★ ★
onexf(金幣+1):謝謝參與
微塵、夢(mèng)想(金幣+2): 謝謝回復(fù)！ 2011-05-01 18:25:08
onexf(金幣+30): xx! 2011-05-02 13:13:59

引用回帖:

I have writen the following code, and the singluar points are disppeared.
Try again, Plz!

C11 = 120;
C12 = 100;
C44 = 60;

n1 = Cos[a] Cos;
n2 = Sin[a] Cos;
n3 = Sin;

m1 = Cos[v] Cos[w];
m2 = Sin[v] Cos[w];
m3 = Sin[w];
AA = - (Cos[a] Cos[v] + Sin[a] Sin[v])/Tan;
Dir = Cos[
   b]^2 (Cos[a]^2  Cos[v]^2 + Sin[a]^2  Sin[v]^2)/(1 +
   AA^2) + (Sin AA)^2/(1 + AA^2);

p = n1^4 + n2^4 + n3^4;

you = 1/(1/3/(C11 + 2*C12) - (1 - 3*p)/3/(C11 - C12) + (1 - p)/2/C44);
((C11 + 2*C12)*(C11 - C12 - 2*C44)*Dir +
   2*C12*C44)/(2*(C11 + 2*C12)*(C11 - C12)*C44)*you;
radi = ((C11 + 2*C12)*(C11 - C12 - 2*C44)*Dir + 2*C12*C44)*
you/(2*(C11 + 2*C12)*(C11 - C12)*C44);
ParametricPlot3D[{Abs[
NMaximize[{radi, 0 + 0.01 < v < \[Pi] - 0.01}, v][[1]]]*n1,
  Abs[NMaximize[{radi, 0 + 0.01 < v < \[Pi] - 0.01}, v][[1]]]*n2,
  Abs[NMaximize[{radi, 0 < v < \[Pi]}, v][[1]]]*n3}, {a, 0,
  2 \[Pi]}, {b, -\[Pi]/2 + 0.01, \[Pi]/2}, PlotPoints -> 40]

贊一下(2人)

回復(fù)此樓

10樓2011-05-01 13:27:47

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

mengxc

木蟲 (著名寫手)

應(yīng)助: 3 (幼兒園)
金幣: 4128.1
帖子: 1342
在線: 491.6小時(shí)
蟲號(hào): 1148508

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

都這么厲害啊，

回復(fù)此樓

11樓2011-06-06 16:50:11

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

huycwork

金蟲 (著名寫手)

程序強(qiáng)帖: 22
應(yīng)助: 0 (幼兒園)
金幣: 953
帖子: 1080
在線: 264.1小時(shí)
蟲號(hào): 1257243

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

這軟件很好玩啊，必須的學(xué)下。頂起先~

贊一下(1人)

回復(fù)此樓

12樓2011-06-09 15:36:51

已閱關(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的回帖

mchen10

金蟲 (正式寫手)

應(yīng)助: 13 (小學(xué)生)
金幣: 1508.2
帖子: 425
在線: 159.5小時(shí)
蟲號(hào): 695973

帖子真精彩!
已經(jīng)收錄到淘貼專輯《大開眼界》

贊一下

回復(fù)此樓

16樓2011-09-06 03:05:50

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

dxwbucea

鐵蟲 (著名寫手)

應(yīng)助: 20 (小學(xué)生)
金幣: 17156.4
帖子: 1722
在線: 80.4小時(shí)
蟲號(hào): 1586258

★ ★
onexf(金幣+1):謝謝參與
jjdg(金幣+1): 感謝參與 2012-02-07 11:05:46

這里還真是藏龍臥虎，有很多的高人啊，應(yīng)該多來這里學(xué)習(xí)學(xué)習(xí)。
本人也知道一些Mathematica軟件，但水平還沒有這么高。
頂一下。

贊一下(2人)

回復(fù)此樓

17樓2012-02-07 11:03:42

已閱關(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 ...

牛

k-li5樓

2011-03-24 08:41

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

zhouyadan27樓

2011-03-24 08:47

lxzhq14樓

2011-09-05 12:30

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

八戒的小弟15樓

2011-09-05 17:23

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 onexf 的主題更新

返回列表

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

[考研] 328求調(diào)劑 +6	嗯滴的基本都 2026-03-27	6/300	2026-03-27 16:49 by sanrepian
[考研] 一志愿北化085600材料專碩275\|有文章專利｜求調(diào)劑 +3	Micky11223 2026-03-25	3/150	2026-03-27 14:52 by 醉在風(fēng)里
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過A區(qū)國家線 +10	哈哈157349 2026-03-21	10/500	2026-03-27 13:06 by zzll406
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-27	3/150	2026-03-27 10:17 by guoweigw
[考研] 304材料求調(diào)劑 +4	鐘llll 2026-03-26	4/200	2026-03-27 03:42 by wxiongid
[考研] 調(diào)劑求收留 +7	果然有我 2026-03-26	7/350	2026-03-27 00:26 by wxiongid
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 0703化學(xué)求調(diào)劑 +3	丹青奶蓋 2026-03-26	5/250	2026-03-26 20:11 by macy2011
[考研] 學(xué)碩274求調(diào)劑 +3	Li李魚 2026-03-26	3/150	2026-03-26 18:32 by 提出方法的提出?/a>
[考研] 22 350 本科985求調(diào)劑，求老登收留 +4	李軼男003 2026-03-20	4/200	2026-03-26 16:05 by 哇啦啦啦xtj
[考研] 290分調(diào)劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] 086003食品工程求調(diào)劑 +6	淼淼111 2026-03-24	6/300	2026-03-25 10:29 by 3Strings
[考研] 0703化學(xué)求調(diào)劑 +6	奶油草莓. 2026-03-22	7/350	2026-03-25 10:00 by shangxh
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +7	天天好運(yùn)來上岸?/a> 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 080500求調(diào)劑 +3	zzzzfan 2026-03-24	3/150	2026-03-24 16:38 by barlinike
[考研] 085404電子信息284分求調(diào)劑 +4	13659058978 2026-03-24	4/200	2026-03-24 12:15 by syl20081243
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過 +6	幸運(yùn)哩哩 2026-03-22	10/500	2026-03-22 20:10 by edmund7
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng

24小時(shí)熱門版塊排行榜

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

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

» 猜你喜歡

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

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

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

查看全部散金貼