查看: 1855  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

onexf

金虫 (小有名气)


[交流] 【求助】帮忙作一个图~



相关的mathematica 公式已经给写好,就是不清楚怎么做,请高手帮忙,给出源代码和结果(最好用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 ]
回复此楼

» 收录本帖的淘帖专辑推荐

大开眼界

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)


★
onexf(金币+1):谢谢参与
这软件很好玩啊,必须的学下。顶起先~
12楼2011-06-09 15:36:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

mshwangg

至尊木虫 (正式写手)


★ ★ ★ ★ ★ ★
onexf(金币+1):谢谢参与
cenwanglai(金币+5): 谢谢,呵呵,很有耐心。等lz给你bb哈! 2011-03-23 23:52:26
onexf(金币+30): 多谢,麻烦继续给力! 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:


相关的mathematica 公式已经给写好,就是不清楚怎么做,请高手帮忙,给出源代码和结果(最好用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[]
计算结果图:

就是算得很慢
2楼2011-03-23 23:45:40
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

circlewu

禁虫 (职业作家)

★
onexf(金币+1):谢谢参与
本帖内容被屏蔽

4楼2011-03-24 08:21:38
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

onexf

金虫 (小有名气)


引用回帖:
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球
最大值这个还可以(如下):

最小会值这个怎么会有些奇点出现呢?(如下)


另我看别人的文献中还可以将这种图型切开(球的里面和外面都不一样,还有立体感,如下),非常直观漂亮 (和文章作者交流过,他说用mathmatica做的,但源码已经丢失了)。请教高手怎么实现?


[ Last edited by onexf on 2011-3-24 at 19:56 ]
8楼2011-03-24 19:52:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
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 ...

牛
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +6 Isak0320 2026-10-12 8/400 2026-10-12 20:09 by caixiaoye1
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见