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

onexf

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

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

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

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

查看全部散金貼

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

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

circlewu

禁蟲(chóng) (職業(yè)作家)

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

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

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

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

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

mshwangg

至尊木蟲(chóng) (正式寫(xiě)手)

程序強(qiáng)帖: 5
應(yīng)助: 206 (大學(xué)生)
金幣: 10702.8
帖子: 597
在線: 195.4小時(shí)
蟲(chóng)號(hào): 576702

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

onexf

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1269.7
帖子: 173
在線: 155.3小時(shí)
蟲(chóng)號(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]做出來(lái)的一個(gè)3D球
最大值這個(gè)還可以（如下）：

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

另我看別人的文獻(xiàn)中還可以將這種圖型切開(kāi)（球的里面和外面都不一樣，還有立體感，如下），非常直觀漂亮 (和文章作者交流過(guò)，他說(shuō)用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

至尊木蟲(chóng) (正式寫(xiě)手)

程序強(qiáng)帖: 5
應(yīng)助: 206 (大學(xué)生)
金幣: 10702.8
帖子: 597
在線: 195.4小時(shí)
蟲(chóng)號(hào): 576702

★ ★
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]做出來(lái)的一個(gè)3D球
最大值這個(gè)還可以（如下）：

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

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

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

贊一下(1人)

回復(fù)此樓

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

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

簡(jiǎn)單回復(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ā)表

[考研] 299求調(diào)劑 +8	shxchem 2026-03-20	10/500	2026-03-27 11:27 by sanrepian
[考研] 材料與化工考研調(diào)劑 +16	孅華 2026-03-22	16/800	2026-03-27 11:23 by 王笑宇一定是研?/a>
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 275求調(diào)劑 +10	Micky11223 2026-03-25	12/600	2026-03-27 09:41 by caszguilin
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語(yǔ)已過(guò)四六級(jí)，六級(jí)521，求調(diào)劑 +3	美味可樂(lè)雞翅 2026-03-26	3/150	2026-03-27 08:46 by hxsm
[考研] 調(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 不吃魚(yú)的貓
[考研] 321求調(diào)劑 +6	wasdssaa 2026-03-26	6/300	2026-03-26 20:57 by sanrepian
[考研] 一志愿南京郵電大學(xué) 288分材料考研求調(diào)劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 打過(guò)很多競(jìng)賽，085406控制工程300分，求調(diào)劑 +3	askeladz 2026-03-26	3/150	2026-03-26 09:08 by 給你你注意休息
[考研] B區(qū)考研調(diào)劑 +4	yqdszhdap－ 2026-03-22	5/250	2026-03-25 08:51 by baoball
[考研] 311求調(diào)劑 +3	冬十三 2026-03-24	3/150	2026-03-24 21:31 by peike
[考研] 292求調(diào)劑 +4	鵝鵝鵝額額額額?/a> 2026-03-24	4/200	2026-03-24 16:41 by peike
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7	jhybd 2026-03-23	12/600	2026-03-24 09:03 by jhybd
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3	開(kāi)開(kāi)心心沒(méi)煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開(kāi)開(kāi)心心沒(méi)煩惱
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 298求調(diào)劑 +8	上岸6666@ 2026-03-20	8/400	2026-03-23 11:02 by laoshidan
[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無(wú)極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan