版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

onexf

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1269.7
帖子: 173
在线: 155.3小时
虫号: 454769

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

相关的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 ]

回复此楼

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

大开眼界

» 猜你喜欢

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

求帮忙作图，急！！！已经有16人回复
液相色谱图主峰本来很好，一个月后突然出现肩峰，求大侠帮忙，附图。已经有49人回复
请帮忙分析两张色谱图已经有4人回复
实验得到一张质谱图，麻烦各位大哥帮忙解析一下这，含有结构式。已经有12人回复
已有核磁图，但有一个峰不知是什么，求高手帮忙，谢谢！已经有17人回复
帮忙用XCrySDen出个图吧. 已经有4人回复
大神来帮忙解释一下，上下两个图是什么意思已经有10人回复
帮忙看一下下面这个谱图是不是我要的？为什么11ppm左右有一个峰？已经有8人回复
求高手帮解一个图谱，谢谢已经有7人回复
帮忙分析一个HR-TEM图！已经有11人回复
以前没做过，求帮忙看下DMA图已经有8人回复
帮忙分析一个XRD图已经有4人回复
急急急，哪位大神帮忙分析一下透镜图~~看不懂~汗一个先已经有20人回复
编辑来信，让提供一个图作为TOC，这是什么意思？已经有8人回复
做一个核磁图谱，有提供原始数据帮忙做成图！已经有36人回复
请问XRD的局部放大的图是怎么做出来的？已经有4人回复
用CASTEP做出来的一个能带图，大家帮忙看看它的导电性已经有21人回复
【请教】帮忙分析一张N2气吸脱附的图，不闭合已经有9人回复
LaTex使用问题，帮忙画两个图标已经有3人回复
【求助】请高手帮忙分析一下这张固化DSC曲线图，非常感谢已经有18人回复
【求助】谁能帮忙分析一张多肽二级质谱图，必有重谢！已经有17人回复

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

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

mengxc

木虫 (著名写手)

应助: 3 (幼儿园)
金币: 4128.1
帖子: 1342
在线: 491.6小时
虫号: 1148508

★
onexf(金币+1):谢谢参与

都这么厉害啊，

回复此楼

11楼2011-06-06 16:50:11

已阅关注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[]
计算结果图:

就是算得很慢

赞一下(3人)

回复此楼

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

金虫 (小有名气)

应助: 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球
最大值这个还可以（如下）：

最小会值这个怎么会有些奇点出现呢？（如下）

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

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

赞一下

回复此楼

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

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

简单回复

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 ...

牛

查看全部 17 个回答