[交流] 【求助】幫忙作一個圖~

相關(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ù)此樓

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

查看全部散金貼

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

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

rachel9610

金蟲 (小有名氣)

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

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

引用回帖:

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

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

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的回帖

查看全部 17 個回答

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[]
計算結(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小時
蟲號: 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]做出來的一個3D球
最大值這個還可以（如下）：

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

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

[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。�！ +3	祁祺祺 2026-03-25	3/150	2026-03-27 13:25 by shangxh
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過A區(qū)國家線 +10	哈哈157349 2026-03-21	10/500	2026-03-27 13:06 by zzll406
[考研] 275求調(diào)劑 +5	jjjjjjjjjjl 2026-03-27	5/250	2026-03-27 12:52 by lature00
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 269專碩求調(diào)劑 +9	金恩貝 2026-03-21	9/450	2026-03-27 11:25 by 王笑宇一定是研?/a>
[考研] 307求調(diào)劑 +8	超級伊昂大王 2026-03-24	8/400	2026-03-27 11:23 by leaiy
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-26	3/150	2026-03-27 07:58 by chemisry
[碩博家園] 招收生物學(xué)/細胞生物學(xué)調(diào)劑 +3	IceGuo 2026-03-26	4/200	2026-03-27 05:35 by user003
[考研] 343求調(diào)劑 +4	贈我一本書 2026-03-23	4/200	2026-03-27 00:40 by wxiongid
[考研] 303求調(diào)劑 +7	元夕元 2026-03-20	8/400	2026-03-26 20:38 by 不吃魚的貓
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner
[考研] 考研調(diào)劑 +6	來好運來來來 2026-03-21	7/350	2026-03-25 22:43 by 418490947
[考研] 材料與化工304求B區(qū)調(diào)劑 +3	邱gl 2026-03-25	3/150	2026-03-25 19:03 by Ainin_
[考研] 材料專碩 335 分求調(diào)劑 +4	拒絕冷暴力 2026-03-25	4/200	2026-03-25 18:45 by haxia
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +7	天天好運來上岸?/a> 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 調(diào)劑 +4	13853210211 2026-03-24	4/200	2026-03-24 19:44 by ms629
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen